package org.simantics.plant3d.scenegraph.controlpoint;

import java.util.ArrayList;
import java.util.Arrays;
import java.util.Collection;
import java.util.List;

import javax.vecmath.AxisAngle4d;
import javax.vecmath.Matrix3d;
import javax.vecmath.Point3d;
import javax.vecmath.Quat4d;
import javax.vecmath.Tuple3d;
import javax.vecmath.Vector3d;

import org.simantics.g3d.math.MathTools;
import org.simantics.g3d.property.annotations.GetPropertyValue;
import org.simantics.g3d.scenegraph.G3DNode;
import org.simantics.plant3d.scenegraph.IP3DNode;
import org.simantics.plant3d.scenegraph.Nozzle;
import org.simantics.plant3d.scenegraph.P3DRootNode;
import org.simantics.plant3d.scenegraph.PipeRun;
import org.simantics.plant3d.scenegraph.PipelineComponent;

import vtk.vtkRenderer;


public class PipeControlPoint extends G3DNode implements IP3DNode {

	private static boolean DEBUG = false;

	public enum PointType{INLINE,TURN,END};
	public enum Direction{NEXT,PREVIOUS};
	public enum PositionType {SPLIT,NEXT,PREVIOUS,PORT}

	private PipelineComponent component;

	private PointType type;
	private boolean isFixed = true;
	private boolean isRotate = false;
	private boolean isReverse = false;
	private boolean isDeletable = true;
	private boolean isSizeChange = false;
	private boolean isSub = false;

	public PipeControlPoint(PipelineComponent component) {
		this.component = component;
		if (component.getPipeRun() != null)
			component.getPipeRun().addChild(this);

	}

	public PipeControlPoint(PipelineComponent component, PipeRun piperun) {
		this.component = component;
		piperun.addChild(this);
	}

	@Override
	public void update(vtkRenderer ren) {
		try {
			PipingRules.requestUpdate(this);
		} catch (Exception e) {
			e.printStackTrace();
		}

	}

	public PipeRun getPipeRun() {
		return (PipeRun)getParent();
	}

	public PipelineComponent getPipelineComponent() {
		return component;
	}

	public PointType getType() {
		return type;
	}

	public void setType(PointType type) {
		this.type = type;
	}

	@GetPropertyValue(name="Fixed",tabId="Debug",value="fixed")
	public boolean isFixed() {
		return isFixed;
	}


	public void setFixed(boolean fixed) {
		this.isFixed = fixed;
	}

	@GetPropertyValue(name="Rotate",tabId="Debug",value="rotate")
	public boolean isRotate() {
		return isRotate;
	}

	public void setRotate(boolean rotate) {
		this.isRotate = rotate;
	}

	@GetPropertyValue(name="Reverse",tabId="Debug",value="reverse")
	public boolean isReverse() {
		return isReverse;
	}

	public void setReverse(boolean reverse) {
		this.isReverse = reverse;
	}

	public void setSub(boolean sub) {
		this.isSub = sub;
	}

	@GetPropertyValue(name="Deletable",tabId="Debug",value="deletable")
	public boolean isDeletable() {
		return isDeletable;
	}

	public void setDeletable(boolean deletable) {
		this.isDeletable = deletable;
	}

	public boolean isPathLegEnd() {
		return type != PointType.INLINE;
	}

	public boolean isEnd() {
		return type == PointType.END;
	}

	public boolean isTurn() {
		return type == PointType.TURN;
	}

	public boolean isInline() {
		return type == PointType.INLINE;
	}

	public boolean isDirected() {
		return isFixed && isEnd();
	}

	public boolean isNonDirected() {
		return !isFixed && isEnd();
	}

	public boolean isVariableLength() {
		return !isFixed && isInline();
	}

	public boolean isVariableAngle() {
		return !isFixed && isTurn();
	}

	public boolean isBranchEnd() {
		return isDeletable && isEnd();
	}

	public boolean isOffset() {
		return offset != null;
	}

	public boolean isDualSub() {
		return parent != null && isSub;
	}

	public boolean isDualInline() {
		return children.size() == 1 && children.get(0).isDualSub();
	}

	public boolean isAxial() {
		return isInline() && !isDualInline();
	}

	public boolean isSizeChange() {
		return isSizeChange;
		//		if (children.size() == 0)
		//			return false;
		//		if (!isDualInline())
		//			return false;
		//		return getPipeRun() != children.get(0).getPipeRun();
	}

	public void setSizeChange(boolean isSizeChange) {
		this.isSizeChange = isSizeChange;
	}


	private PipeControlPoint next;
	private PipeControlPoint previous;

	public PipeControlPoint getNext() {
		return next;
	}

	public PipeControlPoint getPrevious() {
		return previous;
	}

	public void setNext(PipeControlPoint next) {
		if (isEnd() && previous != null && next != null)
			throw new RuntimeException("End control points are allowed to have only one connection");
		if (next == this)
			throw new RuntimeException("Cannot connect to self");
		if (this.next == next)
			return;
		if (DEBUG) System.out.println(this + " next " + next);
		this.next = next;
		if (component != null) {
			if (parent == null || isSub)
				component.setNext(next != null ? next.component : null);
			else
				component.setBranch0(next != null ? next.component : null);
			updateSubPoint();
		}

	}

