I think a spatial indexing scheme should be possible (and much faster), but the details are a bit elusive to me. I'm leaning towards the latter, but even with increased precision (64-digit precision via MATLAB variable precision arithmetic vpa), no intersecting facet is found.Īny suggestions for a more numerically stable projection of a ray onto the hyperplane would be much appreciated as well as any other general comments/suggestions on this approach.ĮDIT: Something else I've also considered is using some kind of spatial indexing scheme to identify the intersecting facet perhaps using linear inequalities in spherical coordinates or some type of projection (e.g. or 5., in that either the convex hull is incorrect due to numerical instability or the projection of the point onto the correct hyperplane is not close enough to the hyperplane to result in all positive barycentric coordinates, respectively. I think I have narrowed down the issue to either 2. identify the hyperplane where all $ \lambda_i \geq 0 $.simplex) and computing barycentric coordinates ( $ \lambda $) projecting the point onto each hyperplane (i.e.identifying all simplices with that nearest neighbor as a vertex.identifying the nearest neighbor in convex hull of a random datapoint MATLAB N-D NN search dsearchn.computing the convex hull of points on the surface of an n-sphere MATLAB N-D convexhull function convhulln.generating points on the surface of an n-sphere (or hyperorthant) based on MATLAB file exchange file hypersphere.In the end, I'm trying to do interpolation in octonion space (7-sphere) using barycentric coordinates which has involved: numerical precision of each other) if I first projected the datapoint onto the hyperplane per 1 and much worse if I used a unit norm datapoint on the surface of the n-sphere (e.g.
I implemented an alternative approach as described in Robust Barycentric Coordinates Computation of the Closest Point to a Hyperplane in E^n (Skala 2013), which I hoped would address the issue, but the results were nearly identical (i.e. at least one of every set of barycentric coordinates was negative). One of the main issues I am experiencing is sometimes no intersecting facet is detected (i.e. Initially, I tried a method to check if a ray intersects a hyperplane, and got it working in 7D Cartesian coordinates, but am running into (I think) numerical instability issues sometimes with 7D and more often in 8D (8D is what I need, also see discussion with author of ).
0 kommentar(er)
