Examples of gov.fnal.eag.healpix.PixTools

gov.fnal.eag.healpix.PixTools
contains methods translated from HEALPix Fortran90 with increased map resolution in comparison to original Fortran code. @author N Kuropatkin

  }


  public void testAngDist() throws Exception {
    Vector3d v1 = new Vector3d(1.5, 1.6, 2.0);
    Vector3d v2 = new Vector3d(1.2, 1.0, 1.5);
    PixTools pt = new PixTools();
    double res1 = pt.AngDist(v1, v2);
    double res2 = v1.angle(v2);
    System.out.println("res1 = " + res1 + " res2=" + res2);
    assertEquals("angular Distance=" + res2, 1.0, res1 / res2, 1e-10);
    /* Check known problem with vecmath for small vector differences */


    Vector3d v3 = new Vector3d(1.5, 1.6, 0.);
    Vector3d v4 = new Vector3d(1.5, 1.601, 0.);
    double res3 = pt.AngDist(v3, v4);
    double res4 = v3.angle(v4);
    double expected = res4 - Math.PI / 2.;
    System.out.println("res3 = " + res3 + " res4=" + res4);
    assertEquals("angular Distance=" + res4, 1., (Math.PI / 2. + res3)
        / res4, 1e-3);
    Vector3d v5 = new Vector3d(1.5, 1.6, 0.);
    Vector3d v6 = new Vector3d(-1.5, -1.75, 0.);
    double res5 = pt.AngDist(v5, v6);
    double res6 = v5.angle(v6);
    System.out.println("res5 = " + res5 + " res6=" + res6);
    assertEquals("angular Distance=" + res6, 2.0, res5 / res6, 1e-10);
    System.out.println(" test of AngDist is done");
  }

View Full Code Here


  public void testSurfaceTriangle() throws Exception {
    Vector3d v1 = new Vector3d(1.0, 0.0, 0.0);
    Vector3d v2 = new Vector3d(0.0, 1.0, 0.0);
    Vector3d v3 = new Vector3d(0.0, 0.0, 1.0);
    PixTools pt = new PixTools();
    double res = pt.SurfaceTriangle(v1, v2, v3);
    System.out.println("Triangle surface is=" + res / Math.PI
        + " steredians");
    assertEquals("Triangle surface=" + res, 0.5, res / Math.PI, 1e-10);
    System.out.println(" test of SurfaceTriangle is done");
  }

View Full Code Here

  }


  public void testNside2Npix() {
    int nside = 1;
    int npix = 0;
    PixTools pt = new PixTools();
    npix = (int) pt.Nside2Npix(nside);
    assertEquals("Npix=" + npix, 12, npix, 1e-10);
    nside = 2;
    npix = (int) pt.Nside2Npix(nside);
    assertEquals("Npix=" + npix, 48, npix, 1e-10);
  }

View Full Code Here

  }


  public void testNpix2Nside() {
    int npix = 12;
    int nside = 0;
    PixTools pt = new PixTools();
    nside = (int) pt.Npix2Nside(npix);
    assertEquals("Nside=" + nside, 1, nside, 1e-10);


  }

View Full Code Here

  }


  public void testVec2Ang() {
    double PI = Math.PI;
    Vector3d v = new Vector3d(0.0, 1.0, 0.0);
    PixTools pt = new PixTools();
    double[] ang_tup = { 0., 0. };
    ang_tup = pt.Vect2Ang(v);
    System.out.println(" Theta=" + ang_tup[0] / PI + " Phi=" + ang_tup[1]
        / PI);
    assertEquals("Theta=" + ang_tup[0], 0.5, ang_tup[0] / PI, 1e-10);
    assertEquals("Phi=" + ang_tup[1], 0.5, ang_tup[1] / PI, 1e-10);
    v = new Vector3d(1.0, 0.0, 0.0);
    ang_tup = pt.Vect2Ang(v);
    assertEquals("phi=" + ang_tup[1], 0., ang_tup[1] / PI, 1e-10);
    System.out.println(" test Vect2Ang is done");
  }

