def computeAMROpoints(B_amp,B_phi_a,B_theta_a):
amroListForPhi = []
for i in range(B_phi_a.shape[0]):
amroListForTheta = []
for j in range(B_theta_a.shape[0]):
dataPoint = Conductivity(band, B_amp, B_phi_a[i], B_theta_a[j], gamma_0=15, gamma_k=0, power=12)
dataPoint.solveMovementFunc()
dataPoint.chambersFunc(i = 2, j = 2)
amroListForTheta.append(dataPoint)
amroListForPhi.append(amroListForTheta)
return amroListForPhi
Bphi=0,
Btheta=0,
gamma_0=1,
gamma_k=0,
power=0)
condObject.Ntime = 2000
## Empty arrays
rhoxx_array = np.empty_like(B_array, dtype=np.float64)
rhoxy_array = np.empty_like(B_array, dtype=np.float64)
RH_array = np.empty_like(B_array, dtype=np.float64)
for i, B in enumerate(B_array):
condObject.Bamp = B
condObject.solveMovementFunc()
condObject.chambersFunc(0, 0)
condObject.chambersFunc(0, 1)
sigma_xx = condObject.sigma[0, 0]
sigma_xy = condObject.sigma[0, 1]
rhoxx = sigma_xx / (sigma_xx**2 + sigma_xy**2) # Ohm.m
rhoxy = sigma_xy / (sigma_xx**2 + sigma_xy**2) # Ohm.m
RH = rhoxy / B # m^3/C
rhoxx_array[i] = rhoxx
rhoxy_array[i] = rhoxy
RH_array[i] = RH
print("{0:.0f}".format(i + 1) + " / " + "{0:.0f}".format(len(B_array)))
gamma_dos_max=0,
factor_arcs=1)
hPocketCondObject.Ntime = 1000 # better for high magnetic field values
## Transport coeffcients -------------
## Empty arrays
rhoxx_array = np.empty_like(B_array, dtype=np.float64)
rhoxy_array = np.empty_like(B_array, dtype=np.float64)
rhozz_array = np.empty_like(B_array, dtype=np.float64)
RH_array = np.empty_like(B_array, dtype=np.float64)
for i, B in enumerate(B_array):
hPocketCondObject.Bamp = B
hPocketCondObject.solveMovementFunc()
hPocketCondObject.chambersFunc(0, 0)
hPocketCondObject.chambersFunc(0, 1)
hPocketCondObject.chambersFunc(2, 2)
sigma_xx = hPocketCondObject.sigma[0, 0]
sigma_xy = hPocketCondObject.sigma[0, 1]
sigma_zz = hPocketCondObject.sigma[2, 2]
rhoxx = sigma_xx / (sigma_xx**2 + sigma_xy**2) # Ohm.m
rhoxy = sigma_xy / (sigma_xx**2 + sigma_xy**2) # Ohm.m
rhozz = 1 / sigma_zz # Ohm.m
RH = rhoxy / B # m^3/C
rhoxx_array[i] = rhoxx
rhoxy_array[i] = rhoxy
rhozz_array[i] = rhozz