[petsc-users] Matrix format mpiaij does not have a built-in PETSc XXX!

[Bojan Niceno]

Dear all,


Max may be correct, but I encounter the same problem as Aron.  Neither 
PCILU nor PCICC work in parallel for me.  Here is the message I get:

[0]PETSC ERROR: --------------------- Error Message 
------------------------------------
[0]PETSC ERROR: No support for this operation for this object type!
[0]PETSC ERROR: Matrix format mpiaij does not have a built-in PETSc ICC!
[0]PETSC ERROR: 
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.2.0, Patch 6, Wed Jan 11 
09:28:45 CST 2012
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.

Portion of the code which I use to set the solver contents follows (I am 
aiming at CG+ICC combination):

/* Create KPS content */
   KSPCreate(PETSC_COMM_WORLD, &ksp);
   KSPSetType(ksp ,KSPCG);

   /* Set operators */
   KSPSetOperators(ksp, A, A, DIFFERENT_NONZERO_PATTERN);

   /* Linear solver defaults (can be ove-ridden) */
   KSPGetPC(ksp, &pc);
   PCSetType(pc, PCICC);
   KSPSetTolerances(ksp, 1.e-5, PETSC_DEFAULT, PETSC_DEFAULT, 
PETSC_DEFAULT);

   /* Run-time options (over-rides above) */
   KSPSetFromOptions(ksp);

What is going wrong here?


     Kind regards,


     Bojan


On 2/26/2012 6:17 PM, Matthew Knepley wrote:
> On Sun, Feb 26, 2012 at 10:48 AM, Max Rudolph <maxwellr at gmail.com 
> <mailto:maxwellr at gmail.com>> wrote:
>
>     MPIAIJ and SEQQIJ matrices are subtypes of the AIJ matrix type.
>     Looking at that table, you should be able to use any of the PCs
>     that supports AIJ and has an X under 'parallel'. 
>
>
> Max is correct. For instance, the most popular general purpose 
> parallel solver is ASM (Additive Schwarz Method), which then
> has a sequential subsolver for each block, which defaults to ILU.
>
>    Matt
>
>     Max
>
>
>     On Sun, Feb 26, 2012 at 8:16 AM, Aron Roland <aaronroland at gmx.de
>     <mailto:aaronroland at gmx.de>> wrote:
>
>         Dear All,
>
>         I hope somebody can help us on this or give at least some
>         clearance.
>
>         We have just included PETSc as an solver for our sparse matrix
>         evolving from an unstructured mesh advection scheme.
>
>         The problem is that we are using the mpiaij matrix type, since
>         our matrix is naturally sparse. However it seems that PETSc
>         has no PC for this, except the PCSOR, which showed to be not
>         very effective for our problem.
>
>         All others give the error msg. of the mail subject, where XXX
>         are the different PC tried.
>
>         The manual is a bit diffuse on this e.g.
>
>         http://www.mcs.anl.gov/petsc/documentation/linearsolvertable.html
>
>         it is claimed that certain PC's are running on aij matrices
>         ... but these are to be defined either as seq. or parallel
>         (mpiaij) matrices. Moreover in the above mentioned list are
>         two columns parallel/seriel, what is the intention of parallel
>         capability when not applicable to matrices stored within the
>         parallel mpiaij framework.
>
>         I guess we just not understanding the concept or have some
>         other difficulties of understanding of all this.
>
>         Any comments help is welcome
>
>         Aron
>
>
>
>
>
> -- 
> 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


-- 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120228/77f1ccd3/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Signature.png
Type: image/png
Size: 6515 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120228/77f1ccd3/attachment.png>