[petsc-users] Calling single-precision MUMPS from PETSC

Evan Um evanum at gmail.com
Fri Jun 5 18:26:36 CDT 2015

Dear Barry and PETSC users,

I am revisiting a problem about how to call a single precision MUMPS from
double precision real/complex PETSC. After taking a look at mumps.c, I feel
that the following changes (attached) can make it possible to always call a
single precision MUMPS from PETSC. The change is basically to replace
double precision MUMPS and their associated types with corresponding single
ones. If someone already has some experience about this work, could you
comment on the change? In advance, I appreciate your help.


On Wed, Oct 22, 2014 at 1:36 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:

>    There is no support for this. You can only call single precision MUMPS
> from single precision PETSc and double precision MUMPS from double
> precision PETSc.
>     You could try to hack the interface that calls MUMPS from PETSc to use
> the single precision version but we don’t support that.
> src/mat/impls/aij/mpi/mumps/mumps.c
>   Barry
> > On Oct 22, 2014, at 3:29 PM, Evan Um <evanum at gmail.com> wrote:
> >
> > Dear PETSC users,
> >
> > When MUMPS is used inside PETSC, The default MUMPS driver seems to be
> double-precision MUMPS (i.e. DMUMPS). To reduce memory costs, I want to
> test a single-precision MUMPS (SMUMPS) from PETSC. Does anyone know how to
> switch from double to single-precision MUMPS inside PETSC? In advance,
> thanks for your comments.
> >
> > Regards,
> > Evan
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150605/d7f7e09e/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: mumps_modified.c
Type: text/x-csrc
Size: 87301 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150605/d7f7e09e/attachment-0001.c>

More information about the petsc-users mailing list