Can I ask you a question about forAll()? Thank you in advance!
In openfoam, there are lots of forAll used, for example (it is from hPsiMixture.C):
const scalarField& hCells = h_.internalField();
const scalarField& pCells = p_.internalField();
scalarField& TCells = T_.internalField();
scalarField& psiCells = psi_.internalField();
scalarField& muCells = mu_.internalField();
scalarField& alphaCells = alpha_.internalField();
const typename MixtureType::thermoType& mixture =
TCells[celli] = mixture.TH(hCells[celli], TCells[celli]);
psiCells[celli] = mixture.psi(pCells[celli], TCells[celli]);
muCells[celli] = mixture.mu(TCells[celli]);
alphaCells[celli] = mixture.alpha(TCells[celli]);
Here, why is the Tcells in forAll() used? I mean why not use other quantities like muCells or alphaCells? I read the source code src/OpenFOAM/containers/Lists/UList/UList.H, but still not got a clear understanding about it.
best regards, h
From what I can see, you can choose either one of the arrays as a reference for the cell ID, as long as you know that all of them have the same IDs. The choice for "TCells" is probably because:
For example, to give another point of view on a similar way of going through an array: in C#, it's possible to go over a complete array with a source code like this:
C# foreach is basically the same thing as OpenFOAM's forAll. But OpenFOAM's is simpler and easier to use! :cool:
Besides all of this, the other advantage is that forAll is generic enough to work in parallel, therefore it won't go through all of the cells of the global mesh in all of the processors, it only goes through the list of cells present in the current processor (at least as far as I can understand).
Thank you so much for the continuous help! About the last point you mentioned, I indeed have some difficulty to understand. For example. I write a loop forAll() in fireFoam. C, which is the top level code we can see. In this case, forAll will be run in a parallel way for a parallel computations?
Thank you !
If I remember this correctly:
forAll is simply:
Here's your example using the new c++11 standard (no need to switch language!):
|All times are GMT -4. The time now is 18:12.|