Hi Jamie, There is no automated way to make a FiPy mesh from Scipy's Delaunay object. The problem is that FiPy needs a face to vertex array and a cell to face array while the Delaunay object just has the cell to vertex array. The functionality to extract the face to vertex array and cell to face array is in FiPy because it must be done for Gmsh, however, it is not abstracted in so that it can be reused.
It is certainly possible to make the correct arrays with Numpy and pass them to the Mesh2D object, but it's a bit of work to write the code. If I find some time I might give it a go, but I don't know when I will get to it. Cheers, Daniel On Tue, Jun 14, 2016 at 4:15 PM, James Pringle <jprin...@unh.edu> wrote: > Daniel et al. -- > > As referenced earlier in this thread, I have a complex domain with holes > in it; I now have it broken up into triangles, based on the Delaunay package > in SciPy. I have > > Locations of all vertex coordinates, > list of vertices that make up faces > list of faces that make cells > list of faces that make up all the internal and external boundaries. > > Can you point me to any code or documentation that would help me understand > how to make this into a Mesh2D object? I am having a devil of a time > figuring out from the manual online. The best would be something that used > the output of either the triangles or scipy.spatial.Delaunay() packaged. My > equation is of the form > > 0=J(Psi,A(x,y)) + \Del(B(x,y)*\Del Psi) > > > and I can get the coefficients A(x,y) and B(x,y) on either the faces or in > the cell centers are needed. > > Thank you. > Jamie Pringle > University of New Hampshire -- Daniel Wheeler _______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
