# 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('================================================================')