[petsc-users] csqrt

Satish Balay balay at mcs.anl.gov
Thu Feb 28 14:04:44 CST 2013


On Thu, 28 Feb 2013, Barry Smith wrote:

> 
> On Feb 28, 2013, at 5:33 AM, Hui Zhang <mike.hui.zhang at hotmail.com> wrote:
> 
> > 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!
> 
>    If you are using the default PETSc builds where it works with real numbers then it will have no effect on.
> 
>    If you are using PETSc with complex numbers you might consider changing the cluster software.

presumably csqrt() is c99 - so c++ std:sqrt() is perhaps ok ok? [so one could use --with-clanguage=cxx]

Satish


More information about the petsc-users mailing list