print subnormals/ inertia

pghysels · Jul 14, 2023 · a081f5f · a081f5f
1 parent 3c7d124
commit a081f5f
Show file tree

Hide file tree

Showing 2 changed files with 23 additions and 23 deletions.
diff --git a/examples/sparse/testMMdouble.cpp b/examples/sparse/testMMdouble.cpp
@@ -60,17 +60,17 @@ test(int argc, char* argv[], CSRMatrix<scalar_t,integer_t>& A) {
   }
   spss.solve(b.data(), x.data());
 
-  // std::size_t subs = 0, zeros = 0;
-  // auto err0 = spss.subnormals(subs, zeros);
-  // std::cout << "# SUBNORMALS = " << subs
-  //           << "   ZEROS = " << zeros
-  //           << " (" << err0 << ")" << std::endl;
+  std::size_t subs = 0, zeros = 0;
+  auto err0 = spss.subnormals(subs, zeros);
+  std::cout << "# SUBNORMALS = " << subs
+            << "   ZEROS = " << zeros
+            << " (" << err0 << ")" << std::endl;
 
-  // integer_t neg, zero, pos;
-  // auto err = spss.inertia(neg, zero, pos);
-  // std::cout << "# INERTIA neg,zero,pos = "
-  //           << neg << ", " << zero << ", " << pos
-  //           <<  " (" << err << ")" << std::endl;
+  integer_t neg, zero, pos;
+  auto err = spss.inertia(neg, zero, pos);
+  std::cout << "# INERTIA neg,zero,pos = "
+            << neg << ", " << zero << ", " << pos
+            <<  " (" << err << ")" << std::endl;
 
   std::cout << "# COMPONENTWISE SCALED RESIDUAL = "
             << A.max_scaled_residual(x.data(), b.data()) << std::endl;

diff --git a/examples/sparse/testMMdoubleMPIDist.cpp b/examples/sparse/testMMdoubleMPIDist.cpp
@@ -75,19 +75,19 @@ test(int argc, char* argv[], CSRMatrix<scalar,integer>& A) {
   }
   spss.solve(b.data(), x.data());
 
-  // std::size_t subs = 0, zeros = 0;
-  // auto err0 = spss.subnormals(subs, zeros);
-  // if (!rank)
-  //   std::cout << "# SUBNORMALS = " << subs
-  //             << "   ZEROS = " << zeros
-  //             << " (" << err0 << ")" << std::endl;
-
-  // integer neg, zero, pos;
-  // auto err = spss.inertia(neg, zero, pos);
-  // if (!rank)
-  //   std::cout << "# INERTIA neg,zero,pos = "
-  //             << neg << ", " << zero << ", " << pos
-  //             <<  " (" << err << ")" << std::endl;
+  std::size_t subs = 0, zeros = 0;
+  auto err0 = spss.subnormals(subs, zeros);
+  if (!rank)
+    std::cout << "# SUBNORMALS = " << subs
+              << "   ZEROS = " << zeros
+              << " (" << err0 << ")" << std::endl;
+
+  integer neg, zero, pos;
+  auto err = spss.inertia(neg, zero, pos);
+  if (!rank)
+    std::cout << "# INERTIA neg,zero,pos = "
+              << neg << ", " << zero << ", " << pos
+              <<  " (" << err << ")" << std::endl;
 
   auto scaled_res = Adist.max_scaled_residual(x.data(), b.data());
   if (!rank)