PorousFlowLumpedKernelBase.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 "GenericKernel.h"
 
template <bool is_ad>
class PorousFlowLumpedKernelBaseTempl : public GenericKernel<is_ad>
{
public:
  PorousFlowLumpedKernelBaseTempl(const InputParameters & parameters);
 
protected:
  virtual void computeJacobian() override;
  virtual void computeOffDiagJacobian(unsigned int jvar) override;
  virtual void jacobianSetup() override;
 
  usingGenericKernelMembers;
 
private:
  // Cache to avoid recomputing AD residuals for each off-diagonal jvar on the same element,
  // mirroring the _my_elem guard in ADKernel::computeOffDiagJacobian.
  const Elem * _my_elem_lma = nullptr;
};
 
typedef PorousFlowLumpedKernelBaseTempl<false> PorousFlowLumpedKernelBase;
typedef PorousFlowLumpedKernelBaseTempl<true> ADPorousFlowLumpedKernelBase;