#include <iostream.h>
#include "sambainit.h"
#include "tod.h"
#include "spheregorski.h"


int main() 
{
  double teta,phi;
  PeidaInit();
  int  m=32;
  SphereGorski<double> sph(m);
  for (int j=0;j<sph.NbPixels();j++) 
    {
      sph.PixThetaPhi(j,teta,phi);
      sph(j)= 0.2* cos(3.*teta)*sin(8*phi);
    }
  // Projection  dans une image

  cout <<  " ===== Fin de test0 ======== " << endl;
  return 0;
}