Hao Zhang <hbcbh1...@gmail.com> writes: > It's 3d incompressible RT simulation. My pressure between serial and > parallel calculation is off by 10^(-14) in relative error. it eventually > build up at later time.
I doubt that is why your solution develops artifacts after many time steps. If you're using __float128, you'll take a big hit on performance so you might as well make everything __float128 (PetscScalar). Mixing precisions is a maintenance burden and exactly the wrong thing when you're trying to rule out potential sources of error.
signature.asc
Description: PGP signature
