My cluster has an old version of glibc 2.3.3, which can not give usual principal value of square root. For example, it gives square root of -i as the value with negative real part. I 'grep' in the source tree of petsc and found csqrt in many places. Is this dangerous? Thank you!