Helper class to define, parameterize and create a level set function used in meshing, often to correct the position of nodes on a surface. More...

#include <LevelSetMeshingHelper.h>

Inheritance diagram for LevelSetMeshingHelper:

[legend]

Public Types
typedef GenericSymFunction< is_ad >	SymFunction
	Shorthand for an autodiff function parser object.

typedef std::shared_ptr< SymFunction >	SymFunctionPtr
	Shorthand for an smart pointer to an autodiff function parser object.

Public Member Functions
	LevelSetMeshingHelper (const InputParameters &parameters)

void	setParserFeatureFlags (SymFunctionPtr &) const
	apply input parameters to internal feature flags of the parser object

Static Public Member Functions
static InputParameters	validParams ()

Protected Types
enum	FailureMethod
	Enum for failure method. More...

Protected Member Functions
Real	levelSetEvaluator (const Point &point)
	Evaluate the level set function at a given point.

void	levelSetCorrection (Node &node)
	Correct the position of a node based on the level set function.

GenericReal< is_ad >	evaluate (SymFunctionPtr &, const std::string &object_name="")
	Evaluate FParser object and check EvalError.

GenericReal< is_ad >	evaluate (SymFunctionPtr &, const std::vector< GenericReal< is_ad > > &, const std::string &object_name="")
	Evaluate FParser object and check EvalError.

void	addFParserConstants (SymFunctionPtr &parser, const std::vector< std::string > &constant_names, const std::vector< std::string > &constant_expressions) const
	add constants (which can be complex expressions) to the parser object

virtual void	functionsOptimize (SymFunctionPtr &parsed_function)
	run FPOptimizer on the parsed function

void	functionsOptimize (SymFunctionPtr &parsed_function)

void	functionsOptimize (SymFunctionPtr &parsed_function)

void	parsedFunctionSetup (SymFunctionPtr &function, const std::string &expression, const std::string &variables, const std::vector< std::string > &constant_names, const std::vector< std::string > &constant_expressions, const libMesh::Parallel::Communicator &comm) const
	Performs setup steps on a SymFunction.

Protected Attributes
const unsigned int	_max_level_set_correction_iterations
	Maximum number of iterations to correct the nodes based on the level set function.

SymFunctionPtr	_func_level_set
	function parser object describing the level set

enum FunctionParserUtils::FailureMethod	_evalerror_behavior

const Real	_quiet_nan
	appropriate not a number value to return

std::vector< GenericReal< is_ad > >	_func_params
	Array to stage the parameters passed to the functions when calling Eval.

const Real	_epsilon
	fuzzy comparison tolerance


bool	_enable_jit
	feature flags

bool	_enable_ad_cache

bool	_disable_fpoptimizer

bool	_enable_auto_optimize

Static Protected Attributes
static const char *	_eval_error_msg []
	table of FParser eval error codes

Detailed Description

Helper class to define, parameterize and create a level set function used in meshing, often to correct the position of nodes on a surface.

Definition at line 19 of file LevelSetMeshingHelper.h.

Member Typedef Documentation

◆ SymFunction

typedef GenericSymFunction<is_ad> FunctionParserUtils< is_ad >::SymFunction

inherited

Shorthand for an autodiff function parser object.

Definition at line 62 of file FunctionParserUtils.h.

◆ SymFunctionPtr

typedef std::shared_ptr<SymFunction> FunctionParserUtils< is_ad >::SymFunctionPtr

inherited

Shorthand for an smart pointer to an autodiff function parser object.

Definition at line 65 of file FunctionParserUtils.h.

Member Enumeration Documentation

◆ FailureMethod

enum FunctionParserUtils::FailureMethod

strongprotectedinherited

Enum for failure method.

Definition at line 116 of file FunctionParserUtils.h.

116{ nan, nan_warning, error, exception } _evalerror_behavior;

FunctionParserUtils< false >::_evalerror_behavior

enum FunctionParserUtils::FailureMethod _evalerror_behavior

Constructor & Destructor Documentation

◆ LevelSetMeshingHelper()

LevelSetMeshingHelper::LevelSetMeshingHelper ( const InputParameters & parameters )

Definition at line 31 of file LevelSetMeshingHelper.C.

  : FunctionParserUtils<false>(parameters),
    _max_level_set_correction_iterations(
        parameters.get<unsigned int>("max_level_set_correction_iterations"))
{
  if (parameters.isParamValid("level_set"))
  {
    _func_level_set = std::make_shared<SymFunction>();
    // set FParser internal feature flags
    setParserFeatureFlags(_func_level_set);
    if (parameters.isParamValid("constant_names") &&
        parameters.isParamValid("constant_expressions"))
      addFParserConstants(_func_level_set,
                          parameters.get<std::vector<std::string>>("constant_names"),
                          parameters.get<std::vector<std::string>>("constant_expressions"));
    if (_func_level_set->Parse(parameters.get<std::string>("level_set"), "x,y,z") >= 0)
      mooseError(
          "Invalid function f(x,y,z)\n", _func_level_set, ".\n", _func_level_set->ErrorMsg());
 
    _func_params.resize(3);
  }
}

