basicMath.h

Go to the documentation of this file.

// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
 
 
#ifndef MATHUTILS_INCLUDE_MATHUTILS_DETAIL_BASICMATH_H_
#define MATHUTILS_INCLUDE_MATHUTILS_DETAIL_BASICMATH_H_
 
#include "GPUCommonDef.h"
#include "GPUCommonMath.h"
#include "CommonConstants/MathConstants.h"
 
#ifndef GPUCA_GPUCODE_DEVICE
#include <cmath>
#include <tuple>
#include <array>
#endif
 
namespace o2
{
namespace math_utils
{
namespace detail
{
 
template <typename T>
GPUhdi() constexpr T copysign(T x, T y)
{
  return o2::gpu::GPUCommonMath::Copysign(x, y);
}
 
template <class T>
GPUhdi() T min(const T x, const T y)
{
  return o2::gpu::GPUCommonMath::Min(x, y);
};
 
template <class T>
GPUhdi() T max(const T x, const T y)
{
  return o2::gpu::GPUCommonMath::Max(x, y);
};
 
template <class T>
GPUhdi() T sqrt(T x)
{
  return o2::gpu::GPUCommonMath::Sqrt(x);
};
 
template <class T>
GPUhdi() T abs(T x)
{
  return o2::gpu::GPUCommonMath::Abs(x);
};
 
template <class T>
GPUdi() int nint(T x)
{
  return o2::gpu::GPUCommonMath::Float2IntRn(x);
};
 
template <class T>
GPUdi() bool finite(T x)
{
  return o2::gpu::GPUCommonMath::Finite(x);
}
 
GPUdi() unsigned int clz(unsigned int val)
{
  return o2::gpu::GPUCommonMath::Clz(val);
};
 
GPUdi() unsigned int popcount(unsigned int val)
{
  return o2::gpu::GPUCommonMath::Popcount(val);
};
 
template <class T>
GPUdi() T log(T x)
{
  return o2::gpu::GPUCommonMath::Log(x);
};
 
template <class T>
GPUdi() T exp(T x)
{
  return o2::gpu::GPUCommonMath::Exp(x);
};
 
#if !defined(GPUCA_GPUCODE_DEVICE) || defined(__CUDACC__) || defined(__HIPCC__) // Double specialisations, not for OpenCL and Metal
template <>
GPUhdi() double copysign(double x, double y)
{
  return o2::gpu::GPUCommonMath::Copysignd(x, y);
}
template <>
GPUhdi() double sqrt(double x)
{
  return o2::gpu::GPUCommonMath::Sqrtd(x);
};
template <>
GPUdi() double log(double x)
{
  return o2::gpu::GPUCommonMath::Logd(x);
};
template <>
GPUdi() double exp(double x)
{
  return o2::gpu::GPUCommonMath::Expd(x);
};
#endif
 
#ifndef GPUCA_GPUCODE_DEVICE // Double overloads using std namespace not visible on GPU
GPUhdi() double abs(double x)
{
  return std::abs(x);
};
template <>
GPUdi() int nint(double x)
{
  return std::nearbyint(x);
};
template <>
GPUdi() bool finite(double x)
{
#ifdef __FAST_MATH__
  return false;
#else
  return std::isfinite(x);
#endif
}
#endif
 
} // namespace detail
} // namespace math_utils
} // namespace o2
 
#endif /* MATHUTILS_INCLUDE_MATHUTILS_DETAIL_BASICMATH_H_ */