KokkosCoupledVarNeumannBC.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 "KokkosIntegratedBCValue.h"
 
class KokkosCoupledVarNeumannBC : public Moose::Kokkos::IntegratedBCValue
{
public:
  static InputParameters validParams();
 
  KokkosCoupledVarNeumannBC(const InputParameters & parameters);
 
  template <typename Derived>
  KOKKOS_FUNCTION Real precomputeQpResidual(const unsigned int qp, AssemblyDatum & datum) const;
  template <typename Derived>
  KOKKOS_FUNCTION Real precomputeQpOffDiagJacobian(const unsigned int j,
                                                   const unsigned int jvar,
                                                   const unsigned int qp,
                                                   AssemblyDatum & datum) const;
 
protected:
  const Moose::Kokkos::VariableValue _coupled_var;
 
  const unsigned int _coupled_num;
 
  const Real _coef;
 
  const Moose::Kokkos::VariableValue _scale_factor;
};
 
template <typename Derived>
KOKKOS_FUNCTION Real
KokkosCoupledVarNeumannBC::precomputeQpResidual(const unsigned int qp, AssemblyDatum & datum) const
{
  return -_scale_factor(datum, qp) * _coef * _coupled_var(datum, qp);
}
 
template <typename Derived>
KOKKOS_FUNCTION Real
KokkosCoupledVarNeumannBC::precomputeQpOffDiagJacobian(const unsigned int j,
                                                       const unsigned int jvar,
                                                       const unsigned int qp,
                                                       AssemblyDatum & datum) const
{
  if (jvar == _coupled_num)
    return -_scale_factor(datum, qp) * _coef * _phi(datum, j, qp);
  else
    return 0;
}