	public void setPrevious(PipeControlPoint previous) {
		if (isEnd() && next != null && previous != null)
			throw new RuntimeException("End control points are allowed to have only one connection");
		if (previous == this)
			throw new RuntimeException("Cannot connect to self");
		if (this.previous == previous)
			return;
		if (DEBUG) System.out.println(this + " previous " + previous);
		this.previous = previous;
		if (component != null) {
			if (parent == null || isSub)
				component.setPrevious(previous != null ? previous.component : null);
			else
				component.setBranch0(previous != null ? previous.component : null);
			updateSubPoint();
		}

	}

	public PipeControlPoint parent;
	public List<PipeControlPoint> children = new ArrayList<PipeControlPoint>();

	public List<PipeControlPoint> getSubPoint() {
		return children;
	}

	public PipeControlPoint getParentPoint() {
		return parent;
	}


	private double length;
	private Double turnAngle;
	private Vector3d turnAxis;

	private Double offset;
	private Double rotationAngle;
	private Boolean reversed;

	@GetPropertyValue(name="Length",tabId="Debug",value="length")
	public double getLength() {
		return length;
	}

	public void setLength(double l) {
		if (Double.isInfinite(l) || Double.isNaN(l)) {
			return;
		}
		if (Math.abs(this.length-l) < MathTools.NEAR_ZERO)
			return;
		this.length = l;
		firePropertyChanged("length");
		if (isDualInline())
			getSubPoint().get(0).setLength(l);
	}

	@GetPropertyValue(name="Turn Angle",tabId="Debug",value="turnAngle")
	public Double getTurnAngle() {
		return turnAngle;
	}

	@GetPropertyValue(name="Turn Axis",tabId="Debug",value="turnAxis")
	public Vector3d getTurnAxis() {
		return turnAxis;
	}

	@GetPropertyValue(name="Offset",tabId="Debug",value="offset")
	public Double getOffset() {
		return offset;
	}

	@GetPropertyValue(name="Rotation Angle",tabId="Debug",value="rotationAngle")
	public Double getRotationAngle() {
		return rotationAngle;
	}

	@GetPropertyValue(name="Reversed",tabId="Debug",value="reversed")
	public Boolean getReversed() {
		return reversed;
	}

	public boolean _getReversed() {
		if (reversed == null)
			return false;
		return reversed;
	}

	public void setTurnAngle(Double turnAngle) {
		if (turnAngle == null || Double.isInfinite(turnAngle) || Double.isNaN(turnAngle)) {
			return;
		}
		if (this.turnAngle != null && Math.abs(this.turnAngle-turnAngle) < MathTools.NEAR_ZERO)
			return;
		this.turnAngle = turnAngle;
		firePropertyChanged("turnAngle");
	}

	public void setTurnAxis(Vector3d turnAxis) {
		if (this.turnAxis != null && MathTools.equals(turnAxis, this.turnAxis))
			return;
		this.turnAxis = turnAxis;
		firePropertyChanged("turnAxis");
	}

	public void setOffset(Double offset) {
		if (Double.isInfinite(offset) || Double.isNaN(offset)) {
			return;
		}
		if (this.offset != null && Math.abs(this.offset-offset) < MathTools.NEAR_ZERO)
			return;
		this.offset = offset;
		firePropertyChanged("offset");
	}

	public void setRotationAngle(Double rotationAngle) {
		if (Double.isInfinite(rotationAngle) || Double.isNaN(rotationAngle)) {
			return;
		}
		if (this.rotationAngle != null && Math.abs(this.rotationAngle-rotationAngle) < MathTools.NEAR_ZERO)
			return;
		this.rotationAngle = rotationAngle;
		firePropertyChanged("rotationAngle");
	}

	public void setReversed(Boolean reversed) {
		this.reversed = reversed;
		firePropertyChanged("reversed");
	}

	public Vector3d getSizeChangeOffsetVector(Vector3d dir) {
		Quat4d q;
		if (rotationAngle == null)
			q = getControlPointOrientationQuat(dir, 0.0);
		else
			q = getControlPointOrientationQuat(dir, rotationAngle);
		Vector3d v = new Vector3d(0.0,offset,0.0);
		Vector3d offset = new Vector3d();
		MathTools.rotate(q, v, offset);
		return offset;
	}

	public Vector3d getSizeChangeOffsetVector() {
		Quat4d q;
		if (rotationAngle == null)
			q = getControlPointOrientationQuat(0.0);
		else
			q = getControlPointOrientationQuat(rotationAngle);
		Vector3d v = new Vector3d(0.0,offset,0.0);
		Vector3d offset = new Vector3d();
		MathTools.rotate(q, v, offset);
		return offset;
	}

	@GetPropertyValue(name="Next",tabId="Debug",value="next")
	private String getNextString() {
		if (next == null)
			return null;
		return next.toString();
	}

	@GetPropertyValue(name="Previous",tabId="Debug",value="previous")
	private String getPrevString() {
		if (previous == null)
			return null;
		return previous.toString();
	}

	@GetPropertyValue(name="Sub",tabId="Debug",value="sub")
	private String getSubString() {
		if (children.size() == 0)
			return "";
		return Arrays.toString(children.toArray());
	}

	@GetPropertyValue(name="Type",tabId="Debug",value="type")
	public String getTypeString() {
		return type.name();
	}

