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

    assertEquals("ipnest = " + ipnest, 123, ipnest, 1e-10);
    System.out.println(" test Ring2Nest is done");
  }


  public void testNext_In_Line_Nest() throws Exception {
    PixTools pt = new PixTools();
    int ipix = 0;
    int nside = 2;
    int ipnext = (int) pt.next_in_line_nest(nside, ipix);
    assertEquals("ipnext=" + ipnext, 23, ipnext, 1e-10);
    ipix = 1;
    nside = 2;
    ipnext = (int) pt.next_in_line_nest(nside, ipix);
    assertEquals("ipnext=" + ipnext, 6, ipnext, 1e-10);
    ipix = 4;
    nside = 2;
    ipnext = (int) pt.next_in_line_nest(nside, ipix);
    assertEquals("ipnext=" + ipnext, 27, ipnext, 1e-10);
    ipix = 27;
    nside = 2;
    ipnext = (int) pt.next_in_line_nest(nside, ipix);
    assertEquals("ipnext=" + ipnext, 8, ipnext, 1e-10);
    ipix = 12;
    nside = 2;
    ipnext = (int) pt.next_in_line_nest(nside, ipix);
    assertEquals("ipnext=" + ipnext, 19, ipnext, 1e-10);
    ipix = 118;
    nside = 4;
    ipnext = (int) pt.next_in_line_nest(nside, ipix);
    assertEquals("ipnext = " + ipnext, 117, ipnext, 1e-10);
    System.out.println(" test next_in_line_nest is done");
  }

View Full Code Here

    System.out.println(" test next_in_line_nest is done");
  }


  public void testInRing() throws Exception {
    System.out.println(" Start test InRing !!!!!!!!!!!!!!!!!!!!!!!!!");
    PixTools pt = new PixTools();
    int[] nestComp = { 19, 0, 23, 4, 27, 8, 31, 12 };
    double PI = Math.PI;
    int nside = 2;
    boolean nest = false;
    int iz = 3;
    double phi = PI;
    double dphi = PI;
    ArrayList ring = pt.InRing(nside, iz, phi, dphi, nest);
    for (int i = 0; i < ring.size(); i++) {
      assertEquals("ipnext = " + ((Long) ring.get(i)).longValue(),
          i + 12, ((Long) ring.get(i)).longValue(), 1e-10);
    }
    nest = true;
    ring = pt.InRing(nside, iz, phi, dphi, nest);
    for (int i = 0; i < ring.size(); i++) {
      assertEquals("ipnext = " + ((Long) ring.get(i)).longValue(),
          nestComp[i], ((Long) ring.get(i)).longValue(), 1e-10);
    }
    nest = false;
    nside = 4;
    phi = 2.1598449493429825;
    iz = 8;
    dphi = 0.5890486225480867;
    //    System.out.println(" iz="+iz+" phi="+phi+" dphi="+dphi);
    ring = pt.InRing(nside, iz, phi, dphi, nest);
    //    for (int i = 0; i<ring.size(); i++) {
    //      System.out.println("ipnext = "+ ((Integer)ring.get(i)).intValue());
    //    }
    nest = false;
    nside = 4;
    dphi = 0. * PI;
    iz = 8;
    phi = 2.1598449493429825;
    //    System.out.println(" iz="+iz+" phi="+phi+" dphi="+dphi);
    ring = pt.InRing(nside, iz, phi, dphi, nest);
    //    for (int i = 0; i<ring.size(); i++) {
    //      System.out.println("ipnext = "+ ((Integer)ring.get(i)).intValue());
    //    }
    System.out.println(" test InRing is done");
  }

