ArrayKernel.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 "KernelBase.h"
#include "MooseVariableInterface.h"
#include "MooseVariableScalar.h"

class ArrayKernel : public KernelBase, public MooseVariableInterface<RealEigenVector>
{
public:
  static InputParameters validParams();

  ArrayKernel(const InputParameters & parameters);

  virtual void computeResidual() override;

  virtual void computeJacobian() override;

  virtual void computeOffDiagJacobian(unsigned int jvar) override;

  virtual void computeOffDiagJacobianScalar(unsigned int jvar) override;

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

protected:
  virtual void computeQpResidual(RealEigenVector & residual) = 0;

  virtual RealEigenVector computeQpJacobian();

  virtual RealEigenMatrix computeQpOffDiagJacobian(const MooseVariableFEBase & jvar);

  virtual RealEigenMatrix computeQpOffDiagJacobianScalar(const MooseVariableScalar & jvar);

  virtual void initQpResidual() {}

  virtual void initQpJacobian() {}

  virtual void initQpOffDiagJacobian(const MooseVariableFEBase & jvar)
  {
    if (jvar.number() == _var.number())
      initQpJacobian();
  }

  ArrayMooseVariable & _var;

  const ArrayVariableTestValue & _test;

  const ArrayVariableTestGradient & _grad_test;
  const MappedArrayVariablePhiGradient & _array_grad_test;

  const ArrayVariablePhiValue & _phi;

  const ArrayVariablePhiGradient & _grad_phi;

  const ArrayVariableValue & _u;

  const ArrayVariableGradient & _grad_u;

  const unsigned int _count;

private:
  RealEigenVector _work_vector;

  RealEigenMatrix _work_matrix;
};