<div dir="ltr"><div><div><div><div><div>hi, all:<br><br></div>I'm developing a CFD project, which all the matrix solver is PETSc based. It works fine with double precision until I need some more, like quadruple precision.<br><br></div><div>Ax = b is the system<br></div>input are double precision<br></div><div><br>shall I use PetscScalar instead of double?<br></div>Or is there a better way to pass double type to PETSc solver? <br></div><div>Or shall I convert double type to quadruple type before passing into the solver?<br><br><br></div>I'm using PETSc 3.5.4. GCC 5.2.0 comiple with float128. I think there should be no compiling error.<br></div>I can provide more code part if necessary. Thanks!<br clear="all"><div><div><div><div><div><div><div><div><div><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">Hao <br></div>
</div></div></div></div></div></div></div></div></div></div>