# Tag Info

10

Randomized linear algebra might be something you'd like. It has direct applications in data analysis and is related to several branches of Mathematics such as Geometry (see the Johnson-Lindenstrauss lemma) and, of course, Probability. A starting point would be https://cacm.acm.org/magazines/2016/6/202647-randnla/abstract and the survey https://epubs.siam.org/...

9

Since your figure is a closed loop, its parametric curves $x(t)$ and $y(t)$ must be periodic functions. This suggests one way to generate such figures, by constructing random smooth periodic functions $x(t)$ and $y(t)$ via summation of sinusoids/harmonics with randomized amplitude and phase. Unfortunately, it would be difficult to guarantee such a figure ...

6

I was surprised for not receiving a satisfactory answer to the question above and my investigations showed me that this, indeed is an unexplored area. Hence, I put some effort in developing solutions to this problem, and published the following manuscripts: T. Birdal, B. Busam, N. Navab, S. Ilic and P. Sturm. "A Minimalist Approach to Type-Agnostic ...

6

Here's another interesting connection: There is a research area that looks into the complexity of computing matrix-matrix products from a geometric perspective. It's a totally bizarre connection, but quite fruitful. I'd start by looking at the papers by J. M. Landsberg at Texas A&M.

5

Choose n random points in the unit square. Compute the minimum point-to-point distance, $D$. Re-scale the space by $1/D$. The minimum point-to-point distance is now equal to 1. Edit: The best method to do step 2. depends on many things. If don't have too many points you can just compute (using matlab syntax): dmin=1; for i=1:n for j=i+1:n dmin = ...

5

