def METHOD_NAME(v, d, k, eta, rho):
    v_abs = np.abs(v)
    v_abs[v_abs < 1e-6] = 1e-6
    re = np.divide(rho * v_abs * d, eta)
    lambda_laminar = np.zeros_like(v)
    lambda_laminar[v != 0] = 64 / re[v != 0]
    lambda_nikuradse = np.divide(1, (-2 * np.log10(k / (3.71 * d))) ** 2)
    return re, lambda_laminar, lambda_nikuradse