ArrayDGKernel.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 "DGKernelBase.h"

#include "MooseVariableScalar.h"

class ArrayDGKernel : public DGKernelBase, public NeighborMooseVariableInterface<RealEigenVector>
{
public:
  static InputParameters validParams();

  ArrayDGKernel(const InputParameters & parameters);

  virtual const MooseVariableFEBase & variable() const override { return _var; }

  virtual void computeOffDiagJacobian(unsigned int jvar) override;

  virtual void computeElemNeighResidual(Moose::DGResidualType type) override;

  virtual void computeElemNeighJacobian(Moose::DGJacobianType type) override;

  virtual void computeOffDiagElemNeighJacobian(Moose::DGJacobianType type,
                                               const MooseVariableFEBase & jvar) override;

protected:
  virtual void computeQpResidual(Moose::DGResidualType type, RealEigenVector & residual) = 0;

  virtual RealEigenVector computeQpJacobian(Moose::DGJacobianType);

  virtual RealEigenMatrix computeQpOffDiagJacobian(Moose::DGJacobianType type,
                                                   const MooseVariableFEBase & jvar);

  virtual void initQpResidual(Moose::DGResidualType) {}

  virtual void initQpJacobian(Moose::DGJacobianType) {}

  virtual void initQpOffDiagJacobian(Moose::DGJacobianType, const MooseVariableFEBase &) {}

  ArrayMooseVariable & _var;
  const ArrayVariableValue & _u;
  const ArrayVariableGradient & _grad_u;
  const ArrayVariablePhiValue & _phi;
  const ArrayVariablePhiGradient & _grad_phi;
  const MappedArrayVariablePhiGradient & _array_grad_phi;
  const ArrayVariableTestValue & _test;
  const ArrayVariableTestGradient & _grad_test;
  const MappedArrayVariablePhiGradient & _array_grad_test;
  const ArrayVariablePhiValue & _phi_neighbor;
  const ArrayVariablePhiGradient & _grad_phi_neighbor;
  const MappedArrayVariablePhiGradient & _array_grad_phi_neighbor;
  const ArrayVariableTestValue & _test_neighbor;
  const ArrayVariableTestGradient & _grad_test_neighbor;
  const MappedArrayVariablePhiGradient & _array_grad_test_neighbor;
  const ArrayVariableValue & _u_neighbor;
  const ArrayVariableGradient & _grad_u_neighbor;
  const std::vector<Eigen::Map<RealDIMValue>> & _array_normals;
  const unsigned int _count;

private:
  RealEigenVector _work_vector;
};