Utility class to use an Exodus mesh to define controllable parameters for optimization problems This class will: More...

#include <ParameterMesh.h>

Public Types
enum class	RegularizationType { L2_GRADIENT }
	Enumerations for regularization computations. More...

Public Member Functions
	ParameterMesh (const libMesh::FEType &param_type, const std::string &exodus_mesh, const bool find_closest=false, const unsigned int kdtree_candidates=5)

dof_id_type	size () const

void	getIndexAndWeight (const Point &pt, std::vector< dof_id_type > &dof_indices, std::vector< Real > &weights) const
	Interpolate parameters onto the computational mesh getIndexAndWeight is only used by ParameterMeshFunction.

void	getIndexAndWeight (const Point &pt, std::vector< dof_id_type > &dof_indices, std::vector< RealGradient > &weights) const
	Performs inner products of parameters with functions on the computational mesh getIndexAndWeight is only used by ParameterMeshFunction.

Real	computeRegularizationObjective (const std::vector< Real > &parameter_values, RegularizationType reg_type) const
	Computes regularization objective value for a given regularization type.

std::vector< Real >	computeRegularizationGradient (const std::vector< Real > &parameter_values, RegularizationType reg_type) const
	Computes regularization gradient for a given regularization type.

Protected Attributes
libMesh::Parallel::Communicator	_communicator

libMesh::ReplicatedMesh	_mesh

const bool	_find_closest
	Find closest projection points.

std::unique_ptr< libMesh::EquationSystems >	_eq

libMesh::System *	_sys

std::unique_ptr< libMesh::PointLocatorBase >	_point_locator

std::unique_ptr< libMesh::ExodusII_IO >	_exodusII_io

dof_id_type	_param_dofs

std::vector< Point >	_mesh_nodes
	Node-based KDTree optimization.

std::unique_ptr< KDTree >	_node_kdtree

std::unordered_map< dof_id_type, std::set< const libMesh::Elem * > >	_node_to_elements

unsigned int	_kdtree_candidates

Private Member Functions
template<typename T >
T	computeRegularizationLoop (const std::vector< Real > &parameter_values, RegularizationType reg_type) const
	Template method containing the element loop for regularization computations.

Point	projectToMesh (const Point &p) const
	Returns the point on the parameter mesh that is projected from the test point.

Point	closestPoint (const Elem &elem, const Point &p) const
	Find closest point on the element to the given point.

Real	computeRegularizationQp (const std::vector< Real > &parameter_values, const std::vector< std::vector< Real > > &phi, const std::vector< std::vector< RealGradient > > &dphi, const unsigned int qp, const std::vector< dof_id_type > &dof_indices, const std::vector< Real > &JxW, RegularizationType reg_type) const
	Compute regularization objective for a single quadrature point This is the main function users should modify to add new regularization types for objectives.

void	computeRegularizationGradientQp (const std::vector< Real > &parameter_values, const std::vector< std::vector< Real > > &phi, const std::vector< std::vector< RealGradient > > &dphi, const unsigned int qp, const std::vector< dof_id_type > &dof_indices, const std::vector< Real > &JxW, RegularizationType reg_type, std::vector< Real > &gradient) const
	Compute regularization gradient for a single quadrature point This is the main function users should modify to add new regularization types for gradients.

Private Attributes
const unsigned short int	_param_var_id

const libMesh::DofMap *	_dof_map

const libMesh::FEType	_fe_type

Detailed Description

Utility class to use an Exodus mesh to define controllable parameters for optimization problems This class will:

Ensure that controllable parameters defined by the mesh are correctly ordered on optimiation main app and forward or adjoint sub-apps
Define the parameter space (nodal/elemental and shape function)
Interpolate parameter values to the forward problem mesh using nodal/elemental shape function

Definition at line 41 of file ParameterMesh.h.

Member Enumeration Documentation

◆ RegularizationType

