sincos.h

Go to the documentation of this file.

// This file's extension implies that it's C, but it's really -*- C++ -*-.
/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
 
#ifndef CXXUTILS_SINCOS_H
#define CXXUTILS_SINCOS_H
 
 
#include <cmath>
 
 
namespace CxxUtils {
 

struct sincos
{
  sincos (double x)
#if defined(__USE_GNU)
  // Version using the GNU sincos() function.
  { ::sincos(x, &sn, &cs); }
#else
  // Generic version.
  : sn (std::sin (x)), cs (std::cos (x)) {}
#endif

  double sn;

  double cs;

  double apply (double a, double b) const { return a*sn + b*cs; }

  double apply2 (double a, double b, double c) const
  { return a*sn*sn + b*sn*cs + c*cs*cs; }
};
 
 
} // namespace CxxUtils
 
 
#endif //not CXXUTILS_SINCOS_H