[estr_sprt]

Hi Foamers,

Please help to eliminate the following error in OF v6:

no match for ‘operator*’ (operand types are ‘Foam::tmp<Foam::GeometricField<Foam::Tensor<doubl e>, Foam::fvPatchField, Foam::volMesh> >’ and ‘Foam::volTensorField {aka Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh>}’)
("rhotau",(mt*dev(twoSymm(tgrad_U))-(scalar(2)/scalar(3))*Ei*mt*tgrad_U_-(scalar(2)/scalar(3))*Ei*rho*(turbulence->k()+dummykk)));

which arises during the compilation. Here Ei is an identity matrix which I define before as

Ei = tensor:ne;

after adding the following to the createfields.H

volTensorField Ei
(
IOobject
(
"Ei",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
mesh,
I
);

Which was suggested in the thread How to access identity matrix in OF?

I also looked through some related topics but can't still deal with it.
Also I tried define the Ei through a constant reference to "I" (which is an edintity tensor) defining is as a volTensor Field, but see the same.

I am only a new user to OF and c++ itself and therefore might be subject to beginner's errors.

"dummykk" is just a scalar, "mt" is a volScalarField and "tgrad_U_" is

const volTensorField& tgradU_ = fvc::grad(U);
volTensorField tgrad_U_ = tgradU_;

Thank you