```--- a/math.c
+++ b/math.c
@@ -601,11 +601,6 @@
lambda = x_dest/distance;
phi = y_dest/distance;

-/*
-    printf("To compute original x %10.4f  y %10.4f distance %10.4f\n", x_dest, y_dest, distance);
-    printf("To compute x %10.4f  y %10.4f\n", lambda, phi);
-    printf("To compute cos(y) %10.4f  cos(x) %10.4f\n", cos(phi), cos(lambda));
-*/
if (abs(cos(phi) * cos(lambda) + 1.0) <= EPSLN) {
*x_src = distance * 2 ;
*y_src = 0;
@@ -614,8 +609,6 @@

k1 = sqrt(2.0 / (1 + cos(phi) * cos(lambda)));

-    printf("To compute k %10.4f  \n", k1);
-
*x_src = distance * k1 * cos(phi) * sin (lambda);
*y_src = distance * k1 * sin(phi);

@@ -630,11 +623,10 @@

x = x_dest/distance;
y = y_dest/distance;
-    //x = x_dest;
-    //y = y_dest;

assert(! isnan(x));
assert(! isnan(y));
+
if (fabs(x) > PI || fabs(y) > PI) {
*y_src = 0;
*x_src = 0;
@@ -645,9 +637,6 @@

if (fabs(ro) <= EPSLN)
{
-      printf("** compute x %10.4f  y %10.4f\n", x, y);
-      printf("** compute c %10.4f\n", ro);
-
*y_src = 0;
*x_src = 0;
return 1;
@@ -659,8 +648,6 @@

if (fabs(ro * cos(c)) <= EPSLN ) {
-      printf("** compute x %10.4f  y %10.4f\n", x, y);
-      printf("** compute cos(y) %10.4f  cos(x) %10.4f\n", cos(c), c);
*x_src = 0;
return 1;
}
```