LevelSetReinitializationMultiApp.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 "LevelSetReinitializationMultiApp.h"
#include "LevelSetReinitializationProblem.h"

#include "Executioner.h"

// libMesh
#include "libmesh/mesh_tools.h"

registerMooseObject("LevelSetApp", LevelSetReinitializationMultiApp);

InputParameters
LevelSetReinitializationMultiApp::validParams()
{
  InputParameters params = MultiApp::validParams();
  params.addClassDescription(
      "MultiApp capable of performing repeated complete solves for level set reinitialization.");
  params.addParam<unsigned int>(
      "interval", 1, "Time step interval when to perform reinitialization.");

  params.suppressParameter<std::vector<Point>>("positions");
  params.suppressParameter<std::vector<FileName>>("positions_file");
  params.suppressParameter<bool>("output_in_position");
  params.suppressParameter<std::vector<Real>>("reset_time");
  params.suppressParameter<std::vector<unsigned int>>("reset_apps");
  params.suppressParameter<Real>("move_time");
  params.suppressParameter<std::vector<unsigned int>>("move_apps");
  params.suppressParameter<std::vector<Point>>("move_positions");

  return params;
}

LevelSetReinitializationMultiApp::LevelSetReinitializationMultiApp(
    const InputParameters & parameters)
  : MultiApp(parameters), _level_set_problem(NULL), _interval(getParam<unsigned int>("interval"))
{
}

void
LevelSetReinitializationMultiApp::initialSetup()
{
  MultiApp::initialSetup();

  if (_has_an_app)
  {
    Executioner * ex = _apps[0]->getExecutioner();

    if (!ex)
      mooseError("Executioner does not exist!");

    ex->init();

    _executioner = ex;

    _level_set_problem = dynamic_cast<LevelSetReinitializationProblem *>(&appProblemBase(0));
    if (!_level_set_problem)
      mooseError("The Problem type must be LevelSetReinitializationProblem.");
  }
}

bool
LevelSetReinitializationMultiApp::solveStep(Real /*dt*/,
                                            Real /*target_time*/,
                                            bool /*auto_advance*/)
{
  // Do nothing if not on interval
  if ((_fe_problem.timeStep() % _interval) != 0)
    return true;

  if (!_has_an_app)
    return true;

  _console << "Solving Reinitialization problem." << std::endl;

  int rank;
  int ierr;
  ierr = MPI_Comm_rank(_orig_comm, &rank);
  mooseCheckMPIErr(ierr);

  bool last_solve_converged = true;

  _level_set_problem->resetTime();
  _executioner->execute();
  if (!_executioner->lastSolveConverged())
    last_solve_converged = false;

  return last_solve_converged;
}