Changeset 27931

Timestamp:

10/01/23 18:46:54 (18 months ago)

Author:

Mathieu Morlighem

Message:

CHG: added print functions for inverse methods

Location:

issm/trunk-jpl/src/c/shared/Elements

Files:

• PrintArrays.cpp (modified) (1 diff)
• elements.h (modified) (1 diff)

issm/trunk-jpl/src/c/shared/Elements/PrintArrays.cpp

@@ -84 +84 @@
         }
 }
+
+void InversionStatsHeader(int NJ){
+        _printf0_("\n");
+        _printf0_("┌────┬─────────────────┬────────────┬────────────────────────\n");
+        _printf0_("│Iter│  Cost function  │ Grad. norm │  List of contributions \n");
+        _printf0_("├────┼─────────────────┼────────────┼────────────────────────\n");
+}
+void InversionStatsIter(int iter,double J, double Gnorm, double* Jlist, int N){
+        _printf0_("│"<<setw(3)<<iter<<" │ ");
+        _printf0_("f(x) = "<<setw(8)<<setprecision(5)<<J<<" │ ");
+        _printf0_("  "<<setw(8)<<setprecision(3)<<Gnorm<<" │");
+        for(int i=0;i<N;i++) _printf0_(" "<<setw(10)<<setprecision(5)<<Jlist[i]);
+        _printf0_("\n");
+}
+void InversionStatsFooter(int NJ){
+        _printf0_("└────┴─────────────────┴────────────┴────────────────────────\n");
+}

issm/trunk-jpl/src/c/shared/Elements/elements.h

@@ -59 +59 @@
 void printsparsity(IssmPDouble* array,int lines,int cols=1);
 void printbinary(int n);
+void InversionStatsHeader(int NJ);
+void InversionStatsIter(int iter,double J, double Gnorm, double* Jlist, int N);
+void InversionStatsFooter(int NJ);
 #endif //ifndef _SHARED_ELEMENTS_H_