	public Quat4d getControlPointOrientationQuat(double angle) {

		if (turnAxis == null) {
			Vector3d dir = getPathLegDirection(Direction.NEXT);
			if (dir.lengthSquared() > MathTools.NEAR_ZERO)
				dir.normalize();
			return getControlPointOrientationQuat(dir, angle);
		} else {
			Vector3d dir = getPathLegDirection(Direction.PREVIOUS);
			dir.negate();
			if (dir.lengthSquared() > MathTools.NEAR_ZERO)
				dir.normalize();
			return getControlPointOrientationQuat(dir, turnAxis, angle);
		}
	}

	public Quat4d getControlPointOrientationQuat(double angle, boolean reversed) {

		if (turnAxis == null) {
			Vector3d dir = getPathLegDirection(Direction.NEXT);
			if (dir.lengthSquared() > MathTools.NEAR_ZERO)
				dir.normalize();
			Quat4d q =  getControlPointOrientationQuat(dir, angle);
			if (reversed) {
				Quat4d q2 = new Quat4d();
				q2.set(new AxisAngle4d(MathTools.Y_AXIS, Math.PI));
				q.mulInverse(q2);
			}
			return q;
		} else {
			Vector3d dir = getPathLegDirection(Direction.PREVIOUS);
			dir.negate();
			if (dir.lengthSquared() > MathTools.NEAR_ZERO)
				dir.normalize();
			return getControlPointOrientationQuat(dir, turnAxis, angle);
		}
	}



	public static Quat4d getControlPointOrientationQuat(Vector3d dir, double angle) {
		if (dir.lengthSquared() < MathTools.NEAR_ZERO)
			return MathTools.getIdentityQuat();


		Vector3d up = new Vector3d(0.0, 1.0, 0.0);
		double a = up.angle(dir);
		if (a < 0.1 || (Math.PI - a) < 0.1) {
			up.set(1.0, 0.0, 0.0);
		}


		return getControlPointOrientationQuat(dir, up, angle);
	}

	public static Quat4d getControlPointOrientationQuat(Vector3d dir, Vector3d up,  double angle) {
		if (dir.lengthSquared() < MathTools.NEAR_ZERO)
			return MathTools.getIdentityQuat();

		final Vector3d front = new Vector3d(1.0,0.0,0.0);

		Quat4d q1 = new Quat4d();


		Vector3d right = new Vector3d();

		right.cross(dir, up);
		up.cross(right, dir);
		right.normalize();
		up.normalize();

		Matrix3d m = new Matrix3d();
		m.m00 = dir.x;
		m.m10 = dir.y;
		m.m20 = dir.z;
		m.m01 = up.x;
		m.m11 = up.y;
		m.m21 = up.z;
		m.m02 = right.x;
		m.m12 = right.y;
		m.m22 = right.z;

		//q1.set(m); MathTools contains more stable conversion
		MathTools.getQuat(m, q1);

		//			if (DEBUG) System.out.println("PipingTools.getPipeComponentOrientationQuat() " + dir+ " " + up + " " + right);

		Quat4d q2 = new Quat4d();
		q2.set(new AxisAngle4d(front, angle));
		q1.mul(q2);
		return q1;
	}

	public void insert(PipeControlPoint previous, PipeControlPoint next) {
		// inserting an offsetpoint is error, 
		if (isDualSub())
			throw new RuntimeException("Dual sub points cannot be inserted.");
		// size change control point cannot be inserted this way, because it ends PipeRun
		if (isSizeChange())
			throw new RuntimeException("Size change points cannot be inserted.");
		PipeRun piperun = previous.getPipeRun();
		// and just to make sure that control point structure is not corrupted
		if (getPipeRun() != null) {
			if (piperun != getPipeRun() || piperun != next.getPipeRun())
				throw new RuntimeException("All controls points must be located on the same pipe run");
		} else {
			piperun.addChild(this);
		}

		// insert new BranchControlPoint between straight's control points
		PipeControlPoint previousNext = previous.getNext();
		PipeControlPoint previousPrevious = previous.getPrevious();

		PipeControlPoint offsetCP = null;
		if (isOffset()) {
			offsetCP = getSubPoint().get(0);
		}
		if (previousNext != null && previousNext == next) {
			if (previous.isDualInline()) {
				throw new RuntimeException();
			}
			if (next.isDualSub()) {
				throw new RuntimeException();
			}
			previous.setNext(this);
			this.setPrevious(previous);
			if (previous.isDualSub()) {
				previous.getParentPoint().setNext(this);
			}
			this.setNext(next);

			if (offsetCP == null) {
				next.setPrevious(this);
			} else {
				next.setPrevious(offsetCP);
				offsetCP.setNext(next);
				offsetCP.setPrevious(previous);
			}

			if (next.isDualInline()) {
				next.getSubPoint().get(0).setPrevious(this);
			}
		} else if (previousPrevious != null && previousPrevious == next) {
			// control point were given in reverse order 
			if (next.isDualInline())
				throw new RuntimeException();
			if (previous.isDualSub())
				throw new RuntimeException();

			this.setNext(previous);
			if (offsetCP == null) {
				previous.setNext(this);
			} else {
				previous.setPrevious(offsetCP);
				offsetCP.setNext(previous);
				offsetCP.setPrevious(next);
			}
			if (previous.isDualInline()) {
				previous.getSubPoint().get(0).setPrevious(this);
			}
			this.setPrevious(next);
			next.setNext(this);
			if (next.isDualSub()) {
				next.getParentPoint().setNext(this);
			}

		} else {
			throw new RuntimeException();
		}	

		PipingRules.validate(piperun);
	}



