TimeKernel.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 "TimeKernel.h"
 
// MOOSE includes
#include "Assembly.h"
#include "MooseVariableFE.h"
#include "SystemBase.h"
 
#include "libmesh/quadrature.h"
 
InputParameters
TimeKernel::validParams()
{
  InputParameters params = Kernel::validParams();
 
  params.set<MultiMooseEnum>("vector_tags") = "time";
  params.set<MultiMooseEnum>("matrix_tags") = "system time";
 
  return params;
}
 
TimeKernel::TimeKernel(const InputParameters & parameters)
  : Kernel(parameters), _u_dot(_var.uDot()), _du_dot_du(_var.duDotDu())
{
}
 
void
TimeKernel::computeResidual()
{
  prepareVectorTag(_assembly, _var.number());
 
  precalculateResidual();
  for (_i = 0; _i < _test.size(); _i++)
    for (_qp = 0; _qp < _qrule->n_points(); _qp++)
      _local_re(_i) += _JxW[_qp] * _coord[_qp] * computeQpResidual();
 
  computeResidualAdditional();
 
  accumulateTaggedLocalResidual();
 
  if (_has_save_in)
    for (unsigned int i = 0; i < _save_in.size(); i++)
      _save_in[i]->sys().solution().add_vector(_local_re, _save_in[i]->dofIndices());
}