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>
