6fd1cef5797fd76e8ea6b68ee9338dce759768ad,trunk/SUAVE/Methods/Aerodynamics/Common/Fidelity_Zero/Lift/compute_induced_velocity_matrix.py,,vortex_leg_from_B_to_inf,#,225

Before Change

    Y1_Y  = Y1-Y
    Z_Z1  = Z-Z1    

    DENUM =  Z_Z1*Z_Z1 + Y1_Y*Y1_Y
    DENUM[DENUM==0] = 1e-32
    XVEC  = -Y1_Y*np.sin(tw)/DENUM
    YVEC  = Z_Z1/DENUM
    ZVEC  = Y1_Y*np.cos(tw)/DENUM
    RVEC   = np.array([XVEC, YVEC, ZVEC])
    BRAC  =  1 + (X_X1 / (np.sqrt(X_X1*X_X1+ Y_Y1*Y_Y1+ Z_Z1*Z_Z1)))    
    COEF  = -(1/(4*np.pi))*RVEC*BRAC      
    
    if np.isnan(COEF).any():
        print("NaN!")   
        
    return COEF

def compute_mach_cone_matrix(XC,YC,ZC,MCM,mach):
    for m_idx in range(len(mach)):

After Change

    Y1_Y  = Y1-Y
    Z_Z1  = Z-Z1    

    DENUM =  np.square(Z_Z1) + np.square(Y1_Y)
    DENUM[DENUM==0] = 1e-32
    XVEC  = -Y1_Y*np.sin(tw)/DENUM
    YVEC  = Z_Z1/DENUM
    ZVEC  = Y1_Y*np.cos(tw)/DENUM
    RVEC  = np.array([XVEC, YVEC, ZVEC])
    BRAC  = 1 + (X_X1 / (np.sqrt(np.square(X_X1)+ np.square(Y_Y1)+ np.square(Z_Z1))))    
    COEF  = -(1/(4*np.pi))*RVEC*BRAC      

    return COEF

Instances