	public void insert(PipeControlPoint pcp, Direction direction) {
		if (isDualSub())
			throw new RuntimeException();
		if (direction == Direction.NEXT) {
			// if direction is next, user must have given OffsetPoint
			if (pcp.isDualInline())
				throw new RuntimeException();
			// basic next/prev links
			pcp.setNext(this);
			this.setPrevious(pcp);
			// and last take care of sizechange / offset points
			if (pcp.isDualSub()) {
				pcp.getParentPoint().setNext(this);
			}
			if (isDualInline()) {
				getSubPoint().get(0).setPrevious(this);
			}
		} else {
			// if direction is previous, user must have given sizechange
			if (pcp.isDualSub())
				throw new RuntimeException();
			// previous direction is more complicated, since if newCP is SizeChangeControlPoint,
			// we must link pcp to newCP's OffsetPoint
			PipeControlPoint nocp = null;
			if (isDualInline()) {
				nocp = getSubPoint().get(0);
				nocp.setNext(pcp);
			}
			if (nocp == null) {
				pcp.setPrevious(this);
			} else {
				pcp.setPrevious(nocp);
			}
			this.setNext(pcp);
			if (pcp.isDualInline()) {
				PipeControlPoint ocp = pcp.getSubPoint().get(0);
				if (nocp == null)
					ocp.setPrevious(this);
				else
					ocp.setPrevious(nocp);
			}

		}
		PipingRules.validate(getPipeRun());
	}

	public Vector3d getDirectedControlPointDirection() {
		assert (isDirected());
		Vector3d dir = new Vector3d();
		MathTools.rotate(getWorldOrientation(), new Vector3d(1.0, 0.0, 0.0), dir);
		dir.normalize();
		return dir;
	}
	
	public Vector3d getDirection(Direction direction) {
        if (isDirected())
            return getDirectedControlPointDirection();
        if (isTurn() && isFixed()) {
            if (direction == Direction.NEXT) {
                if (previous != null) {
                    PipeControlPoint pcp = this;
                    Vector3d dir = new Vector3d();
                    dir.sub(pcp.getWorldPosition(),previous.getWorldPosition());
                    if (dir.lengthSquared() > MathTools.NEAR_ZERO)
                        dir.normalize();
                    else
                        return null;
                    Quat4d q = getControlPointOrientationQuat(dir, pcp.getRotationAngle() != null ? pcp.getRotationAngle() : 0.0);
                    AxisAngle4d aa = new AxisAngle4d(MathTools.Y_AXIS,pcp.getTurnAngle() == null ? 0.0 : pcp.getTurnAngle());
                    Quat4d q2 = MathTools.getQuat(aa);
                    Vector3d v = new Vector3d(1.,0.,0.);
                    Vector3d offset = new Vector3d();
                    MathTools.rotate(q2, v, offset);
                    MathTools.rotate(q, offset, dir);
                    return dir;
                }
            } else {
                if (next != null) {
                    PipeControlPoint pcp = this;
                    Vector3d dir = new Vector3d();
                    dir.sub(next.getWorldPosition(),pcp.getWorldPosition());
                    if (dir.lengthSquared() > MathTools.NEAR_ZERO)
                        dir.normalize();
                    else
                        return null;
                    Quat4d q = getControlPointOrientationQuat(dir, pcp.getRotationAngle() != null ? pcp.getRotationAngle() : 0.0);
                    AxisAngle4d aa = new AxisAngle4d(MathTools.Y_AXIS,pcp.getTurnAngle() == null ? 0.0 : pcp.getTurnAngle());
                    Quat4d q2 = MathTools.getQuat(aa);
                    Vector3d v = new Vector3d(1.,0.,0.);
                    Vector3d offset = new Vector3d();
                    MathTools.rotate(q2, v, offset);
                    MathTools.rotate(q, offset, dir);
                    return dir;
                }
            }
        }
        return null;
    }

	public Vector3d getPathLegDirection(Direction direction) {
		if (direction == Direction.NEXT) {
			if (next != null) {
				PipeControlPoint pcp = this;
				if (pcp.isDualInline()) {
					pcp = pcp.getSubPoint().get(0);
				}
				Vector3d v = new Vector3d();
				v.sub(next.getWorldPosition(),pcp.getWorldPosition());
				return v;
			} else {
				if (previous == null) {
					if (!isDirected())
						throw new RuntimeException("Cannot calculate path leg direction for unconnected control point " + this);
					return getDirectedControlPointDirection();

				} else {
					if (isVariableAngle())
						throw new RuntimeException("Cannot calculate path leg direction for unconnected variable angle control point " + this);
					if (isInline()) {
						PipeControlPoint pcp = this;
						if (pcp.isDualSub()) {
							pcp = pcp.getParentPoint();
						}
						Vector3d v = new Vector3d();
						v.sub(pcp.getWorldPosition(),previous.getWorldPosition());
						return v;
					} else if (isDirected()) {
						return getDirectedControlPointDirection();
					} else if (isEnd()) {
						Vector3d v = new Vector3d();
						v.sub(getWorldPosition(),previous.getWorldPosition());
						return v;
					} else if (isTurn() && isFixed() && !_getReversed()) {
						return getDirection(Direction.NEXT);
					}
					throw new RuntimeException("Missing implementation " + this);
				}
			}
		} else {
			if (previous != null) {
				PipeControlPoint pcp = this;
				if (isDualSub()) 
					pcp = getParentPoint();
				Vector3d v = new Vector3d();
				v.sub(previous.getWorldPosition(),pcp.getWorldPosition());
				return v;
			} else {
				if (next == null)  {
					if (!isDirected())
						throw new RuntimeException("Cannot calculate path leg direction for unconnected control point " + this);
					Vector3d v = getDirectedControlPointDirection();
					v.negate();
					return v;
				} else {
					if (isVariableAngle())
						throw new RuntimeException("Cannot calculate path leg direction for unconnected variable angle control point " + this);
					if (isInline()) {
						PipeControlPoint pcp = this;
						if (pcp.isDualInline()) {
							pcp = pcp.getSubPoint().get(0);
						}
						Vector3d v = new Vector3d();
						v.sub(pcp.getWorldPosition(),next.getWorldPosition());
						return v;
					} else if (isDirected()) {
						Vector3d v = getDirectedControlPointDirection();
						v.negate();
						return v;
					} else if (isEnd()) {
						Vector3d v = new Vector3d();
						v.sub(getWorldPosition(),next.getWorldPosition());
						return v;
					} else if (isTurn() && isFixed() && _getReversed()) {
						return getDirection(Direction.PREVIOUS);
					}
					throw new RuntimeException("Missing implementation " + this);
				}
			}
		}
	}

