Tricrystal2CircleGrainsICAction.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* 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 "Tricrystal2CircleGrainsICAction.h"
#include "Factory.h"
#include "Parser.h"
#include "FEProblem.h"
#include "Conversion.h"
 
#include <sstream>
#include <stdexcept>
 
#include "libmesh/libmesh.h"
#include "libmesh/exodusII_io.h"
#include "libmesh/equation_systems.h"
#include "libmesh/nonlinear_implicit_system.h"
#include "libmesh/explicit_system.h"
#include "libmesh/string_to_enum.h"
 
const Real Tricrystal2CircleGrainsICAction::_abs_zero_tol = 1e-12;
 
registerMooseAction("PhaseFieldApp", Tricrystal2CircleGrainsICAction, "add_ic");
 
template <>
InputParameters
validParams<Tricrystal2CircleGrainsICAction>()
{
  InputParameters params = validParams<Action>();
  params.addRequiredParam<unsigned int>("op_num", "number of order parameters to create");
  params.addRequiredParam<std::string>("var_name_base", "specifies the base name of the variables");
 
  return params;
}
 
Tricrystal2CircleGrainsICAction::Tricrystal2CircleGrainsICAction(const InputParameters & params)
  : Action(params),
    _var_name_base(getParam<std::string>("var_name_base")),
    _op_num(getParam<unsigned int>("op_num"))
{
}
 
void
Tricrystal2CircleGrainsICAction::act()
{
#ifdef DEBUG
  Moose::err << "Inside the Tricrystal2CircleGrainsICAction Object\n";
#endif
 
  // Loop through the number of order parameters
  for (unsigned int op = 0; op < _op_num; op++)
  {
    // Create variable names
    std::string var_name = _var_name_base;
    std::stringstream out;
    out << op;
    var_name.append(out.str());
 
    // Set parameters for BoundingBoxIC
    InputParameters poly_params = _factory.getValidParams("Tricrystal2CircleGrainsIC");
    poly_params.set<VariableName>("variable") = var_name;
    poly_params.set<unsigned int>("op_num") = _op_num;
    poly_params.set<unsigned int>("op_index") = op;
 
    // Add initial condition
    _problem->addInitialCondition("Tricrystal2CircleGrainsIC",
                                  "Tricrystal2CircleGrainsIC_" + Moose::stringify(op),
                                  poly_params);
  }
}