Hi Hong, Thank you for quick reply.
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 Regards, Alexander 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