enum class ParameterMesh::RegularizationType

strong

Enumerations for regularization computations.

Enumerator
L2_GRADIENT

Definition at line 50 of file ParameterMesh.h.

  {
    L2_GRADIENT,
    // Future regularization types can be added here:
    // L1,
    // H1,
    // TV (Total Variation)
  };

Constructor & Destructor Documentation

◆ ParameterMesh()

ParameterMesh::ParameterMesh	(	const libMesh::FEType &	param_type,
		const std::string &	exodus_mesh,
		const bool	find_closest = `false`,
		const unsigned int	kdtree_candidates = `5`
	)

Definition at line 38 of file ParameterMesh.C.

  : _communicator(MPI_COMM_SELF),
    _mesh(_communicator),
    _find_closest(find_closest),
    _kdtree_candidates(kdtree_candidates),
    _param_var_id(0),
    _dof_map(nullptr),
    _fe_type(param_type)
{
  _mesh.allow_renumbering(false);
  _mesh.prepare_for_use();
  _exodusII_io = std::make_unique<ExodusII_IO>(_mesh);
  _exodusII_io->read(exodus_mesh);
  _mesh.read(exodus_mesh);
  // Create system to store parameter values
  _eq = std::make_unique<libMesh::EquationSystems>(_mesh);
  _sys = &_eq->add_system<ExplicitSystem>("_parameter_mesh_sys");
  _sys->add_variable("_parameter_mesh_var", param_type);
 
  // Create point locator
  _point_locator = PointLocatorBase::build(TREE_LOCAL_ELEMENTS, _mesh);
  _point_locator->enable_out_of_mesh_mode();
 
  // Initialize the equations systems
  _eq->init();
 
  // getting number of parameter dofs for size() function
  const unsigned short int var_id = _sys->variable_number("_parameter_mesh_var");
  std::set<dof_id_type> var_indices;
  _sys->local_dof_indices(var_id, var_indices);
  _param_dofs = var_indices.size();
 
  if (_find_closest)
  {
    for (const auto & elem : _mesh.element_ptr_range())
      if (elem->default_order() != FIRST)
        mooseError("Closet point projection currently does not support second order elements.");
  }
 
  // Initialize node-based KDTree optimization
  _mesh_nodes.clear();
  _node_to_elements.clear();
 
  // Extract all node coordinates
  for (const auto & node : _mesh.node_ptr_range())
    _mesh_nodes.push_back(*node);
 
  // Build node-to-elements connectivity map
  for (const auto & elem : _mesh.element_ptr_range())
  {
    for (const auto n : make_range(elem->n_nodes()))
    {
      dof_id_type node_id = elem->node_id(n);
      _node_to_elements[node_id].insert(elem);
    }
  }
 
  // Create KDTree from node coordinates
  if (!_mesh_nodes.empty())
    _node_kdtree = std::make_unique<KDTree>(_mesh_nodes, 10);
  // Update cached values for gradient computations
  const_cast<unsigned short int &>(_param_var_id) = var_id;
  const_cast<const libMesh::DofMap *&>(_dof_map) = &_sys->get_dof_map();
  const_cast<FEType &>(_fe_type) = _dof_map->variable_type(_param_var_id);
}

Member Function Documentation

◆ closestPoint()

Point ParameterMesh::closestPoint	(	const Elem &	elem,
		const Point &	p
	)		const

private

Find closest point on the element to the given point.

Parameters

elem
p

Returns: Point

Definition at line 235 of file ParameterMesh.C.

