Add abs error printout to TacsPrintErrorComponents

A-CGray · Jun 1, 2023 · eaf712c · eaf712c
1 parent 90b1a72
commit eaf712c
Showing 1 changed file with 5 additions and 4 deletions.
diff --git a/src/elements/TACSElementVerification.cpp b/src/elements/TACSElementVerification.cpp
@@ -84,13 +84,14 @@ double TacsGetMaxRelError(TacsScalar *a, TacsScalar *b, int size,
 */
 void TacsPrintErrorComponents(FILE *fp, const char *descript, TacsScalar *a,
                               TacsScalar *b, int size) {
-  fprintf(fp, "%*s[   ] %15s %15s %15s\n", (int)strlen(descript), "Val",
-          "Analytic", "Approximate", "Rel. Error");
+  fprintf(fp, "%*s[   ] %15s %15s %15s %15s\n", (int)strlen(descript), "Val",
+          "Analytic", "Approximate", "Rel. Error", "Abs. Error");
   for (int i = 0; i < size; i++) {
     if (a[i] != 0.0) {
-      fprintf(fp, "%s[%3d] %15.6e %15.6e %15.4e\n", descript, i,
+      fprintf(fp, "%s[%3d] %15.6e %15.6e %15.4e %15.4e\n", descript, i,
               TacsRealPart(a[i]), TacsRealPart(b[i]),
-              fabs(TacsRealPart((a[i] - b[i]) / a[i])));
+              fabs(TacsRealPart((a[i] - b[i]) / a[i])),
+              fabs(TacsRealPart(a[i] - b[i])));
     } else {
       fprintf(fp, "%s[%3d] %15.6e %15.6e\n", descript, i, TacsRealPart(a[i]),
               TacsRealPart(b[i]));