Fix LLT solve with compressed hmat RHS

jeromerobert · Dec 5, 2024 · ce0c73e · ce0c73e
1 parent 368b88d
commit ce0c73e
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/h_matrix.cpp b/src/h_matrix.cpp
@@ -2670,7 +2670,7 @@ template<typename T> void HMatrix<T>::solve(
         /* Solve LX=B, result in B */
         this->solveLowerTriangularLeft(b, algo, Diag::NONUNIT, Uplo::LOWER);
         /* Solve L^tX=B, result in B */
-        this->solveUpperTriangularLeft(b, algo, Diag::NONUNIT, Uplo::UPPER);
+        this->solveUpperTriangularLeft(b, algo, Diag::NONUNIT, Uplo::LOWER);
         break;
     default:
         HMAT_ASSERT(false);