FastMath.cc

Go to the documentation of this file.

#include "DataFormats/Math/interface/FastMath.h"
namespace fastmath_details {
  float atanbuf_[257 * 2];
  double datanbuf_[513 * 2];
 
  namespace {
    // ====================================================================
    // arctan initialization
    // =====================================================================
    struct Initatan {
      Initatan() {
        unsigned int ind;
        for (ind = 0; ind <= 256; ind++) {
          double v = ind / 256.0;
          double asinv = ::asin(v);
          atanbuf_[ind * 2] = ::cos(asinv);
          atanbuf_[ind * 2 + 1] = asinv;
        }
        for (ind = 0; ind <= 512; ind++) {
          double v = ind / 512.0;
          double asinv = ::asin(v);
          datanbuf_[ind * 2] = ::cos(asinv);
          datanbuf_[ind * 2 + 1] = asinv;
        }
      }
    };
    Initatan initAtan;
  }  // namespace
}  // namespace fastmath_details