bad mesh generated.

[duncan]

This is from ticket #9, which was overloaded.

here's an error from least-squares.py

v 6
self.mesh.vertexlist[v] None
Traceback (most recent call last):
  File "I:\inundation\debug\mesh_error_reporting\show_mesh_bug.py", li
ne 123, in ?
    verbose=True)
  File "I:\inundation\pyvolution\least_squares.py", line 211, in fit_t
o_mesh
    precrop = precrop)
  File "I:\inundation\pyvolution\least_squares.py", line 360, in __ini
t__
    precrop = precrop)
  File "I:\inundation\pyvolution\least_squares.py", line 545, in build
_interpolation_matrix_A
    element_found, sigma0, sigma1, sigma2, k = \
  File "I:\inundation\pyvolution\least_squares.py", line 608, in searc
h_triangles_of_vertices
    for k, _ in self.mesh.vertexlist[v]:
TypeError: iteration over non-sequence

It's caused since vertex 6 is not part of any triangle.

I:\inundation\debug\mesh_error_reporting\show_bad_mesh.py (changeset:1954) produces a different mesh that shows the same error. Vertex 10 is not associated with any triangle. This is evident in the .tsh file.

See show_mesh_bug.py (changeset:1878) for 1st report of bug.

Currently least squares assumes a mesh is ok. Should this be changed? probably.

I just had a look at the mesh.py, at the triang.genMesh call. The error seems to be coming from triangle.

[ole]

Duncan said: "Currently least squares assumes a mesh is ok. Should this be changed? probably."

Perhaps that test should be part of check_integrity() in mesh ?

[duncan]

All user points are turned into triangulation vertices in triangle. added code in pmesh.mesh and triagle so vertices without a triangle are removed.