View Full Code Here

    System.out.println(" test InRing is done");
  }


  public void testNeighbours_Nest() throws Exception {
    System.out.println(" Start test Neighbours_Nest !!!!!!!!!!!!!!!!!");
    PixTools pt = new PixTools();
    long nside = 2;
    long ipix = 25;
    long[] test17 = { 34, 16, 18, 19, 2, 0, 22, 35 };
    long[] test32 = { 40, 44, 45, 34, 35, 33, 38, 36 };
    long[] test3 = { 0, 2, 13, 15, 11, 7, 6, 1 };
    long[] test25 = { 42, 24, 26, 27, 10, 8, 30, 43 };
    //
    ArrayList npixList = new ArrayList();
    npixList = pt.neighbours_nest(nside, ipix);
    for (int i = 0; i < npixList.size(); i++) {
      assertEquals("ip = " + ((Long) npixList.get( i)).longValue(),
          test25[ i], ((Long) npixList.get(  i)).longValue(), 1e-10);
    }
    ipix = 17;


    npixList = pt.neighbours_nest(nside, ipix);
    for (int i = 0; i < npixList.size(); i++) {
      assertEquals("ip = " + ((Long) npixList.get(i)).longValue(),
          test17[i], ((Long) npixList.get(i)).longValue(), 1e-10);
    }
    ipix = 32;


    npixList = pt.neighbours_nest(nside, ipix);
    for (int i = 0; i < npixList.size(); i++) {
      assertEquals("ip = " + ((Long) npixList.get(i)).longValue(),
          test32[i], ((Long) npixList.get(i)).longValue(), 1e-10);
    }
    ipix = 3;


    npixList = pt.neighbours_nest(nside, ipix);
    for (int i = 0; i < npixList.size(); i++) {
      assertEquals("ip = " + ((Long) npixList.get(i)).longValue(),
          test3[i], ((Long) npixList.get(i)).longValue(), 1e-10);
    }
    System.out.println(" test NeighboursNest is done");

View Full Code Here

    System.out.println(" test NeighboursNest is done");
  }


  public void testIntrs_Intrv() throws Exception {
    System.out.println(" test intrs_intrv !!!!!!!!!!!!!!!!!!!!!!!!!!!");
    PixTools pt = new PixTools();
    double[] d1 = { 1.0, 9.0 };
    double[] d2 = { 3.0, 16.0 };
    double[] di;
    //    System.out.println("Case "+d1[0]+" "+d1[1]+" | "+d2[0]+" "+d2[1]);
    di = pt.intrs_intrv(d1, d2);
    //    System.out.println("Result "+di[0]+" - "+di[1]);
    int n12 = di.length / 2;
    assertEquals("n12 = " + n12, 1, n12, 1e-6);
    assertEquals("di[0] = " + di[0], 3.0, di[0], 1e-6);
    assertEquals("di[1] = " + di[1], 9.0, di[1], 1e-6);
    d1 = new double[] { 0.537, 4.356 };
    d2 = new double[] { 3.356, 0.8 };
    //    System.out.println("Case "+d1[0]+" "+d1[1]+" | "+d2[0]+" "+d2[1]);
    di = pt.intrs_intrv(d1, d2);
    n12 = di.length / 2;
    assertEquals("n12 = " + n12, 2, n12, 1e-6);
    assertEquals("di[0] = " + di[0], 0.537, di[0], 1e-6);
    assertEquals("di[1] = " + di[1], 0.8, di[1], 1e-6);
    assertEquals("di[2] = " + di[2], 3.356, di[2], 1e-6);
    assertEquals("di[1] = " + di[3], 4.356, di[3], 1e-6);


    d1 = new double[] { 2.356194490092345, 2.356194490292345 };
    d2 = new double[] { 1.251567, 4.17 };
    //    System.out.println("Case "+d1[0]+" "+d1[1]+" | "+d2[0]+" "+d2[1]);
    di = pt.intrs_intrv(d1, d2);
    n12 = di.length / 2;
    assertEquals("n12 = " + n12, 1, n12, 1e-6);
    assertEquals("di[0] = " + di[0], 2.35619449009, di[0], 1e-6);
    assertEquals("di[1] = " + di[1], 2.35619449029, di[1], 1e-6);

View Full Code Here

  public void testPix2Vect_ring() throws Exception {
    System.out.println("Start test Pix2Vect_ring !!!!!!!!!!!!!!!!!!!");
    double TWOPI = 2.0 * Math.PI;
    int nside = 2;
    int ipix = 0;
    PixTools pt = new PixTools();
    Vector3d v1 = new Vector3d(0., 0., 0.);
    v1 = pt.pix2vect_ring(nside, ipix);
    assertEquals("v1.z = " + v1.z, 1.0, v1.z, 1e-1);


    ipix = 20;
    Vector3d v2 = new Vector3d(0., 0., 0.);
    v2 = pt.pix2vect_ring(nside, ipix);
    assertEquals("v2.x = " + v2.x, 1.0, v2.x, 1e-1);
    assertEquals("v2.z = " + v2.z, 0.0, v2.z, 1e-1);
    ipix = 22;
    Vector3d v3 = new Vector3d();
    v3 = pt.pix2vect_ring(nside, ipix);
    assertEquals("v3.y = " + v3.y, 1.0, v3.y, 1e-1);
    assertEquals("v3.z = " + v3.z, 0.0, v3.z, 1e-1);
    //    System.out.println("Vector3 x="+v3.x+" y="+v3.y+" z="+v3.z);
    ipix = 95;
    nside = 4;
    v1 = pt.pix2vect_ring(nside, ipix);
    v1.normalize();
    double phi1 = Math.atan2(v1.y, v1.x);
    double[] tetphi = new double[2];
    tetphi = pt.pix2ang_ring(nside, ipix);
    assertEquals("phi = " + phi1, 0.0, Math.abs(phi1 - tetphi[1]), 1e-10);
    ipix = 26;
    nside = 4;
    v1 = pt.pix2vect_ring(nside, ipix);
    v1.normalize();
    phi1 = Math.atan2(v1.y, v1.x);
    if (phi1 < 0.)
      phi1 += TWOPI;
    tetphi = new double[2];
    tetphi = pt.pix2ang_ring(nside, ipix);
    assertEquals("phi = " + phi1, 0.0, Math.abs(phi1 - tetphi[1]), 1e-10);
    System.out.println("------------------------------------------");
    System.out.println(" test pix2vect_ring is done");
  }

View Full Code Here

  public void testPix2Vect_nest() throws Exception {
    double TWOPI = 2.0 * Math.PI;
    int nside = 2;
    int ipix = 3;
    System.out.println(" Start test Pix2Vect_nest !!!!!!!!!!!!!!");
    PixTools pt = new PixTools();
    Vector3d v1 = new Vector3d(0., 0., 0.);
    v1 = pt.pix2vect_nest(nside, ipix);
    assertEquals("v1.z = " + v1.z, 1.0, v1.z, 1e-1);
    ipix = 17;
    Vector3d v2 = new Vector3d(0., 0., 0.);
    v2 = pt.pix2vect_nest(nside, ipix);
    assertEquals("v2.x = " + v2.x, 1.0, v2.x, 1e-1);
    assertEquals("v2.z = " + v2.z, 0.0, v2.z, 1e-1);


    ipix = 21;
    Vector3d v3 = new Vector3d();
    v3 = pt.pix2vect_nest(nside, ipix);
    assertEquals("v3.y = " + v3.y, 1.0, v3.y, 1e-1);
    assertEquals("v3.z = " + v3.z, 0.0, v3.z, 1e-1);
    nside = 4;
    ipix = 105;
    v1 = pt.pix2vect_nest(nside, ipix);
    v1.normalize();
    double phi1 = Math.atan2(v1.y, v1.x);
    if (phi1 < 0.)
      phi1 += TWOPI;
    double[] tetphi = new double[2];
    tetphi = pt.pix2ang_nest(nside, ipix);
    assertEquals("phi = " + phi1, 0.0, Math.abs(phi1 - tetphi[1]), 1e-10);
    nside = 4;
    ipix = 74;
    v1 = pt.pix2vect_nest(nside, ipix);
    v1.normalize();
    phi1 = Math.atan2(v1.y, v1.x);
    if (phi1 < 0.)
      phi1 += TWOPI;
    tetphi = new double[2];
    tetphi = pt.pix2ang_nest(nside, ipix);
    assertEquals("phi = " + phi1, 0.0, Math.abs(phi1 - tetphi[1]), 1e-10);


    System.out.println(" test pix2vect_nest is done");
    System.out.println("-------------------------------------------");
  }

View Full Code Here

    System.out.println("Start test Vect2Pix_ring !!!!!!!!!!!!!!!!!!!");
    double TWOPI = 2.0 * Math.PI;
    int nside = 4;
    int ipix = 83;
    int respix = 0;
    PixTools pt = new PixTools();
    Vector3d v1 = new Vector3d(0., 0., 0.);
    v1 = pt.pix2vect_ring(nside, ipix);
    respix = (int) pt.vect2pix_ring(nside, v1);
    assertEquals("respix = " + respix, 83, respix, 1e-10);


    System.out.println("------------------------------------------");
    System.out.println(" test vect2pix_ring is done");
  }

