def init_value(self):
rhoA = init_value['C42A']
rhoB = init_value['C42B']
rhoC = init_value['C42C']
print('rhoA:', rhoA)
print('rhoB:', rhoB)
print('rhoC:', rhoC)
box = np.array([[4.1, 0], [0, 4.1]], dtype=np.float)
NS = 8
space = FourierSpace(box, NS)
rhoA = space.fourier_interpolation(rhoA)
rhoB = space.fourier_interpolation(rhoB)
rhoC = space.fourier_interpolation(rhoC)
print('rhoA:', rhoA)
print('rhoB:', rhoB)
print('rhoC:', rhoC)
def init_value(self):
rhoA = init_value['LAM_A']
rhoB = init_value['LAM_B']
rhoC = init_value['LAM_C']
print('rhoA:', rhoA)
print('rhoB:', rhoB)
print('rhoC:', rhoC)
box = np.array([[8, 0, 0], [0, 8, 0], [0, 0, 8]], dtype=np.float)
NS = 64
space = FourierSpace(box, NS)
rhoA = space.fourier_interpolation(rhoA)
rhoB = space.fourier_interpolation(rhoB)
rhoC = space.fourier_interpolation(rhoC)
print('rhoA:', rhoA)
print('rhoB:', rhoB)
print('rhoC:', rhoC)
