/*
GeoGebra - Dynamic Mathematics for Everyone
http://www.geogebra.org
This file is part of GeoGebra.
This program is free software; you can redistribute it and/or modify it
under the terms of the GNU General Public License as published by
the Free Software Foundation.
*/
/*
* AlgoPolarLine.java
*
* Created on 30. August 2001, 21:37
*/
package org.geogebra.common.geogebra3D.kernel3D.algos;
import org.geogebra.common.geogebra3D.kernel3D.geos.GeoLine3D;
import org.geogebra.common.kernel.Construction;
import org.geogebra.common.kernel.Kernel;
import org.geogebra.common.kernel.Matrix.Coords;
import org.geogebra.common.kernel.algos.AlgoPolarLineND;
import org.geogebra.common.kernel.geos.GeoLine;
import org.geogebra.common.kernel.kernelND.GeoConicND;
import org.geogebra.common.kernel.kernelND.GeoLineND;
import org.geogebra.common.kernel.kernelND.GeoPointND;
/**
*
* @author mathieu
*/
public class AlgoPolarLine3D extends AlgoPolarLineND {
private GeoLine polar2D;
private Coords coords2D;
private double[] polarCoords;
private Coords polarOrigin, polarDirection;
/**
* Creates new AlgoPolarLine
*
* @param cons
* construction
* @param label
* label
* @param c
* conic
* @param P
* polar point
*/
public AlgoPolarLine3D(Construction cons, String label, GeoConicND c,
GeoPointND P) {
super(cons, label, c, P);
}
@Override
protected GeoLineND newGeoLine(Construction cons1) {
// we also need this
polar2D = new GeoLine(cons1);
polarCoords = new double[3];
// create 3D line
return new GeoLine3D(cons1);
}
// calc polar line of P relativ to c
@Override
public final void compute() {
// check if point lies on conic coord sys
coords2D = c.getCoordSys()
.getNormalProjection(P.getInhomCoordsInD3())[1];
if (!Kernel.isZero(coords2D.getZ())) {
polar.setUndefined();
return;
}
// now it's a 2D point in coord sys
coords2D.setZ(1);
// update polar line in conic coord sys
c.polarLine(coords2D, polar2D);
// update 3D polar
polar2D.getCoords(polarCoords);
polarDirection = c.getCoordSys().getVector(-polarCoords[1],
polarCoords[0]);
if (Kernel.isZero(polarCoords[0])) {
polarOrigin = c.getCoordSys().getPoint(0,
-polarCoords[2] / polarCoords[1]);
} else {
polarOrigin = c.getCoordSys()
.getPoint(-polarCoords[2] / polarCoords[0], 0);
}
((GeoLine3D) polar).setCoord(polarOrigin, polarDirection);
}
}