songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

ENH: lowerCSR: ignore lowerCSR for matrix type.

Browse Source 
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 ...
This commit is contained in:
mattijs 2025-01-08 11:53:16 +00:00
parent b86ba60d35
commit 0a99702f9e
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
src/OpenFOAM/matrices/lduMatrix/lduMatrix/lduMatrix.H
View File

@ -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_);
}