Dear Petscs,<br><br>Use the uniprocessor complex-value based version petsc, I recently successfully make a FETI_DP domain <br>decomposition approach working for 3D electromagnetic induction (earth) problem.  The number of iteration of <br>
the interface problem seems to be scalable with regard to the number of sub-domains.  <br><br>To do this, I had two subroutines for petsc<br><br>(1) int mat_vec_product_interface_problem(Mat A, Vec X, Vec Y) for matrix-free GMRES solver<br>
(2) int preconditioner_mat_vec(PC pc,Vec X,Vec Y) for shell preconditioner. <br clear="all"><br>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. <br>
<br><br>Thanks in advance.<br>Have a nice weekeed<br>Zhengyong <br><br><br>-- <br>Zhengyong Ren<br>AUG Group, Institute of Geophysics<br>Department of Geosciences, ETH Zurich<br>NO H 47 Sonneggstrasse 5<br>CH-8092, Zürich, Switzerland<br>
Tel: +41 44 633 37561<br>e-mail: <a href="mailto:zhengyong.ren@aug.ig.erdw.ethz.ch">zhengyong.ren@aug.ig.erdw.ethz.ch</a><br>Gmail: <a href="mailto:renzhengyong@gmail.com">renzhengyong@gmail.com</a><br>