	public void getInlineControlPointEnds(Tuple3d p1, Tuple3d p2) {
		assert (isInline());

		PipeControlPoint sub = isAxial() ? this : getSubPoint().get(0);
		Vector3d pos = getWorldPosition(), pos2 = sub == this ? pos : sub.getWorldPosition();
		Vector3d dir = sub.getPathLegDirection(Direction.NEXT);
		
		dir.normalize();
		dir.scale(length * 0.5);
		p1.set(pos);
		p2.set(pos2);
		p1.sub(dir);
		p2.add(dir);
	}

	public void getControlPointEnds(Tuple3d p1, Tuple3d p2) {
		PipeControlPoint sub = isAxial() ? this : getSubPoint().get(0);
		Vector3d pos = getWorldPosition(), pos2 = sub == this ? pos : sub.getWorldPosition();
		
		Vector3d dir1 = getPathLegDirection(Direction.PREVIOUS);
		dir1.normalize();
		Vector3d dir2 = sub.getPathLegDirection(Direction.NEXT);
		dir2.normalize();
		if (isInline()) {
			dir1.scale(length * 0.5);
			dir2.scale(length * 0.5);
		} else {
			dir1.scale(length);
			dir2.scale(length);
		}
		p1.set(pos);
		p2.set(pos2);
		p1.add(dir1);
		p2.add(dir2);
	}

	public void getEndDirections(Tuple3d v1, Tuple3d v2) {
		PipeControlPoint sub = isAxial() ? this : getSubPoint().get(0);
		
		Vector3d dir1 = getPathLegDirection(Direction.PREVIOUS);
		dir1.normalize();
		Vector3d dir2 = sub.getPathLegDirection(Direction.NEXT);
		dir2.normalize();
		v1.set(dir1);
		v2.set(dir2);
	}

	public void getInlineControlPointEnds(Tuple3d p1, Tuple3d p2, Vector3d dir) {
		assert (isInline());

		Vector3d pos = getWorldPosition();
		dir.set(getPathLegDirection(Direction.NEXT));
		dir.normalize();
		dir.scale(length * 0.5);
		p1.set(pos);
		p2.set(pos);
		p1.sub(dir);
		p2.add(dir);
	}

	public void getInlineControlPointEnds(Tuple3d center, Tuple3d p1, Tuple3d p2, Vector3d dir) {
		assert (isInline());

		Vector3d pos = getWorldPosition();
		center.set(pos);
		dir.set(getPathLegDirection(Direction.NEXT));
		dir.normalize();
		dir.scale(length * 0.5);
		p1.set(pos);
		p2.set(pos);
		p1.sub(dir);
		p2.add(dir);
	}

	public double getInlineLength() {
		if (type == PointType.TURN)
			return length;
		else if (type == PointType.INLINE)
			return length * 0.5;
		return 0;
	}

	public Vector3d getRealPosition(PositionType type) {
		Vector3d pos = getWorldPosition();
		switch (type) {
		case NEXT: {
			Vector3d dir = getPathLegDirection(Direction.NEXT);
			double length = getInlineLength();
			dir.normalize();
			dir.scale(length);
			pos.add(dir);
			break;
		}
		case PREVIOUS: {
			Vector3d dir = getPathLegDirection(Direction.PREVIOUS);
			double length = getInlineLength();
			dir.normalize();
			dir.scale(length);
			pos.add(dir);
			break;
		}
		case PORT:
			// IEntity portDir = pcp.getSingleRelatedObject(ProcessResource.plant3Dresource.HasDirection);
			// TODO : how we calculated needed space for a port; does it has an offset from control point's position or not?
			break;
		case SPLIT:
			// do nothing
			break;

		}
		return pos;
	}

	public void getInlineMovement(Tuple3d start, Tuple3d end) {
		// FIXME : check type of neighbor components and allow movement on top of variable length components,
		//         find proper range for movement (pcp's position is not)
		PipeControlPoint p = previous.getPrevious();
		PipeControlPoint n = next.getNext();
		start.set(p.getWorldPosition());
		end.set(n.getWorldPosition());
	}