{
  mooseAssert(!elem.contains_point(p),
              "Points inside of elements shouldn't need to find closestPoint.");
 
  // Lambda to find closest point from range without storing temporary vectors
  auto findClosest = [&p](auto range, auto point_func) -> Point
  {
    Real min_distance = std::numeric_limits<Real>::max();
    Point min_point = p;
 
    for (const auto & item : range)
    {
      Point candidate = point_func(item);
      Real distance = (candidate - p).norm();
      if (distance < min_distance)
      {
        min_distance = distance;
        min_point = candidate;
      }
    }
    return min_point;
  };
 
  switch (elem.type())
  {
    case EDGE2:
    {
      LineSegment ls(*(elem.node_ptr(0)), *(elem.node_ptr(1)));
      return ls.closest_point(p);
    }
 
    case TRI3:
    {
      Point a = *(elem.node_ptr(0));
      Point b = *(elem.node_ptr(1));
      Point c = *(elem.node_ptr(2));
      libMesh::Plane pl(a, b, c);
      Point trial = pl.closest_point(p);
      if (elem.contains_point(trial))
        return trial;
 
      return findClosest(make_range(elem.n_edges()),
                         [&](dof_id_type i) { return closestPoint(*elem.build_edge_ptr(i), p); });
    }
    case QUAD4:
    {
      Point a = *(elem.node_ptr(0));
      Point b = *(elem.node_ptr(1));
      Point c = *(elem.node_ptr(2));
      Point d = *(elem.node_ptr(3));
      libMesh::Plane pl1(a, b, c);
      libMesh::Plane pl2(b, c, d);
      Point trial1 = pl1.closest_point(p);
      Point trial2 = pl2.closest_point(p);
      if (!trial1.absolute_fuzzy_equals(trial2, TOLERANCE * TOLERANCE))
        mooseError("Quad4 element is not coplanar");
 
      if (elem.contains_point(trial1))
        return trial1;
 
      return findClosest(make_range(elem.n_edges()),
                         [&](dof_id_type i) { return closestPoint(*elem.build_edge_ptr(i), p); });
    }
 
    default:
    {
      if (elem.dim() == 3)
      {
        return findClosest(make_range(elem.n_sides()),
                           [&](dof_id_type i) { return closestPoint(*elem.build_side_ptr(i), p); });
      }
      else
      {
        mooseError("Unsupported element type ",
                   Utility::enum_to_string(elem.type()),
                   " for projection of parameter mesh.");
      }
    }
  }
}

Referenced by projectToMesh().

◆ computeRegularizationGradient()

std::vector< Real > ParameterMesh::computeRegularizationGradient	(	const std::vector< Real > &	parameter_values,
		RegularizationType	reg_type
	)		const

Computes regularization gradient for a given regularization type.

Parameters

parameter_values	vector of parameter values to compute gradient for
reg_type	type of regularization (L2_GRADIENT, etc.)

Returns: vector of gradient values (same size as parameter_values)

Definition at line 380 of file ParameterMesh.C.

{
  return computeRegularizationLoop<std::vector<Real>>(parameter_values, reg_type);
}

◆ computeRegularizationGradientQp()

void ParameterMesh::computeRegularizationGradientQp	(	const std::vector< Real > &	parameter_values,
		const std::vector< std::vector< Real > > &	phi,
		const std::vector< std::vector< RealGradient > > &	dphi,
		const unsigned int	qp,
		const std::vector< dof_id_type > &	dof_indices,
		const std::vector< Real > &	JxW,
		RegularizationType	reg_type,
		std::vector< Real > &	gradient
	)		const

private

Compute regularization gradient for a single quadrature point This is the main function users should modify to add new regularization types for gradients.

Parameters

parameter_values	all parameter values
phi	shape function values (full array)
dphi	shape function gradients (full array)
qp	quadrature point index
dof_indices	element DOF indices
JxW	quadrature weights array
reg_type	type of regularization to compute
gradient	gradient vector to update

Definition at line 419 of file ParameterMesh.C.

