cbm755 / cp_matrices Goto Github PK

In cpCutHole, there is a little nudge given in some rare conditions to prevent an infinite loop when certain normal vectors line up. This is detected with if dist2 < tol2....

This if statement is a vector because dist2 is generally all the points we're processing. This is wrong. IMHO, Matlab should raise an error here but it does not (Octave might, I argued that it should but I cannot recall the outcome).

anyway the right thing to do is to find if there are any points that are not converging and perturb only those points...

CC @KellyHu

If we're going to refine anyway, we should be able to hand in a full meshgrid to the refining code.

I tried running the heat_triangulated_sphere.py module, but this failed at some different points in the code. It appears the index shape used in build_diff_matrix and build_interp_matrix needs to be converted to an integer list.
I've just ran the module with shape = shape.astype(np.int64) and it runs with no errors now.