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

SVD: Improved overflow/division protection.

Browse Source 
Resolves bug report https://bugs.openfoam.org/view.php?id=2486
This commit is contained in:
Will Bainbridge
2017-04-20 10:54:08 +01:00
parent ea257737df
commit e4ac3fba77
1 changed files with 3 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
src/OpenFOAM/matrices/scalarMatrices/SVD/SVD.C
View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -156,7 +156,7 @@ Foam::SVD::SVD(const scalarRectangularMatrix& A, const scalar minCondition)
{ {
if (i < Un-1) if (i < Un-1)
{ {
if (g != 0) if (g*U_(i, l) != 0)
{ {
for (label j=l; j<Un; j++) for (label j=l; j<Un; j++)
{ {
@ -199,7 +199,7 @@ Foam::SVD::SVD(const scalarRectangularMatrix& A, const scalar minCondition)
U_(i, j) = 0; U_(i, j) = 0;
} }
if (g != 0) if (g*U_(i, i) != 0)
{ {
g = 1.0/g; g = 1.0/g;