{
  // Switch on regularization type
  switch (reg_type)
  {
    case RegularizationType::L2_GRADIENT:
    {
      // Compute parameter gradient at this quadrature point
      RealGradient param_grad;
      for (const auto i : index_range(dof_indices))
        param_grad += parameter_values[dof_indices[i]] * dphi[i][qp];
 
      // Compute gradient contribution: 2 * grad(p) * dphi_j
      for (const auto j : index_range(dof_indices))
        gradient[dof_indices[j]] += 2.0 * param_grad * dphi[j][qp] * JxW[qp];
      break;
    }
 
    default:
      mooseError("Unknown Regularization Type");
  }
}

Referenced by computeRegularizationLoop().

◆ computeRegularizationLoop()

template<typename T >

T ParameterMesh::computeRegularizationLoop	(	const std::vector< Real > &	parameter_values,
		RegularizationType	reg_type
	)		const

private

Template method containing the element loop for regularization computations.

Parameters

parameter_values	vector of parameter values
reg_type	type of regularization

Returns: result of type T (Real for objective, std::vector<Real> for gradient)

Definition at line 319 of file ParameterMesh.C.

{
  if (parameter_values.size() != _param_dofs)
    mooseError("Parameter values size (",
               parameter_values.size(),
               ") does not match mesh DOFs (",
               _param_dofs,
               ")");
 
  T result;
  if constexpr (std::is_same_v<T, Real>)
    result = 0.0;
  else if constexpr (std::is_same_v<T, std::vector<Real>>)
    result.resize(_param_dofs, 0.0);
 
  // Iterate over all elements in the mesh
  for (const auto & elem : _mesh.element_ptr_range())
  {
    // Get DOF indices for this element
    std::vector<dof_id_type> dof_indices;
    _dof_map->dof_indices(elem, dof_indices, _param_var_id);
 
    // Get quadrature rule for this element
    const unsigned int dim = elem->dim();
    QGauss qrule(dim, _fe_type.default_quadrature_order());
 
    // Create finite element objects
    std::unique_ptr<FEBase> fe(FEBase::build(dim, _fe_type));
    fe->attach_quadrature_rule(&qrule);
 
    // Request shape functions and derivatives before reinit
    const std::vector<Real> & JxW = fe->get_JxW();
    const std::vector<std::vector<Real>> & phi = fe->get_phi();
    const std::vector<std::vector<RealGradient>> & dphi = fe->get_dphi();
 
    // Reinitialize for current element
    fe->reinit(elem);
 
    for (const auto qp : make_range(qrule.n_points()))
    {
      if constexpr (std::is_same_v<T, Real>)
        result +=
            computeRegularizationQp(parameter_values, phi, dphi, qp, dof_indices, JxW, reg_type);
      else if constexpr (std::is_same_v<T, std::vector<Real>>)
        computeRegularizationGradientQp(
            parameter_values, phi, dphi, qp, dof_indices, JxW, reg_type, result);
    }
  }
 
  return result;
}

◆ computeRegularizationObjective()

Real ParameterMesh::computeRegularizationObjective	(	const std::vector< Real > &	parameter_values,
		RegularizationType	reg_type
	)		const

Computes regularization objective value for a given regularization type.

Parameters

parameter_values	vector of parameter values to compute regularization for
reg_type	type of regularization (L2_GRADIENT, etc.)

Returns: scalar objective value

Definition at line 373 of file ParameterMesh.C.

{
  return computeRegularizationLoop<Real>(parameter_values, reg_type);
}

◆ computeRegularizationQp()

Real ParameterMesh::computeRegularizationQp	(	const std::vector< Real > &	parameter_values,
		const std::vector< std::vector< Real > > &	phi,
		const std::vector< std::vector< RealGradient > > &	dphi,
		const unsigned int	qp,
		const std::vector< dof_id_type > &	dof_indices,
		const std::vector< Real > &	JxW,
		RegularizationType	reg_type
	)		const

private

Compute regularization objective for a single quadrature point This is the main function users should modify to add new regularization types for objectives.

Parameters

