On Wed, Nov 30, 2011 at 2:37 PM, Hong Zhang <span dir="ltr"><<a href="mailto:hzhang@mcs.anl.gov">hzhang@mcs.anl.gov</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Alexander :<br>
<br>
><br>
> I just rewrote code concerning mumps from this example (lines 150-170):<br>
> <a href="http://www.mcs.anl.gov/petsc/petsc-dev/src/ksp/ksp/examples/tutorials/ex52.c.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-dev/src/ksp/ksp/examples/tutorials/ex52.c.html</a><br>
<br>
Where is your Fortran code? ex52.c works fine.<br></blockquote><div><br></div><div>Its in the version he attached at line 150-170.</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hong<br>
<br>
><br>
> On 30.11.2011 17:40, Hong Zhang wrote:<br>
>><br>
>> Alexander:<br>
>><br>
>>> Has anybody tried to use MatMumpsSetIcntl from Fortran?<br>
>><br>
>> We are not aware of it.<br>
>><br>
>>> Because when I try to call it I fall into infinite recursion in function:<br>
>><br>
>> Can you give me a short Fortran code that repeats this error for<br>
>> investigating?<br>
>> Meanwhile, you can use runtime option '-mat_mumps_icntl_xxx<>' to get<br>
>> your code run.<br>
>><br>
>> Hong<br>
>>><br>
>>> PetscErrorCode MatMumpsSetIcntl(Mat F,PetscInt icntl,PetscInt ival)<br>
>>> {<br>
>>> PetscErrorCode ierr;<br>
>>><br>
>>> PetscFunctionBegin;<br>
>>> PetscValidLogicalCollectiveInt(F,icntl,2);<br>
>>> PetscValidLogicalCollectiveInt(F,ival,3);<br>
>>> ierr =<br>
>>><br>
>>> PetscTryMethod(F,"MatMumpsSetIcntl_C",(Mat,PetscInt,PetscInt),(F,icntl,ival));CHKERRQ(ierr);<br>
>>> PetscFunctionReturn(0);<br>
>>> }<br>
>>><br>
>>> At the moment when program crashes call stack looks like:<br>
>>><br>
>>> __libc_memalign,<br>
>>> FP=7fff342ca010<br>
>>> PetscMallocAlign, FP=7fff342ca080<br>
>>> PetscTrMallocDefault, FP=7fff342ca180<br>
>>> PetscStrallocpy, FP=7fff342ca230<br>
>>> PetscFListGetPathAndFunction, FP=7fff342cb2e0<br>
>>> PetscFListFind, FP=7fff342cb520<br>
>>> PetscObjectQueryFunction_Petsc, FP=7fff342cb590<br>
>>> PetscObjectQueryFunction, FP=7fff342cb620<br>
>>> MatMumpsSetIcntl, FP=7fff342cb720<br>
>>> MatMumpsSetIcntl, FP=7fff342cb820<br>
>>> MatMumpsSetIcntl, FP=7fff342cb920<br>
>>> MatMumpsSetIcntl, FP=7fff342cba20<br>
>>> MatMumpsSetIcntl, FP=7fff342cbb20<br>
>>> MatMumpsSetIcntl, FP=7fff342cbc20<br>
>>> MatMumpsSetIcntl, FP=7fff342cbd20<br>
>>> MatMumpsSetIcntl, FP=7fff342cbe20<br>
>>> MatMumpsSetIcntl, FP=7fff342cbf20<br>
>>> ... (Hundreds of MatMumpsSetIcntl) ...<br>
>>><br>
>>> What can I do about that?<br>
>>><br>
>>> Regards,<br>
>>> Alexander<br>
><br>
><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>