[adhiraj]

For a chemically reacting flow if we use ODE solver as RK for instance, does each thread get its own copy of the RK object?
The code for RK4 seems thread unsafe otherwise--am I missing something?

$FOAM_SRC/ODE/ODESolvers/RK/RK.C

[wyldckat]

Greetings Adhiraj,

Originally Posted by adhiraj

For a chemically reacting flow if we use ODE solver as RK for instance, does each thread get its own copy of the RK object?
The code for RK4 seems thread unsafe otherwise--am I missing something?

OpenFOAM (the official one, not the variant Extend) doesn't use threads. It relies on MPI doing all of the communication magic

The Extend variant on the other hand, does have at least some modifications for running with threads instead of MPI, but I don't know if it's already working as intended or not.

Best regards,
Bruno