<div dir="ltr"><font size="3" face="Calibri">Thank you Matthew and Hong.<br>
</font>
<br><font size="3" face="Calibri">That's too bad...</font>
<br><font size="3" face="Calibri">I was wondering if it would be possible
to add this alternative only in the MUMPS interface. <br>
Indeed, in the field of solid mechanics (where 3D, beams and shell are
often mixed thanks to Lagrange multipliers), it's very usefull to factorize
a matrix in</font> <font size="3" face="Calibri">single precision thanks to MUMPS, and use
it to precondition a sequence of systems solved in double precision. It
turns to be very robust and quite fast.<br>
</font>
<br><font size="3" face="Calibri">In  </font><a href="http://www.mcs.anl.gov/petsc/petsc-current/src/mat/impls/aij/mpi/mumps/mumps.c.html#MATSOLVERMUMPS" target="_blank"><font size="3" color="blue" face="Calibri"><u>src/mat/impls/aij/mpi/mumps/mumps.c,</u></font></a><font size="3" color="blue" face="Calibri">
</font><font size="3" face="Calibri">a new type, say PetscMumpsSingleScalar,
could be defined in order to pass the required information to the
single precision</font>
<br><font size="3" face="Calibri">MUMPS solver.</font>
<br><font size="3" face="Calibri"><br>
Would the Petsc team accept to integrate such a feature?</font>
<br>
<br><font size="3" face="Calibri">Thanks,</font>
<br><font size="3" face="Calibri">Sylvain</font></div><div class="gmail_extra"><br><br><div class="gmail_quote">2014-08-20 17:57 GMT+02:00 Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span>:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div><div class="h5">On Wed, Aug 20, 2014 at 7:53 AM, Sylvain Mercier <span dir="ltr"><<a href="mailto:sylvainmercier85@gmail.com" target="_blank">sylvainmercier85@gmail.com</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Hi everyone,<br><br></div>I'm using PETSc to solve saddle point systems. I have found a block preconditionner which requires to solve a smaller symmetric positive definite system. I do it with  MUMPS which is fast an robust.<br>


</div><div>However, I would like to use this solver in single precision in order to compare the performances.<br></div><div><br></div><div>So is it possible to use MUMPS in single precision inside a double precision solving with PETSc?<br>

</div></div></blockquote><div><br></div></div></div><div>No, right now PETSc must have all the same precision. We have investigated several methods for</div><div>using multi-precision, but none has satisfied all our needs.</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"><div dir="ltr"><div></div><div>Thanks,<br>

</div><div>Sylvain<span class="HOEnZb"><font color="#888888"><br></font></span></div></div><span class="HOEnZb"><font color="#888888">
</font></span></blockquote></div><span class="HOEnZb"><font color="#888888"><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
</font></span></div></div>
</blockquote></div><br></div>