parameter_values	all parameter values
phi	shape function values (full array)
dphi	shape function gradients (full array)
qp	quadrature point index
dof_indices	element DOF indices
JxW	quadrature weights array
reg_type	type of regularization to compute

Returns: contribution to objective function

Definition at line 387 of file ParameterMesh.C.

{
  Real objective_contribution = 0.0;
 
  // Switch on regularization type
  switch (reg_type)
  {
    case RegularizationType::L2_GRADIENT:
    {
      // Compute parameter gradient at this quadrature point
      RealGradient param_grad;
      for (const auto i : index_range(dof_indices))
        param_grad += parameter_values[dof_indices[i]] * dphi[i][qp];
 
      // Add L2 norm squared of gradient for regularization
      objective_contribution = param_grad.norm_sq() * JxW[qp];
      break;
    }
    default:
      mooseError("Unknown Regularization Type");
  }
 
  return objective_contribution;
}

Referenced by computeRegularizationLoop().

◆ getIndexAndWeight() [1/2]

void ParameterMesh::getIndexAndWeight	(	const Point &	pt,
		std::vector< dof_id_type > &	dof_indices,
		std::vector< Real > &	weights
	)		const

Interpolate parameters onto the computational mesh getIndexAndWeight is only used by ParameterMeshFunction.

Parameters

pt	location to compute elemnent dof_indices weights
dof_indices	return dof indices for element containing pt
weights	returns element shape function weights at pt

Definition at line 108 of file ParameterMesh.C.

{
  Point test_point = (_find_closest ? projectToMesh(pt) : pt);
 
  const Elem * elem = (*_point_locator)(test_point);
  if (!elem)
    mooseError("No element was found to contain point ", test_point);
 
  // Get the dof_indices for our element
  // variable id is hard coded to _param_var_id
  // this is probably the only variable in the ParameterMesh system used by ParameterMeshFunction
  _dof_map->dof_indices(elem, dof_indices, _param_var_id);
 
  // Map the physical co-ordinates to the reference co-ordinates
  Point coor = FEMap::inverse_map(elem->dim(), elem, test_point);
  // get the shape function value via the FEInterface
  libMesh::FEComputeData fe_data(*_eq, coor);
  FEInterface::compute_data(elem->dim(), _fe_type, elem, fe_data);
  // Set weights to the value of the shape functions
  weights = fe_data.shape;
 
  if (dof_indices.size() != weights.size())
    mooseError("Internal error: weights and DoF indices do not have the same size.");
}

Referenced by ParameterMeshFunction::gradient(), ParameterMeshFunction::parameterGradient(), ParameterMeshFunction::timeDerivative(), and ParameterMeshFunction::value().

◆ getIndexAndWeight() [2/2]

void ParameterMesh::getIndexAndWeight	(	const Point &	pt,
		std::vector< dof_id_type > &	dof_indices,
		std::vector< RealGradient > &	weights
	)		const

Performs inner products of parameters with functions on the computational mesh getIndexAndWeight is only used by ParameterMeshFunction.

Parameters

pt	location to compute elemnent dof_indices weights
dof_indices	return dof indices for element containing pt
weights	returns element shape function gradient weights at pt

Definition at line 136 of file ParameterMesh.C.

{
  if (!_sys->has_variable("_parameter_mesh_var"))
    mooseError("Internal error: System being read does not contain _parameter_mesh_var.");
  Point test_point = (_find_closest ? projectToMesh(pt) : pt);
  // Locate the element the point is in
  const Elem * elem = (*_point_locator)(test_point);
 
  // Get the dof_indices for our element
  // variable id is hard coded to _param_var_id
  // this is probably the only variable in the ParameterMesh system used by ParameterMeshFunction
  _dof_map->dof_indices(elem, dof_indices, _param_var_id);
 
  // Map the physical co-ordinates to the reference co-ordinates
  Point coor = FEMap::inverse_map(elem->dim(), elem, test_point);
  // get the shape function value via the FEInterface
  libMesh::FEComputeData fe_data(*_eq, coor);
  fe_data.enable_derivative();
  FEInterface::compute_data(elem->dim(), _fe_type, elem, fe_data);
  // Set weights to the value of the shape functions
  weights = fe_data.dshape;
 
  if (dof_indices.size() != weights.size())
    mooseError("Internal error: weights and DoF indices do not have the same size.");
}

