zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

Improved docs

Browse Source
This commit is contained in:
Henry
2013-11-03 23:26:32 +00:00
parent 734d88a875
commit 29b3e9adfe
1 changed files with 4 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
src/ODE/ODESolvers/Rosenbrock43/Rosenbrock43.C
View File

@ -113,7 +113,7 @@ Foam::scalar Foam::Rosenbrock43::solve
LUDecompose(a_, pivotIndices_); LUDecompose(a_, pivotIndices_);
// compute k1: // Compute k1:
forAll(k1_, i) forAll(k1_, i)
{ {
k1_[i] = dydx0[i] + dx*d1*dfdx_[i]; k1_[i] = dydx0[i] + dx*d1*dfdx_[i];
@ -121,7 +121,7 @@ Foam::scalar Foam::Rosenbrock43::solve
LUBacksubstitute(a_, pivotIndices_, k1_); LUBacksubstitute(a_, pivotIndices_, k1_);
// compute k2: // Compute k2:
forAll(y, i) forAll(y, i)
{ {
y[i] = y0[i] + a21*k1_[i]; y[i] = y0[i] + a21*k1_[i];
@ -136,7 +136,7 @@ Foam::scalar Foam::Rosenbrock43::solve
LUBacksubstitute(a_, pivotIndices_, k2_); LUBacksubstitute(a_, pivotIndices_, k2_);
// compute k3: // Compute k3:
forAll(y, i) forAll(y, i)
{ {
y[i] = y0[i] + a31*k1_[i] + a32*k2_[i]; y[i] = y0[i] + a31*k1_[i] + a32*k2_[i];
@ -151,6 +151,7 @@ Foam::scalar Foam::Rosenbrock43::solve
LUBacksubstitute(a_, pivotIndices_, k3_); LUBacksubstitute(a_, pivotIndices_, k3_);
// Compute k4:
forAll(k4_, i) forAll(k4_, i)
{ {
k4_[i] = dydx_[i] + dx*d4*dfdx_[i] k4_[i] = dydx_[i] + dx*d4*dfdx_[i]