[petsc-users] petsc-3.5 quadruple precision mixed with double

[Matthew Knepley]

On Wed, Jun 28, 2017 at 8:22 PM, Hao Zhang <hbcbh1999 at gmail.com> wrote:

> It's 3d incompressible RT simulation. My pressure between serial and
> parallel calculation is off by 10^(-14) in relative error.
>

This could just be reordering of the calculation.


> it eventually build up at later time. I want to rule out the possibilities
> that PETSc give me bad solution. pressure scale is 10^(-2).
>
> I use PetscScalar. thanks @Jed Brown for confirming that but I have
> Segmentation Violation when retrieving x. I allocated memory for the array
> x (PetscScalar type). if not for quadruple precision, there is no error.
>

It sounds like maybe you are passing a double where a PetscScalar is
expected, or vice versa. Run under valgrind.

  Matt


> thanks @Satish Balay. I will update code petsc-3.7 later.
>



-- 
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

http://www.caam.rice.edu/~mk51/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20170628/361d83db/attachment-0001.html>