[petsc-users] MatMumpsSetIcntl from Fortran
Hong Zhang
hzhang at mcs.anl.gov
Wed Nov 30 14:37:50 CST 2011
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
>
>
More information about the petsc-users
mailing list