#ifndef physics_constants
#define physics_constants

#include "Double.h"
#include <complex>
const std::complex<LDouble_t> I_(0,1);
const LDouble_t PI_=2*atan2(1,0);
#endif