<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Oct 25, 2016 at 12:19 PM, Stefano Zampini <span dir="ltr"><<a href="mailto:stefano.zampini@gmail.com" target="_blank">stefano.zampini@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I have a working conversion from HypreParCSR to PETSc MPIAIJ format.<br>
I could add this code to PETSc, maybe in the contrib folder. Barry, what do you think?<br></blockquote><div><br></div><div>No, no one looks there. Add it to src/mat/utils and make an interface function like MatCreateFromHypreParCSR().</div><div><br></div><div>  Thanks,</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">
> We tried a similar approach to get MFEM objects to PETSc and the real<br>
> problem is that all other convenient functions like creating<br>
> gridfunctions and projections, you have to convert them every time<br>
> which is basically nothing more than a bad workaround.<br>
<br>
So far, my interface covers matrices and Krylov solvers (PCFieldSplit and PCBDDC are explicitly supported).<br>
<br>
Can you tell me how would you like to use these objects with PETSc? What you would like to achieve?<br>
<br>
So far, my work on the PETSc interface to MFEM originated from a wishlist for solvers, but I could expand it.<br>
<br>
<br>
> Stefano Zampini<br>
> replied to the issue on github and was stating that there is some<br>
> intention to get MFEM working with PETSc but there is no specific<br>
> timeframe.<br>
><br>
<br>
There’s currently an open pull request for PETSc solvers inside the private MFEM repo.<br>
However, I don’t know when the code, if merged, will be released.<br>
<br>
<br>
> Regards<br>
> Julian Andrej<br>
><br>
> On Tue, Oct 25, 2016 at 6:30 PM, Abdullah Ali Sivas<br>
> <<a href="mailto:abdullahasivas@gmail.com">abdullahasivas@gmail.com</a>> wrote:<br>
>> I will check that. I am preallocating but it may be that I am not allocating<br>
>> big enough. I still have to figure out nuance differences between these<br>
>> formats to solve the bugs. I appreciate your answer and hope Satish knows<br>
>> it.<br>
>><br>
>> Thank you,<br>
>> Abdullah Ali Sivas<br>
>><br>
>><br>
>> On 2016-10-25 12:15 PM, Matthew Knepley wrote:<br>
>><br>
>> On Tue, Oct 25, 2016 at 10:54 AM, Abdullah Ali Sivas<br>
>> <<a href="mailto:abdullahasivas@gmail.com">abdullahasivas@gmail.com</a>> wrote:<br>
>>><br>
>>> Hello,<br>
>>><br>
>>> I want to use PETSc with mfem and I know that mfem people will figure out<br>
>>> a way to do it in few months. But for now as a temporary solution I just<br>
>>> thought of converting hypre PARCSR matrices (that is what mfem uses as<br>
>>> linear solver package) into PETSc MPIAIJ matrices and I have a semi-working<br>
>>> code with some bugs. Also my code is dauntingly slow and seems like not<br>
>>> scaling. I have used MatHYPRE_IJMatrixCopy from myhp.c of PETSc and<br>
>>> hypre_ParCSRMatrixPrintIJ from par_csr_matrix.c of hypre as starting points.<br>
>>> Before starting I checked whether there was anything done similar to this, I<br>
>>> could not find anything.<br>
>>><br>
>>> My question is, are you aware of such a conversion code (i.e. something<br>
>>> like hypre_ParCSRtoPETScMPIAIJ( hypre_ParCSRMatrix *matrix, Mat *A)?<br>
>><br>
>> No, but maybe Satish knows. Slow running times most likely come from lack of<br>
>> preallocation for the target matrix.<br>
>><br>
>>  Thanks,<br>
>><br>
>>     Matt<br>
>>><br>
>>> Thanks in advance,<br>
>>> Abdullah Ali Sivas<br>
>><br>
>><br>
>><br>
<span class="HOEnZb"><font color="#888888">>><br>
>> --<br>
>> What most experimenters take for granted before they begin their experiments<br>
>> is infinitely more interesting than any results to which their experiments<br>
>> lead.<br>
>> -- Norbert Wiener<br>
>><br>
>><br>
<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">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</div>
</div></div>