[vcvedant]

Hello,
I am modifying a turbulence model and for a part of the code I get an error:

error: no match for 'operator[]' in 'FST[cellI]'

for the function

tmp<volScalarField> myModel::rCm
(
     const volScalarField& tke
)
{
     tmp<volScalarField> rC = tke/tke;
     tmp<volScalarField> FST = tke;
     const fvPatchScalarField* topTKE(tke.boundaryField()[2]);
     forAll(tke, cellI)
     {
          FST[cellI] = topTKE[3];
     }
     return rC;
}

I get an additional

error: 'class Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >' has no member named 'size'

when I tried

     forAll(FST,cellI)
     {
         FST[cellI] = topTKE[3];
     }

What should I change for FST such that I can access individual value of FST for each cell in the mesh and assign it some value?

Thanks,

[Zeppo]

FST()[cellI] = topTKE[3];

[DoQuocVu]

Hi Vedamt,

As far as i know, the forAll(FST, cellI) means: for(int i=0; i< A; i++), with A is a integer which is the size of FST.

However, you defined FST as a volScalarField which is a 2D or 3D matrix according to your mesh, and size of a matrix is mxn but not a scalar so you cannot use it for the "for" loop.

The only solution I can come up with is using another scalarField:

scalarField& FSTin_ = FST.internalField();
forAll(FSTin_, cellI)
     FSTin_[cellI] = topTKE[3];

Good luck,
Vu

[Zeppo]

scalarField& FSTin_ = FST.internalField();

->

scalarField& FSTin_ = FST().internalField();

[Sugajen]

Hi all,

I did a similar expression as below but got an error.
error:

error: invalid initialization of reference of type ‘Foam::scalarField& {aka Foam::Field<double>&}’ from expression of type ‘const Internal {aka const Foam::DimensionedField<double, Foam::volMesh>}’
     scalarField& varNuIn_ = varNu_().internalField();

code:

    Foam::tmp<Foam::volScalarField> varNu_;
   
    scalarField& varNuIn_ = varNu_().internalField();

Any ideas?

[Zeppo]

Originally Posted by Sugajen

scalarField& varNuIn_ = varNu_().internalField();

If you are not going to modify the field later on:

const scalarField& varNuIn_ = varNu_();

otherwise:

scalarField& varNuIn_ = varNu_.ref();

[sibo]

Hi Sergel,

I have a similar error message when compiling a customized solver in OF_4.1.

It is a dynamic solver and in order to assign my calculated displacement value "dispVals" to cellDisplacement boundary, I use this line:

PointDisplacement.boundaryField()[patchID] == dispVals;

and the error is:

error: no match for ‘operator==’ (operand types are ‘const Foam:ointPatchField<Foam::Vector<double> >’ and ‘Foam::vectorField {aka Foam::Field<Foam::Vector<double> >}’)
PointDisplacement.boundaryField()[patchID] == dispVals;
^

here I want to use this "==" to force assignment.
Do you know how to solve this problem? Any suggestion is welcome!
Thanks a lot!

[Zeppo]

PointDisplacement.boundaryField()[patchID] == dispVals;

->

PointDisplacement.boundaryFieldRef()[patchID] == dispVals;

[sibo]

Thanks a lot! I solved this problem a few days ago, so I think in the new version they name the const and non-const function differently.

[NICKTAVOUKT]

Hi guys I also have a similar problem... more specifically:

When I try to compile the following:

inline tmp<volVectorField> solidParticleCloud::momentumSource() const
{
tmp<volVectorField> tsource
(
new volVectorField
(
IOobject
(
"smom",
mesh_.time().timeName(),
mesh_,
IOobject::NO_READ,
IOobject::NO_WRITE
),
mesh_,
dimensionedVector
(
"zero",
dimensionSet(1, -2, -2, 0, 0),
vector::zero
)
)
);
tsource().internalField() = smom_/mesh_.cellVolumes();
return tsource;
}

and I get the following error:

solidParticleCloudI.H:91:57: error: no match for 'operator=' (operand types are 'const Internal {aka const Foam:imensionedField<Foam::Vector<double>, Foam::volMesh>}' and 'Foam::tmp<Foam::Field<Foam::Vector<double> > >')
tsource().internalField() = smom_/mesh_.cellVolumes();

Interestingly enough when I compiled the same code in OF3 I did not have any problems... The problem occurs only when compiled n OF5
I would really appreciate any help!

Thanks!

[sibo]

HI Nick,

Try

tsource().internalFieldRef() = smom_/mesh_.cellVolumes();

I believe in the version 4.0 and after that, they seperate the const and non-const function with different names. InternalField() provides the const reference to the internal field, internalFieldRef() provides a non-const reference.

[NICKTAVOUKT]

Thanks Sibo,

I'll give it a try and see what comes out of it!

[NICKTAVOUKT]

Hi all,

I found this link which describes how to update older code versions to the newest... which also gives an answer to most of the questions posed in this thread. I hope you find it helpful!

https://www.openfoam.com/documentati...rade-guide.php

Regards,
Nick

[JulianS]

Hey everyone,

This thread already helped me a lot. Thanks for this guys!
Now I would like to do something quite similar, just the other way round as the issue solved before here:

Originally Posted by Zeppo

PointDisplacement.boundaryField()[patchID] == dispVals;

