Fix memory leak for direction_dependent BCs (#147)

Author: cgilet

src/projection/incflo_apply_nodal_projection.cpp

@@ -166,16 +166,16 @@ void incflo::ApplyNodalProjection (Vector<MultiFab const*> const& density,
     // Enforce solvability by matching outflow to inflow.
     if (has_inout_bndry && set_inflow_bc)
     {
-        Vector<Array<MultiFab*, AMREX_SPACEDIM>> vel_vec(finest_level+1);
+        Vector<Array<MultiFab, AMREX_SPACEDIM>> vel_vec(finest_level+1);
 
         for (int lev = 0; lev <= finest_level; lev++) {
             auto& ld = *m_leveldata[lev];
-            AMREX_D_TERM(vel_vec[lev][0] = new MultiFab(ld.velocity, amrex::make_alias, 0, 1);,
-                         vel_vec[lev][1] = new MultiFab(ld.velocity, amrex::make_alias, 1, 1);,
-                         vel_vec[lev][2] = new MultiFab(ld.velocity, amrex::make_alias, 2, 1););
+            AMREX_D_TERM(vel_vec[lev][0] = MultiFab(ld.velocity, amrex::make_alias, 0, 1);,
+                         vel_vec[lev][1] = MultiFab(ld.velocity, amrex::make_alias, 1, 1);,
+                         vel_vec[lev][2] = MultiFab(ld.velocity, amrex::make_alias, 2, 1););
         }
 
-        HydroUtils::enforceInOutSolvability(vel_vec, get_velocity_bcrec().data(), geom, true);
+        HydroUtils::enforceInOutSolvability(GetVecOfArrOfPtrs(vel_vec), get_velocity_bcrec().data(), geom, true);
     }
 
     LPInfo info;