Missing third index in the write to file script

In the commented out code for writing the solution to a file, the third index is missing. Just adding it here.
MattBolitho · Nov 22, 2022 · 5de3177 · 5de3177
1 parent a666652
commit 5de3177
Show file tree

Hide file tree

Showing 3 changed files with 9 additions and 3 deletions.
diff --git a/examples/ScalableProblems/MittelmannBndryCntrlDiri3D.cpp b/examples/ScalableProblems/MittelmannBndryCntrlDiri3D.cpp
@@ -906,7 +906,9 @@ void MittelmannBndryCntrlDiriBase3D::finalize_solution(
 
     for (Index i=0; i<=N_+1; i++) {
     for (Index j=0; j<=N_+1; j++) {
-    fprintf(fp, "y[%6d,%6d] = %15.8e\n", i, j, x[y_index(i,j)]);
+    for (Index k=0; k<=N_+1; k++) {
+    fprintf(fp, "y[%6d,%6d,%6d] = %15.8e\n", i, j, k, x[y_index(i,j,k)]);
+    }
     }
     }
 

diff --git a/examples/ScalableProblems/MittelmannBndryCntrlDiri3D_27.cpp b/examples/ScalableProblems/MittelmannBndryCntrlDiri3D_27.cpp
@@ -1074,7 +1074,9 @@ void MittelmannBndryCntrlDiriBase3D_27::finalize_solution(
 
     for (Index i=0; i<=N_+1; i++) {
     for (Index j=0; j<=N_+1; j++) {
-    fprintf(fp, "y[%6d,%6d] = %15.8e\n", i, j, x[y_index(i,j)]);
+    for (Index k=0; k<=N_+1; k++) {
+    fprintf(fp, "y[%6d,%6d,%6d] = %15.8e\n", i, j, k, x[y_index(i,j,k)]);
+    }
     }
     }
 

diff --git a/examples/ScalableProblems/MittelmannBndryCntrlDiri3Dsin.cpp b/examples/ScalableProblems/MittelmannBndryCntrlDiri3Dsin.cpp
@@ -1001,7 +1001,9 @@ void MittelmannBndryCntrlDiriBase3Dsin::finalize_solution(
 
     for (Index i=0; i<=N_+1; i++) {
     for (Index j=0; j<=N_+1; j++) {
-    fprintf(fp, "y[%6d,%6d] = %15.8e\n", i, j, x[y_index(i,j)]);
+    for (Index k=0; k<=N_+1; k++) {
+    fprintf(fp, "y[%6d,%6d,%6d] = %15.8e\n", i, j, k, x[y_index(i,j,k)]);
+    }
     }
     }