albiremo / dualGPy

Python library to generate the dual graph of a generic mesh
MIT License
2 stars 0 forks source link

Use arrays instead of flattened lists for points #17

Closed RiMillo closed 2 years ago

RiMillo commented 2 years ago

So far in ComputeGeometry points are flattened in one list, passed to Face2D (resp. 3D) where they are put back into matrix form of shape (n_vertices, 2)' (resp.3`). This avoids doing unnecessary works and keeps the matrix form.

This PR passes the test of the present repository. It passes 2D coupling tests with CoMMA (with and without anisotropy). However, unfortunately, when running this kind of tests in 3D I got a segmentation fault.