[petsc-users] Problems about parallel version in SNES_MF

[Yingjie Wu]

Dear Petsc developers:
Hi,
I am currently using Petsc to solve a two-dimensional system of non-linear
PDEs, which is a thermal problem, including pressure field, temperature
field and velocity field. Since it is not easy to construct Jacobian matrix
explicitly, I adopt snes_mf method and provide a preconditioning matrix.
Now I want to adapt the program to parallel version.
The main steps of the current procedure are as fellows:

   - A solution vector U is constructed. The dimension is the sum of the
   meshed of each physical field.
   - FormFunction is provided to calculate residual.
   - Preconditioning matrix Pmat is assembled using FormJacibian, but Amat
   (Jacobian matrix) is not assembled.
   - Compute with -snes_mf_operator.

There are following questions about the development of parallel version:

   1. If one processor is assigned to each physical field, such as pressure
   field P,velocity field V, temperature field T, belongs to three processors.
   Because the problem is nonlinear, the information of temperature T may be
   used in the calculation of pressure P(in FormFunction). How to transfer the
   information?
   2. Since the snes_mf method needs a good preconditioning matrix to
   ensure convergence. How to use preconditioning matrix in parallel? As far
   as I know, BJACOBI can run in parallel. How can it be used in snes_mf
   method?

Because there are few snes_mf examples with parallel computation, I need
some advice from you.

Thanks,
Yingjie
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20181225/b85d4f40/attachment.html>