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();
};