The question really boils down to how far two points $x,y$ are from each other in direction $u$. This is easily answered: You need to compute the component of $y-x$ onto $u$, i.e., their (signed) distance in direction $u$ is $\frac{(y-x)\cdot u}{\|u\|}$. For your point cloud, you can then compute for each $i$ the (unsigned) distance $$d_i = \frac{(\bar x-... 4 I'm in the middle of doing this for myself. I think the most appropriate analogue to the Gaussian would be the heat kernel in hyperbolic space. Fortunately, this has been figured out before: https://www.math.uni-bielefeld.de/~grigor/nog.pdf (also available in a Bulletin of the London Mathematical Society). If you use the standard decay (e^{-dist^2/constant}... 4 How about something like radial integral channel features? Affordable person detection in omnidirectional cameras using radial integral channel features, Barış Evrim Demiröz, Albert Ali Salah, Yalin Bastanlar, Lale Akarun, Machine Vision and Appliactions 2019 https://www.cmpe.boun.edu.tr/~salah/demiroz19omnidirectional They would allow you to: Specify ... 4 You could also rewrite the expression to avoid square roots entirely. Reorder a,b so that a>b>0, and let the inequality be$$ \sqrt{a}-\sqrt{b}<\tau. $$This is equivalent to$$ a < \tau^2 + 2\tau\sqrt{b} + b \qquad\Leftrightarrow\qquad 2\tau\sqrt{b} > a-\tau^2-b, $$which is in turn equivalent to$$ \Big( a < \tau^2+b \quad\text{or}\...

4

As @Nicoguaro and @Paul have said in the comments to the question post, there are a great many ways to do this kind of thing, and I'm not sure if there is a single "best" approach. From a review study of Jonathan Richard Shewchuck at Berkley, an answer is: Please refer to the original document (version 31/12/2002) for symbology, terminology, special ...

4

I do not think that there exists an answer to this question in general, because it all depends on the intended use for the mesh. For instance, if you are doing computational fluid dynamics, you may want to have a mesh that is extremely anisotropic near the boundary layer. Now if you are doing computational electromagnetics, the best mesh will be probably ...

3

Here is a suggestion. Partition your problem into two parts: (1) Construct the intersection, (2) Find the smallest cube. (1) The intersection is a polytope defined by the union of your inequalities for both shapes. This is known as an H-representation of the polytope. There is software for constructing a polytope from its H-representation. polymake is one ...

3

If $n$ is too large and the side length of your squarer $a$ is fixed then there may not be a solution to your problem. If you are willing to increase the size of the square to fit the points then the solution by @DougLipinski will work. If the square size if fixed (which is why I assume you asked the question in the first place), then the re-scale step in ...

3

If you multiply the length of every side by a constant then all the angles will be preserved, this just produces a similar polygon. Multiply by a very small constant and you get a very small length, as small as you want, all the way down to zero. You can do this easily if you have the $(x,y)$ coordinates of the vertices of the polygon with perimeter $L$ and ...

3

An algorithm for this computation is described in Computational Geometry in C, Chapter 7, Section 6. Code is available at that link. Much of the tricky code is concerned with "degenerate" cases. Here is a nice Univ Montreal web page (Eric Plante) that describes a different algorithm: link.

3

The Haversine formula is really the best way to calculate distances on the Earth unless you have outstanding accuracy concerns. The next step passed the Haversine involves taking into account that the Earth is not a sphere but ellipsoidal and which involves a much more expensive computation. There's a good page with some embedded javascript at http://www....

3

What you're describing is a very specific way to solve the Laplace equation: using the Jacobi iteration to solve the five-point stencil when using finite differences. But there are much more general methods that would be difficult to derive from your starting point. A more promising approach is to start with the Poisson equation and ask how it can be ...

3

For me, a classic problem in linear algebra that got me climbing down a rabbit hole of complex analysis, conformal mappings, and polynomials, was in trying to prove converge rate bounds and iteration bounds for Krylov subspace methods. At a high level, it asks one of the following questions: What is the fastest way of computing $A^{-1}b$, under the ...

3

What I ended up doing is to download thousands of random images, smooth them with a Gaussian filter, and extract contours at different levels. I took 2-3 closed contours from each blurred image and packed them all into a collection of curves. My MATLAB code is available here: https://github.com/royvelich/smooth-curves-generator However, you will have to ...

2

You can use FEniCS: from fenics import ( UnitSquareMesh, FunctionSpace, Expression, interpolate, assemble, sqrt, inner, grad, dx, TrialFunction, TestFunction, Function, solve, DirichletBC, DomainBoundary, MPI, XDMFFile, ) # Create mesh and define function space mesh = UnitSquareMesh(100, ...

2

The problem of finding the smallest $a$-square that can contain $n$ points whose distance is not less than $1$ is equivalent to the circle packing in a square problem, see e.g Equivalence between circle packing and point packing. Info about this problem can be found also in Circle Packing and a list of optimal or best known solution in The best known ...

2

Here I devise a novel strategy, based on only 3D points, that I think, would work. I will parametrize a 3D plane by a point $\mathbf{p}$ and its normal $\mathbf{n}$. Imaging you take a pair of oriented points $\mathbf{p}_1$ and $\mathbf{p}_2$, on the point cloud $\mathbf{P}$ with corresponding normals $\mathbf{n}_1$ and $\mathbf{n}_2$. Let $\mathbf{d}$ ...

2

Generally speaking, the $\mathtt{sqrt}$ function is going to be the slowest part of that. Fortunately, if their lengths are the same, then the squares of their lengths are also the same. Therefore, you can just omit that part and compare the squares of the norm. I don't know of any way to make that go too much faster using cross-platform code. However, if ...

2

Given a small (infinitesimal) 2D displacement field $u = (u_x, u_y)$ the infinitesimal counterclockwise rotation angle $\theta$ is simply $$\theta = \frac 12 \left ( \frac{\partial u_y}{\partial x} - \frac{\partial u_x}{\partial y} \right )$$ In general, for $n$ dimensional displacement fields $u\in\mathbb{R}^n$ it can be shown that the infinitesimal ...

2

You can use the Bentley–Ottmann algorithm for this. Given a set of $n$ line segments with $k$ intersections, the algorithm can identify all intersections in $O((n+k)\log n)$ time and $O(n)$ space. In cases where $k=o(\frac{n^2}{\log n})$ (that is, cases in which $k$ has an appropriate upper bound) this offers time savings versus a naive $O(n^2)$ algorithm ...

2

Currently there is a GMSH API in the works: https://gitlab.onelab.info/gmsh/gmsh/tree/master/api Also, there are rumors that there will be a fully documented API by version 4.0. In short, there is no documentation on the functions for GMSH (except for what is in the source all ready). But, they are planning on creating a documented version by GMSH 4.0. ...

2

I want to submit a separate answer to this question, because there is a way to automatically simplify such formulas to remove square roots. I will use sage and QEPCAD (both free). The tutorial is here. The inequality $$|\sqrt{a}-\sqrt{b}|<t$$ is equivalent to the statement that there exist $u, v$ such that  \exists u,v: \qquad (u-v)^2 < t^2, \...

2

I don't fully understand where the problem is. Where the bonds are supposed to be, S-Se bonds? Like this, those distances don't seem much longer than Cd-S. In any case, you should not be troubled by the "sticks" connecting atoms in the visual representation: there are thresholds telling the program when to show them, and when not to show them. Actually, in ...

Only top voted, non community-wiki answers of a minimum length are eligible