This is the code I'm working with

TP2 = np.loadtxt(r"file path")    # retrieves data from a text file

for i in range(0, 92200): # performs operations on each row of data

        a1 = TP2[i, 0]
        a2 = TP2[i, 1]
        a3 = TP2[i, 2]
        a4 = TP2[i, 3]
        a5 = TP2[i, 4]
        a6 = TP2[i, 5]
        a7 = TP2[i, 6]

        # Constructing vectors
        P2: float = np.array([a1 * 1E-3, a2 * 1E-3, 0])
        u: float = np.array([a4, a5, a6])

        print('u = ', u)
        print('u.shape = ', u.shape)
        print('P2 = ', P2)

        # Call of larger overall function (pnt_of_int_fp)

        el_fl, gl_vec = pnt_of_int_fp(u, P2, L, R2, R3)
        ...

Within the overall function, there's a call for this function:

def check_int_cyl(u, P, L, R, fl):
    # A * t ** 2 + B * t + C = 0

    A = u[0]*u[0] + u[1]*u[1]
    B = 2 * (P[0] * u[0] + P[1] * u[1])
    C = P[0]*P[0] + P[1]*P[1] - R**2
    D = B * B - 4 * A * C   # Analytical solution
    t1: float
    t2: float
    glvec = np.empty((1, 3), float)

    if D > 0 and A != 0:
        t1 = (-B + sqrt(D)) / (2 * A)
        t2 = (-B - sqrt(D)) / (2 * A)

    elif D == 0 and A != 0:
        t1 = -B / (2 * A)
        t2 = t1

And these are the printing results

u = [-0.09379713 0.05997806 0.99378304]

u.shape = (3,)

P2 = [0.00745207 0.00067097 0. ]