Another way to do this is to construct a matrix from the point normal and your two vectors, then find the determinant of the matrix and use it and the dot product to find the angle.
vector N = point(0, "N", P0);
float dot = dot(v0, v1);
matrix3 mat = set(v0.x, v1.x, N.x,
v0.y, v1.y, N.y,
v0.z, v1.z, N.z);
float det = determinant(mat);
return (360 + degrees(atan2(det, dot))) % 360;
getHoleAngles.hipnc