[petsc-users] MatMumpsSetIcntl from Fortran

[Matthew Knepley]

On Wed, Nov 30, 2011 at 2:37 PM, Hong Zhang <hzhang at mcs.anl.gov> wrote:

> Alexander :
>
> >
> > I just rewrote code concerning mumps from this example (lines 150-170):
> >
> http://www.mcs.anl.gov/petsc/petsc-dev/src/ksp/ksp/examples/tutorials/ex52.c.html
>
> Where is your Fortran code? ex52.c works fine.
>

Its in the version he attached at line 150-170.

   Matt


> Hong
>
> >
> > On 30.11.2011 17:40, Hong Zhang wrote:
> >>
> >> Alexander:
> >>
> >>> Has anybody tried to use MatMumpsSetIcntl from Fortran?
> >>
> >> We are not aware of it.
> >>
> >>> Because when I try to call it I fall into infinite recursion in
> function:
> >>
> >> Can you give me a short Fortran code that repeats this error for
> >> investigating?
> >> Meanwhile, you can use runtime option '-mat_mumps_icntl_xxx<>' to get
> >> your code run.
> >>
> >> Hong
> >>>
> >>>  PetscErrorCode MatMumpsSetIcntl(Mat F,PetscInt icntl,PetscInt ival)
> >>>  {
> >>>   PetscErrorCode ierr;
> >>>
> >>>   PetscFunctionBegin;
> >>>   PetscValidLogicalCollectiveInt(F,icntl,2);
> >>>   PetscValidLogicalCollectiveInt(F,ival,3);
> >>>   ierr =
> >>>
> >>>
> PetscTryMethod(F,"MatMumpsSetIcntl_C",(Mat,PetscInt,PetscInt),(F,icntl,ival));CHKERRQ(ierr);
> >>>   PetscFunctionReturn(0);
> >>>  }
> >>>
> >>> At the moment when program crashes call stack looks like:
> >>>
> >>>     __libc_memalign,
> >>> FP=7fff342ca010
> >>> PetscMallocAlign,
>  FP=7fff342ca080
> >>> PetscTrMallocDefault,
>  FP=7fff342ca180
> >>> PetscStrallocpy,
> FP=7fff342ca230
> >>> PetscFListGetPathAndFunction,
>  FP=7fff342cb2e0
> >>> PetscFListFind,
>  FP=7fff342cb520
> >>> PetscObjectQueryFunction_Petsc,
>  FP=7fff342cb590
> >>> PetscObjectQueryFunction,
>  FP=7fff342cb620
> >>> MatMumpsSetIcntl,
>  FP=7fff342cb720
> >>> MatMumpsSetIcntl,
>  FP=7fff342cb820
> >>> MatMumpsSetIcntl,
>  FP=7fff342cb920
> >>> MatMumpsSetIcntl,
>  FP=7fff342cba20
> >>> MatMumpsSetIcntl,
>  FP=7fff342cbb20
> >>> MatMumpsSetIcntl,
>  FP=7fff342cbc20
> >>> MatMumpsSetIcntl,
>  FP=7fff342cbd20
> >>> MatMumpsSetIcntl,
>  FP=7fff342cbe20
> >>> MatMumpsSetIcntl,
>  FP=7fff342cbf20
> >>> ... (Hundreds of MatMumpsSetIcntl) ...
> >>>
> >>> What can I do about that?
> >>>
> >>> Regards,
> >>> Alexander
> >
> >
>



-- 
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/20111130/19bdf094/attachment.htm>