THMParsedFunctionWrapper.h

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
 
#pragma once
 
#include "MooseError.h"
#include "MooseTypes.h"
#include "ControlData.h"
#include "libmesh/parsed_function.h"
#include <string>
#include <vector>
 
class Simulation;
class FEProblemBase;
class Function;
 
class THMParsedFunctionWrapper
{
public:
  THMParsedFunctionWrapper(Simulation & sim,
                           FEProblemBase & feproblem,
                           const std::string & function_str,
                           const std::vector<std::string> & vars,
                           const std::vector<std::string> & vals,
                           const THREAD_ID tid = 0);
 
  Real evaluate(Real t, const Point & p);
 
  const std::vector<ControlData<Real> *> getRealControlData() { return _cd_real_vals; }
 
  const std::vector<ControlData<bool> *> getBoolControlData() { return _cd_bool_vals; }
 
private:
  Simulation & _sim;
 
  FEProblemBase & _feproblem;
 
  const std::string & _function_str;
 
  const std::vector<std::string> & _vars;
 
  const std::vector<std::string> & _vals_input;
 
  std::vector<Real> _initial_vals;
 
  std::unique_ptr<libMesh::ParsedFunction<Real>> _function_ptr;
 
  std::vector<unsigned int> _cd_real_index;
 
  std::vector<ControlData<Real> *> _cd_real_vals;
 
  std::vector<unsigned int> _cd_bool_index;
 
  std::vector<ControlData<bool> *> _cd_bool_vals;
 
  std::vector<unsigned int> _scalar_index;
 
  std::vector<const VariableValue *> _scalar_vals;
 
  std::vector<unsigned int> _function_index;
 
  std::vector<const Function *> _functions;
 
  std::vector<Real *> _addr;
 
  const THREAD_ID _tid;
 
  void initialize();
 
  void update();
 
  void updateFunctionValues(Real t, const Point & pt);
 
  void updateControlDataValues();
};