Hello!
I just pushed some fixes to the jit interface of DOLFIN. Now one can jit on
different mpi groups. Previously jiting was only done on rank 1 of the
mpi_comm_world. Now it is done on rank 1 of any passed group communicator.
There is no demo atm showing this but a test has been added:
test/unit/python/jit/test_jit_with_mpi_groups.py
Here an expression, a subdomain, and a form is constructed on different
ranks using group. It is somewhat tedious as one need to initialize PETSc
with the same group, otherwise PETSc will deadlock during initialization
(the moment a PETSc la object is constructed).
The procedure in Python for this is:
1) Construct mpi groups using mpi4py
2) Initalize petscy4py using the groups
3) Wrap groups to petsc4py comm (dolfin only support petsc4py not mpi4py)
4) import dolfin
5) Do group specific stuff:
a) Function and forms no change needed as communicator
is passed via mesh
b) domain = CompiledSubDomain("...", mpi_comm=group_comm)
c) e = Expression("...", mpi_comm=group_comm)
Please try it out and report any sharp edges. A demo would also be fun to
include :)
Johan
_______________________________________________
fenics mailing list
[email protected]
http://fenicsproject.org/mailman/listinfo/fenics
