/******************************** fmod_m90_p90 *******************************/ /* Fold input angle in the [-90,+90[ domain. */ double fmod_m90_p90(double angle) { return angle>0.0? fmod(angle+90.0,180.0)-90.0 : fmod(angle-90.0,180.0)+90.0; }