From 9fbe2ef0b34eac9ccf77fb8631d1376888961efc Mon Sep 17 00:00:00 2001 From: Kutalmis Bercin Date: Tue, 6 Dec 2022 15:30:36 +0000 Subject: [PATCH] ENH: leastSquaresFaVectors: use pseudo inverse instead of direct inverse --- .../gradSchemes/leastSquaresFaGrad/leastSquaresFaVectors.C | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/finiteArea/finiteArea/gradSchemes/leastSquaresFaGrad/leastSquaresFaVectors.C b/src/finiteArea/finiteArea/gradSchemes/leastSquaresFaGrad/leastSquaresFaVectors.C index 730aeaa699..2370dcac31 100644 --- a/src/finiteArea/finiteArea/gradSchemes/leastSquaresFaGrad/leastSquaresFaVectors.C +++ b/src/finiteArea/finiteArea/gradSchemes/leastSquaresFaGrad/leastSquaresFaVectors.C @@ -6,7 +6,7 @@ \\/ M anipulation | ------------------------------------------------------------------------------- Copyright (C) 2016-2017 Wikki Ltd - Copyright (C) 2020 OpenCFD Ltd. + Copyright (C) 2020-2022 OpenCFD Ltd. ------------------------------------------------------------------------------- License This file is part of OpenFOAM. @@ -159,7 +159,7 @@ void Foam::leastSquaresFaVectors::makeLeastSquaresVectors() const // Invert the dd tensor - const symmTensorField invDd(inv(dd)); + const symmTensorField invDd(pinv(dd)); // Revisit all faces and calculate the lsP and lsN vectors