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

COMP: EigenMatrix: remove unused iter variable

Browse Source
This commit is contained in:
Kutalmis Bercin
2022-11-24 13:47:08 +00:00
parent 9d212dfd32
commit 8c02820db6
2 changed files with 2 additions and 20 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
applications/test/matrices/EigenMatrix/Test-EigenMatrix.C
View File

@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2020 OpenCFD Ltd.
Copyright (C) 2020-2022 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is derivative work of OpenFOAM.
@ -49,20 +49,6 @@ using namespace Foam;
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// Return the absolute tolerance value for bitwise comparisons of floatScalars
floatScalar getTol(floatScalar)
{
return 1e-2;
}
// Return the absolute tolerance value for bitwise comparisons of doubleScalars
doubleScalar getTol(doubleScalar)
{
return 1e-10;
}
// Create each constructor of EigenMatrix<Type>, and print output
template<class Type>
void test_constructors(Type)

6
src/OpenFOAM/matrices/EigenMatrix/EigenMatrix.C
View File

@ -9,7 +9,7 @@
Contributed 2018-07-31 to the OpenFOAM Foundation
Copyright (C) 2018 OpenFOAM Foundation
Copyright (C) 2019-2020 Alberto Passalacqua
Copyright (C) 2020 OpenCFD Ltd.
Copyright (C) 2020-2022 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -221,12 +221,8 @@ void Foam::EigenMatrix<cmptType>::symmTridiagonalQL()
// If m == l, EValsRe_[l] is an eigenvalue, otherwise, iterate
if (m > l)
{
label iter = 0;
do
{
iter += 1;
// Compute implicit shift
cmptType g = EValsRe_[l];
cmptType p = (EValsRe_[l + 1] - g)/(2.0*EValsIm_[l]);