[saumavadey]

Hello Foamers!

I am using Modified Bessel Functions 'cyl_bessel_k' from the boost C++ library in my OpenFOAM program.

I need to get the Modified Bessel Function values for all the values of a scalarField. Do I need to run a loop for that? I seek suggestions in this regard.

[olesen]

Something like this?

scalarField  fld =... ;  // some field

for (scalar& val : fld)
{
    val = your_func(val);
}

[saumavadey]

The functional value of a field should be another field. But in the above example, the value of the function has been defined as a scalar only.

Is it correctly done?

[olesen]

Originally Posted by saumavadey

The functional value of a field should be another field. But in the above example, the value of the function has been defined as a scalar only.

Is it correctly done?

Wasn't clear from your question, but I don't see that the problem is OpenFOAM-specific, but looks much like any C++ solution. For example,

scalarField input =.... ;

scalarField output(input.size());
auto iter = output.begin();

for (const scalar val : input) 
{
    *iter = my_func(val);
    ++iter;
}

Can also use anything from the standard <algorithm> header too. I don't see where the problem is.

[dlahaye]

You can create the output scalarField by extending the file createField.H as AndreaCFD shows in his thread "Linearized Navier Stokes Solver". Once the solver has terminated, you should be able to see the file with the computed output scalar field. This file can visualized using ParaView.

Can you please provide more details on what exactly your intention is? What does it mean to apply the modified Bessel function to a scalarField in your case?