◆ projectToMesh()

Point ParameterMesh::projectToMesh ( const Point & p ) const

private

Returns the point on the parameter mesh that is projected from the test point.

Parameters

p	test point

Returns: Point

Definition at line 165 of file ParameterMesh.C.

{
  // quick path: p already inside an element
  if ((*_point_locator)(p))
    return p;
 
  // Lambda to find closest point from elements using squared distance for efficiency
  auto findClosestElement = [&p, this](const auto & elements) -> Point
  {
    Real best_d2 = std::numeric_limits<Real>::max();
    Point best_point = p;
 
    for (const auto * elem : elements)
    {
      Point trial = closestPoint(*elem, p);
      Real d2 = (trial - p).norm_sq();
      if (d2 < best_d2)
      {
        best_d2 = d2;
        best_point = trial;
      }
    }
 
    if (best_d2 == std::numeric_limits<Real>::max())
      mooseError("project_to_mesh failed - no candidate elements.");
    return best_point;
  };
 
  // Use KDTree optimization if available
  if (_node_kdtree && !_mesh_nodes.empty())
  {
    // Find K nearest nodes using KDTree
    std::vector<std::size_t> nearest_node_indices;
    _node_kdtree->neighborSearch(p, _kdtree_candidates, nearest_node_indices);
 
    // Collect all elements connected to these nodes
    std::set<const Elem *> candidate_elements;
    for (auto node_idx : nearest_node_indices)
    {
      // Get the actual node from the mesh using the index
      if (node_idx < _mesh.n_nodes())
      {
        const Node * node = _mesh.node_ptr(node_idx);
        dof_id_type node_id = node->id();
        auto it = _node_to_elements.find(node_id);
        if (it != _node_to_elements.end())
        {
          const auto & connected_elems = it->second;
          candidate_elements.insert(connected_elems.begin(), connected_elems.end());
        }
      }
    }
 
    // Convert set to vector for consistent type
    std::vector<const Elem *> candidate_vector(candidate_elements.begin(),
                                               candidate_elements.end());
    return findClosestElement(candidate_vector);
  }
  else
  {
    // Fallback to original O(n) method if KDTree not available
    std::vector<const Elem *> all_elements;
    for (const auto & elem : _mesh.element_ptr_range())
      all_elements.push_back(elem);
 
    return findClosestElement(all_elements);
  }
}

Referenced by getIndexAndWeight(), and getIndexAndWeight().

◆ size()

dof_id_type ParameterMesh::size ( ) const

inline

Returns: the number of parameters read from the mesh for a single timestep

Definition at line 62 of file ParameterMesh.h.

62{ return _param_dofs; }

Referenced by ParameterMeshFunction::checkSize(), ParameterMeshFunction::gradient(), ParameterMeshFunction::parameterGradient(), ParameterMeshFunction::timeDerivative(), and ParameterMeshFunction::value().

Member Data Documentation

◆ _communicator

libMesh::Parallel::Communicator ParameterMesh::_communicator

protected

Definition at line 113 of file ParameterMesh.h.

◆ _dof_map

const libMesh::DofMap* ParameterMesh::_dof_map

private

Definition at line 188 of file ParameterMesh.h.

Referenced by computeRegularizationLoop(), getIndexAndWeight(), getIndexAndWeight(), and ParameterMesh().

◆ _eq

std::unique_ptr<libMesh::EquationSystems> ParameterMesh::_eq

protected

Definition at line 117 of file ParameterMesh.h.