	public PipeControlPoint findNextEnd() {
		ArrayList<PipeControlPoint> t = new ArrayList<PipeControlPoint>();
		return findNextEnd( t);
	}

	public PipeControlPoint findPreviousEnd() {
		ArrayList<PipeControlPoint> t = new ArrayList<PipeControlPoint>();
		return findPreviousEnd(t);
	}

	public PipeControlPoint findNextEnd(List<PipeControlPoint> nextList) {
		while (true) {
			PipeControlPoint pcp = null;
			PipeControlPoint p = null;
			if (nextList.size() == 0)
				p = this;

			else
				p = nextList.get(nextList.size() - 1);

			pcp = p.getNext();
			if (pcp == null) {
				pcp = p;
				if (nextList.size() > 0)
					nextList.remove(nextList.size() - 1);
				//            	if (DEBUG) System.out.println(" " + pcp.getResource() + " not full");
				return pcp;
				//break;
			}
			if (pcp.isPathLegEnd()) {
				//if (DEBUG) System.out.println(" " + pcp.getResource());
				return pcp;
			} else {
				nextList.add(pcp);
				// if (DEBUG) System.out.print(" " + pcp.getResource());
			}
		}
	}

	public PipeControlPoint findPreviousEnd(List<PipeControlPoint> prevList) {
		while (true) {
			PipeControlPoint pcp = null;
			PipeControlPoint p = null;
			if (prevList.size() == 0)
				p = this;

			else
				p = prevList.get(prevList.size() - 1);

			pcp = p.getPrevious();
			if (pcp == null) {
				pcp = p;
				if (prevList.size() > 0)
					prevList.remove(prevList.size() - 1);
				//				if (DEBUG) System.out.println(" " + pcp.getResource() + " not full");
				return pcp;
			}
			if (pcp.isPathLegEnd()) {
				//				if (DEBUG)	System.out.println(" " + pcp.getResource());
				return pcp;
			} else {
				prevList.add(pcp);
				//				if (DEBUG)System.out.print(" " + pcp.getResource());
			}
		}
	}
	
	public void _remove() {
	    _remove(true);
	}

	public void _remove(boolean renconnect) {
		if (component == null && next == null && previous == null)
			return;
		if (DEBUG) System.out.println(this + " Remove " + renconnect);
		if (isDualInline() || isDualSub()) {
			removeDualPoint();
			return;
		}
		if (getParentPoint() != null) {
		    getParentPoint()._remove(renconnect);
		}
		PipeRun pipeRun = getPipeRun();
		if (pipeRun == null)
			return;

		PipeControlPoint additionalRemove = null;
		if (!PipingRules.isEnabled()) {
			setPrevious(null);
			setNext(null);
		} else {

			PipeControlPoint currentPrev = previous;
			PipeControlPoint currentNext = next;
			if (currentNext == null && currentPrev == null) {
				removeComponent();
				pipeRun.remChild(this);
				return;
			}
			if (currentNext != null && currentPrev != null) {
				boolean link = renconnect;
				if (currentNext.isBranchEnd()) {
					link = false;
					//					currentNext.setPrevious(null);
					//					currentNext.setNext(null);
					currentNext.remove();
					currentNext = null;
					setNext(null);
				}
				if (currentPrev.isBranchEnd()) {
					link = false;
					//					currentPrev.setPrevious(null);
					//					currentPrev.setNext(null);
					currentPrev.remove();
					currentPrev = null;
					setPrevious(null);
				}
				if (link && currentPrev.isDirected() && currentNext.isDirected()) {
					link = false;
				}
				if (currentNext == null) {
					// Nothing to do
				} else if (currentNext.isDualInline()) {
					PipeControlPoint sccp = currentNext;
					PipeControlPoint ocp = sccp.getSubPoint().get(0);
					if (ocp == null) {
						throw new RuntimeException("Removing PipeControlPoint " + this+ " structure damaged, no offset control point");
					}
					if (link) {
						sccp.setPrevious(currentPrev);
						ocp.setPrevious(currentPrev);
					} else {
						sccp.setPrevious(null);
						ocp.setPrevious(null);
					}
					setNext(null);
				} else if (currentNext.isDualSub()) {
					throw new RuntimeException("Removing PipeControlPoint " + this + " structure damaged, next control point is offset control point");
				} else if (currentNext.previous == this) {
					if (link) {
						currentNext.setPrevious(currentPrev);
					} else {
						currentNext.setPrevious(null);
					}
					setNext(null);
				} else {
					throw new RuntimeException("Removing PipeControlPoint "+ this+ " structure damaged");
				}
				if (currentPrev == null) {
					// Nothing to do
				} else if (currentPrev.isDualInline()) {
					throw new RuntimeException("Removing PipeControlPoint " + this + " structure damaged, previous control point is size change control point");
				} else if (currentPrev.isDualSub()) {
					PipeControlPoint ocp = currentPrev;
					PipeControlPoint sccp = ocp.getParentPoint();
					if (sccp == null)
						throw new RuntimeException("Removing PipeControlPoint " + this + " structure damaged, no size change control point");
					if (link) {
						ocp.setNext(currentNext);
						sccp.setNext(currentNext);
					} else {
						ocp.setNext(null);
						sccp.setNext(null);
					}
					setPrevious(null);
				} else if (currentPrev.next == this) {
					if (link)
						currentPrev.setNext(currentNext);
					else
						currentPrev.setNext(null);

					setPrevious(null);
				} else {
					throw new RuntimeException("Removing PipeControlPoint " + this + " structure damaged");
				}
				if (link) {
					if (currentNext.isVariableLength() && currentPrev.isVariableLength()) {
						// we have to join them into single variable length component.
						additionalRemove = currentPrev;
						// combine lengths and set the location of remaining control point to the center.
						Point3d ps = new Point3d();
						Point3d pe = new Point3d();
						Point3d ns = new Point3d();
						Point3d ne = new Point3d();
						currentPrev.getInlineControlPointEnds(ps, pe);
						currentNext.getInlineControlPointEnds(ns, ne);
						double l = currentPrev.getLength() + currentNext.getLength();
						Vector3d cp = new Vector3d();
						cp.add(ps, ne);
						cp.scale(0.5);
						currentNext.setLength(l);
						currentNext.setWorldPosition(cp);
					}
				} else {
					// FIXME : pipe run must be split into two parts, since the control point structure is no more continuous. 
				}
			} else if (next != null) {
				if (next.isDualInline()) {
					PipeControlPoint sccp = next;
					PipeControlPoint ocp = sccp.getSubPoint().get(0);
					if (ocp == null) {
						throw new RuntimeException("Removing PipeControlPoint " + this+ " structure damaged, no offset control point");
					}
					sccp.setPrevious(null);
					ocp.setPrevious(null);
				} else if (next.isDualSub()) {
					throw new RuntimeException("Removing PipeControlPoint " + this + " structure damaged, next control point is offset control point");
				} else if (next.previous == this) {
					next.setPrevious(null);
				} else {
					throw new RuntimeException("Removing PipeControlPoint "+ this+ " structure damaged");
				}
				setNext(null);
			} else {  //(previous != null)
				if(previous.isDualInline()) {
					throw new RuntimeException("Removing PipeControlPoint " + this + " structure damaged, previous control point is size change control point");
				} else if (previous.isDualSub()) {
					PipeControlPoint ocp = previous;
					PipeControlPoint sccp = ocp.getParentPoint();
					if (sccp == null) {
						throw new RuntimeException("Removing PipeControlPoint " + this + " structure damaged, no size change control point");
					}
					ocp.setNext(null);
					sccp.setNext(null);
				} else if (previous.next == this) {
					previous.setNext(null);
				} else {
					throw new RuntimeException("Removing PipeControlPoint "+ this+ " structure damaged");
				}
				setPrevious(null);
			}
			if (children.size() > 0 ) {
				removeSubPoints();
			} else if (parent!= null) {
				removeParentPoint();
			}

		}

		removeComponent();
		pipeRun.remChild(this);
		checkRemove(pipeRun);
		if (PipingRules.isEnabled() && pipeRun.getParent() != null && pipeRun.getControlPoints().size() > 0)
			PipingRules.validate(pipeRun);
		if (additionalRemove != null)
			additionalRemove.remove();
	}

