The TS constants are weird in the code. Some of have enough digits for double, some for quad, and some for quad double. I suspect mathematica output some of the constants.
Todd. > On Sep 1, 2016, at 10:29 PM, Jed Brown <j...@jedbrown.org> wrote: > > "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.