elementaryop.h

Go to the documentation of this file.

/*
 MCL
 Copyright (c) 2012-18, Enzo De Sena
 All rights reserved.
 
 Authors: Enzo De Sena, enzodesena@gmail.com
 */

#ifndef MCL_ELEMENTARYOP_H
#define MCL_ELEMENTARYOP_H

#include "mcltypes.h"

namespace mcl {
  
template<class T>
T Max(const T& scalar_a, const T& scalar_b) noexcept {
  if (scalar_a >= scalar_b) { return scalar_a; }
  else { return scalar_b; }
}
  
template<class T>
T Min(const T& scalar_a, const T& scalar_b) noexcept {
  if (scalar_a < scalar_b) { return scalar_a; }
  else { return scalar_b; }
}

Real Rem(const Real& scalar_a, const Real& scalar_b);

Real Mod(const Real& scalar_a, const Real& scalar_b);

Int Mod(const Int& scalar_a, const Int& scalar_b);

Int Fix(const Real scalar);

Real Abs(Real input);
  
Real Abs(Complex input);
  
Real Pow(Real input, Real exponent);
  
Real Sqrt(Real input);
  
inline Int RoundToInt(Real input) {
  Int output = static_cast<int>(input);
  output += (input-output >= 0.5) - (input-output <= -0.5);
  return output;
}
    
Int Sign(const Real scalar);
  
Complex Conj(Complex scalar);

Real RealPart(Complex scalar);
  
Real ImagPart(Complex scalar);
  
Int NextPow2(Real input);
  
Int Next2(Int input);
  
double StringToDouble(const std::string& s);
  
Int Factorial(const Int input);
  
Real LinearInterpolation(Real x1, Real y1, Real x2, Real y2, Real x);
  
bool IsReal(const std::vector<Complex>& input);
  
Real Entropy(std::vector<Real> pdf, Real base);
  
bool ElementaryOpTest();
  
#if MCL_LOAD_BOOST

Real AssociatedLegendreP(Int n, Int m, Real x);
  
Complex SphericalHarmonic(Int n, Int m, Real theta, Real phi);
#endif
  
  
} 
#endif