// GeometricAlgebra Java package file
   // Current version
   // Go to ... http://sinai.mech.fukui-u.ac.jp/gcj/software/KamiWaAi/index.html
   // Copyright (C) 2003, Eckhard M.S. Hitzer
   //
   // This library is free software; you can redistribute it and/or
   // modify it under the terms of the GNU Lesser General Public
   // License as published by the Free Software Foundation; either
   // version 2.1 of the License, or any later version.
  //
  // This library is distributed in the hope that it will be useful,
  // but WITHOUT ANY WARRANTY; without even the implied warranty of
  // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
  // Lesser General Public License for more details.
  //
  // You should have received a copy of the GNU Lesser General Public
  // License along with this library; if not, write to the Free Software
  // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
  // Go to ... http://www.gnu.org/copyleft/lesser.html
  // 
  // hitzer@mech.fukui-u.ac.jp
  // 
  // Department of Mechanical Engineering, Fukui University
  // 3-9-1 Bunkyo, 910-8507 Fukui, Japan
  // 
  //
  //
  // It the personal wish of E.M.S. Hitzer,
  // that you never use this software for infringing human dignity:
  //
  // "God created man in his own image,
  // in the image of God he created him;
  // male and female he created them.
  // [The Bible, Genesis chapter 1 verse 27]
  //
  //
  //
  // 
  // File PointC.java
  // 2003-3-4
  //
  // 
  //
  
  package net.sourceforge.kamiwaai.geometricalgebra;
  
  import java.awt.Color;
  import java.awt.Graphics;
  import java.awt.Graphics2D;
  import java.awt.geom.Ellipse2D;
  
  public class PointC extends GeometricObject {
  	// The default color of points will be blue in the front view
  	// and green in the right side view.
  	private Color color = Color.blue;
  
  	private Color rightcolor = Color.green;
  
  	private int pointsize = 10;
  
  	public PointC() {
  		L = nbar; // This sets the point at the origin.
  	}
  
  	public PointC(double e1c, double e2c, double e3c) {
  		L = nbar;
  		//    	MultiVector Tc1, Tc2, Tc3, Tcr1, Tcr2, Tcr3;
  		//		Tc1 = One.add(n.mult(e1.multSc(e1c - L.ScProd(e1))).multSc(0.5));
  		//		Tcr1 = Tc1.reverse();
  		//		Tc2 = One.add(n.mult(e2.multSc(e2c - L.ScProd(e2))).multSc(0.5));
  		//		Tcr2 = Tc2.reverse();
  		//		Tc3 = One.add(n.mult(e3.multSc(e3c - L.ScProd(e3))).multSc(0.5));
  		//		Tcr3 = Tc3.reverse();
  		//		L = Tc1.mult(Tc2).mult(Tc3).mult(L).mult(Tcr3).mult(Tcr2).mult(Tcr1);
  		MultiVector Tc = One.add(n.mult(e1.multSc(e1c - L.ScProd(e1))).multSc(0.5)).add(
  				n.mult(e2.multSc(e2c - L.ScProd(e2))).multSc(0.5)).add(
  				n.mult(e3.multSc(e3c - L.ScProd(e3))).multSc(0.5));
  		L = Tc.mult(L).mult(Tc.reverse());
  	}
  
  	public PointC(MultiVector L1) {
  		L = L1;
  	}
  
  	public void sete1(double e1c) {
  		MultiVector Tc, Tcr;
  		Tc = One.add(n.mult(e1.multSc(e1c - L.ScProd(e1))).multSc(0.5));
  		Tcr = Tc.reverse();
  		L = Tc.mult(L).mult(Tcr);
  	}
  
  	public void sete2(double e2c) {
  		MultiVector Tc, Tcr;
  		Tc = One.add(n.mult(e2.multSc(e2c - L.ScProd(e2))).multSc(0.5));
  		Tcr = Tc.reverse();
  		L = Tc.mult(L).mult(Tcr);
  	}
  
  	public void sete3(double e3c) {
 		MultiVector Tc, Tcr;
 		Tc = One.add(n.mult(e3.multSc(e3c - L.ScProd(e3))).multSc(0.5));
 		Tcr = Tc.reverse();
 		L = Tc.mult(L).mult(Tcr);
 	}
 
 	public void setColor(Color c) {
 		color = c;
 	}
 
 	public void setRightColor(Color rc) {
 		rightcolor = rc;
 	}
 
 	public void setSize(int s) {
 		pointsize = s;
 	}
 
 	double[] coord() {
 		double[] lps = new double[3];
 		lps[0] = L.ScProd(e1);
 		lps[1] = L.ScProd(e2);
 		lps[2] = L.ScProd(e3);
 		return lps;
 	}
 
 	public void draw(Graphics gg) {
 		Graphics2D g = (Graphics2D) gg;
 		// Draw the disk in front panel
 		Ellipse2D.Double disk = new Ellipse2D.Double(L.ScProd(e1) - pointsize
 				/ 2, L.ScProd(e2) - pointsize / 2, pointsize, pointsize);
 		g.setColor(color);
 		g.fill(disk);
 		// and drawing an outer black ring around the disk
 		g.setColor(Color.black);
 		g.draw(disk);
 	}
 
 	public void erase(Graphics gg) {
 		Graphics2D g = (Graphics2D) gg;
 		// Erase the disk in front panel
 		Ellipse2D.Double disk = new Ellipse2D.Double(L.ScProd(e1) - pointsize
 				/ 2, L.ScProd(e2) - pointsize / 2, pointsize, pointsize);
 		g.setColor(g.getBackground());
 		g.fill(disk);
 		// and erase the outer black ring around the disk
 		g.draw(disk);
 	}
 
 	public void eraseZ(Graphics gg) {
 		Graphics2D g = (Graphics2D) gg;
 		// Erase the disk in front panel
 		Ellipse2D.Double disk = new Ellipse2D.Double(L.ScProd(e3) - pointsize
 				/ 2, L.ScProd(e2) - pointsize / 2, pointsize, pointsize);
 		g.setColor(g.getBackground());
 		g.fill(disk);
 		// and erase outer black ring around the disk
 		g.draw(disk);
 	}
 
 	public void drawZ(Graphics gg) {
 		Graphics2D g = (Graphics2D) gg;
 		// Draw the disk in front panel
 		Ellipse2D.Double disk = new Ellipse2D.Double(L.ScProd(e3) - pointsize
 				/ 2, L.ScProd(e2) - pointsize / 2, pointsize, pointsize);
 		g.setColor(rightcolor);
 		g.fill(disk);
 		// and drawing an outer black ring around the disk
 		g.setColor(Color.black);
 		g.draw(disk);
 	}
 
 	double ppdistance(PointC pc2) {
 		return Math.sqrt(-2.0 * L.ScProd(pc2.getMultiVector()));
 	}
 
 	public boolean isXYpoint(double x, double y) {
 		boolean checkpoint = false;
 		double ddxy2 = (L.ScProd(e1) - x) * (L.ScProd(e1) - x)
 				+ (L.ScProd(e2) - y) * (L.ScProd(e2) - y);
 		if (ddxy2 <= pointsize / 2 * pointsize / 2)
 			checkpoint = true;
 		return checkpoint;
 	}
 
 	public boolean isYZpoint(double y, double z) {
 		boolean checkpoint = false;
 		double ddyz2 = (L.ScProd(e2) - y) * (L.ScProd(e2) - y)
 				+ (L.ScProd(e3) - z) * (L.ScProd(e3) - z);
 		if (ddyz2 <= pointsize / 2 * pointsize / 2)
 			checkpoint = true;
 		return checkpoint;
 	}
 
 }