From 0a99702f9e2ca166a35f687d3d0a6292bf07fbc0 Mon Sep 17 00:00:00 2001 From: mattijs Date: Wed, 8 Jan 2025 11:53:16 +0000 Subject: [PATCH] ENH: lowerCSR: ignore lowerCSR for matrix type. Reverting 267bd1af36 since makes it hard to run cell-based algorithm on symmetric matrix. - might become a problem if all discretisation becomes lowerCSR based and lower() is never constructed ... --- src/OpenFOAM/matrices/lduMatrix/lduMatrix/lduMatrix.H | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/OpenFOAM/matrices/lduMatrix/lduMatrix/lduMatrix.H b/src/OpenFOAM/matrices/lduMatrix/lduMatrix/lduMatrix.H index ee6b8d5bd2..40179089fd 100644 --- a/src/OpenFOAM/matrices/lduMatrix/lduMatrix/lduMatrix.H +++ b/src/OpenFOAM/matrices/lduMatrix/lduMatrix/lduMatrix.H @@ -708,19 +708,19 @@ public: //- Matrix has diagonal only bool diagonal() const noexcept { - return (diagPtr_ && !lowerPtr_ && !lowerCSRPtr_ && !upperPtr_); + return (diagPtr_ && !lowerPtr_ && !upperPtr_); } //- Matrix is symmetric bool symmetric() const noexcept { - return (diagPtr_ && !(lowerPtr_ || lowerCSRPtr_) && upperPtr_); + return (diagPtr_ && !lowerPtr_ && upperPtr_); } //- Matrix is asymmetric (ie, full) bool asymmetric() const noexcept { - return (diagPtr_ && (lowerPtr_ || lowerCSRPtr_) && upperPtr_); + return (diagPtr_ && lowerPtr_ && upperPtr_); }