[petsc-users] Fill-in of factorization

Matthew Knepley knepley at gmail.com
Fri Feb 28 17:23:10 CST 2014


On Fri, Feb 28, 2014 at 2:49 PM, Jan Blechta <blechta at karlin.mff.cuni.cz>wrote:

> On Fri, 28 Feb 2014 14:05:35 -0800
> Matthew Knepley <knepley at gmail.com> wrote:
>
> > On Fri, Feb 28, 2014 at 1:59 PM, Jan Blechta
> > <blechta at karlin.mff.cuni.cz>wrote:
> >
> > > On Fri, 28 Feb 2014 12:30:05 -0800
> > > Matthew Knepley <knepley at gmail.com> wrote:
> > >
> > > > On Fri, Feb 28, 2014 at 12:13 PM, Jan Blechta
> > > > <blechta at karlin.mff.cuni.cz>wrote:
> > > >
> > > > > Hi,
> > > > >
> > > > > is it possible to check an amount of fill-in produced during
> > > > > factorization with PCLU or PCCHOLESKY? Or should I ask the
> > > > > library, say MUMPS, directly?
> > > >
> > > >
> > > > This is from the PCFactorView():
> > > >
> > > >       MatInfo info;
> > > >       ierr =
> > > > MatGetInfo(factor->fact,MAT_LOCAL,&info);CHKERRQ(ierr);
> > >
> > > MatGetInfo does not seem to be available in petsc4py. Am I grepping
> > > correctly?
> > >
> >
> > It appears to be missing, probably because it passes a struct around.
>
> Is it a bug in petsc4py?
>

Its the lack of a feature.

   Matt


> Jan
>
> >
> >   Matt
> >
> >
> > > Jan
> > >
> > > >       ierr = PetscViewerASCIIPrintf(viewer,"  factor fill ratio
> > > > given %g, needed
> > > >
> > >
> %g\n",(double)info.fill_ratio_given,(double)info.fill_ratio_needed);CHKERRQ(ierr);
> > > >
> > > >   THanks,
> > > >
> > > >      MAtt
> > > >
> > > >
> > > > >
> > > > > Jan
> > > > >
> > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>


-- 
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/20140228/c6b4769b/attachment-0001.html>


More information about the petsc-users mailing list