	/**
	 * Removes control point and attempts to reconnect next/prev
	 * 
	 * If this point is size change (PipeRuns are different on both sides), then reconnection cannot be made.
	 */
	public void remove() {
		PipeControlPoint currentPrev = previous;
		PipeControlPoint currentNext = next;
		_remove();
		try {
			if (currentNext != null)
			    if (!currentNext.checkRemove())
			        PipingRules.requestUpdate(currentNext);
			if (currentPrev != null)
			    if (!currentPrev.checkRemove())
			        PipingRules.requestUpdate(currentPrev);
		} catch (Exception e) {
			e.printStackTrace();
		}
	}
	
	
	/**
	 * Removes control point without attempting to reconnect next/prev.
	 * This usually leads to creation of another PipeRun for the control points after this point. 
	 */
	public void removeAndSplit() {
        PipeControlPoint currentPrev = previous;
        PipeControlPoint currentNext = next;
        
        if (next != null && previous != null) {
            P3DRootNode root = (P3DRootNode)getPipelineComponent().getRootNode();
            PipeRun nextPipeRun = new PipeRun();
            nextPipeRun.setName(root.getUniqueName("PipeRun"));
            root.addChild(nextPipeRun);
            
            PipeRun previousRun = previous.getPipeRun();
            nextPipeRun.setPipeDiameter(previousRun.getPipeDiameter());
            nextPipeRun.setTurnRadius(previousRun.getTurnRadius());
            
            PipelineComponent n = next.getPipelineComponent();
            while (n != null) {
                if (! (n instanceof Nozzle)) {
                    n.deattach();
                    nextPipeRun.addChild(n);
                } else
                    n.setPipeRun(nextPipeRun);
                n = n.getNext();
            }
        }
        _remove(false);
        try {
            if (currentNext != null)
                if (!currentNext.checkRemove())
                    PipingRules.requestUpdate(currentNext);
            if (currentPrev != null)
                if (!currentPrev.checkRemove())
                    PipingRules.requestUpdate(currentPrev);
        } catch (Exception e) {
            e.printStackTrace();
        }
    }
	