Member Function Documentation

◆ addFParserConstants()

void FunctionParserUtils< is_ad >::addFParserConstants	(	SymFunctionPtr &	parser,
		const std::vector< std::string > &	constant_names,
		const std::vector< std::string > &	constant_expressions
	)		const

protectedinherited

add constants (which can be complex expressions) to the parser object

◆ evaluate() [1/2]

GenericReal< is_ad > FunctionParserUtils< is_ad >::evaluate	(	SymFunctionPtr &	,
		const std::string &	object_name = `""`
	)

protectedinherited

Evaluate FParser object and check EvalError.

◆ evaluate() [2/2]

GenericReal< is_ad > FunctionParserUtils< is_ad >::evaluate	(	SymFunctionPtr &	,
		const std::vector< GenericReal< is_ad > > &	,
		const std::string &	object_name = `""`
	)

protectedinherited

Evaluate FParser object and check EvalError.

This version uses a supplied vector of function parameters, which is useful if an object uses more than one parsed function, which may have different function parameter values.

◆ functionsOptimize() [1/3]

virtual void FunctionParserUtils< is_ad >::functionsOptimize ( SymFunctionPtr & parsed_function )

protectedvirtualinherited

run FPOptimizer on the parsed function

◆ functionsOptimize() [2/3]

void FunctionParserUtils< true >::functionsOptimize ( SymFunctionPtr & parsed_function )

protectedinherited

◆ functionsOptimize() [3/3]

void FunctionParserUtils< true >::functionsOptimize ( SymFunctionPtr & parsed_function )

protectedinherited

Definition at line 211 of file FunctionParserUtils.C.

{
  // set desired epsilon for optimization!
  auto tmp_eps = parsed_function->epsilon();
  parsed_function->setEpsilon(_epsilon);
 
  // base function
  if (!_disable_fpoptimizer)
    parsed_function->Optimize();
  if (!_enable_jit || !parsed_function->JITCompile())
    mooseError("AD parsed objects require JIT compilation to be enabled and working.");
 
  parsed_function->setEpsilon(tmp_eps);
}

◆ levelSetCorrection()

void LevelSetMeshingHelper::levelSetCorrection ( Node & node )

protected

Correct the position of a node based on the level set function.

Parameters

node	The node to be corrected

Definition at line 61 of file LevelSetMeshingHelper.C.

{
  // Based on the given level set, we try to move the node in its normal direction
  const Real diff = libMesh::TOLERANCE * 10.0; // A small value to perturb the node
  const Real original_eval = levelSetEvaluator(node);
  const Real xp_eval = levelSetEvaluator(node + Point(diff, 0.0, 0.0));
  const Real yp_eval = levelSetEvaluator(node + Point(0.0, diff, 0.0));
  const Real zp_eval = levelSetEvaluator(node + Point(0.0, 0.0, diff));
  const Real xm_eval = levelSetEvaluator(node - Point(diff, 0.0, 0.0));
  const Real ym_eval = levelSetEvaluator(node - Point(0.0, diff, 0.0));
  const Real zm_eval = levelSetEvaluator(node - Point(0.0, 0.0, diff));
  const Point grad = Point((xp_eval - xm_eval) / (2.0 * diff),
                           (yp_eval - ym_eval) / (2.0 * diff),
                           (zp_eval - zm_eval) / (2.0 * diff));
  const Real xyz_diff = -original_eval / grad.contract(grad);
  node(0) += xyz_diff * grad(0);
  node(1) += xyz_diff * grad(1);
  node(2) += xyz_diff * grad(2);
}

Referenced by Boundary2DDelaunayGenerator::General2DDelaunay(), and SurfaceSubdomainsDelaunayRemesher::General2DDelaunay().

◆ levelSetEvaluator()

Real LevelSetMeshingHelper::levelSetEvaluator ( const Point & point )

protected

Evaluate the level set function at a given point.

Parameters

point The point at which the level set function is to be evaluated

Returns: the value of the level set function at the given point

Definition at line 55 of file LevelSetMeshingHelper.C.

{
  return evaluate(_func_level_set, std::vector<Real>({point(0), point(1), point(2), 0}));
}

Referenced by Boundary2DDelaunayGenerator::General2DDelaunay(), SurfaceSubdomainsDelaunayRemesher::General2DDelaunay(), and levelSetCorrection().

◆ parsedFunctionSetup()

void FunctionParserUtils< is_ad >::parsedFunctionSetup	(	SymFunctionPtr &	function,
		const std::string &	expression,
		const std::string &	variables,
		const std::vector< std::string > &	constant_names,
		const std::vector< std::string > &	constant_expressions,
		const libMesh::Parallel::Communicator &	comm
	)		const

protectedinherited

Performs setup steps on a SymFunction.

Parameters

