ADDiracKernel.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 "DiracKernelBase.h"
#include "MooseVariableInterface.h"
#include "ADFunctorInterface.h"

class ADDiracKernel : public DiracKernelBase,
                      public MooseVariableInterface<Real>,
                      public ADFunctorInterface
{
public:
  static InputParameters validParams();

  ADDiracKernel(const InputParameters & parameters);

  virtual void jacobianSetup() override;

  virtual void computeResidual() override;
  virtual void computeJacobian() override;
  virtual void computeOffDiagJacobian(unsigned int jvar) override;
  virtual void computeResidualAndJacobian() override;

  virtual const MooseVariableField<Real> & variable() const override { return _var; }

protected:
  virtual ADReal computeQpResidual() = 0;

  MooseVariableField<Real> & _var;

  const ADTemplateVariablePhiValue<Real> & _phi;

  const ADTemplateVariableTestValue<Real> & _test;

  const ADTemplateVariableValue<Real> & _u;
  const ADTemplateVariableGradient<Real> & _grad_u;

private:
  void computeADResiduals();
  void computeFullJacobian();

  std::vector<ADReal> _ad_residuals;
  const Elem * _last_jacobian_elem;
};