View Full Code Here

    double PI = Math.PI;
    long pix = -1;
    double theta = PI / 2. - 0.2;
    double phi = PI / 2. ; 
    long nside = 4;
    PixTools pt = new PixTools();
    try {
      pix =  pt.ang2pix_ring(nside,theta, phi);
    } catch (Exception e) {
      e.printStackTrace();
    }
    Vector3d v = pt.Ang2Vec(theta,phi);
    long pix1 =  pt.vect2pix_ring(nside,v);
    assertEquals("pix=" + pix, pix1, pix, 1e-10);
    assertEquals("pix=" + pix, 76, pix, 1e-10);


    long pix2 =  pt.ang2pix_nest(nside,theta,phi);
    long pix3 =  pt.vect2pix_nest(nside,v);
//    assertEquals("pix2=" + pix2, pix3, pix2, 1e-10);
//    assertEquals("pix2=" + pix2, 92, pix2, 1e-10);




    double[] radec = pt.pix2ang_ring(nside,76);
    assertEquals("theta=" + theta, theta, radec[0], 4e-2);
    assertEquals("phi=" + phi, radec[1], phi, 1e-2);
    double[] radec1 = pt.pix2ang_nest(nside,92);
    System.out.println("theta="+radec1[0]+" phi="+radec1[1]);
    assertEquals("theta=" + theta, theta, radec1[0], 4e-2);
    assertEquals("phi=" + phi, radec1[1], phi, 1e-2);  
    System.out.println(" test Ang2Pix is done");
  }

View Full Code Here

  public void testAng2Vect() throws Exception {
    System.out.println(" Start test Ang2Vect----------------");
    double PI = Math.PI;
    double theta = PI / 2.;
    double phi = PI / 2;
    PixTools pt = new PixTools();
    Vector3d v = pt.Ang2Vec(theta, phi);
    System.out.println("Vector x=" + v.x + " y=" + v.y + " z=" + v.z);
    assertEquals("x=" + v.x, 0., v.x, 1e-10);
    assertEquals("y=" + v.y, 1., v.y, 1e-10);
    assertEquals("z=" + v.z, 0., v.z, 1e-10);
    System.out.println(" test Ang2Vect is done");

View Full Code Here


  public void testRingNum() throws Exception {
    double z = 0.25;
    int nside = 1;
    System.out.println("Start test RingNum !!!!!!!!!!!!!!!!!!!!");
    PixTools pt = new PixTools();
    int nring = (int) pt.RingNum(nside, z);
    System.out.println("z=" + z + " ring number =" + nring);
    assertEquals("z=" + z, 2, nring, 1e-10);
    z = -0.25;
    nring = (int) pt.RingNum(nside, z);
    assertEquals("z=" + z, 2, nring, 1e-10);
    z = 0.8;
    nring = (int) pt.RingNum(nside, z);
    assertEquals("z=" + z, 1, nring, 1e-10);
    z = -0.8;
    nring = (int) pt.RingNum(nside, z);
    assertEquals("z=" + z, 3, nring, 1e-10);
    System.out.println(" test RingNum is done");
    nside = 4;
    int pix = 3;
    Vector3d v = pt.pix2vect_ring(nside, pix);
    z = v.z;
    nring = (int) pt.RingNum(nside, z);
    assertEquals("z=" + z, 1, nring, 1e-10);
    pix = 11;
    v = pt.pix2vect_ring(nside, pix);
    z = v.z;
    nring = (int) pt.RingNum(nside, z);
    assertEquals("z=" + z, 2, nring, 1e-10);
    pix = 23;
    v = pt.pix2vect_ring(nside, pix);
    z = v.z;
    nring = (int) pt.RingNum(nside, z);
    assertEquals("z=" + z, 3, nring, 1e-10);
    pix = 39;
    v = pt.pix2vect_ring(nside, pix);
    z = v.z;
    nring = (int) pt.RingNum(nside, z);
    assertEquals("z=" + z, 4, nring, 1e-10);
    pix = 55;
    v = pt.pix2vect_ring(nside, pix);
    z = v.z;
    nring = (int) pt.RingNum(nside, z);
    assertEquals("z=" + z, 5, nring, 1e-10);
    pix = 71;
    v = pt.pix2vect_ring(nside, pix);
    z = v.z;
    nring = (int) pt.RingNum(nside, z);
    assertEquals("z=" + z, 6, nring, 1e-10);
    pix = 87;
    v = pt.pix2vect_ring(nside, pix);
    z = v.z;
    nring = (int) pt.RingNum(nside, z);
    assertEquals("z=" + z, 7, nring, 1e-10);
    pix = 103;
    v = pt.pix2vect_ring(nside, pix);
    z = v.z;
    nring = (int) pt.RingNum(nside, z);
    assertEquals("z=" + z, 8, nring, 1e-10);
    pix = 119;
    v = pt.pix2vect_ring(nside, pix);
    z = v.z;
    nring = (int) pt.RingNum(nside, z);
    assertEquals("z=" + z, 9, nring, 1e-10);
    pix = 135;
    v = pt.pix2vect_ring(nside, pix);
    z = v.z;
    nring = (int) pt.RingNum(nside, z);
    assertEquals("z=" + z, 10, nring, 1e-10);
    pix = 151;
    v = pt.pix2vect_ring(nside, pix);
    z = v.z;
    nring = (int) pt.RingNum(nside, z);
    assertEquals("z=" + z, 11, nring, 1e-10);
    pix = 167;
    v = pt.pix2vect_ring(nside, pix);
    z = v.z;
    nring = (int) pt.RingNum(nside, z);
    assertEquals("z=" + z, 12, nring, 1e-10);
    pix = 169;
    v = pt.pix2vect_ring(nside, pix);
    z = v.z;
    nring = (int) pt.RingNum(nside, z);
    assertEquals("z=" + z, 13, nring, 1e-10);
    pix = 180;
    v = pt.pix2vect_ring(nside, pix);
    z = v.z;
    nring = (int) pt.RingNum(nside, z);
    assertEquals("z=" + z, 14, nring, 1e-10);
    System.out.println("End test RingNum");
  }

View Full Code Here

    assertEquals("z=" + z, 14, nring, 1e-10);
    System.out.println("End test RingNum");
  }


  public void testNest2Ring() throws Exception {
    PixTools pt = new PixTools();
    int ipnest = 3;
    int nside = 2;
    int ipring = (int) pt.nest2ring(nside, ipnest);
    assertEquals("ipring=" + ipring, 0, ipring, 1e-10);
    ipnest = 0;
    nside = 2;
    ipring = (int) pt.nest2ring(nside, ipnest);
    assertEquals("ipring=" + ipring, 13, ipring, 1e-10);
    ipnest = 18;
    nside = 2;
    ipring = (int) pt.nest2ring(nside, ipnest);
    assertEquals("ipring=" + ipring, 27, ipring, 1e-10);
    ipnest = 23;
    nside = 2;
    ipring = (int) pt.nest2ring(nside, ipnest);
    assertEquals("ipring=" + ipring, 14, ipring, 1e-10);
    ipnest = 5;
    nside = 4;
    ipring = (int) pt.nest2ring(nside, ipnest);
    assertEquals("ipring = " + ipring, 27, ipring, 1e-10);
    System.out.println(" test Nest2Ring is done");
  }

