Skip to content
Prev 37984 / 63424 Next

long double, C, fortran

I'm using .Call() to call C code from R under Windows (on an Intel 
Core 2 duo). The C code involves some very small numbers, and I think 
I'm losing precision using doubles. I thought I might use long doubles 
to see if I can get that precision back. I have a few questions:

1. Does this affect the portability to other OSs or processors?
2. I'm returning the results in a matrix. Will a matrix of REALs be 
sufficient for holding long doubles, or will it be cast back to doubles?
3. Will calls to FORTRAN BLAS (like dsymv, dpotrf, dpotri) still work 
with long doubles?

Thanks for any help you can provide.