MooseParsedFunction.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#include "MooseError.h"

// MOOSE includes
#include "InputParameters.h"
#include "MooseParsedFunction.h"
#include "MooseParsedFunctionWrapper.h"
#include "FEProblemBase.h"

registerMooseObjectAliased("MooseApp", MooseParsedFunction, "ParsedFunction");

InputParameters
MooseParsedFunction::validParams()
{
  InputParameters params = Function::validParams();
  params += MooseParsedFunctionBase::validParams();
  params.addRequiredCustomTypeParam<std::string>(
      "expression", "FunctionExpression", "The user defined function.");
  params.addClassDescription("Function created by parsing a string");
  return params;
}

MooseParsedFunction::MooseParsedFunction(const InputParameters & parameters)
  : Function(parameters),
    MooseParsedFunctionBase(parameters),
    _value(verifyFunction(this->template getRenamedParam<std::string>("value", "expression")))
{
}

Real
MooseParsedFunction::value(Real t, const Point & p) const
{
  mooseAssert(_function_ptr, "ParsedFunction should have been initialized");
  return _function_ptr->evaluate<Real>(t, p);
}

RealGradient
MooseParsedFunction::gradient(Real t, const Point & p) const
{
  mooseAssert(_function_ptr, "ParsedFunction should have been initialized");
  return _function_ptr->evaluateGradient(t, p);
}

Real
MooseParsedFunction::timeDerivative(Real t, const Point & p) const
{
  mooseAssert(_function_ptr, "ParsedFunction should have been initialized");
  return _function_ptr->evaluateDot(t, p);
}

RealVectorValue
MooseParsedFunction::vectorValue(Real /*t*/, const Point & /*p*/) const
{
  mooseError("The vectorValue method is not defined in ParsedFunction");
}

void
MooseParsedFunction::initialSetup()
{
  // Check for non-scalar variables.
  for (const auto i : index_range(_vars))
    if (_pfb_feproblem.hasVariable(_vals[i]) && !_pfb_feproblem.hasScalarVariable(_vals[i]) &&
        !_pfb_feproblem.hasFunction(_vals[i]) &&
        !_pfb_feproblem.hasPostprocessorValueByName(_vals[i]))
      mooseError("The only variables supported by ParsedFunction are scalar variables, and var '" +
                 _vals[i] + "' is not scalar.");

  if (!_function_ptr)
  {
    THREAD_ID tid = this->isParamValid("_tid") ? this->template getParam<THREAD_ID>("_tid") : 0;

    _function_ptr =
        std::make_unique<MooseParsedFunctionWrapper>(_pfb_feproblem, _value, _vars, _vals, tid);
  }
}