[david_flo1]

Hello, I am trying to understand the implementation of ODESolver

When reading the line
yStart[0] = ::Foam::j0(xStart);
in ODETest.C, I was completely lost after some digging, all I found was

transFunc(j0)

I can't seem to be able to find its definition or body

The thing is I want to be able to understand all the mechanism so I can try and integrate ODE solving as input for a Foam solver

Any pointers are mostly welcomed, I tried google with no good results.

Regards

[olwi]

Hi,

If you're after the ODE features, I think your following the wrong lead if you track j0. I think j0 is just a Bessel function. (j0 and a lot of other standard math functions are defined in a generic way, using the transFunc macro you found.)

I will follow this thread, hoping someone helps you out with the ODE problem. It would be useful for me, too. If you solve your problem, it would be nice to see some feedback and recommendations here.

Good luck!

/Ola