KokkosMatCoupledForce.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 "KokkosKernelValue.h"
#include "KokkosMap.h"

class KokkosMatCoupledForce : public Moose::Kokkos::KernelValue
{
public:
  static InputParameters validParams();

  KokkosMatCoupledForce(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 j,
                                                const unsigned int jvar,
                                                const unsigned int qp,
                                                AssemblyDatum & datum) const;

private:
  const unsigned int _n_coupled;
  const bool _coupled_props;
  std::vector<unsigned int> _v_var;
  const Moose::Kokkos::VariableValue _v;
  Moose::Kokkos::Array<Real> _coef;
  Moose::Kokkos::Map<unsigned int, unsigned int> _v_var_to_index;
  Moose::Kokkos::Array<Moose::Kokkos::MaterialProperty<Real>> _mat_props;
};

template <typename Derived>
KOKKOS_FUNCTION Real
KokkosMatCoupledForce::computeQpResidual(const unsigned int qp, AssemblyDatum & datum) const
{
  Real r = 0;
  if (_coupled_props)
    for (unsigned int j = 0; j < _n_coupled; ++j)
      r += -_coef[j] * _mat_props[j](datum, qp) * _v(datum, qp, j);
  else
    for (unsigned int j = 0; j < _n_coupled; ++j)
      r += -_coef[j] * _v(datum, qp, j);
  return r;
}

template <typename Derived>
KOKKOS_FUNCTION Real
KokkosMatCoupledForce::computeQpOffDiagJacobian(const unsigned int j,
                                                const unsigned int jvar,
                                                const unsigned int qp,
                                                AssemblyDatum & datum) const
{
  if (!_v_var_to_index.exists(jvar))
    return 0;

  unsigned int p = _v_var_to_index[jvar];

  if (_coupled_props)
    return -_coef[p] * _mat_props[p](datum, qp) * _phi(datum, j, qp);
  return -_coef[p] * _phi(datum, j, qp);
}