#include "rqmc.hpp" double frac_part(double x){ return x - floor(x); }