Merge pull request #3650 from dknez/preevaluate_parametrized_function_cleanup

Added RBParametrizedFunction::preevaluate_parametrized_function_cleanup()

Author: jwpeterson

include/reduced_basis/rb_parametrized_function.h

@@ -336,6 +336,14 @@ class RBParametrizedFunction
                                           const std::vector<subdomain_id_type> & sbd_ids,
                                           const std::vector<boundary_id_type> & boundary_ids);
 
+  /**
+   * Virtual function that performs cleanup after each
+   * "preevaluate parametrized function" evaluation. This function
+   * is a no-op by default, but it can be overridden in subclasses
+   * in order to do necessary cleanup, such as clearing cached data.
+   */
+  virtual void preevaluate_parametrized_function_cleanup();
+
   /**
    * Storage for pre-evaluated values. The indexing is given by:
    *   parameter index --> point index --> component index --> value.

src/reduced_basis/rb_parametrized_function.C

@@ -484,6 +484,8 @@ void RBParametrizedFunction::preevaluate_parametrized_function_on_mesh(const RBP
                       all_xyz_perturb_vec,
                       phi_i_qp_vec,
                       preevaluated_values);
+
+  preevaluate_parametrized_function_cleanup();
 }
 
 void RBParametrizedFunction::preevaluate_parametrized_function_on_mesh_sides(const RBParameters & mu,
@@ -639,6 +641,8 @@ void RBParametrizedFunction::preevaluate_parametrized_function_on_mesh_sides(con
                            all_xyz_perturb_vec,
                            phi_i_qp_vec,
                            preevaluated_values);
+
+  preevaluate_parametrized_function_cleanup();
 }
 
 void RBParametrizedFunction::preevaluate_parametrized_function_on_mesh_nodes(const RBParameters & mu,
@@ -677,6 +681,8 @@ void RBParametrizedFunction::preevaluate_parametrized_function_on_mesh_nodes(con
                            node_ids_vec,
                            boundary_ids_vec,
                            preevaluated_values);
+
+  preevaluate_parametrized_function_cleanup();
 }
 
 Number RBParametrizedFunction::lookup_preevaluated_value_on_mesh(unsigned int comp,
@@ -823,6 +829,11 @@ void RBParametrizedFunction::initialize_spatial_indices(const std::vector<std::v
   // No-op by default
 }
 
+void RBParametrizedFunction::preevaluate_parametrized_function_cleanup()
+{
+  // No-op by default
+}
+
 const std::set<boundary_id_type> & RBParametrizedFunction::get_parametrized_function_boundary_ids() const
 {
   return _parametrized_function_boundary_ids;