KokkosMatchedValueBC.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 "KokkosNodalBC.h"

class KokkosMatchedValueBC : public Moose::Kokkos::NodalBC
{
public:
  static InputParameters validParams();

  KokkosMatchedValueBC(const InputParameters & parameters);

  template <typename Derived>
  KOKKOS_FUNCTION Real computeQpResidual(const unsigned int qp, AssemblyDatum & datum) const;
  template <typename Derived>
  KOKKOS_FUNCTION Real computeQpOffDiagJacobian(const unsigned int jvar,
                                                const unsigned int qp,
                                                AssemblyDatum & datum) const;

protected:
  const Moose::Kokkos::VariableValue _v;

  const unsigned int _v_num;
  const Real _u_coeff;
  const Real _v_coeff;
};

template <typename Derived>
KOKKOS_FUNCTION Real
KokkosMatchedValueBC::computeQpResidual(const unsigned int qp, AssemblyDatum & datum) const
{
  return _u_coeff * _u(datum, qp) - _v_coeff * _v(datum, qp);
}

template <typename Derived>
KOKKOS_FUNCTION Real
KokkosMatchedValueBC::computeQpOffDiagJacobian(const unsigned int jvar,
                                               const unsigned int /* qp */,
                                               AssemblyDatum & /* datum */) const
{
  if (jvar == _v_num)
    return -_v_coeff;
  else
    return 0;
}