View Full Code Here

    System.out.println("Start test Vect2Pix_nest !!!!!!!!!!!!!!!!!!!");
    double TWOPI = 2.0 * Math.PI;
    int nside = 4;
    int ipix = 83;
    int respix = 0;
    PixTools pt = new PixTools();
    Vector3d v1 = new Vector3d(0., 0., 0.);
    v1 = pt.pix2vect_ring(nside, ipix);
    respix = (int) pt.vect2pix_nest(nside, v1);
    assertEquals("respix = " + respix, 123, respix, 1e-10);


    System.out.println("------------------------------------------");
    System.out.println(" test vect2pix_nest is done");
  }

View Full Code Here

    System.out.println(" test vect2pix_nest is done");
  }


  public void testQuery_Strip() throws Exception {
    System.out.println(" Start test query Strip !!!!!!!!!!!!!!!!");
    PixTools pt = new PixTools();
    int[] pixel1 = { 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15,
        16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31,
        32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47,
        48, 49, 50, 51, 52, 53, 54, 55 };
    int nside = 4;
    int nest = 0;
    double theta1 = 0.0;
    double theta2 = Math.PI / 4.0 + 0.2;
    ArrayList pixlist;
    pixlist = pt.query_strip(nside, theta1, theta2, nest);
    int nlist = pixlist.size();
    for (int i = 0; i < nlist; i++) {
      long ipix = ((Long) pixlist.get(i)).longValue();
      assertEquals("pixel = " + ipix, pixel1[i], ipix, 1e-10);
    }

View Full Code Here


  }


  public void testQuery_Disc() throws Exception {
    System.out.println(" Start test query_disc !!!!!!!!!!!!!!!!!!!!!");
    PixTools pt = new PixTools();
    long nside = 4;
    int nest = 0;
    long ipix = 0;
    int[] pixel1 = { 45, 46, 60, 61, 62, 76, 77, 78, 79, 91, 92, 93, 94,
        95, 108, 109, 110, 111, 124, 125, 126, 141, 142 };
    int[] pixel2 = { 24, 19, 93, 18, 17, 92, 87, 16, 107, 89, 86, 85, 106,
        105, 83, 84, 159, 104, 81, 158, 157, 155, 156 };
    int inclusive = 1;
    double radius = Math.PI / 8.0;
    Vector3d v = pt.pix2vect_ring(nside, 93);
    ArrayList pixlist;
    pixlist = pt.query_disc(nside, v, radius, nest, inclusive);


    int nlist = pixlist.size();


    for (int i = 0; i < nlist; i++) {
      ipix = ((Long) pixlist.get(i)).longValue();
      assertEquals("pixel = " + ipix, pixel1[i], ipix, 1e-10);
    }
/*    
    v = pt.pix2vect_ring(nside, 103);
    
    pixlist = pt.query_disc(nside, v, radius, nest, inclusive);


    nlist = pixlist.size();


    for (int i = 0; i < nlist; i++) {
      ipix = ((Integer) pixlist.get(i)).intValue();
//      assertEquals("pixel = " + ipix, pixel1[i], ipix, 1e-10);
      System.out.println(" i="+i+" pixel="+ipix);
    }
*/  
    //
    v = pt.pix2vect_ring(nside, 93);
    nest = 1;
    inclusive = 1;
    pixlist = pt.query_disc(nside, v, radius, nest, inclusive);
    nlist = pixlist.size();


    for (int i = 0; i < nlist; i++) {
      ipix = ((Long) pixlist.get(i)).longValue();
//            System.out.println("ipix="+ ipix);

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.