[petsc-users] Using PETSc solvers and preconditioners with mfem

[Matthew Knepley]

On Tue, Oct 25, 2016 at 10:54 AM, Abdullah Ali Sivas <
abdullahasivas at gmail.com> wrote:

> Hello,
>
> I want to use PETSc with mfem and I know that mfem people will figure out
> a way to do it in few months. But for now as a temporary solution I just
> thought of converting hypre PARCSR matrices (that is what mfem uses as
> linear solver package) into PETSc MPIAIJ matrices and I have a semi-working
> code with some bugs. Also my code is dauntingly slow and seems like not
> scaling. I have used *MatHYPRE_IJMatrixCopy *from myhp.c of PETSc and
> *hypre_ParCSRMatrixPrintIJ* from par_csr_matrix.c of hypre as starting
> points. Before starting I checked whether there was anything done similar
> to this, I could not find anything.
>
> My question is, are you aware of such a conversion code (i.e. something
> like *hypre_ParCSRtoPETScMPIAIJ( hypre_ParCSRMatrix *matrix, **Mat *A*)?
>
No, but maybe Satish knows. Slow running times most likely come from lack
of preallocation for the target matrix.

  Thanks,

     Matt

> Thanks in advance,
> Abdullah Ali Sivas
>



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20161025/1c59f094/attachment.html>