[massive_turbulence]

Hello All,

I'm trying to get this code to compile so that I can see what it does but for some reason it's hung up over this

solver.C:34:65: error: passing ‘const Internal {aka const Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh>}’ as ‘this’ argument discards qualifiers [-fpermissive]
     U.internalField() = axis ^ (centres.internalField() - origin);

Here's the full code

#include "fvCFD.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

int main(int argc, char *argv[])
{

	#include "setRootCase.H"
	#include "createTime.H"
	#include "createMesh.H"

    Info<< "Time = " << runTime.value() << endl;

    Info<< "    Reading U" << endl;
    volVectorField U
    (
        IOobject
        (
            "U",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ,
            IOobject::NO_WRITE
        ),
        mesh
    );

    // Do cells
    const volVectorField& centres = mesh.C();

    point origin(1, 1, 0.05);
    vector axis(0, 0, -1);

    U.internalField() = axis ^ (centres.internalField() - origin);

    U.boundaryField()[0] == (axis ^ (centres.boundaryField()[0] - origin));
    U.boundaryField()[1] == (axis ^ (centres.boundaryField()[1] - origin));
    U.boundaryField()[2] == (axis ^ (centres.boundaryField()[2] - origin));

    U.write();

    Info<< "End\n" << endl;

    return(0);
}

and the link to the original Github

Please let me know if there are other ways to do the same thing or something similar.

[Daniel_Khazaei]

you haven't said what OpenFOAM version you are trying to compile this piece of code on, but if the latest version is involved, you may try this:

U.primitiveFieldRef() = axis ^ (centres.primitiveField() - origin);

U.boundaryFieldRef()[0] == (axis ^ (centres.boundaryField()[0] - origin));
U.boundaryFieldRef()[1] == (axis ^ (centres.boundaryField()[1] - origin));
U.boundaryFieldRef()[2] == (axis ^ (centres.boundaryField()[2] - origin));

In recent versions, internalField() and boundaryField() functions both return a "const" reference to the internal field and boundary field of an object, respectively! So you can't modify them and need to use the above method.

[massive_turbulence]

Originally Posted by Daniel_Khazaei

you haven't said what OpenFOAM version you are trying to compile this piece of code on, but if the latest version is involved, you may try this:

U.primitiveFieldRef() = axis ^ (centres.primitiveField() - origin);

U.boundaryFieldRef()[0] == (axis ^ (centres.boundaryField()[0] - origin));
U.boundaryFieldRef()[1] == (axis ^ (centres.boundaryField()[1] - origin));
U.boundaryFieldRef()[2] == (axis ^ (centres.boundaryField()[2] - origin));

In recent versions, internalField() and boundaryField() functions both return a "const" reference to the internal field and boundary field of an object, respectively! So you can't modify them and need to use the above method.

OpenFOAM 6 v1806. Thank you sir!