Hi, Matthew and Hong, I will use PETSc for my 4 years project. :) Thanks Zhengyong Ren
On Fri, 2009-03-13 at 08:46 -0500, Hong Zhang wrote: > Zhengyong, > > > I known your excellent codes several months ago. I want to use Mumps as > > direct solver for multi-sources problem and BoomerAMG for real-value based > > iterative solver. I went though the doc of the petsc. It showed that petsc > > offered a easy and top interface to these two packages. So, I think if I > > chose the petsc as my tools, the coding work should be easier for me as only > > data structure of petsc should be leaned, not one for Mumps and one for > > Hypre. > > > > Am I right? can pestc do the work Mumps and BloomAMG do? I want to get your > > Yes, you can use Mumps, BloomAMG and other packages > (see http://www.mcs.anl.gov/petsc/petsc-as/miscellaneous/external.html) > without changing your application code. > > > confirm on my decision. As you said, sometimes it is important to make a > > correct decision and also to learn petsc should not be a short way. If petsc > > really do the fast direct solving and fast algebra multigrid algorithms by > > Mumps and BloomAMG, respectively, I think PESTC should definitely be the > > first choice for my following PhD project. > > Here is how to get started: > > 1. install petsc with configure options > '--download-scalapack --download-superlu --download-superlu_dist > --download-mumps --download-blacs --download-hypre' > > (you need F90 compiler for mumps. Mumps requires > scalapack and blacs. I also suggest you install > superlu and superlu_dist - commonly used sparse direct sovlers) > > 2. build petsc and test the installation > > 3. test mumps and hypre, e.g., > cd ~petsc/src/ksp/ksp/examples/tutorials > make ex2 > ./ex2 -ksp_monitor -ksp_type preonly -pc_type lu > -pc_factor_mat_solver_package mumps > > Norm of error < 1.e-12 iterations 1 > > ./ex2 -ksp_monitor -pc_type hypre -pc_hypre_type boomeramg > > 0 KSP Residual norm 7.372609060992e+00 > 1 KSP Residual norm 5.824128070306e-02 > 2 KSP Residual norm 8.370364383637e-05 > Norm of error 8.39541e-05 iterations 2 > > Let us know if you encounter difficulty, > > Hong > > > > > > -- > > Zhengyong Ren > > AUG Group, Institute of Geophysics > > Department of Geoscience > > NO H 47 Sonneggstrasse 5 > > CH-8092, Zrich, Switzerland > > Tel: +41 44 633 37561 > > e-mail: renzh at ethz.ch > > Gmail: renzhengyong at gmail.com > >