Update linear algebra function templates.cpp

SNEWS2 · Jul 2, 2024 · 97b09ce · 97b09ce
1 parent 350e30d
commit 97b09ce
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/BEMEWS/_ext/mstl/math2/algebra/linear algebra function templates.cpp b/src/BEMEWS/_ext/mstl/math2/algebra/linear algebra function templates.cpp
@@ -2076,7 +2076,7 @@ template <std::size_t N> CVECTOR<double,N> TridiagonalSolve(MVECTOR<double,N> &L
      { CVECTOR<double,N> X, Y0(Y);
 
        try{ // effectively zero all the subdiagonal terms
-            int i,imax=static_cast<int>(D.N())-2;
+            int imax=static_cast<int>(D.N())-2;
             for(int i=0;i<=imax;i++){ D[i+1]-=L[i]/D[i]*U[i]; Y[i+1]-=L[i]/D[i]*Y[i]; }
 
             // back substitution