function	reference to pointer to the function to set up
expression	expression to parse
variables	comma separated string holding all the variables of the expression
constant_names	vector of names (symbols) of constants in the expression
constant_expressions	vectors of expressions (=values) of constants in the expression
comm	communicator used to stagger JIT file creations

◆ setParserFeatureFlags()

void FunctionParserUtils< is_ad >::setParserFeatureFlags ( SymFunctionPtr & ) const

inherited

apply input parameters to internal feature flags of the parser object

◆ validParams()

InputParameters LevelSetMeshingHelper::validParams ( )

static

Definition at line 13 of file LevelSetMeshingHelper.C.

{
  InputParameters params = FunctionParserUtils<false>::validParams();
  params.renameParameterGroup("Parsed expression advanced", "Level set shape parsed expression");
 
  params.addParam<std::string>(
      "level_set",
      "Level set used to achieve more accurate reverse projection compared to interpolation.");
  params.addParam<unsigned int>(
      "max_level_set_correction_iterations",
      3,
      "Maximum number of iterations to correct the nodes based on the level set function.");
  params.addParamNamesToGroup("level_set max_level_set_correction_iterations",
                              "Level set shape correction");
 
  return params;
}

Referenced by Boundary2DDelaunayGenerator::validParams(), and SurfaceSubdomainsDelaunayRemesher::validParams().

Member Data Documentation

◆ _disable_fpoptimizer

bool FunctionParserUtils< is_ad >::_disable_fpoptimizer

protectedinherited

Definition at line 111 of file FunctionParserUtils.h.

◆ _enable_ad_cache

bool FunctionParserUtils< is_ad >::_enable_ad_cache

protectedinherited

Definition at line 110 of file FunctionParserUtils.h.

◆ _enable_auto_optimize

bool FunctionParserUtils< is_ad >::_enable_auto_optimize

protectedinherited

Definition at line 112 of file FunctionParserUtils.h.

◆ _enable_jit

bool FunctionParserUtils< is_ad >::_enable_jit

protectedinherited

feature flags

Definition at line 109 of file FunctionParserUtils.h.

◆ _epsilon

const Real FunctionParserUtils< is_ad >::_epsilon

protectedinherited

fuzzy comparison tolerance

Definition at line 128 of file FunctionParserUtils.h.

◆ _eval_error_msg

const char* FunctionParserUtils< is_ad >::_eval_error_msg[]

staticprotectedinherited

table of FParser eval error codes

Definition at line 122 of file FunctionParserUtils.h.

◆ _evalerror_behavior

enum FunctionParserUtils::FailureMethod FunctionParserUtils< is_ad >::_evalerror_behavior

protectedinherited

◆ _func_level_set

SymFunctionPtr LevelSetMeshingHelper::_func_level_set

protected

function parser object describing the level set

Definition at line 31 of file LevelSetMeshingHelper.h.

Referenced by Boundary2DDelaunayGenerator::General2DDelaunay(), SurfaceSubdomainsDelaunayRemesher::General2DDelaunay(), levelSetEvaluator(), and LevelSetMeshingHelper().

◆ _func_params

std::vector<GenericReal<is_ad> > FunctionParserUtils< is_ad >::_func_params

protectedinherited

Array to stage the parameters passed to the functions when calling Eval.

Definition at line 125 of file FunctionParserUtils.h.

◆ _max_level_set_correction_iterations

const unsigned int LevelSetMeshingHelper::_max_level_set_correction_iterations

protected

Maximum number of iterations to correct the nodes based on the level set function.

Definition at line 28 of file LevelSetMeshingHelper.h.

Referenced by Boundary2DDelaunayGenerator::General2DDelaunay(), and SurfaceSubdomainsDelaunayRemesher::General2DDelaunay().

◆ _quiet_nan

const Real FunctionParserUtils< is_ad >::_quiet_nan

protectedinherited

appropriate not a number value to return

Definition at line 119 of file FunctionParserUtils.h.

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

include/meshgenerators/LevelSetMeshingHelper.h
src/meshgenerators/LevelSetMeshingHelper.C

Public Types

Public Member Functions

Static Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Static Protected Attributes

Detailed Description

Member Typedef Documentation

◆ SymFunction

◆ SymFunctionPtr

Member Enumeration Documentation

◆ FailureMethod

Constructor & Destructor Documentation

◆ LevelSetMeshingHelper()

Member Function Documentation

◆ addFParserConstants()

◆ evaluate() [1/2]

◆ evaluate() [2/2]

◆ functionsOptimize() [1/3]

◆ functionsOptimize() [2/3]

◆ functionsOptimize() [3/3]

◆ levelSetCorrection()

◆ levelSetEvaluator()

◆ parsedFunctionSetup()

◆ setParserFeatureFlags()

◆ validParams()

Member Data Documentation

◆ _disable_fpoptimizer

◆ _enable_ad_cache

◆ _enable_auto_optimize

◆ _enable_jit

◆ _epsilon

◆ _eval_error_msg

◆ _evalerror_behavior

◆ _func_level_set

◆ _func_params

◆ _max_level_set_correction_iterations

◆ _quiet_nan