Multilevel solver

[Matthew Knepley]

On Wed, Apr 23, 2008 at 1:32 PM,  <Amit.Itagi at seagate.com> wrote:
> Barry,
>
>  Is the installation of petsc-dev different from the installation of the
>  2.3.3 release ? I ran the config. But the folder tree seems to be
>  different. Hence, make is giving problems.

1) Always always send the error log. I cannot tell anything from the
description "problems".

2) Some things have moved, but of course, make will work with the new
organization.

  Matt

>  Amit
>
>              Barry Smith
>              <bsmith at mcs.anl.g
>              ov>                                                        To
>              Sent by:                  petsc-users at mcs.anl.gov
>              owner-petsc-users                                          cc
>              @mcs.anl.gov
>              No Phone Info                                         Subject
>              Available                 Re: Multilevel solver
>
>
>              04/22/2008 10:08
>              PM
>
>
>              Please respond to
>              petsc-users at mcs.a
>                   nl.gov
>
>
>
>
>
>
>   Amit,
>
>      Using a a PCSHELL should be fine (it can be used with GMRES),
>  my guess is there is a memory corruption error somewhere that is
>  causing the crash. This could be tracked down with www.valgrind.com
>
>     Another way to you could implement this is with some very recent
>  additions I made to PCFIELDSPLIT that are in petsc-dev
>  (http://www-unix.mcs.anl.gov/petsc/petsc-as/developers/index.html)
>  With this you would chose
>  PCSetType(pc,PCFIELDSPLIT
>  PCFieldSplitSetIS(pc,is1
>  PCFieldSplitSetIS(pc,is2
>  PCFieldSplitSetType(pc,PC_COMPOSITE_SYMMETRIC_MULTIPLICATIVE
>  to use LU on A11 use the command line options
>  -fieldsplit_0_pc_type lu -fieldsplit_0_ksp_type preonly
>  and SOR on A22
>  -fieldsplit_1_pc_type sor -fieldsplit_1_ksp_type preonly -
>  fieldsplit_1_pc_sor_lits <lits> where
>     <its> is the number of iterations you want to use block A22
>
>  is1 is the IS that contains the indices for all the vector entries in
>  the 1 block while is2 is all indices in the
>  vector for the 2 block. You can use ISCreateGeneral() to create these.
>
>    Probably it is easiest just to try this out.
>
>    Barry
>
>
>  On Apr 22, 2008, at 8:45 PM, Amit.Itagi at seagate.com wrote:
>
>  >
>  > Hi,
>  >
>  > I am trying to implement a multilevel method for an EM problem. The
>  > reference is : "Comparison of hierarchical basis functions for
>  > efficient
>  > multilevel solvers", P. Ingelstrom, V. Hill and R. Dyczij-Edlinger,
>  > IET
>  > Sci. Meas. Technol. 2007, 1(1), pp 48-52.
>  >
>  > Here is the summary:
>  >
>  > The matrix equation Ax=b is solved using GMRES with a multilevel
>  > pre-conditioner. A has a block structure.
>  >
>  > A11    A12       *         x1  =  b1
>  > A21    A22                  x2       b2
>  >
>  > A11 is mxm and A33 is nxn, where m is not equal to n.
>  >
>  > Step 1  :      Solve  A11 *  e1   = b1     (parallel LU using
>  > superLU or
>  > MUMPS)
>  >
>  > Step 2:        Solve   A22 * e2    =b2-A21*e1    (might either user
>  > a SOR
>  > solver or a parallel LU)
>  >
>  > Step 3:        Solve   A11* e1 = b1-A12*e2   (parallel LU)
>  >
>  > This gives the approximate solution to
>  >
>  > A11     A12     *      e1   =  b1
>  > A21     A22             e2       b2
>  >
>  > and is used as the pre-conditioner for the GMRES.
>  >
>  >
>  > Which PetSc method can implement this pre-conditioner ? I tried a
>  > PCSHELL
>  > type PC. With Hong's help, I also got the parallel LU to work
>  > withSuperLU/MUMPS. My program runs successfully on multiple
>  > processes on a
>  > single machine. But when I submit the program over multiple
>  > machines, I get
>  > a crash in the PCApply routine after several GMRES iterations. I
>  > think this
>  > has to do with using PCSHELL with GMRES (which is not a good idea). Is
>  > there a different way to implement this ? Does this resemble the usage
>  > pattern of one of the AMG preconditioners ?
>  >
>  >
>  > Thanks
>  >
>  > Rgds,
>  > Amit
>  >
>
>
>
>



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which
their experiments lead.
-- Norbert Wiener