Leaks in block matrices.

Author: prj-

plugin/mpi/PETSc-code.hpp

@@ -1354,20 +1354,7 @@ namespace PETSc {
         }
       }
     }
-    ~initCSRfromBlockMatrix( ) {
-      if (e_M) {
-        for (int i = 0; i < N; ++i) {
-          delete[] e_M[i];
-          delete[] t_M[i];
-        }
-        delete[] e_M;
-        delete[] t_M;
-        N = 0;
-        M = 0;
-        e_M = nullptr;
-        t_M = nullptr;
-      }
-    }
+    ~initCSRfromBlockMatrix( ) { }
     static ArrayOfaType typeargs( ) {
       return ArrayOfaType(atype< Result >( ), atype< E_Array >( ));
     }
@@ -1503,6 +1490,14 @@ namespace PETSc {
           MatDestroy(a + p.first * M + p.second);
       sparse_mat->_exchange = reinterpret_cast<HPDDM::Subdomain<PetscScalar>**>(exchange);
       delete[] a;
+      if (e_M) {
+        for (int i = 0; i < N; ++i) {
+          delete[] e_M[i];
+          delete[] t_M[i];
+        }
+        delete[] e_M;
+        delete[] t_M;
+      }
       return sparse_mat;
     }
   };

src/bamglib/write_xdmf.cpp

@@ -1,3 +1,4 @@
+extern long verbosity;
 #include "write_xdmf.hpp"
 #ifdef DEBUG
 int debug_w = 1;