"Munson, Todd" <tmun...@mcs.anl.gov> writes: > Any idea how to do the same thing in C?
#define PetscRealLiteral(x) (x ## q) PetscReal a = PetscRealLiteral(0.7); Seems like we should do that for the precomputed methods in TS. Actually, it would be better to do all the setup arithmetic in the most accurate supported type. Unfortunately, those methods are often computed using tools that only support double precision or are taken From papers that only report that many digits.
signature.asc
Description: PGP signature