Matthew Knepley wrote: > On 8/27/07, Anders Logg <[EMAIL PROTECTED]> wrote: >> jesperc wrote: >>> Hi, >>> >> 4. The first thing we need is to be able to communicate Mesh and >> MeshFunction between processes. There is a sketch in >> src/kernel/mesh/MPIMeshCommunicator.cpp (empty). One process reads the >> mesh, partitions it and sends it to all processes. (To begin with, each >> process will have a copy of the mesh.) > > I am attaching a preprint which details our current scheme for this.
Excellent, but the attachment seems to have gone missing. /Anders >> 5. The tricky part will be how to order degrees of freedom (class >> DofMap). And this relates to ongoing work on merging the linear algebra >> interfaces of DOLFIN and PyCC (Simula in-house code). We need to find >> a design that works well for communicating both with PETSc and Epetra >> (Trilinos). > > You can see how we are doing it in PETSc by looking at > src/dm/mesh/sieve/Mesh.hh > in setupField(). However, this is still not completely correct since > it does not handle > vector elements correctly. I am fixing this now. > > Thanks, > > Matt > >> I'm sure Garth and Martin will have further opinions. >> >> /Anders >> >> _______________________________________________ >> DOLFIN-dev mailing list >> [email protected] >> http://www.fenics.org/mailman/listinfo/dolfin-dev >> > > _______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
