Add if guards to verify deal.II version

chaos-polymtl · Aug 20, 2024 · d21a346 · d21a346
1 parent 0a75bfd
commit d21a346
Show file tree

Hide file tree

Showing 2 changed files with 24 additions and 0 deletions.
diff --git a/prototypes/direct_gls_navier_stokes/direct_gls_navier_stokes.cc b/prototypes/direct_gls_navier_stokes/direct_gls_navier_stokes.cc
@@ -535,7 +535,13 @@ DirectSteadyGLSNavierStokes<dim>::refine_mesh()
   triangulation.execute_coarsening_and_refinement();
   setup_dofs();
   BlockVector<double> tmp(dofs_per_block);
+
+#if DEAL_II_VERSION_GTE(9, 7, 0)
   solution_transfer.interpolate(tmp);
+#else
+  solution_transfer.interpolate(tmp, present_solution);
+#endif
+
   nonzero_constraints.distribute(tmp);
   initialize_system();
   present_solution = tmp;
@@ -550,7 +556,13 @@ DirectSteadyGLSNavierStokes<dim>::refine_mesh_uniform()
   triangulation.refine_global(1);
   setup_dofs();
   BlockVector<double> tmp(dofs_per_block);
+
+#if DEAL_II_VERSION_GTE(9, 7, 0)
   solution_transfer.interpolate(tmp);
+#else
+  solution_transfer.interpolate(tmp, present_solution);
+#endif
+
   nonzero_constraints.distribute(tmp);
   initialize_system();
   present_solution = tmp;

diff --git a/prototypes/direct_steady_navier_stokes/direct_steady_navier_stokes.cc b/prototypes/direct_steady_navier_stokes/direct_steady_navier_stokes.cc
@@ -429,7 +429,13 @@ DirectSteadyNavierStokes<dim>::refine_mesh()
   triangulation.execute_coarsening_and_refinement();
   setup_dofs();
   BlockVector<double> tmp(dofs_per_block);
+
+#if DEAL_II_VERSION_GTE(9, 7, 0)
   solution_transfer.interpolate(tmp);
+#else
+  solution_transfer.interpolate(tmp, present_solution);
+#endif
+
   nonzero_constraints.distribute(tmp);
   initialize_system();
   present_solution = tmp;
@@ -444,7 +450,13 @@ DirectSteadyNavierStokes<dim>::refine_mesh_uniform()
   triangulation.refine_global(1);
   setup_dofs();
   BlockVector<double> tmp(dofs_per_block);
+
+#if DEAL_II_VERSION_GTE(9, 7, 0)
   solution_transfer.interpolate(tmp);
+#else
+  solution_transfer.interpolate(tmp, present_solution);
+#endif
+
   nonzero_constraints.distribute(tmp);
   initialize_system();
   present_solution = tmp;