	/**
	 * This is called when adjacent control point is removed.
	 * 
	 * This call should remove the give point, if the point cannot exist alone. 
	 * At the moment there is one such case: branch.
	 * 
	 * @return
	 */
	protected boolean checkRemove() {
	    if (getParentPoint() != null) {
	        return getParentPoint().checkRemove();
	    } else {
	        if (getPipelineComponent() == null)
	            return true; // already removed
    	    if (getPipelineComponent().getType().equals("Plant3D.URIs.Builtin_BranchSplitComponent")) {
    	        if (getSubPoint().get(0).getNext() == null && getSubPoint().get(0).getPrevious() == null) {
        	        remove();
        	        return true;
    	        }
    	    }
            return false;
	    }
    }

	private void checkRemove(PipeRun pipeRun) {
		Collection<PipeControlPoint> points = pipeRun.getControlPoints();
		if (points.size() == 0) {
			pipeRun.remove();
		} else if (points.size() == 1) {
			PipeControlPoint pcp = points.iterator().next();
			if (pcp.isDeletable())
				pcp._remove();
		}
	}

	private void removeDualPoint() {
		if (previous != null)
			previous.setNext(null);
		if (next != null)
			next.setPrevious(null);
		PipeControlPoint ocp;
		PipeControlPoint sccp;
		if (isDualInline()) {
			sccp = this;
			ocp = getSubPoint().get(0);
		} else {
			ocp = this;
			sccp = getParentPoint();
		}
		PipeRun p1 = ocp.getPipeRun();
		PipeRun p2 = sccp.getPipeRun();

		ocp.removeComponent();
		sccp.removeComponent();

		if (p1 != null)
			p1.remChild(ocp);
		if (p2 != null)
			p2.remChild(sccp);

		// TODO : now we assume that this is size change, and we do
		if (ocp.next != null)
			ocp.next.setPrevious(null);
		if (ocp.previous != null)
			ocp.previous.setNext(null);
		if (sccp.next != null)
			sccp.next.setPrevious(null);
		if (sccp.previous != null)
			sccp.previous.setNext(null);
		ocp.setNext(null);
		ocp.setPrevious(null);
		sccp.setNext(null);
		sccp.setPrevious(null);

		if (p1 != null)
			checkRemove(p1);
		if (p2 != null)
			checkRemove(p2);
	}

	private void removeSubPoints() {
		for (PipeControlPoint p : children) {
			// TODO : this may affect delete routine, since classification of the point changes.
			p.parent = null;
			p._remove();
		}
		children.clear();
	}

	private void removeParentPoint() {
		throw new RuntimeException("Child points cannot be removed directly");
	}

	private void removeComponent() {
		if (component == null)
			return;
		PipelineComponent next = component.getNext();
		PipelineComponent prev = component.getPrevious();
		PipelineComponent br0 = component.getBranch0();
		component.setNext(null);
		component.setPrevious(null);
		component.setBranch0(null);
		if (next != null) {
			if (next.getNext() == component)
				next.setNext(null);
			else if (next.getPrevious() == component)
				next.setPrevious(null);
			else if (next.getBranch0() == component)
				next.setBranch0(null);
		}
		if (prev != null) {
			if (prev.getNext() == component)
				prev.setNext(null);
			else if (prev.getPrevious() == component)
				prev.setPrevious(null);
			else if (prev.getBranch0() == component)
				prev.setBranch0(null);
		}
		if (br0 != null) {
			if (br0.getNext() == component)
				br0.setNext(null);
			else if (br0.getPrevious() == component)
				br0.setPrevious(null);
			else if (br0.getBranch0() == component)
				br0.setBranch0(null);
		}
		PipelineComponent comp = component;
		component = null;

		comp.remove();
	}

	@Override
	public void setOrientation(Quat4d orientation) {
		if (MathTools.equals(orientation, getOrientation()))
			return;
		super.setOrientation(orientation);
		if (getParentPoint() == null && component != null)
			component._setWorldOrientation(getWorldOrientation());
		updateSubPoint();
	}

	@Override
	public void setPosition(Vector3d position) {
		if (MathTools.equals(position, getPosition()))
			return;
		if (Double.isNaN(position.x) || Double.isNaN(position.y) || Double.isNaN(position.z))
			throw new IllegalArgumentException("NaN is not supported");
		super.setPosition(position);
		if (getParentPoint() == null && component != null)
			component._setWorldPosition(getWorldPosition());
		updateSubPoint();
	}

	private void updateSubPoint() {
		if (isOffset()) {
			if (next == null && previous == null) {
				for (PipeControlPoint sub : getSubPoint()) {
					sub.setWorldPosition(getWorldPosition());
					sub.setWorldOrientation(getWorldOrientation());
				}
				return;
			}
			for (PipeControlPoint sub : getSubPoint()) {
				Vector3d wp = getWorldPosition();
				wp.add(getSizeChangeOffsetVector());
				sub.setWorldPosition(wp);
				sub.setWorldOrientation(getWorldOrientation());
			}
		} else {
			for (PipeControlPoint sub : getSubPoint()) {
				sub.setWorldPosition(getWorldPosition());
				sub.setWorldOrientation(getWorldOrientation());
			}
		}
	}


	public void _setWorldPosition(Vector3d position) {
		Vector3d localPos = getLocalPosition(position);
		super.setPosition(localPos);
		updateSubPoint();
	}

	public void _setWorldOrientation(Quat4d orientation) {
		Quat4d localOr = getLocalOrientation(orientation);
		super.setOrientation(localOr);
		updateSubPoint();
	}

	@Override
	public String toString() {
		return getClass().getName() + "@" + Integer.toHexString(hashCode());
	}

}