Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #include "NonlinearTimeIntegratorInterface.h"
      11             : #include "FEProblem.h"
      12             : #include "NonlinearSystemBase.h"
      13             : 
      14             : #include "libmesh/nonlinear_implicit_system.h"
      15             : #include "libmesh/nonlinear_solver.h"
      16             : #include "libmesh/dof_map.h"
      17             : 
      18       61486 : NonlinearTimeIntegratorInterface::NonlinearTimeIntegratorInterface(FEProblemBase & problem,
      19       61486 :                                                                    SystemBase & system)
      20       61486 :   : _nl(dynamic_cast<NonlinearSystemBase *>(&system)),
      21       61486 :     _nonlinear_implicit_system(_nl ? dynamic_cast<NonlinearImplicitSystem *>(&_nl->system())
      22             :                                    : nullptr),
      23       61486 :     _Re_time(_nl ? &_nl->getResidualTimeVector() : nullptr),
      24       61486 :     _Re_non_time(_nl ? &_nl->getResidualNonTimeVector() : nullptr),
      25       61486 :     _u_dot_factor_tag(problem.addVectorTag("u_dot_factor", Moose::VECTOR_TAG_SOLUTION)),
      26      122972 :     _u_dotdot_factor_tag(problem.addVectorTag("u_dotdot_factor", Moose::VECTOR_TAG_SOLUTION))
      27             : {
      28       61486 : }
      29             : 
      30             : NumericVector<Number> *
      31        5566 : NonlinearTimeIntegratorInterface::addVector(const std::string & name,
      32             :                                             const bool project,
      33             :                                             const libMesh::ParallelType type)
      34             : {
      35        5566 :   if (_nl)
      36        2783 :     return &_nl->addVector(name, project, type);
      37             :   else
      38        2783 :     return nullptr;
      39             : }