Improved printing of the residual

labmec · Jul 3, 2024 · 749b548 · 749b548
1 parent 6d900c6
commit 749b548
Showing 1 changed file with 8 additions and 1 deletion.
diff --git a/Analysis/TPZAnalysis.cpp b/Analysis/TPZAnalysis.cpp
@@ -1394,6 +1394,7 @@ void TPZAnalysis::PrintVectorByElement(std::ostream &out, TPZFMatrix<STATE> &vec
         if (hasgeometry) {
             out << " Gel " << gel->Index() << " matid " << gel->MaterialId() << " Center " << xcenter << std::endl;
         }
+        TPZManVector<REAL,3> xco(3);
         for (ic = 0; ic<nc; ic++) {
             TPZManVector<STATE> connectsol;
             int64_t cindex = cel->ConnectIndex(ic);
@@ -1408,7 +1409,13 @@ void TPZAnalysis::PrintVectorByElement(std::ostream &out, TPZFMatrix<STATE> &vec
                     connectsol[i] = 0.;
                 }
             }
-            out << ic << " index " << cindex << " values " << connectsol << std::endl;
+            if(gel && ic < gel->NCornerNodes()) {
+                gel->NodePtr(ic)->GetCoordinates(xco);
+                out << "co " << xco << " ic ";
+            }
+            if(connectsol.size()) {
+                out << ic << " index " << cindex << " values " << connectsol << std::endl;
+            }
         }
     }