[petsc-users] MatMumpsSetIcntl from Fortran

[agrayver@gfz-potsdam.de]

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" <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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111130/12e9139d/attachment-0001.htm>