From 4df93b61768b47d3c4527a00f9f4196e32696527 Mon Sep 17 00:00:00 2001
From: Francisco Orlandini <francisco.orlandini@gmail.com>
Date: Mon, 27 May 2024 18:41:18 -0300
Subject: [PATCH] test sparse mat

---
 Matrix/TPZYSMPMatrix.cpp | 8 ++++++--
 1 file changed, 6 insertions(+), 2 deletions(-)

diff --git a/Matrix/TPZYSMPMatrix.cpp b/Matrix/TPZYSMPMatrix.cpp
index fc798145a..9d6a2096d 100644
--- 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: "