source: Sophya/trunk/Poubelle/archTOI.old/cmvplgalc.cc@ 649

#include <iostream.h>
#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <math.h>

#include "plgalcross.h"

#define PSCAL2(_v1_,_v2_) (_v1_[0]*_v2_[0]+_v1_[1]*_v2_[1]+_v1_[2]*_v2_[2])
#define TOCART(_a_,_d_,_v_) \
        {_v_[0]=cos(_d_)*cos(_a_); _v_[1]=cos(_d_)*sin(_a_); _v_[2]=sin(_d_);}

////////////////////////////////////////////////////////////////////////////////////////
int main(int narg,char *arg[])
{
{
// test FindPerpEquat
/*
cmvplgalc >! toto
cat toto | awk -v f=1e-10 '{if($1==0&&($9>f||$9<-f||$10>f||$10<-f)) print $0}'
*/
double ap,dp,v1[3],v2[3],vp[3],a1,d1,a2,d2;
int rc;
for(int aa1=0;aa1<=24;aa1+=1) for(int dd1=-90;dd1<=90;dd1+=15)
for(int aa2=0;aa2<=24;aa2+=1) for(int dd2=-90;dd2<=90;dd2+=15) {
  a1 = (double) aa1; d1 = (double) dd1;
  a2 = (double) aa2; d2 = (double) dd2;
  rc = FindPerpEquat(a1,d1,a2,d2,ap,dp);
  printf("%d   %g %g   %g %g   %g %g",rc,a1,d1,a2,d2,ap,dp);
  a1 /=  HRS_IN_RADIAN; d1 /=  DEG_IN_RADIAN;
  a2 /=  HRS_IN_RADIAN; d2 /=  DEG_IN_RADIAN;
  ap /=  HRS_IN_RADIAN; dp /=  DEG_IN_RADIAN;
  TOCART(a1,d1,v1)
  TOCART(a2,d2,v2)
  TOCART(ap,dp,vp)
  printf("   %g   %g %g\n"
        ,PSCAL2(vp,vp),PSCAL2(v1,vp),PSCAL2(v2,vp));
}
}

exit(0);
}