Initialize LU to zero, probably not needed.

modified: Solvers_Mod/Dense/Lu_Factorization_Gauss.f90
Niceno · Feb 18, 2024 · a0b38e2 · a0b38e2
1 parent 520705d
commit a0b38e2
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/Solvers_Mod/Dense/Lu_Factorization_Gauss.f90 b/Solvers_Mod/Dense/Lu_Factorization_Gauss.f90
@@ -39,10 +39,10 @@ subroutine Solvers_Mod_Dense_Lu_Factorization_Gauss(LU, A)
   LU % text_u ="L"
   LU % text_l ="U"
 
-
   !----------------------------------------------------------------------!
   !   Initialize the values by copying the original matrix to LU first   !
   !----------------------------------------------------------------------!
+  LU % val(:,:) = 0.0
   do i = 1, n  ! <-A
     do j = 1, n
       LU % val(i,j) = A % val(i,j)