Computational Physics, Python, Math, Physics.
celer(X[i],Y[j]) ############## CALCULATION CONSTANTS ############### Cx2 = (dt/dx)**2 Cy2 = (dt/dy)**2 CFL_1 = dt/dy*c[:,0] CFL_2 = dt/dy*c[:,N_y] CFL_3 = dt/dx*c[0,:] CFL_4 =dt/dx*c[N_x,:] ############## PROCESSING LOOP ############### if loop_exec: # $\forall i \in {0,...,N_x}$ U =