> 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. This would save my time :-) Appreciate.
Hong > > Regards, > Alexander > > > ----- Reply message ----- > From: "Hong Zhang" <hzhang at mcs.anl.gov> > 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 >> >>
