print ' SOLVE THE LINEARS EQUATIONS:' print ' ---------------------------- \n' print ' Saving simulation in %s \n' %folderResults solution_start_time = time() os.chdir(initial_path) # ------------------------ Export VTK File --------------------------------------- # Linear and Mini Elements if polynomial_option == 0 or polynomial_option == 1 or polynomial_option == 2: save = exportVTK.Linear2D(x,y,IEN,numNodes,numElements,w,w,psi,vx,vy) save.create_dir(folderResults) save.saveVTK(folderResults + str(0)) # Quad Element elif polynomial_option == 3: save = exportVTK.Quad2D(x,y,IEN,numNodes,numElements,w,w,psi,vx,vy) save.create_dir(folderResults) save.saveVTK(folderResults + str(0)) # --------------------------------------------------------------------------------- vorticityAux1BC = np.zeros([numNodes,1], dtype = float) x_old = np.zeros([numNodes,1], dtype = float) y_old = np.zeros([numNodes,1], dtype = float)
print ' SOLVE THE LINEARS EQUATIONS:' print ' ---------------------------- \n' print ' Saving simulation in %s \n' %folderResults solution_start_time = time() os.chdir(initial_path) # ------------------------ Export VTK File --------------------------------------- # Linear and Mini Elements if polynomial_option == 0 or polynomial_option == 1 or polynomial_option == 2: save = exportVTK.Linear2D(x,y,IEN,numVerts,numElements,p,p,c,vx,vy) save.create_dir(folderResults) save.saveVTK(folderResults + str(0)) # Quad Element elif polynomial_option == 3: save = exportVTK.Quad2D(x,y,IEN,numNodes,numElements,w,psi,c,vx,vy) save.create_dir(folderResults) save.saveVTK(folderResults + str(0)) # --------------------------------------------------------------------------------- x_old = np.zeros([numNodes,1], dtype = float) y_old = np.zeros([numNodes,1], dtype = float) vx_old = np.zeros([numNodes,1], dtype = float)
print ' Schmidt number: %s' % Sc
print ""
print ' ----------------------------'
print ' SOLVE THE LINEARS EQUATIONS:'
print ' ---------------------------- \n'
print ' Saving simulation in %s \n' % folderResults
solution_start_time = time()
os.chdir(initial_path)
# ------------------------ Export VTK File ---------------------------------------
# Linear and Mini Elements
if polynomial_option == 1 or polynomial_option == 2:
save = exportVTK.Linear2D(z, r, IEN, numNodes, numElements, w, w, psi, vz,
vr)
save.create_dir(folderResults)
save.saveVTK(folderResults + str(0))
# Quad Element
elif polynomial_option == 3:
save = exportVTK.Quad2D(z, r, IEN, numNodes, numElements, w, w, psi, vz,
vr)
save.create_dir(folderResults)
save.saveVTK(folderResults + str(0))
# ---------------------------------------------------------------------------------
vorticityAux1BC = np.zeros([numNodes, 1], dtype=float)
vz_old = np.zeros([numNodes, 1], dtype=float)
vr_old = np.zeros([numNodes, 1], dtype=float)
end_type = 0