import java.util.Vector;

public class HVatlas extends PointLayout {

    public String getFrame(int which){
		if(which == 1) return "RA---";
		else           return "DEC--";
	}

	public Vector getPoints (String atlasParam){
		double dec, dech, c, ra, rainc;
		Vector3d v;
		Vector oldSet, newSet;
		int level, lev, d, p, q, m, n, i;

		try {
			level = Integer.parseInt(atlasParam);
        } catch (NumberFormatException e) {
            System.out.println("Cannot parse Level from "
                + atlasParam);
            level = 4;
            System.out.println("Using Level 4");
        }

		newSet = new Vector();
		newSet.addElement(new Vector3d(0., 0., 1.));
		newSet.addElement(new Vector3d(0., 1., 0.));
		newSet.addElement(new Vector3d(1., 0., 0.));

		for (lev=0; lev<level; lev++){
			oldSet = newSet;
			newSet = new Vector();

			d = 1<<lev;
			n = 0;
			for (p=0; p<=d; p++){
				m = n;
				n += p;
				if(p > 0){
					for(q=0; q<p; q++){
						avg(oldSet, n+q, m+q, newSet);
						avg(oldSet, n+q+1, m+q, newSet);
//						System.out.println((n+q) + "--" + (m+q));
//						System.out.println((n+q+1) + "--" + (m+q));
					}
				}
				for(q=0; q<=p; q++){
					avg(oldSet, n+q, n+q, newSet);
//					System.out.println(q + " & " + (n+q) + "--" + (n+q));
					if(q != p){
						avg(oldSet, n+q, n+q+1, newSet);
//						System.out.println((n+q) + "--" + (n+q+1));
					}
				}
			}
		}

		double Epsilon = 1.0E-8;

		oldSet = newSet;
		newSet = new Vector();
		for (i=0; i<oldSet.size(); i++){
			v =(Vector3d) oldSet.elementAt(i);
			newSet.addElement(v);
			if(Math.abs(v.x()) > Epsilon) 
				newSet.addElement(new Vector3d(-v.x(),  v.y(),  v.z()));
		}

		oldSet = newSet;
		newSet = new Vector();
		for (i=0; i<oldSet.size(); i++){
			v =(Vector3d) oldSet.elementAt(i);
			newSet.addElement(v);
			if(Math.abs(v.y()) > Epsilon) 
				newSet.addElement(new Vector3d( v.x(), -v.y(),  v.z()));
		}

		oldSet = newSet;
		newSet = new Vector();
		for (i=0; i<oldSet.size(); i++){
			v =(Vector3d) oldSet.elementAt(i);
			newSet.addElement(v);
			if(Math.abs(v.z()) > Epsilon) 
				newSet.addElement(new Vector3d( v.x(),  v.y(), -v.z()));
		}
		return newSet;
	}

	private void avg(Vector oldSet, int i, int j, Vector newSet){
		Vector3d vi = (Vector3d)oldSet.elementAt(i);
		Vector3d vj = (Vector3d)oldSet.elementAt(j);
		Vector3d v = vi.add(vj);
		v.normalize();
		newSet.addElement(v);
	}
}