->

PointDisplacement.boundaryFieldRef()[patchID] == dispVals;

When I try to write the PointDisplacement entries to a vectorField, like this:

dispVals=PointDisplacement.boundaryFieldRef()[patchID];

I get the error message:

error: no match for ‘operator=’ (operand types are ‘Foam::vectorField {aka Foam::Field<Foam::Vector<double> >}’ and ‘Foam::pointPatchField<Foam::Vector<double> >’)

Does anyone know what I'm doing wrong or how to solve this?

Thanks!

[sujata]

hello everyone.
I am using multiphaseEulerfoam solver. I am new to OpenFOAM and using version 4.1 . I am trying to print u velocity, Reynolds No in the phasepair.C file
Foam::tmp<Foam::volScalarField> Foam:hasePair::magUr() const
{
return mag(phase1().U() - phase2().U());
Info<< "magUr =" <<magUr()[0] << endl;
}


Foam::tmp<Foam::volVectorField> Foam:hasePair::Ur() const
{
return dispersed().U() - continuous().U();
}


Foam::tmp<Foam::volScalarField> Foam:hasePair::Re() const
{
return magUr()*dispersed().d()/continuous().nu();
Info<< "Re =" <<Re()[0] << endl;
}



The error I am getting is as follows

phasePair/phasePair/phasePair.C(108): error: no operator "[]" matches these operands
operand types are: Foam::tmp<Foam::GeometricField<Foam::scalar, Foam::fvPatchField, Foam::volMesh>> [ int ]
Info<< "magUr =" <<magUr()[0] << endl;
^

phasePair/phasePair/phasePair.C(121): error: no operator "[]" matches these operands
operand types are: Foam::tmp<Foam::GeometricField<Foam::scalar, Foam::fvPatchField, Foam::volMesh>> [ int ]
Info<< "Re =" <<Re()[0] << endl;
^

I want to print it for 0th cell. Can anyone help me with this.

[DoQuocVu]

Hi Sujata,

the function mag() return a scalar, thus magUr() is a scalar and it doesn't have [] operator

[Zeppo]

Originally Posted by sujata

hello everyone.
I am using multiphaseEulerfoam solver. I am new to OpenFOAM and using version 4.1 . I am trying to print u velocity, Reynolds No in the phasepair.C file
Foam::tmp<Foam::volScalarField> Foam:hasePair::magUr() const
{
return mag(phase1().U() - phase2().U());
Info<< "magUr =" <<magUr()[0] << endl;
}


Foam::tmp<Foam::volVectorField> Foam:hasePair::Ur() const
{
return dispersed().U() - continuous().U();
}


Foam::tmp<Foam::volScalarField> Foam:hasePair::Re() const
{
return magUr()*dispersed().d()/continuous().nu();
Info<< "Re =" <<Re()[0] << endl;
}



The error I am getting is as follows

phasePair/phasePair/phasePair.C(108): error: no operator "[]" matches these operands
operand types are: Foam::tmp<Foam::GeometricField<Foam::scalar, Foam::fvPatchField, Foam::volMesh>> [ int ]
Info<< "magUr =" <<magUr()[0] << endl;
^

phasePair/phasePair/phasePair.C(121): error: no operator "[]" matches these operands
operand types are: Foam::tmp<Foam::GeometricField<Foam::scalar, Foam::fvPatchField, Foam::volMesh>> [ int ]
Info<< "Re =" <<Re()[0] << endl;
^

I want to print it for 0th cell. Can anyone help me with this.

Foam::tmp<Foam::volScalarField> Foam:hasePair::magUr() const
{
    tmp<volScalarField> ret = mag(phase1().U() - phase2().U());
    Info<< "magUr =" << ret()[0] << endl;
    return ret;
}

Function mag here returns tmp<volScalarField>. You can access its first element through [0];

[Zeppo]

Originally Posted by JulianS

When I try to write the PointDisplacement entries to a vectorField, like this:

dispVals=PointDisplacement.boundaryFieldRef()[patchID];

I get the error message:

error: no match for ‘operator=’ (operand types are ‘Foam::vectorField {aka Foam::Field<Foam::Vector<double> >}’ and ‘Foam::pointPatchField<Foam::Vector<double> >’)

Does anyone know what I'm doing wrong or how to solve this?

Thanks!

dispVals=PointDisplacement.boundaryField()[patchID].primitiveField();

This way you get a copy of internal patch field. When you modify any value of dispVals no changes will be made to PointDisplacement.boundaryField()[patchID]!

[wind_]

Originally Posted by NICKTAVOUKT

Hi guys I also have a similar problem... more specifically:

When I try to compile the following:

and I get the following error:



Interestingly enough when I compiled the same code in OF3 I did not have any problems... The problem occurs only when compiled n OF5
I would really appreciate any help!

Thanks!

Hi，Nick
I also have this problem in OF5x, and I changed
tsource().internalField() =...
->
tsource().primitiveField() =...
but it doesn't work, how to fix it?
Thanks
Feng

[wind_]

Originally Posted by wind_

Hi，Nick
I also have this problem in OF5x, and I changed
tsource().internalField() =...
->
tsource().primitiveField() =...
but it doesn't work, how to fix it?
Thanks
Feng

This thread is the detail.Error while compiling a solidParticle library based on interFoam in OF5x