From 29b3e9adfec20221ea40c1b619310e182c5a39e2 Mon Sep 17 00:00:00 2001 From: Henry Date: Sun, 3 Nov 2013 23:26:32 +0000 Subject: [PATCH] Improved docs --- src/ODE/ODESolvers/Rosenbrock43/Rosenbrock43.C | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/src/ODE/ODESolvers/Rosenbrock43/Rosenbrock43.C b/src/ODE/ODESolvers/Rosenbrock43/Rosenbrock43.C index 73df3c4da6..df1680aa1e 100644 --- a/src/ODE/ODESolvers/Rosenbrock43/Rosenbrock43.C +++ b/src/ODE/ODESolvers/Rosenbrock43/Rosenbrock43.C @@ -113,7 +113,7 @@ Foam::scalar Foam::Rosenbrock43::solve LUDecompose(a_, pivotIndices_); - // compute k1: + // Compute k1: forAll(k1_, i) { k1_[i] = dydx0[i] + dx*d1*dfdx_[i]; @@ -121,7 +121,7 @@ Foam::scalar Foam::Rosenbrock43::solve LUBacksubstitute(a_, pivotIndices_, k1_); - // compute k2: + // Compute k2: forAll(y, i) { y[i] = y0[i] + a21*k1_[i]; @@ -136,7 +136,7 @@ Foam::scalar Foam::Rosenbrock43::solve LUBacksubstitute(a_, pivotIndices_, k2_); - // compute k3: + // Compute k3: forAll(y, i) { y[i] = y0[i] + a31*k1_[i] + a32*k2_[i]; @@ -151,6 +151,7 @@ Foam::scalar Foam::Rosenbrock43::solve LUBacksubstitute(a_, pivotIndices_, k3_); + // Compute k4: forAll(k4_, i) { k4_[i] = dydx_[i] + dx*d4*dfdx_[i]