Vector projection onto another vector

#1

Given two splines, I’ve managed to calculate the length of the vector projection of a onto b (see attached image).
The result is represented by the green square.
The question is: how do I calculate the same result but ALONG the vector b?

getProjection
getProjection.jpg1823x893 110 KB

Here’s my code
–get the vectors
firstPoint1 = getKnotPoint selection[1] 1 1
secondPoint1 = getKnotPoint selection[1] 1 2
newVec1 = secondPoint1 - firstPoint1
normalizedVector1 = normalize newVec1
firstPoint2 = getKnotPoint selection[2] 1 1
secondPoint2 = getKnotPoint selection[2] 1 2
newVec2 = secondPoint2 - firstPoint2
normalizedVector2 = normalize newVec2

–get the angle between the splines
result = acos (dot normalizedVector1 normalizedVector2)

–get the length of the projection
result2 = (length newVec1)*(cos result)

b = box width:2 length:2 height:2
b.position = [result2,0,0]

I’ve found the solution myself: I needed to multiply that result by the normalized b vector

#2 
/*********************** MATH *********************/
fn projectPoint2Vec a b p = 
(
	t = dot (normalize (b - a)) (normalize (p - a))
	c += normalize (b - a) * length (p - a) * t
)

/******************* TEST SCENE *******************/
delete objects 
gc()

ab = splineshape name:"AB" wirecolor:orange
addnewspline ab
addKnot ab 1 #corner #line [0,0,0]
addKnot ab 1 #corner #line [100,0,0]
updateshape ab

cd = splineshape name:"CD" wirecolor:green
addnewspline cd
addKnot cd 1 #corner #line [0,0,0]
addKnot cd 1 #corner #line [150,0,-50]
updateshape cd


/********************* CALC ***********************/
a = getknotpoint ab 1 1
b = getknotpoint ab 1 2

c = getknotpoint cd 1 1
d = getknotpoint cd 1 2

pt = projectPoint2Vec c d b

pr = point name:"PROJ" pos:pt wirecolor:red cross:on