Error in calculation of tensor eigenvalues
Using OpenFOAM 1706, with pimpleFoam, and the turbulentDFSEMInlet boundary condition, defined for velocity as below.
inlet
{
type turbulentDFSEMInlet;
value uniform (0 18 0);
delta 1;
R uniform (3.2 0.1 0.1 3.2 0.1 3.2);
L uniform 10;
U uniform (0 10 0);
mapMethod nearestCell;
}
Results in warning: --> FOAM Warning : From function Foam::vector Foam::eigenValues(const tensor&) in file primitives/Tensor/tensor/tensor.C at line 100 Complex eigenvalues detected for tensor: (3.2 0.1 0.1 0.1 3.2 0.1 0.1 0.1 3.2)
However, eigenvalues of this tensor are real (3.1, 3.1, 3.4).
The code in tensor.C appears correct, so I suspect the problem may lie in the cubicEqn code, which is significantly more opaque.