scalarwrapper.hh

Go to the documentation of this file.

// SPDX-FileCopyrightText: 2021-2024 The Ikarus Developers mueller@ibb.uni-stuttgart.de
// SPDX-License-Identifier: LGPL-3.0-or-later
 
#pragma once
 
#include <dune/common/referencehelper.hh>
 
template <typename T>
struct ScalarWrapper
{
  ScalarWrapper(T val)
      : val{val} {}
 
  using RawScalarType = Dune::ResolveRef_t<T>;
 
  operator RawScalarType&() const { return Dune::resolveRef(val); }
 
  const RawScalarType& value() const { return Dune::resolveRef(val); }
 
  RawScalarType& value() { return Dune::resolveRef(val); }
 
  RawScalarType operator+(const ScalarWrapper& v) const {
    return RawScalarType{Dune::resolveRef(val) + Dune::resolveRef(v.val)};
  }
  RawScalarType operator-(const ScalarWrapper& v) const {
    return RawScalarType{Dune::resolveRef(val) - Dune::resolveRef(v.val)};
  }
 
  RawScalarType operator-() const { return RawScalarType{-Dune::resolveRef(val)}; }
 
  RawScalarType operator*(RawScalarType value) const { return RawScalarType{val * value}; }
 
  ScalarWrapper& operator+=(const ScalarWrapper& v) {
    Dune::resolveRef(val) += Dune::resolveRef(v.val);
    return *this;
  }
 
  ScalarWrapper& operator-=(const ScalarWrapper& v) {
    Dune::resolveRef(val) -= Dune::resolveRef(v.val);
    return *this;
  }
 
  ScalarWrapper& operator*=(RawScalarType v) {
    Dune::resolveRef(val) *= v;
    return *this;
  }
 
  ScalarWrapper& operator/=(RawScalarType v) {
    Dune::resolveRef(val) /= v;
    return *this;
  }
 
private:
  T val;
 
  friend RawScalarType operator*(RawScalarType f, const ScalarWrapper& v) {
    return RawScalarType{f * Dune::resolveRef(v.val)};
  }
};