KokkosNodalSum.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 "KokkosNodalVariablePostprocessor.h"

class KokkosNodalSum : public KokkosNodalVariablePostprocessor
{
public:
  static InputParameters validParams();

  KokkosNodalSum(const InputParameters & parameters);

  virtual void initialize() override;
  virtual void finalize() override;
  virtual Real getValue() const override;

  template <typename Derived>
  KOKKOS_FUNCTION void reduce(Datum & datum, Real * result) const;

  KOKKOS_FUNCTION Real computeValue(const unsigned int qp, Datum & datum) const
  {
    return _u(datum, qp);
  }

  template <typename Derived>
  KOKKOS_FUNCTION void join(Real * result, const Real * source) const;
  template <typename Derived>
  KOKKOS_FUNCTION void init(Real * result) const;
};

template <typename Derived>
KOKKOS_FUNCTION void
KokkosNodalSum::reduce(Datum & datum, Real * result) const
{
  if (datum.isNodalDefined(_u.variable()))
    result[0] += static_cast<const Derived *>(this)->computeValue(0, datum);
}

template <typename Derived>
KOKKOS_FUNCTION void
KokkosNodalSum::join(Real * result, const Real * source) const
{
  result[0] += source[0];
}

template <typename Derived>
KOKKOS_FUNCTION void
KokkosNodalSum::init(Real * result) const
{
  result[0] = 0;
}