[petsc-users] MatMumpsSetIcntl from Fortran
Alexander Grayver
agrayver at gfz-potsdam.de
Wed Nov 30 11:03:13 CST 2011
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
More information about the petsc-users
mailing list