Referenced by getIndexAndWeight(), getIndexAndWeight(), and ParameterMesh().

◆ _exodusII_io

std::unique_ptr<libMesh::ExodusII_IO> ParameterMesh::_exodusII_io

protected

Definition at line 120 of file ParameterMesh.h.

Referenced by ParameterMesh().

◆ _fe_type

const libMesh::FEType ParameterMesh::_fe_type

private

Definition at line 189 of file ParameterMesh.h.

Referenced by computeRegularizationLoop(), getIndexAndWeight(), getIndexAndWeight(), and ParameterMesh().

◆ _find_closest

const bool ParameterMesh::_find_closest

protected

Find closest projection points.

Definition at line 116 of file ParameterMesh.h.

Referenced by getIndexAndWeight(), getIndexAndWeight(), and ParameterMesh().

◆ _kdtree_candidates

unsigned int ParameterMesh::_kdtree_candidates

protected

Definition at line 128 of file ParameterMesh.h.

Referenced by projectToMesh().

◆ _mesh

libMesh::ReplicatedMesh ParameterMesh::_mesh

protected

Definition at line 114 of file ParameterMesh.h.

Referenced by computeRegularizationLoop(), ParameterMesh(), and projectToMesh().

◆ _mesh_nodes

std::vector<Point> ParameterMesh::_mesh_nodes

protected

Node-based KDTree optimization.

Definition at line 125 of file ParameterMesh.h.

Referenced by ParameterMesh(), and projectToMesh().

◆ _node_kdtree

std::unique_ptr<KDTree> ParameterMesh::_node_kdtree

protected

Definition at line 126 of file ParameterMesh.h.

Referenced by ParameterMesh(), and projectToMesh().

◆ _node_to_elements

std::unordered_map<dof_id_type, std::set<const libMesh::Elem *> > ParameterMesh::_node_to_elements

protected

Definition at line 127 of file ParameterMesh.h.

Referenced by ParameterMesh(), and projectToMesh().

◆ _param_dofs

dof_id_type ParameterMesh::_param_dofs

protected

Definition at line 122 of file ParameterMesh.h.

Referenced by computeRegularizationLoop(), ParameterMesh(), and size().

◆ _param_var_id

const unsigned short int ParameterMesh::_param_var_id

private

Definition at line 187 of file ParameterMesh.h.

Referenced by computeRegularizationLoop(), getIndexAndWeight(), getIndexAndWeight(), and ParameterMesh().

◆ _point_locator

std::unique_ptr<libMesh::PointLocatorBase> ParameterMesh::_point_locator

protected

Definition at line 119 of file ParameterMesh.h.

Referenced by ParameterMesh(), and projectToMesh().

◆ _sys

libMesh::System* ParameterMesh::_sys

protected

Definition at line 118 of file ParameterMesh.h.

Referenced by getIndexAndWeight(), and ParameterMesh().

The documentation for this class was generated from the following files:

optimization/include/utils/ParameterMesh.h
optimization/src/utils/ParameterMesh.C

Public Types

Public Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Member Enumeration Documentation

◆ RegularizationType

Constructor & Destructor Documentation

◆ ParameterMesh()

Member Function Documentation

◆ closestPoint()

◆ computeRegularizationGradient()

◆ computeRegularizationGradientQp()

◆ computeRegularizationLoop()

◆ computeRegularizationObjective()

◆ computeRegularizationQp()

◆ getIndexAndWeight() [1/2]

◆ getIndexAndWeight() [2/2]

◆ projectToMesh()

◆ size()

Member Data Documentation

◆ _communicator

◆ _dof_map

◆ _eq

◆ _exodusII_io

◆ _fe_type

◆ _find_closest

◆ _kdtree_candidates

◆ _mesh

◆ _mesh_nodes

◆ _node_kdtree

◆ _node_to_elements

◆ _param_dofs

◆ _param_var_id

◆ _point_locator

◆ _sys