[Users] REAL_PRECISION 16

Roland Haas rhaas at illinois.edu
Tue Jan 25 10:53:04 CST 2022


Hello David,

> It still looks like the Vectors Thorn is required to Carpet to work.
> Is there a way to not compile this Thorn and still have Carpet work?
Hmm, I have not tried this out actually. 

VECTORISE=no

is an option to the vector's thorn which tells it to not use any simd
vectorization but just plain scalars. So the thorn is present, making
CarpetLib happy. But Vectors (see below) may not compile.

When using VECTORISE=no it should be possible to compile even
with CCTK_REAL being extended precision (16 bytes), though it is
possible that there are explicit checks for the CCTK_REAL width even in
that case in the code.

The good news is that, in that case only, one should be able to change
the source code to make things compile by (worst case) add a specific
sizeof(CCTK_REAL) == 16 case and using "long double" instead of
"double" for the type of variable.

Looking at the source code of Vectors one probably has to duplicate
vectors-8-default.h into vectors-16-default.h and add a branch for
CCTK_REAL_PRECISION_16 to vectors.h.

Does this explanation make sense to you and gives you an idea how you
could change the code?

Yours,
Roland

-- 
My email is as private as my paper mail. I therefore support encrypting
and signing email messages. Get my PGP key from http://pgp.mit.edu .
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 833 bytes
Desc: OpenPGP digital signature
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20220125/841960e2/attachment.bin 


More information about the Users mailing list