Hong, Sorry if I wasn't clear. "c" example works fine, I know, what I meant is that if you try to implement lines 150-171 from it on FORTRAN you will see the problem I reported. If you need particularly my FORTRAN code I can send it tomorrow.
Regards, Alexander ----- Reply message ----- From: "Hong Zhang" <[email protected]> To: "PETSc users list" <petsc-users at mcs.anl.gov> Subject: [petsc-users] MatMumpsSetIcntl from Fortran Date: Wed, Nov 30, 2011 9:37 pm 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. 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 > > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111130/12e9139d/attachment-0001.htm>
