[lindstroem]

Hi Foamers,

I was wondering about a defintion of a little function within a solver.
If I, for example, want to calculate the square of x in a function within a solver it would look like this:

        scalar xsquare
        (
                const scalar& x
        ) const
        {
                return pow(x,2);
        }

But my problem is, that I don't know where to place it. Using it in a BC I would put it in the header, but I cannot put it in createFields.H of my solver (in my case interFoam).

Any advises?
Thanks in advance!

[Bernhard]

Did you try putting them in front of "int main()"? If I remember correctly, functions are always defined outside the main program. createFields.H is included in the main program.

[lindstroem]

Good hint, but then I need to set the namespace correctly, right?

interFoam.C:57:4: warning: type qualifiers ignored on function return type [-Wignored-qualifiers]
interFoam.C:57:4: error: non-member function 'const scalar xsquare(const scalar&)' cannot have cv-qualifier

You know how to handle that?
thanks!

[Bernhard]

The only other hint I can give you, is to put it below other includes like fvCFD.H, but probably you tried that.

[lindstroem]

Yes, what I did was

#include "fvCFD.H"
#include "MULES.H"
#include "subCycle.H"
#include "interfaceProperties.H"
#include "twoPhaseMixture.H"
#include "turbulenceModel.H"

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

scalar xsquare
    (
        const scalar& x
    ) const
  {
      return pow(x,2);
     }

int main(int argc, char *argv[])
{
    #include "setRootCase.H"
    ...

But thanks anyway, I'll keep on digging for the namespace stuff, I am pretty sure thats the reason!

Greetings

[Bernhard]

In that case you can try Foam::scalar instead of just scalar, but I am absolutely not sure about that. Good luck!

[lindstroem]

Yeah, tried several of these combinations, but unsuccessfully. I got to understand that, otherwise it's just trial and error..

[Lieven]

Hey all,

This is how I would do it. Add

scalar xsquare(const scalar& x);

right underneath the #include statements and before int main (...

This way you can simply add your function at the bottom of the file so after the } of int main. But you should remove the 'const' after the (const scalar& x). So the function becomes:

scalar xsquare(const scalar& x)
 {
     return pow(x,2);
 }

I'm not a C++ expert, but this will make it work...

[lindstroem]

Thank you! Just works

Greetings