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>
