source: Sophya/trunk/SophyaExt/XephemAstroLib/sphcart.c@ 2631

#include <math.h>
#include <stdio.h>

#include "astro.h"

/* transformation from spherical to cartesian coordinates */
void
sphcart (
double l, double b, double r,           /* source: spherical coordinates */
double *x, double *y, double *z)        /* result: rectangular coordinates */
{
        double rcb = r * cos(b);

        *x = rcb * cos(l);
        *y = rcb * sin(l);
        *z = r * sin(b);
}

/* transformation from cartesian to spherical coordinates */
void
cartsph (
double x, double y, double z,           /* source: rectangular coordinates */
double *l, double *b, double *r)        /* result: spherical coordinates */
{
        double rho = x*x + y*y;

        if (rho > 1e-35) {      /* standard case: off axis */
            *l = atan2(y, x);
            range (l, 2*PI);
            *b = atan2(z, sqrt(rho));
            *r = sqrt(rho + z*z);
        } else {                /* degenerate case; avoid math error */
            *l = 0.0;
            if (z == 0.0)
                *b = 0.0;
            else
                *b = (z > 0.0) ? PI/2. : -PI/2.;
            *r = fabs(z);
        }
}

/* For RCS Only -- Do Not Edit */
static char *rcsid[2] = {(char *)rcsid, "@(#) $RCSfile: sphcart.c,v $ $Date: 2004-06-15 16:52:40 $ $Revision: 1.3 $ $Name: not supported by cvs2svn $"};