Line.py

from Vector3D import *
class Line:
    """Vectorial data structure representing a line"""
    def __init__(self,point,direction):
        self.position= point
        self.direction = direction.normalize()
        
    def __repr__(self):
        return "Line({},{})".format(self.position,self.direction.normalize())
    
    def origin_point(self):
        return self.position - self.position.x*self.direction
    
    def point(self,t):
        return self.position + t*self.direction
    
    def plane_intersect(self,plane):
        """Calculates the intersection of a plane of the plane class and itself, requires the use of the plane module! if it does not find an intercept it retuns None"""
        #print("{},{}".format((plane.distance_to_origin() - plane.normal.x*self.position.x - plane.normal.y*self.position.y - plane.normal.z*self.position.z),(plane.normal.x*self.direction.x + plane.normal.y*self.direction.y + plane.normal.z*self.direction.z)))   
        dp = self.direction.dot_product(plane.normal)
        if dp == 0:
            raise Exception("Line {} and plane {} are parallel".format(self, plane))
            return None
        
        tVal=((plane.position-self.position).dot_product(plane.normal))/ dp # ACCORDING to the formula on wikipedia vector form: https://en.wikipedia.org/wiki/Line%E2%80%93plane_intersection
        return Vector(tVal*self.direction.x + self.position.x,tVal*self.direction.y + self.position.y,tVal*self.direction.z + + self.position.z)