test sparse mat

labmec · May 27, 2024 · 4df93b6 · 4df93b6
1 parent cbcc1d0
commit 4df93b6
Showing 1 changed file with 6 additions and 2 deletions.
diff --git a/Matrix/TPZYSMPMatrix.cpp b/Matrix/TPZYSMPMatrix.cpp
@@ -213,12 +213,14 @@ void TPZFYsmpMatrix<TVar>::AddKel(TPZFMatrix<TVar>&elmat, TPZVec<int64_t> &sourc
     auto k = fIA[ipos];
     const auto maxj = fIA[ipos+1];
     const auto source_i = sourceindex[i];
+		int64_t *ja_ptr = &fJA[k];
     for(auto dummy_j=0;dummy_j<neq;dummy_j++){
       const auto j = idx[dummy_j];
       const auto jpos=destinationindex[j];
       const auto source_j = sourceindex[j];
       const auto &value=elmat.GetVal(source_i,source_j);
-      while(fJA[k]!=jpos){
+      while(*ja_ptr!=jpos){
+				ja_ptr++;
 				k++;
 				if(k==maxj){
 					std::cout << "TPZFYsmpMatrix::AddKelAtomic: "
@@ -257,12 +259,14 @@ void TPZFYsmpMatrix<TVar>::AddKelAtomic(TPZFMatrix<TVar>&elmat, TPZVec<int64_t>
     auto k = fIA[ipos];
     const auto maxj = fIA[ipos+1];
     const auto source_i = sourceindex[i];
+		int64_t *ja_ptr = &fJA[k];
     for(auto dummy_j=0;dummy_j<neq;dummy_j++){
       const auto j = idx[dummy_j];
       const auto jpos=destinationindex[j];
       const auto source_j = sourceindex[j];
       const auto &value=elmat.GetVal(source_i,source_j);
-      while(fJA[k]!=jpos){
+      while(*ja_ptr!=jpos){
+				ja_ptr++;
 				k++;
 				if(k==maxj){
 					std::cout << "TPZFYsmpMatrix::AddKelAtomic: "