Contents:
1 2 3 4 5 6 7 8
double sawsin(double x) { double t = fmod(x/(2*M_PI),(double)1.0); if (t>0.5) return -sin(x); if (t<=0.5) return (double)2.0*t-1.0; }