View Full Code Here

    assertEquals("ipring = " + ipring, 27, ipring, 1e-10);
    System.out.println(" test Nest2Ring is done");
  }


  public void testRing2Nest() throws Exception {
    PixTools pt = new PixTools();
    System.out.println(" start test Ring2Nest !!!!!!!!!!!!!!!!!!!!!!");
    int ipring = 0;
    int nside = 2;


    int ipnest = (int) pt.ring2nest(nside, ipring);
    assertEquals("ipnest=" + ipnest, 3, ipnest, 1e-10);
    ipring = 13;
    nside = 2;
    ipnest = (int) pt.ring2nest(nside, ipring);
    assertEquals("ipnest=" + ipnest, 0, ipnest, 1e-10);
    ipring = 27;
    nside = 2;
    ipnest = (int) pt.ring2nest(nside, ipring);
    assertEquals("ipnest=" + ipnest, 18, ipnest, 1e-10);
    ipring = 14;
    nside = 2;
    ipnest = (int) pt.ring2nest(nside, ipring);
    assertEquals("ipnest=" + ipnest, 23, ipnest, 1e-10);
    ipring = 27;
    nside = 4;
    ipnest = (int) pt.ring2nest(nside, ipring);
    assertEquals("ipnest = " + ipnest, 5, ipnest, 1e-10);
    ipring = 83;
    nside = 4;
    ipnest = (int) pt.ring2nest(nside, ipring);
    assertEquals("ipnest = " + ipnest, 123, ipnest, 1e-10);
    System.out.println(" test Ring2Nest is done");
  }

View Full Code Here

0 1 2

TOP

Related Classes of gov.fnal.eag.healpix.PixTools

gov.fnal.eag.healpix.test.PixToolsTest

java.util.Vector

java.util.ArrayList

javax.vecmath.Vector3d

All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.