[petsc-users] about parallel preconditioned matrix-free gmres

[RenZhengYong]

Dear Petscs,

Use the uniprocessor complex-value based version petsc, I recently
successfully make a FETI_DP domain
decomposition approach working for 3D electromagnetic induction (earth)
problem.  The number of iteration of
the interface problem seems to be scalable with regard to the number of
sub-domains.

To do this, I had two subroutines for petsc

(1) int mat_vec_product_interface_problem(Mat A, Vec X, Vec Y) for
matrix-free GMRES solver
(2) int preconditioner_mat_vec(PC pc,Vec X,Vec Y) for shell preconditioner.

Now, I want to solve the interface problem by paralleled GMRES solver so
that I can solve real large-scale problems.  Could you please tell me the
easiest way to accomplish it. Which specific data structures of petsc
should be used.  I have been using Petsc for 3.5 years, I really want to
have a try the real MPI-based Petsc.


Thanks in advance.
Have a nice weekeed
Zhengyong


-- 
Zhengyong Ren
AUG Group, Institute of Geophysics
Department of Geosciences, ETH Zurich
NO H 47 Sonneggstrasse 5
CH-8092, Zürich, Switzerland
Tel: +41 44 633 37561
e-mail: zhengyong.ren at aug.ig.erdw.ethz.ch
Gmail: renzhengyong at gmail.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120630/319220d7/attachment.html>