# SET PARAMETERS AND LIMITS OF INTEGRATION
kxi = -math.pi / a
kxf = math.pi / a
kyi = -math.pi / a
kyf = math.pi / a
resultArr = np.zeros(spacing)
errorArr = np.zeros(spacing)
timeArr = np.zeros(spacing)
j = 0
for i in np.linspace(qinitial, qfinal, spacing):
helpers.setqx(i)
MC = ZMCIntegral.MCintegral(modDsN2, [[kxi, kxf], [kyi, kyf]])
# Setting the zmcintegral parameters
MC.depth = depths
MC.sigma_multiplication = sigmults
MC.num_trials = trials
start = time.time()
result = MC.evaluate()
print('Result for qx = ', i, ': ', result[0], ' with error: ', result[1])
print('================================================================')
end = time.time()
print('Computed in ', end - start, ' seconds.')
print('================================================================')
resultArr[j] = result[0]
errorArr[j] = result[1]
timeArr[j] = end - start
# SET PARAMETERS AND LIMITS OF INTEGRATION
kxi = -math.pi / a
kxf = math.pi / a
kyi = -math.pi / a
kyf = math.pi / a
# qx = [0.001,6pi/a]
resultArr = np.zeros(100)
errorArr = np.zeros(100)
timeArr = np.zeros(100)
j = 0
for i in np.linspace(.01, .785, 100):
cudahelpers.setqx(i)
MC = ZMCIntegral.MCintegral(modDsN2, [[kxi, kxf], [kyi, kyf]])
# Setting the zmcintegral parameters
MC.depth = depths
MC.sigma_multiplication = sigmults
MC.num_trials = trials
start = time.time()
result = MC.evaluate()
print('Result for qx = ', i, ': ', result[0], ' with error: ', result[1])
print('================================================================')
end = time.time()
print('Computed in ', end - start, ' seconds.')
print('================================================================')
resultArr[j] = result[0]
errorArr[j] = result[1]
timeArr[j] = end - start
print('kxi = - math.pi / a')
print('kxf = math.pi / a')
print('kyi = - math.pi / a')
print('kyf = math.pi / a')
print('\n========================================================')
# qx = getqx()
# SET PARAMETERS AND LIMITS OF INTEGRATION
kxi = -math.pi / a
kxf = math.pi / a
kyi = -math.pi / a
kyf = math.pi / a
# qx = [0.001,6pi/a]
cudahelpers.setqx(math.pi / 4 - 0.05)
MC = ZMCIntegral.MCintegral(modDsN2, [[kxi, kxf], [kyi, kyf]])
# Setting the zmcintegral parameters
MC.depth = 2
MC.sigma_multiplication = 1E6
MC.num_trials = 5
start = time.time()
result = MC.evaluate()
print('Result for qx: ', result[0], ' with error: ', result[1])
print('================================================================')
end = time.time()
print('Computed in ', end - start, ' seconds.')
print('================================================================')