use constant PI => 4 * atan2(1, 1);