From 72e07ed568aff5690506e26b0fa49f58697d997a Mon Sep 17 00:00:00 2001 From: sergio Date: Wed, 2 Apr 2014 10:03:13 +0100 Subject: [PATCH 1/2] BUG: Correct index for wideBandAbsorptionEmission for coeff array --- .../wideBandAbsorptionEmission/wideBandAbsorptionEmission.C | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/thermophysicalModels/radiationModels/submodels/absorptionEmissionModel/wideBandAbsorptionEmission/wideBandAbsorptionEmission.C b/src/thermophysicalModels/radiationModels/submodels/absorptionEmissionModel/wideBandAbsorptionEmission/wideBandAbsorptionEmission.C index a1bf1cc0..b2294fe9 100644 --- a/src/thermophysicalModels/radiationModels/submodels/absorptionEmissionModel/wideBandAbsorptionEmission/wideBandAbsorptionEmission.C +++ b/src/thermophysicalModels/radiationModels/submodels/absorptionEmissionModel/wideBandAbsorptionEmission/wideBandAbsorptionEmission.C @@ -102,7 +102,7 @@ Foam::radiation::wideBandAbsorptionEmission::wideBandAbsorptionEmission << nl << exit(FatalError); } } - coeffs_[nSpec][nBand].initialise(specDicts.subDict(key)); + coeffs_[nBand][nSpec].initialise(specDicts.subDict(key)); nSpec++; } nBand++; From 072e93f46359e8e67db0c0d525a8c9b928062024 Mon Sep 17 00:00:00 2001 From: mattijs Date: Thu, 3 Apr 2014 09:11:21 +0100 Subject: [PATCH 2/2] ENH: tensor: limit calculation of small eigenvalues --- src/OpenFOAM/primitives/Tensor/tensor/tensor.C | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/OpenFOAM/primitives/Tensor/tensor/tensor.C b/src/OpenFOAM/primitives/Tensor/tensor/tensor.C index 8ddf9d29..d0d3b30b 100644 --- a/src/OpenFOAM/primitives/Tensor/tensor/tensor.C +++ b/src/OpenFOAM/primitives/Tensor/tensor/tensor.C @@ -119,7 +119,7 @@ Foam::vector Foam::eigenValues(const tensor& t) + t.xy()*t.yx()*t.zz() + t.xx()*t.yz()*t.zy(); // If there is a zero root - if (mag(c) < ROOTVSMALL) + if (mag(c) < 1e-100) { scalar disc = sqr(a) - 4*b; @@ -313,7 +313,7 @@ Foam::vector Foam::eigenValues(const symmTensor& t) + t.xy()*t.xy()*t.zz() + t.xx()*t.yz()*t.yz(); // If there is a zero root - if (mag(c) < ROOTVSMALL) + if (mag(c) < 1e-100) { scalar disc = sqr(a) - 4*b;