[doctorWho]

Hi,

I want to do some calculation with point coordinate with pow function.

The code is roughly as follows:

const pointMesh& pMesh = pointMesh::New(mesh);
        vectorField coord = this->points0();
        scalarField sf(this->points0().size());

        forAll(sf, pointi)
        {
             sf[pointi] = pow(coord[pointi].x(), 0.6);
        }

I have searched several threads about pow function of OpenFOAM. I have tried Foam:: pow instead pow without any success.

when the exponent is integer type e.g., 2 instead of 0.6, then it works.
The problem is always when the exponent is scalar type or float type.
Any advice?

[wyldckat]

Quick answer: My guess is that you have "x" values that are negative. For example:

pow(a,0.5)

is the same as:

sqrt(a)

so if "a" is negative, the answer is a complex number, not a real number, hence the crash.

[doctorWho]

Yes, you are right. Thanks alot.