Changeset 22194b1 in sasview for src/sas/sascalc/corfunc

Timestamp:

Jul 17, 2017 7:53:16 AM (7 years ago)

Author:

lewis

Branches:

master, ESS_GUI, ESS_GUI_Docs, ESS_GUI_batch_fitting, ESS_GUI_bumps_abstraction, ESS_GUI_iss1116, ESS_GUI_iss879, ESS_GUI_iss959, ESS_GUI_opencl, ESS_GUI_ordering, ESS_GUI_sync_sascalc, magnetic_scatt, release-4.2.2, ticket-1009, ticket-1094-headless, ticket-1242-2d-resolution, ticket-1243, ticket-1249, ticket885, unittest-saveload

Children:

7cde638

Parents:

2a54ba5

Message:

Only calculate gamma3 up to x=200, as this is all that's plotted

File:

• src/sas/sascalc/corfunc/transform_thread.py (modified) (2 diffs)

src/sas/sascalc/corfunc/transform_thread.py

@@ -35 +35 @@
             # gamma3(R) = 1/R int_{0}^{R} gamma1(x) dx
             # trapz uses the trapezium rule to calculate the integral
-            gamma3 = [trapz(gamma1[:n], xs[:n])/xs[n-1] for n in range(1, len(xs+1))]
+            mask = xs <= 200.0 # Only calculate gamma3 up to x=200 (as this is all that's plotted)
+            gamma3 = [trapz(gamma1[:n], xs[:n])/xs[n-1] for n in range(2, len(xs[mask]) + 1)]
+            gamma3.insert(0, 1.0) # Gamma_3(0) is defined as 1
             gamma3 = np.array(gamma3)
         except Exception as e:
@@ -50 +52 @@
 
         transform1 = Data1D(xs, gamma1)
-        transform3 = Data1D(xs, gamma3)
+        transform3 = Data1D(xs[xs <= 200], gamma3)
 
         transforms = (transform1, transform3)