run_piece_GdM.py

# piece_GdM (guepardo00)

from Testes import *
from multiprocessing import Pool

# degree=np.arange(3,31,2)
degree=np.arange(4,31,2)
# degree=np.arange(3,31)

# pool = Pool(processes=22) # start 23 worker processes
pool = Pool(processes=40) # start 23 worker processes

T=3
T_frac_snapshot=1/2
methods=np.array(['RK7','RK2','RK4','2MS','FA','HORK','KRY'])

# pool.apply_async(wave_eq, args=(0.005,'scalar_dx2',2,1,0.8,30,'8',T,T_frac_snapshot,np.array([1]),0,'piece_GdM_c','RK7','H_amplified',))

for l in range(50):
# for l in range(1):
    for i in range(4):
    # for i in range(1):
        pool.apply_async(wave_eq, args=(0.01,'scalar_dx2',2,1,0.8,30,'8',T,T_frac_snapshot,np.array([l+1]),0,'piece_GdM_b',methods[i],'H_amplified',))

    for i in range(len(degree)):
        for j in range(3):
                pool.apply_async(wave_eq, args=(0.01,'scalar_dx2',2,1,0.8,30,'8',T,T_frac_snapshot,np.array([l+1]),np.array([degree[i]]),'piece_GdM_b',methods[4+j],'H_amplified',))
pool.close()
pool.join()