Improve GPU performance of adaptive integration on Mac M2

[pkienzle]

GPU performance on mac M2 is bad compared to cpu for small models, though large models see some benefit.

Speed check cpu vs gpu on mac ocl, after setting gpu_speed_check and speed_only in explore/check_adaptive.py:

$ python explore/check_adaptive.py          
== Speed and accuracy tests for all adaptive integration models ==
* target evaluation time is 2 s (running on a mac M2 chip)
*     q in [1e-5, 1] with 40 points per decade for 201 points total
*     warns if the adaptive model is 2x slower than a 76-point gaussian
* large models tested against 5000 point gaussian integration
*     q=[5e-4, 1e-3, 2e-3] with tol=1e-5 relative (measured q)
*     q=[0.01, 0.1] with tol=0.2 relative (slit resolution limits)
* small models tested against 5000 point gaussian integration
*     q in [1e-3, 1] with 1 points per decade
!!!! These tests run very slowly --- don't use as part of CI !!!!



=== small rods: a=20 b=40 c=200 ===
core_shell_bicelle background=0 radius=18.0 thick_rim=2.0 thick_face=2.0 length=196.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0
  cpu double: 0.5 ms  gpu single: 3.9 ms  [*** 6.1x slow down]
core_shell_cylinder background=0 radius=18.0 thickness=2.0 length=196.0 sld_core=0 sld_shell=1 sld_solvent=0
  cpu double: 0.3 ms  gpu single: 3.9 ms  [*** 12.0x slow down]
core_shell_ellipsoid background=0 radius_equat_core=18.0 x_core=5.444444444444445 thick_shell=2.0 x_polar_shell=1 sld_core=0 sld_shell=1 sld_solvent=0
  cpu double: 0.2 ms  gpu single: 2.2 ms  [*** 11.2x slow down]
cylinder background=0 sld=1 sld_solvent=0 radius=20.0 length=200
  cpu double: 0.2 ms  gpu single: 2.6 ms  [*** 11.1x slow down]
ellipsoid background=0 sld=1 sld_solvent=0 radius_equatorial=20.0 radius_polar=100.0
  cpu double: 0.1 ms  gpu single: 1.6 ms  [*** 12.8x slow down]
flexible_cylinder background=0 length=2000 kuhn_length=200 radius=40 sld=1 sld_solvent=0
  cpu double: 0.1 ms  gpu single: 0.7 ms  [*** 7.1x slow down]
hollow_cylinder background=0 radius=18.0 thickness=2.0 length=200 sld=1 sld_solvent=0
  cpu double: 0.3 ms  gpu single: 3.1 ms  [*** 10.1x slow down]
core_shell_bicelle_elliptical_belt_rough background=0 radius=9.0 x_core=2.111111111111111 thick_rim=1.0 thick_face=1.0 length=198.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0 sigma=0.1
  cpu double: 6.8 ms  gpu single: 45.9 ms  [*** 5.7x slow down]
core_shell_bicelle_elliptical background=0 radius=9.0 x_core=2.111111111111111 thick_rim=1.0 thick_face=1.0 length=198.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0
  cpu double: 6.9 ms  gpu single: 45.8 ms  [*** 5.6x slow down]
core_shell_parallelepiped background=0 sld_core=0 sld_a=1 sld_b=1 sld_c=1 sld_solvent=0 length_a=18.0 length_b=38.0 length_c=198.0 thick_rim_a=1.0 thick_rim_b=1.0 thick_rim_c=1.0
  cpu double: 4.8 ms  gpu single: 41.4 ms  [*** 7.7x slow down]
elliptical_cylinder background=0 radius_minor=10.0 axis_ratio=0.2 length=200 sld=1 sld_solvent=0
  cpu double: 2.5 ms  gpu single: 29.2 ms  [*** 10.5x slow down]
hollow_rectangular_prism_thin_walls background=0 sld=1 sld_solvent=0 length_a=20 b2a_ratio=2.0 c2a_ratio=10.0
  cpu double: 3.0 ms  gpu single: 24.5 ms  [*** 7.1x slow down]
hollow_rectangular_prism background=0 sld=1 sld_solvent=0 length_a=20 b2a_ratio=2.0 c2a_ratio=10.0 thickness=1.0
  cpu double: 4.2 ms  gpu single: 39.5 ms  [*** 8.3x slow down]
parallelepiped background=0 sld=1 sld_solvent=0 length_a=20 length_b=40 length_c=200
  cpu double: 2.9 ms  gpu single: 25.8 ms  [*** 7.8x slow down]
rectangular_prism background=0 sld=1 sld_solvent=0 length_a=20 b2a_ratio=2.0 c2a_ratio=10.0
  cpu double: 2.8 ms  gpu single: 24.4 ms  [*** 7.6x slow down]


=== small disks: a=180 b=200 c=40 ===
core_shell_ellipsoid background=0 radius_equat_core=90.0 x_core=0.1111111111111111 thick_shell=10.0 x_polar_shell=1 sld_core=0 sld_shell=1 sld_solvent=0
  cpu double: 0.2 ms  gpu single: 1.2 ms  [*** 6.1x slow down]
cylinder background=0 sld=1 sld_solvent=0 radius=100.0 length=40
  cpu double: 0.2 ms  gpu single: 1.3 ms  [*** 5.6x slow down]
ellipsoid background=0 sld=1 sld_solvent=0 radius_equatorial=100.0 radius_polar=20.0
  cpu double: 0.1 ms  gpu single: 0.9 ms  [*** 8.3x slow down]
flexible_cylinder background=0 length=400 kuhn_length=40 radius=200 sld=1 sld_solvent=0
  cpu double: 0.1 ms  gpu single: 0.6 ms  [*** 7.3x slow down]
triaxial_ellipsoid background=0 sld=1 sld_solvent=0 radius_equat_minor=180 radius_equat_major=200 radius_polar=40
  cpu double: 29.0 ms  gpu single: 110.6 ms  [*** 2.8x slow down]


=== small cubes: a=200 b=200 c=200 ===
core_shell_ellipsoid background=0 radius_equat_core=90.0 x_core=1.0 thick_shell=10.0 x_polar_shell=1 sld_core=0 sld_shell=1 sld_solvent=0
  cpu double: 0.2 ms  gpu single: 1.2 ms  [*** 6.6x slow down]
cylinder background=0 sld=1 sld_solvent=0 radius=100.0 length=200
  cpu double: 0.2 ms  gpu single: 1.3 ms  [*** 5.5x slow down]
ellipsoid background=0 sld=1 sld_solvent=0 radius_equatorial=100.0 radius_polar=100.0
  cpu double: 0.1 ms  gpu single: 0.9 ms  [*** 8.8x slow down]
flexible_cylinder background=0 length=2000 kuhn_length=200 radius=200 sld=1 sld_solvent=0
  cpu double: 0.1 ms  gpu single: 0.5 ms  [*** 7.2x slow down]
elliptical_cylinder background=0 radius_minor=100.0 axis_ratio=1.0 length=200 sld=1 sld_solvent=0
  cpu double: 11.0 ms  gpu single: 214.6 ms  [*** 18.5x slow down]
triaxial_ellipsoid background=0 sld=1 sld_solvent=0 radius_equat_minor=200 radius_equat_major=200 radius_polar=200
  cpu double: 30.2 ms  gpu single: 110.0 ms  [*** 2.6x slow down]


=== big rods: a=1000 b=2000 c=200000 ===
flexible_cylinder background=0 length=2000000 kuhn_length=200000 radius=2000 sld=1 sld_solvent=0
  cpu double: 0.1 ms  gpu single: 0.6 ms  [*** 8.1x slow down]
barbell background=0 radius=500.0 radius_bell=1000.0 length=196267.94919243112 sld=1 sld_solvent=0
  cpu double: 778.3 ms  gpu single: 334.8 ms  [56% speed up]
capped_cylinder background=0 radius=500.0 radius_cap=1000.0 length=199732.05080756888 sld=1 sld_solvent=0
  cpu double: 771.3 ms  gpu single: 319.3 ms  [58% speed up]
core_shell_bicelle_elliptical_belt_rough background=0 radius=450.0 x_core=2.111111111111111 thick_rim=50.0 thick_face=50.0 length=199900.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0 sigma=5.0
  cpu double: 1297.7 ms  gpu single: 549.2 ms  [57% speed up]
core_shell_bicelle_elliptical background=0 radius=450.0 x_core=2.111111111111111 thick_rim=50.0 thick_face=50.0 length=199900.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0
  cpu double: 1302.1 ms  gpu single: 549.4 ms  [57% speed up]
core_shell_parallelepiped background=0 sld_core=0 sld_a=1 sld_b=1 sld_c=1 sld_solvent=0 length_a=900.0 length_b=1900.0 length_c=199900.0 thick_rim_a=50.0 thick_rim_b=50.0 thick_rim_c=50.0
  cpu double: 891.7 ms  gpu single: 497.2 ms  [44% speed up]
elliptical_cylinder background=0 radius_minor=500.0 axis_ratio=0.01 length=200000 sld=1 sld_solvent=0
  cpu double: 599.2 ms  gpu single: 335.3 ms  [44% speed up]
hollow_rectangular_prism background=0 sld=1 sld_solvent=0 length_a=1000 b2a_ratio=2.0 c2a_ratio=200.0 thickness=50.0
  cpu double: 778.6 ms  gpu single: 485.6 ms  [37% speed up]
parallelepiped background=0 sld=1 sld_solvent=0 length_a=1000 length_b=2000 length_c=200000
  cpu double: 546.8 ms  gpu single: 309.3 ms  [43% speed up]
rectangular_prism background=0 sld=1 sld_solvent=0 length_a=1000 b2a_ratio=2.0 c2a_ratio=200.0
  cpu double: 512.8 ms  gpu single: 296.2 ms  [42% speed up]


=== big disks: a=180000 b=200000 c=1000 ===
flexible_cylinder background=0 length=10000 kuhn_length=1000 radius=200000 sld=1 sld_solvent=0
  cpu double: 0.1 ms  gpu single: 0.6 ms  [*** 8.1x slow down]
barbell background=0 radius=90000.0 radius_bell=100000.0 length=0 sld=1 sld_solvent=0
  cpu double: 1046.0 ms  gpu single: 337.2 ms  [67% speed up]
capped_cylinder background=0 radius=90000.0 radius_cap=100000.0 length=0 sld=1 sld_solvent=0
  cpu double: 936.5 ms  gpu single: 337.3 ms  [63% speed up]
core_shell_bicelle_elliptical_belt_rough background=0 radius=81000.0 x_core=1.123456790123457 thick_rim=9000.0 thick_face=9000.0 length=-17000.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0 sigma=900.0
  cpu double: 1874.4 ms  gpu single: 836.5 ms  [55% speed up]
core_shell_bicelle_elliptical background=0 radius=81000.0 x_core=1.123456790123457 thick_rim=9000.0 thick_face=9000.0 length=-17000.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0
  cpu double: 1874.6 ms  gpu single: 839.9 ms  [55% speed up]
core_shell_parallelepiped background=0 sld_core=0 sld_a=1 sld_b=1 sld_c=1 sld_solvent=0 length_a=162000.0 length_b=182000.0 length_c=-17000.0 thick_rim_a=9000.0 thick_rim_b=9000.0 thick_rim_c=9000.0
  cpu double: 1721.5 ms  gpu single: 749.8 ms  [56% speed up]
elliptical_cylinder background=0 radius_minor=90000.0 axis_ratio=200.0 length=1000 sld=1 sld_solvent=0
! ** elliptical_cylinder is slow: 2.2 s for 201 points in [1e-05, 1.0]
  cpu double: 2225.0 ms  gpu single: 523.9 ms  [76% speed up]
hollow_rectangular_prism_thin_walls background=0 sld=1 sld_solvent=0 length_a=180000 b2a_ratio=1.1111111111111112 c2a_ratio=0.005555555555555556
  cpu double: 883.0 ms  gpu single: 453.4 ms  [48% speed up]
hollow_rectangular_prism background=0 sld=1 sld_solvent=0 length_a=180000 b2a_ratio=1.1111111111111112 c2a_ratio=0.005555555555555556 thickness=9000.0
  cpu double: 1633.8 ms  gpu single: 715.9 ms  [56% speed up]
parallelepiped background=0 sld=1 sld_solvent=0 length_a=180000 length_b=200000 length_c=1000
  cpu double: 1113.5 ms  gpu single: 463.5 ms  [58% speed up]
rectangular_prism background=0 sld=1 sld_solvent=0 length_a=180000 b2a_ratio=1.1111111111111112 c2a_ratio=0.005555555555555556
  cpu double: 1075.5 ms  gpu single: 437.0 ms  [59% speed up]
triaxial_ellipsoid background=0 sld=1 sld_solvent=0 radius_equat_minor=180000 radius_equat_major=200000 radius_polar=1000
  cpu double: 600.6 ms  gpu single: 259.9 ms  [56% speed up]


=== big cubes: a=200000 b=200000 c=200000 ===
flexible_cylinder background=0 length=2000000 kuhn_length=200000 radius=200000 sld=1 sld_solvent=0
  cpu double: 0.1 ms  gpu single: 0.6 ms  [*** 8.2x slow down]
barbell background=0 radius=100000.0 radius_bell=100000.0 length=0.0 sld=1 sld_solvent=0
  cpu double: 1000.5 ms  gpu single: 336.3 ms  [66% speed up]
capped_cylinder background=0 radius=100000.0 radius_cap=100000.0 length=0.0 sld=1 sld_solvent=0
  cpu double: 1000.8 ms  gpu single: 335.6 ms  [66% speed up]
core_shell_bicelle_elliptical_belt_rough background=0 radius=90000.0 x_core=1.0 thick_rim=10000.0 thick_face=10000.0 length=180000.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0 sigma=1000.0
  cpu double: 1765.2 ms  gpu single: 835.1 ms  [52% speed up]
core_shell_bicelle_elliptical background=0 radius=90000.0 x_core=1.0 thick_rim=10000.0 thick_face=10000.0 length=180000.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0
  cpu double: 1764.4 ms  gpu single: 839.7 ms  [52% speed up]
core_shell_parallelepiped background=0 sld_core=0 sld_a=1 sld_b=1 sld_c=1 sld_solvent=0 length_a=180000.0 length_b=180000.0 length_c=180000.0 thick_rim_a=10000.0 thick_rim_b=10000.0 thick_rim_c=10000.0
  cpu double: 1726.7 ms  gpu single: 751.9 ms  [56% speed up]
elliptical_cylinder background=0 radius_minor=100000.0 axis_ratio=1.0 length=200000 sld=1 sld_solvent=0
  cpu double: 912.3 ms  gpu single: 523.2 ms  [42% speed up]
hollow_rectangular_prism_thin_walls background=0 sld=1 sld_solvent=0 length_a=200000 b2a_ratio=1.0 c2a_ratio=1.0
  cpu double: 888.5 ms  gpu single: 453.8 ms  [48% speed up]
hollow_rectangular_prism background=0 sld=1 sld_solvent=0 length_a=200000 b2a_ratio=1.0 c2a_ratio=1.0 thickness=10000.0
  cpu double: 1641.9 ms  gpu single: 714.1 ms  [56% speed up]
parallelepiped background=0 sld=1 sld_solvent=0 length_a=200000 length_b=200000 length_c=200000
  cpu double: 1112.4 ms  gpu single: 464.0 ms  [58% speed up]
rectangular_prism background=0 sld=1 sld_solvent=0 length_a=200000 b2a_ratio=1.0 c2a_ratio=1.0
  cpu double: 1076.6 ms  gpu single: 436.8 ms  [59% speed up]

Originally posted by @pkienzle in #658 (comment)

[pkienzle]

A speed comparison of gpu vs dll for fixed gauss-76 quadrature also shows a performance penalty for fast models.

[krzywon]

If I'm reading this correctly, the vast majority of the models that are slower are models that were taking tenths of milliseconds to calculate, but now take a few milliseconds. Is my understanding correct here? My big question - will anyone notice this?

Side note: I saw these two errors while running check_adaptive.py. I'm guessing these are divide by 0 errors, but not certain. I was dumping the output to a text file, but these showed up in the prompt, so I don't know exactly what models and what loops they happened.

\sasmodels\explore\check_adaptive.py:325: RuntimeWarning: invalid value encountered in divide
relerr = abs(Iq_test_adaptive - Iq_test_accurate)/Iq_test_accurate
\sasmodels\explore\check_adaptive.py:337: RuntimeWarning: invalid value encountered in divide
relerr_76 = abs(Iq_test_76 - Iq_test_accurate)/Iq_test_accurate

================================================
No GPU, and no options given to check_adaptive.py:

== Speed and accuracy tests for all adaptive integration models ==

• target evaluation time is 2 s (running on a mac M2 chip)
• q in [1e-5, 1] with 40 points per decade for 201 points total
• large models tested against 5000 point gaussian integration
• q=[5e-4, 1e-3, 2e-3] with tol=1e-5 relative (measured q)
• q=[0.01, 0.1] with tol=0.2 relative (slit resolution limits)
• small models tested against 5000 point gaussian integration
• q in [1e-3, 1] with 1 points per decade
  !!!! These tests run very slowly --- don't use as part of CI !!!!

=== small rods: a=20 b=40 c=200 ===
core_shell_bicelle background=0 radius=18.0 thick_rim=2.0 thick_face=2.0 length=196.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [1.19946526073276e-11]
rel-76 [2.07980488165952e-06]
core_shell_cylinder background=0 radius=18.0 thickness=2.0 length=196.0 sld_core=0 sld_shell=1 sld_solvent=0
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [1.19983558035632e-11]
rel-76 [2.07980488240016e-06]
cylinder background=0 sld=1 sld_solvent=0 radius=20.0 length=200
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [8.76342261326569e-12]
rel-76 [5.10725132142934e-04]
stacked_disks background=0 thick_core=18.0 thick_layer=1.0 radius=20.0 n_stacking=10 sigma_d=1.0 sld_core=0 sld_layer=1 sld_solvent=0
cpu double: 1.8 ms gauss-76: 5.3 ms [gauss-76 2.9x slower]
qvalue [1.00000000000000e-01]
target [1.86796449661117e-02]
actual [1.86796449683030e-02]
relerr [1.17313926260878e-10]
barbell background=0 radius=10.0 radius_bell=20.0 length=125.35898384862244 sld=1 sld_solvent=0
cpu double: 12.4 ms gauss-76: 48.7 ms [gauss-76 3.9x slower]
capped_cylinder background=0 radius=10.0 radius_cap=20.0 length=194.64101615137756 sld=1 sld_solvent=0
cpu double: 17.2 ms gauss-76: 80.2 ms [gauss-76 4.7x slower]
qvalue [1.00000000000000e+00]
target [1.60757256812769e-05]
actual [1.60759828022906e-05]
relerr [1.59943643466140e-05]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [1.59943643466140e-05]
rel-76 [1.59943643466140e-05]
core_shell_bicelle_elliptical_belt_rough background=0 radius=9.0 x_core=2.111111111111111 thick_rim=1.0 thick_face=1.0 length=198.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0 sigma=0.1
cpu double: 22.3 ms gauss-76: 87.8 ms [gauss-76 3.9x slower]
core_shell_bicelle_elliptical background=0 radius=9.0 x_core=2.111111111111111 thick_rim=1.0 thick_face=1.0 length=198.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0
cpu double: 19.1 ms gauss-76: 92.7 ms [gauss-76 4.8x slower]
core_shell_parallelepiped background=0 sld_core=0 sld_a=1 sld_b=1 sld_c=1 sld_solvent=0 length_a=18.0 length_b=38.0 length_c=198.0 thick_rim_a=1.0 thick_rim_b=1.0 thick_rim_c=1.0
cpu double: 10.2 ms gauss-76: 46.1 ms [gauss-76 4.5x slower]
elliptical_cylinder background=0 radius_minor=10.0 axis_ratio=0.2 length=200 sld=1 sld_solvent=0
cpu double: 9.1 ms gauss-76: 55.5 ms [gauss-76 6.1x slower]
flexible_cylinder_elliptical background=0 length=2000 kuhn_length=200 radius=20 axis_ratio=2.0 sld=1 sld_solvent=0
cpu double: 0.5 ms gauss-76: 1.3 ms [gauss-76 2.5x slower]
hollow_rectangular_prism_thin_walls background=0 sld=1 sld_solvent=0 length_a=20 b2a_ratio=2.0 c2a_ratio=10.0
cpu double: 8.6 ms gauss-76: 40.9 ms [gauss-76 4.8x slower]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [1.36395284619967e-11]
rel-76 [1.28327955419546e-05]
hollow_rectangular_prism background=0 sld=1 sld_solvent=0 length_a=20 b2a_ratio=2.0 c2a_ratio=10.0 thickness=1.0
cpu double: 11.0 ms gauss-76: 46.4 ms [gauss-76 4.2x slower]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [1.38398876615951e-11]
rel-76 [9.42506049129388e-06]
parallelepiped background=0 sld=1 sld_solvent=0 length_a=20 length_b=40 length_c=200
cpu double: 6.6 ms gauss-76: 30.9 ms [gauss-76 4.7x slower]
pringle background=0 radius=20.0 thickness=200 alpha=0.5 beta=0.5 sld=1 sld_solvent=0
cpu double: 596.6 ms gauss-76: 1377.9 ms [gauss-76 2.3x slower]
qvalue [1.00000000000000e+00]
target [8.94852694171988e-05]
actual [8.90702124575545e-05]
relerr [4.63827133054966e-03]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [4.63827133054966e-03]
rel-76 [4.17405432356878e-03]
rectangular_prism background=0 sld=1 sld_solvent=0 length_a=20 b2a_ratio=2.0 c2a_ratio=10.0
cpu double: 6.2 ms gauss-76: 30.8 ms [gauss-76 4.9x slower]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [1.28765055306283e-11]
rel-76 [4.83452757762925e-05]
triaxial_ellipsoid background=0 sld=1 sld_solvent=0 radius_equat_minor=20 radius_equat_major=40 radius_polar=200
qvalue [1.00000000000000e+00]
target [5.55750011411866e-05]
actual [5.52135920592964e-05]
relerr [6.50308725990077e-03]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [6.50308725990077e-03]
rel-76 [6.50308725990430e-03]

=== small disks: a=180 b=200 c=40 ===
core_shell_bicelle background=0 radius=90.0 thick_rim=10.0 thick_face=10.0 length=20.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [6.46895700428258e-13]
rel-76 [8.21234112102640e-05]
core_shell_cylinder background=0 radius=90.0 thickness=10.0 length=20.0 sld_core=0 sld_shell=1 sld_solvent=0
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [6.46733286242863e-13]
rel-76 [8.21234112089647e-05]
core_shell_ellipsoid background=0 radius_equat_core=90.0 x_core=0.1111111111111111 thick_shell=10.0 x_polar_shell=1 sld_core=0 sld_shell=1 sld_solvent=0
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [1.41044582771822e-11]
rel-76 [2.47407320458531e-08]
cylinder background=0 sld=1 sld_solvent=0 radius=100.0 length=40
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [6.93315921756591e-13]
rel-76 [1.77239587494715e-04]
ellipsoid background=0 sld=1 sld_solvent=0 radius_equatorial=100.0 radius_polar=20.0
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [9.98283109190283e-13]
rel-76 [1.26418472150758e-08]
hollow_cylinder background=0 radius=90.0 thickness=10.0 length=40 sld=1 sld_solvent=0
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [1.26431292541489e-11]
rel-76 [3.48312044020518e-05]
stacked_disks background=0 thick_core=3.6 thick_layer=0.2 radius=100.0 n_stacking=10 sigma_d=0.2 sld_core=0 sld_layer=1 sld_solvent=0
cpu double: 2.6 ms gauss-76: 5.3 ms [gauss-76 2.0x slower]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [5.81201116980365e-13]
rel-76 [1.90161439457918e-06]
barbell background=0 radius=90.0 radius_bell=100.0 length=0 sld=1 sld_solvent=0
cpu double: 25.4 ms gauss-76: 52.1 ms [gauss-76 2.1x slower]
capped_cylinder background=0 radius=90.0 radius_cap=100.0 length=0 sld=1 sld_solvent=0
cpu double: 15.7 ms gauss-76: 47.1 ms [gauss-76 3.0x slower]
qvalue [1.00000000000000e+00]
target [2.75569192104124e-05]
actual [2.75569192194923e-05]
relerr [3.29498227142533e-10]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [3.29498227142533e-10]
rel-76 [3.29498227142533e-10]
core_shell_bicelle_elliptical_belt_rough background=0 radius=81.0 x_core=1.123456790123457 thick_rim=9.0 thick_face=9.0 length=22.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0 sigma=0.9
cpu double: 38.6 ms gauss-76: 97.2 ms [gauss-76 2.5x slower]
qvalue [1.00000000000000e+00]
target [1.57161805207409e-04]
actual [1.57161805250852e-04]
relerr [2.76420605290949e-10]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [2.76420605290949e-10]
rel-76 [9.77355925821890e-07]
core_shell_bicelle_elliptical background=0 radius=81.0 x_core=1.123456790123457 thick_rim=9.0 thick_face=9.0 length=22.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0
cpu double: 38.1 ms gauss-76: 96.0 ms [gauss-76 2.5x slower]
qvalue [1.00000000000000e+00]
target [2.35631493977476e-04]
actual [2.35631494069685e-04]
relerr [3.91324540503391e-10]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [3.91324540503391e-10]
rel-76 [3.43046688236032e-06]
core_shell_parallelepiped background=0 sld_core=0 sld_a=1 sld_b=1 sld_c=1 sld_solvent=0 length_a=162.0 length_b=182.0 length_c=22.0 thick_rim_a=9.0 thick_rim_b=9.0 thick_rim_c=9.0
qvalue [1.00000000000000e+00]
target [2.42326923165195e-04]
actual [2.42333667936020e-04]
relerr [2.78333531259064e-05]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [2.78333531259064e-05]
rel-76 [3.47658398048087e-03]
elliptical_cylinder background=0 radius_minor=90.0 axis_ratio=5.0 length=40 sld=1 sld_solvent=0
qvalue [1.00000000000000e+00]
target [6.12554808965977e-05]
actual [6.02643927068993e-05]
relerr [1.61795838542423e-02]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [1.61795838542423e-02]
rel-76 [3.73519475863190e-01]
hollow_rectangular_prism_thin_walls background=0 sld=1 sld_solvent=0 length_a=180 b2a_ratio=1.1111111111111112 c2a_ratio=0.2222222222222222
cpu double: 22.6 ms gauss-76: 47.8 ms [gauss-76 2.1x slower]
qvalue [1.00000000000000e-01 1.00000000000000e+00]
target [2.75787141463231e-02 3.07578043101305e-04]
actual [2.75787141752361e-02 3.07450050065018e-04]
relerr [1.04838366179869e-09 4.16131902642557e-04]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [4.16131902642557e-04]
rel-76 [4.81433600979508e-04]
hollow_rectangular_prism background=0 sld=1 sld_solvent=0 length_a=180 b2a_ratio=1.1111111111111112 c2a_ratio=0.2222222222222222 thickness=9.0
qvalue [1.00000000000000e-01 1.00000000000000e+00]
target [1.60102657069035e-01 4.21190511154586e-04]
actual [1.60102657153591e-01 4.21207536067096e-04]
relerr [5.28137364655840e-10 4.04209308117173e-05]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [4.04209308117173e-05]
rel-76 [6.76748214026063e-05]
parallelepiped background=0 sld=1 sld_solvent=0 length_a=180 length_b=200 length_c=40
qvalue [1.00000000000000e-01 1.00000000000000e+00]
target [6.54471658075910e-01 6.73870249389426e-05]
actual [6.54471658146499e-01 6.74084672293596e-05]
relerr [1.07856264229397e-10 3.18196128059957e-04]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [3.18196128059957e-04]
rel-76 [7.23380119752346e-02]
pringle background=0 radius=100.0 thickness=40 alpha=0.09999999999999998 beta=0.09999999999999998 sld=1 sld_solvent=0
cpu double: 1720.4 ms gauss-76: 1374.0 ms [gauss-76 20% faster] ***
qvalue [1.00000000000000e-01 1.00000000000000e+00]
target [2.18733343077074e-01 6.83106606401077e-05]
actual [2.18755844943364e-01 6.85148995490768e-05]
relerr [1.02873507870985e-04 2.98985410264337e-03]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e-01 1.00000000000000e+00]
relerr [1.02873507870985e-04 2.98985410264337e-03]
rel-76 [1.02873484433851e-04 1.83251623778961e-03]
rectangular_prism background=0 sld=1 sld_solvent=0 length_a=180 b2a_ratio=1.1111111111111112 c2a_ratio=0.2222222222222222
qvalue [1.00000000000000e-01 1.00000000000000e+00]
target [6.54471658079303e-01 6.73870249688472e-05]
actual [6.54471658146502e-01 6.74084672274739e-05]
relerr [1.02675734173720e-10 3.18195656160577e-04]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [3.18195656160577e-04]
rel-76 [3.82482537489621e-04]
triaxial_ellipsoid background=0 sld=1 sld_solvent=0 radius_equat_minor=180 radius_equat_major=200 radius_polar=40
qvalue [1.00000000000000e+00]
target [2.65781861041188e-05]
actual [2.63440364974738e-05]
relerr [8.80984148909772e-03]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [8.80984148909772e-03]
rel-76 [6.91751920168694e-03]

=== small cubes: a=200 b=200 c=200 ===
core_shell_bicelle background=0 radius=90.0 thick_rim=10.0 thick_face=10.0 length=180.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [9.54561141906767e-12]
rel-76 [2.55391950531571e-06]
core_shell_cylinder background=0 radius=90.0 thickness=10.0 length=180.0 sld_core=0 sld_shell=1 sld_solvent=0
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [9.54891714403618e-12]
rel-76 [2.55391950491902e-06]
cylinder background=0 sld=1 sld_solvent=0 radius=100.0 length=200
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [1.06394803378813e-11]
rel-76 [2.26887671016852e-05]
hollow_cylinder background=0 radius=90.0 thickness=10.0 length=200 sld=1 sld_solvent=0
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [8.51956596890793e-12]
rel-76 [4.21961053386673e-03]
stacked_disks background=0 thick_core=18.0 thick_layer=1.0 radius=100.0 n_stacking=10 sigma_d=1.0 sld_core=0 sld_layer=1 sld_solvent=0
cpu double: 2.9 ms gauss-76: 5.8 ms [gauss-76 2.0x slower]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [3.07851592739912e-14]
rel-76 [1.19994354638996e-05]
barbell background=0 radius=100.0 radius_bell=100.0 length=0.0 sld=1 sld_solvent=0
cpu double: 18.5 ms gauss-76: 50.1 ms [gauss-76 2.7x slower]
capped_cylinder background=0 radius=100.0 radius_cap=100.0 length=0.0 sld=1 sld_solvent=0
cpu double: 17.9 ms gauss-76: 49.0 ms [gauss-76 2.7x slower]
core_shell_bicelle_elliptical_belt_rough background=0 radius=90.0 x_core=1.0 thick_rim=10.0 thick_face=10.0 length=180.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0 sigma=1.0
cpu double: 38.8 ms gauss-76: 97.5 ms [gauss-76 2.5x slower]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [6.98446614112596e-11]
rel-76 [1.16443555406110e-04]
core_shell_bicelle_elliptical background=0 radius=90.0 x_core=1.0 thick_rim=10.0 thick_face=10.0 length=180.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0
cpu double: 40.0 ms gauss-76: 105.2 ms [gauss-76 2.6x slower]
qvalue [1.00000000000000e+00]
target [1.02492851673480e-04]
actual [1.02492851684422e-04]
relerr [1.06759181241761e-10]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [1.06759181241761e-10]
rel-76 [9.31522529271066e-04]
core_shell_parallelepiped background=0 sld_core=0 sld_a=1 sld_b=1 sld_c=1 sld_solvent=0 length_a=180.0 length_b=180.0 length_c=180.0 thick_rim_a=10.0 thick_rim_b=10.0 thick_rim_c=10.0
qvalue [1.00000000000000e+00]
target [6.98345673054341e-05]
actual [7.00984945279586e-05]
relerr [3.77932065319739e-03]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [3.77932065319739e-03]
rel-76 [7.17822024108545e-03]
elliptical_cylinder background=0 radius_minor=100.0 axis_ratio=1.0 length=200 sld=1 sld_solvent=0
cpu double: 25.6 ms gauss-76: 61.1 ms [gauss-76 2.4x slower]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [7.06608580851663e-11]
rel-76 [3.80598188626683e-03]
hollow_rectangular_prism_thin_walls background=0 sld=1 sld_solvent=0 length_a=200 b2a_ratio=1.0 c2a_ratio=1.0
cpu double: 23.0 ms gauss-76: 48.9 ms [gauss-76 2.1x slower]
qvalue [1.00000000000000e-01 1.00000000000000e+00]
target [1.04904095454156e-01 8.56430787049371e-04]
actual [1.04904094452217e-01 8.56144929699485e-04]
relerr [9.55099076047659e-09 3.33777526694336e-04]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [3.33777526694336e-04]
rel-76 [9.36321746796359e-04]
hollow_rectangular_prism background=0 sld=1 sld_solvent=0 length_a=200 b2a_ratio=1.0 c2a_ratio=1.0 thickness=10.0
qvalue [1.00000000000000e-01 1.00000000000000e+00]
target [9.82008425165988e-01 2.42978141719551e-04]
actual [9.82008421794487e-01 2.42955329124303e-04]
relerr [3.43327069183032e-09 9.38874381322809e-05]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [9.38874381322809e-05]
rel-76 [2.99875418191674e-04]
parallelepiped background=0 sld=1 sld_solvent=0 length_a=200 length_b=200 length_c=200
qvalue [1.00000000000000e-01 1.00000000000000e+00]
target [8.65241732031059e-02 1.16766869284293e-05]
actual [8.65241716425500e-02 1.16947856496101e-05]
relerr [1.80360680858781e-08 1.54998770557011e-03]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [1.54998770557011e-03]
rel-76 [7.44548186547528e-03]
pringle background=0 radius=100.0 thickness=200 alpha=0.001 beta=0.001 sld=1 sld_solvent=0
cpu double: 308.6 ms gauss-76: 1304.1 ms [gauss-76 4.2x slower]
qvalue [1.00000000000000e+00]
target [2.87134261111240e-05]
actual [2.86855009494169e-05]
relerr [9.72547183990834e-04]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [9.72547183990834e-04]
rel-76 [9.72547204537502e-04]
rectangular_prism background=0 sld=1 sld_solvent=0 length_a=200 b2a_ratio=1.0 c2a_ratio=1.0
qvalue [1.00000000000000e-01 1.00000000000000e+00]
target [8.65241732033492e-02 1.16766869304387e-05]
actual [8.65241716425503e-02 1.16947856493809e-05]
relerr [1.80388768710616e-08 1.54998751358618e-03]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e+00]
relerr [1.54998751358618e-03]
rel-76 [4.41276355454061e-03]

=== big rods: a=1000 b=2000 c=200000 ===
core_shell_bicelle background=0 radius=900.0 thick_rim=100.0 thick_face=100.0 length=199800.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0
cpu double: 75.9 ms gauss-76: 2.1 ms [gauss-76 36.0x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [2.46280714784500e-11 2.66379987083313e-11]
rel-76 [5.54931903193811e-03 7.65455256606515e-01]
core_shell_cylinder background=0 radius=900.0 thickness=100.0 length=199800.0 sld_core=0 sld_shell=1 sld_solvent=0
cpu double: 68.6 ms gauss-76: 2.2 ms [gauss-76 31.5x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [2.46372023738518e-11 2.66383146538426e-11]
rel-76 [5.54931903193891e-03 7.65455256607306e-01]
core_shell_ellipsoid background=0 radius_equat_core=900.0 x_core=111.0 thick_shell=100.0 x_polar_shell=1 sld_core=0 sld_shell=1 sld_solvent=0
cpu double: 35.9 ms gauss-76: 1.3 ms [gauss-76 27.0x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03]
relerr [1.47794477540076e-11]
rel-76 [1.26042332853879e-04]
cylinder background=0 sld=1 sld_solvent=0 radius=1000.0 length=200000
cpu double: 44.8 ms gauss-76: 1.5 ms [gauss-76 30.2x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [2.48213959403765e-11 2.66750444957834e-11]
rel-76 [4.07974950245269e-03 7.64448118101773e-01]
ellipsoid background=0 sld=1 sld_solvent=0 radius_equatorial=1000.0 radius_polar=100000.0
cpu double: 15.8 ms gauss-76: 0.6 ms [gauss-76 24.7x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-02]
relerr [1.39260259171204e-11 2.40774440171201e-13]
rel-76 [5.57349012497014e-05 2.22995907445804e-01]
hollow_cylinder background=0 radius=900.0 thickness=100.0 length=200000 sld=1 sld_solvent=0
cpu double: 67.9 ms gauss-76: 1.9 ms [gauss-76 35.2x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [1.56084825521937e-11 1.67574931989118e-12]
rel-76 [4.08003973033052e-03 7.26051566017397e-01]
stacked_disks background=0 thick_core=18000.0 thick_layer=1000.0 radius=1000.0 n_stacking=10 sigma_d=1000.0 sld_core=0 sld_layer=1 sld_solvent=0
cpu double: 127.3 ms gauss-76: 8.2 ms [gauss-76 15.5x faster] ***
qvalue [5.00000000000000e-04 1.00000000000000e-03 2.00000000000000e-03]
target [2.96917938719083e+04 2.60914421213130e+04 1.42434747458768e+04]
actual [2.96507625109331e+04 2.60745206784249e+04 1.42398416282128e+04]
relerr [1.38190912789591e-03 6.48543794915104e-04 2.55072426417611e-04]
! ** gauss-76 is out of tolerance
qvalue [5.00000000000000e-04 1.00000000000000e-03 2.00000000000000e-03]
relerr [1.38190912789591e-03 6.48543794915104e-04 2.55072426417611e-04]
rel-76 [1.37999315295319e-03 6.48543794915104e-04 2.55072426417611e-04]
barbell background=0 radius=500.0 radius_bell=1000.0 length=196267.94919243112 sld=1 sld_solvent=0
! ** barbell is slow: 2.2 s for 201 points in [1e-05, 1.0]
cpu double: 2237.6 ms gauss-76: 64.4 ms [gauss-76 34.7x faster] ***
qvalue [2.00000000000000e-03 1.00000000000000e-01]
target [9.96690490148787e+04 3.69725105766636e-02]
actual [9.84699740235029e+04 8.97299419744593e-03]
relerr [1.20305651877625e-02 7.57306332258932e-01]
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [1.20305651877625e-02 7.57306332258932e-01]
rel-76 [1.20305651877610e-02 7.57306332258939e-01]
capped_cylinder background=0 radius=500.0 radius_cap=1000.0 length=199732.05080756888 sld=1 sld_solvent=0
cpu double: 1618.4 ms gauss-76: 57.1 ms [gauss-76 28.3x faster] ***
qvalue [2.00000000000000e-03 1.00000000000000e-01]
target [9.54509605532282e+04 3.75843968780453e-02]
actual [9.50819794058316e+04 8.75212961204768e-03]
relerr [3.86566196147295e-03 7.67133961456218e-01]
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [3.86566196147295e-03 7.67133961456218e-01]
rel-76 [3.86566196147295e-03 7.67133961456219e-01]
core_shell_bicelle_elliptical_belt_rough background=0 radius=450.0 x_core=2.111111111111111 thick_rim=50.0 thick_face=50.0 length=199900.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0 sigma=5.0
cpu double: 450.7 ms gauss-76: 120.3 ms [gauss-76 3.7x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [1.57765903724379e-11 4.69577776780822e-03]
rel-76 [4.49968767675039e-03 7.23613501072581e-01]
core_shell_bicelle_elliptical background=0 radius=450.0 x_core=2.111111111111111 thick_rim=50.0 thick_face=50.0 length=199900.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0
cpu double: 483.2 ms gauss-76: 126.6 ms [gauss-76 3.8x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [1.57862531629417e-11 4.86470905932948e-03]
rel-76 [4.56043282562372e-03 7.25405995437354e-01]
core_shell_parallelepiped background=0 sld_core=0 sld_a=1 sld_b=1 sld_c=1 sld_solvent=0 length_a=900.0 length_b=1900.0 length_c=199900.0 thick_rim_a=50.0 thick_rim_b=50.0 thick_rim_c=50.0
cpu double: 224.4 ms gauss-76: 64.3 ms [gauss-76 3.5x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [1.57715719685174e-11 1.39750606963823e-02]
rel-76 [4.60832083635136e-03 7.26804261235388e-01]
elliptical_cylinder background=0 radius_minor=500.0 axis_ratio=0.01 length=200000 sld=1 sld_solvent=0
cpu double: 230.8 ms gauss-76: 66.4 ms [gauss-76 3.5x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [1.57530244837658e-11 4.94371271269406e-03]
rel-76 [4.14476738914578e-03 7.25949013513855e-01]
flexible_cylinder_elliptical background=0 length=2000000 kuhn_length=200000 radius=1000 axis_ratio=2.0 sld=1 sld_solvent=0
cpu double: 10.3 ms gauss-76: 1.3 ms [gauss-76 7.8x faster] ***
hollow_rectangular_prism_thin_walls background=0 sld=1 sld_solvent=0 length_a=1000 b2a_ratio=2.0 c2a_ratio=200.0
cpu double: 160.1 ms gauss-76: 51.0 ms [gauss-76 3.1x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [2.48826284685329e-11 3.92623384057090e-03]
rel-76 [5.16259775023319e-03 7.65090746035799e-01]
hollow_rectangular_prism background=0 sld=1 sld_solvent=0 length_a=1000 b2a_ratio=2.0 c2a_ratio=200.0 thickness=50.0
cpu double: 213.2 ms gauss-76: 70.0 ms [gauss-76 3.0x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [2.48853464897084e-11 3.54831835887896e-03]
rel-76 [5.06685086364162e-03 7.65012913489626e-01]
parallelepiped background=0 sld=1 sld_solvent=0 length_a=1000 length_b=2000 length_c=200000
cpu double: 130.7 ms gauss-76: 38.3 ms [gauss-76 3.4x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [1.57426411813081e-11 3.33713279321374e-04]
rel-76 [4.16822479463203e-03 7.15187352718665e-01]
pringle background=0 radius=1000.0 thickness=200000 alpha=0.5 beta=0.5 sld=1 sld_solvent=0
! ** pringle is slow: 7.6 s for 201 points in [1e-05, 1.0]
cpu double: 7620.8 ms gauss-76: 1620.5 ms [gauss-76 4.7x faster] ***
qvalue [5.00000000000000e-04 1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
target [5.29285175436895e+05 2.25418078593827e+05 6.00092800236802e+04 5.04753936368984e+02 8.73018584972186e-02]
actual [5.29442866649060e+05 2.26829405624547e+05 6.03425922786133e+04 2.53341781736512e+02 6.97860160994661e-03]
relerr [2.97932418067012e-04 6.26093097556351e-03 5.55434517463874e-03 4.98088546750204e-01 9.20063538966138e-01]
! ** gauss-76 is out of tolerance
qvalue [5.00000000000000e-04 1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [2.97932418067012e-04 6.26093097556351e-03 5.55434517463874e-03 4.98088546750204e-01 9.20063538966138e-01]
rel-76 [2.89931312847352e-04 6.26609257914980e-03 5.54864494400059e-03 4.97985975463520e-01 9.12616592934624e-01]
rectangular_prism background=0 sld=1 sld_solvent=0 length_a=1000 b2a_ratio=2.0 c2a_ratio=200.0
cpu double: 125.5 ms gauss-76: 39.9 ms [gauss-76 3.1x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [2.48538847559773e-11 9.83902881017714e-04]
rel-76 [4.06561455635949e-03 7.55363841653083e-01]
triaxial_ellipsoid background=0 sld=1 sld_solvent=0 radius_equat_minor=1000 radius_equat_major=2000 radius_polar=200000
cpu double: 137.7 ms gauss-76: 39.1 ms [gauss-76 3.5x faster] ***
qvalue [5.00000000000000e-04 1.00000000000000e-03 2.00000000000000e-03]
target [2.76304345549839e+06 9.28680416872372e+05 1.11806895058799e+05]
actual [2.77380253296110e+06 8.56140969167502e+05 1.35838503298750e+05]
relerr [3.89392263856951e-03 7.81102372645801e-02 2.14938517229313e-01]
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03]
relerr [2.14938517229313e-01]
rel-76 [1.16121876382577e-03]

=== big disks: a=180000 b=200000 c=1000 ===
core_shell_bicelle background=0 radius=90000.0 thick_rim=10000.0 thick_face=10000.0 length=-19000.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0
cpu double: 82.4 ms gauss-76: 2.9 ms [gauss-76 28.2x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [2.66653665528962e-13 3.59113672634819e-11]
rel-76 [1.14111846847644e-02 3.59732048904305e-01]
core_shell_cylinder background=0 radius=90000.0 thickness=10000.0 length=-19000.0 sld_core=0 sld_shell=1 sld_solvent=0
cpu double: 84.0 ms gauss-76: 2.9 ms [gauss-76 29.0x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [2.68578962752998e-13 3.59110035696926e-11]
rel-76 [1.14111846847627e-02 3.59732048904482e-01]
core_shell_ellipsoid background=0 radius_equat_core=90000.0 x_core=-0.10555555555555556 thick_shell=10000.0 x_polar_shell=1 sld_core=0 sld_shell=1 sld_solvent=0
cpu double: 41.6 ms gauss-76: 1.5 ms [gauss-76 27.6x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [7.58386063268560e-10 2.84749621286725e-10 3.45258756102195e-10]
rel-76 [3.94733854946654e-01 4.33595180465061e-01 6.35201112226380e-01]
cylinder background=0 sld=1 sld_solvent=0 radius=100000.0 length=1000
cpu double: 48.0 ms gauss-76: 1.7 ms [gauss-76 27.6x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [2.71334409023091e-13 7.56825181748260e-13 6.07109088371037e-11]
rel-76 [1.57071805281326e-02 2.09757544751596e-01 8.66468363440380e-01]
ellipsoid background=0 sld=1 sld_solvent=0 radius_equatorial=100000.0 radius_polar=500.0
cpu double: 15.6 ms gauss-76: 0.7 ms [gauss-76 21.6x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [1.75173203814238e-09 4.59174526370547e-10 5.99005632332623e-10]
rel-76 [9.78273741627898e-01 8.99057443691645e-01 9.78176397361045e-01]
hollow_cylinder background=0 radius=90000.0 thickness=10000.0 length=1000 sld=1 sld_solvent=0
cpu double: 84.8 ms gauss-76: 2.4 ms [gauss-76 34.9x faster] ***
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [9.41245557634718e-11 3.35609592027369e-10 3.25633611491467e-10 2.67460126643722e-09]
rel-76 [4.47580343066644e-04 5.21355724513471e-02 7.64648514663852e-01 6.70228206782561e-01]
stacked_disks background=0 thick_core=90.0 thick_layer=5.0 radius=100000.0 n_stacking=10 sigma_d=5.0 sld_core=0 sld_layer=1 sld_solvent=0
cpu double: 184.9 ms gauss-76: 6.2 ms [gauss-76 29.8x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [2.71205675968554e-13 7.63063707781760e-13 6.18197793784475e-11]
rel-76 [1.52716728379616e-02 2.10610449137554e-01 8.82791614540920e-01]
barbell background=0 radius=90000.0 radius_bell=100000.0 length=0 sld=1 sld_solvent=0
! ** barbell is slow: 2.9 s for 201 points in [1e-05, 1.0]
cpu double: 2920.3 ms gauss-76: 87.8 ms [gauss-76 33.2x faster] ***
qvalue [2.00000000000000e-03]
target [7.27129821730916e+02]
actual [7.41190530349740e+02]
relerr [1.93372740308644e-02]
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [1.93372740308644e-02 1.49955750185846e-02 1.56140457879098e-01]
rel-76 [3.79228920960551e+01 3.12170071541252e+03 6.31373613309942e+04]
capped_cylinder background=0 radius=90000.0 radius_cap=100000.0 length=0 sld=1 sld_solvent=0
! ** capped_cylinder is slow: 2.8 s for 201 points in [1e-05, 1.0]
cpu double: 2810.2 ms gauss-76: 90.7 ms [gauss-76 31.0x faster] ***
qvalue [2.00000000000000e-03]
target [1.72595023131100e+03]
actual [1.75785915885211e+03]
relerr [1.84877448736591e-02]
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [1.84877448736591e-02 8.36372279095532e-02 6.28846172470313e-02]
rel-76 [1.84877448736932e-02 6.26437363392380e+02 4.78989358648669e+05]
core_shell_bicelle_elliptical_belt_rough background=0 radius=81000.0 x_core=1.123456790123457 thick_rim=9000.0 thick_face=9000.0 length=-17000.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0 sigma=900.0
cpu double: 750.0 ms gauss-76: 159.4 ms [gauss-76 4.7x faster] ***
core_shell_bicelle_elliptical background=0 radius=81000.0 x_core=1.123456790123457 thick_rim=9000.0 thick_face=9000.0 length=-17000.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0
cpu double: 769.0 ms gauss-76: 160.8 ms [gauss-76 4.8x faster] ***
qvalue [1.00000000000000e-02 1.00000000000000e-01]
target [3.28417524346727e+02 8.39942404743102e-03]
actual [1.07636695503978e+02 3.33410038077655e-03]
relerr [6.72256540761388e-01 6.03056071231897e-01]
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [2.44624933387515e-11 6.72256540761388e-01 6.03056071231897e-01]
rel-76 [4.85736303495077e-02 8.84249794867423e-01 7.59229788601636e-01]
core_shell_parallelepiped background=0 sld_core=0 sld_a=1 sld_b=1 sld_c=1 sld_solvent=0 length_a=162000.0 length_b=182000.0 length_c=-17000.0 thick_rim_a=9000.0 thick_rim_b=9000.0 thick_rim_c=9000.0
cpu double: 446.6 ms gauss-76: 89.9 ms [gauss-76 5.0x faster] ***
elliptical_cylinder background=0 radius_minor=90000.0 axis_ratio=200.0 length=1000 sld=1 sld_solvent=0
cpu double: 798.0 ms gauss-76: 127.3 ms [gauss-76 6.3x faster] ***
qvalue [5.00000000000000e-04 1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
target [2.45889265678694e+06 5.75220887632304e+05 1.10114203175099e+05 2.18146592678313e+02 4.28525676717862e-04]
actual [2.38720219223509e+06 5.44187896382958e+05 9.77620682209403e+04 3.98314862747898e+01 1.52817775624962e-04]
relerr [2.91555893479042e-02 5.39496946591810e-02 1.12175673963846e-01 8.17409541970125e-01 6.43387120241161e-01]
! ** gauss-76 is out of tolerance
qvalue [5.00000000000000e-04 1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [2.91555893479042e-02 5.39496946591810e-02 1.12175673963846e-01 8.17409541970125e-01 6.43387120241161e-01]
rel-76 [2.02742538886734e-01 4.53268497520969e-01 8.23315432452319e-01 9.49553348296445e-01 7.81474747412623e-01]
flexible_cylinder_elliptical background=0 length=10000 kuhn_length=1000 radius=180000 axis_ratio=1.1111111111111112 sld=1 sld_solvent=0
cpu double: 56.3 ms gauss-76: 1.7 ms [gauss-76 34.1x faster] ***
hollow_rectangular_prism_thin_walls background=0 sld=1 sld_solvent=0 length_a=180000 b2a_ratio=1.1111111111111112 c2a_ratio=0.005555555555555556
cpu double: 330.6 ms gauss-76: 70.9 ms [gauss-76 4.7x faster] ***
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-01]
relerr [3.72925751437142e-06 1.11591897606304e-05 8.24512257918917e-03]
rel-76 [5.12431584806911e-05 1.01563057441760e-02 7.74169625577148e-01]
hollow_rectangular_prism background=0 sld=1 sld_solvent=0 length_a=180000 b2a_ratio=1.1111111111111112 c2a_ratio=0.005555555555555556 thickness=9000.0
cpu double: 435.5 ms gauss-76: 105.5 ms [gauss-76 4.1x faster] ***
qvalue [2.00000000000000e-03 1.00000000000000e-01]
target [2.45447847519954e+03 9.50759234490639e-04]
actual [2.45671044003829e+03 7.54232382882304e-04]
relerr [9.09343822445768e-04 2.06705172538895e-01]
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [9.09343822445768e-04 2.06705172538895e-01]
rel-76 [1.01976462198123e-01 8.48723865063934e-01]
parallelepiped background=0 sld=1 sld_solvent=0 length_a=180000 length_b=200000 length_c=1000
cpu double: 252.8 ms gauss-76: 57.1 ms [gauss-76 4.4x faster] ***
qvalue [1.00000000000000e-02 1.00000000000000e-01]
target [2.32420758030332e+02 4.41200772547863e-04]
actual [4.22830124104198e+01 9.11548370845696e-05]
relerr [8.18075576515839e-01 7.93393750064931e-01]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [3.36770260687829e-06 7.28446596159316e-08 8.18075576515839e-01 7.93393750064931e-01]
rel-76 [9.14356578529087e-02 6.91990808578044e-01 9.38664534996578e-01 9.00154508762172e-01]
pringle background=0 radius=100000.0 thickness=1000 alpha=0.09999999999999998 beta=0.09999999999999998 sld=1 sld_solvent=0
! ** pringle is slow: 17.2 s for 201 points in [1e-05, 1.0]
cpu double: 17238.0 ms gauss-76: 2204.6 ms [gauss-76 7.8x faster] ***
qvalue [5.00000000000000e-04 1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
target [1.37488167339623e+02 1.98707910687414e+01 2.57324811715918e+00 1.64260746490172e-02 2.02003678328499e-05]
actual [1.29985658342895e+02 1.90252329036105e+01 3.10973392740184e+00 3.37300437763607e-02 3.96386850904321e-04]
relerr [5.45683977166929e-02 4.25528184663491e-02 2.08485845832442e-01 1.05344517768758e+00 1.86227541094433e+01]
! ** gauss-76 is out of tolerance
qvalue [5.00000000000000e-04 1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [5.45683977166929e-02 4.25528184663491e-02 2.08485845832442e-01 1.05344517768758e+00 1.86227541094433e+01]
rel-76 [5.77647800680464e-01 3.75281899950231e-01 1.36186239413146e+00 2.11405987369260e+01 8.05667916660927e+01]
rectangular_prism background=0 sld=1 sld_solvent=0 length_a=180000 b2a_ratio=1.1111111111111112 c2a_ratio=0.005555555555555556
cpu double: 278.3 ms gauss-76: 59.8 ms [gauss-76 4.7x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [7.49616000450599e-08 5.81249935225696e-02]
rel-76 [9.58009001051876e-03 7.69647326091659e-01]
triaxial_ellipsoid background=0 sld=1 sld_solvent=0 radius_equat_minor=180000 radius_equat_major=200000 radius_polar=1000
cpu double: 259.6 ms gauss-76: 46.7 ms [gauss-76 5.6x faster] ***
qvalue [2.00000000000000e-03]
target [9.33814656984402e+04]
actual [9.32551999397901e+04]
relerr [1.35215010500996e-03]
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03]
relerr [1.35215010500996e-03]
rel-76 [7.26075765197704e-04]

=== big cubes: a=200000 b=200000 c=200000 ===
core_shell_bicelle background=0 radius=90000.0 thick_rim=10000.0 thick_face=10000.0 length=180000.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0
cpu double: 115.0 ms gauss-76: 4.1 ms [gauss-76 28.3x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [2.27343634549065e-11 3.65789998841924e-11]
rel-76 [1.44777807487966e-02 8.71561027296069e-01]
core_shell_cylinder background=0 radius=90000.0 thickness=10000.0 length=180000.0 sld_core=0 sld_shell=1 sld_solvent=0
cpu double: 91.4 ms gauss-76: 3.0 ms [gauss-76 30.9x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [2.27420947563974e-11 3.65735042411808e-11]
rel-76 [1.44777807487972e-02 8.71561027296168e-01]
core_shell_ellipsoid background=0 radius_equat_core=90000.0 x_core=1.0 thick_shell=10000.0 x_polar_shell=1 sld_core=0 sld_shell=1 sld_solvent=0
cpu double: 33.9 ms gauss-76: 1.4 ms [gauss-76 25.0x faster] ***
cylinder background=0 sld=1 sld_solvent=0 radius=100000.0 length=200000
cpu double: 54.2 ms gauss-76: 1.9 ms [gauss-76 28.7x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-01]
relerr [1.78808472554716e-11 9.73879242072866e-12]
rel-76 [2.77587645138200e-03 7.83687518915097e-01]
ellipsoid background=0 sld=1 sld_solvent=0 radius_equatorial=100000.0 radius_polar=100000.0
cpu double: 13.2 ms gauss-76: 0.6 ms [gauss-76 22.5x faster] ***
hollow_cylinder background=0 radius=90000.0 thickness=10000.0 length=200000 sld=1 sld_solvent=0
cpu double: 75.8 ms gauss-76: 3.5 ms [gauss-76 21.7x faster] ***
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [8.51837753916576e-12 4.68087865456433e-11 9.10379787605037e-11 1.17935182258835e-10]
rel-76 [4.21961053386734e-03 5.04260898233641e-03 2.05684509620272e-01 6.97868607596239e-01]
stacked_disks background=0 thick_core=18000.0 thick_layer=1000.0 radius=100000.0 n_stacking=10 sigma_d=1000.0 sld_core=0 sld_layer=1 sld_solvent=0
cpu double: 238.5 ms gauss-76: 8.1 ms [gauss-76 29.6x faster] ***
qvalue [5.00000000000000e-04]
target [9.02716549915350e+03]
actual [9.02052754443698e+03]
relerr [7.35331009179351e-04]
! ** gauss-76 is out of tolerance
qvalue [5.00000000000000e-04 2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [7.35331009179351e-04 9.17521234507877e-08 7.55197224910177e-13 4.69463238360040e-11]
rel-76 [7.35331009179351e-04 1.63987317375929e-02 2.10337159071940e-01 6.77470703873726e-01]
barbell background=0 radius=100000.0 radius_bell=100000.0 length=0.0 sld=1 sld_solvent=0
! ** barbell is slow: 3.0 s for 201 points in [1e-05, 1.0]
cpu double: 3046.6 ms gauss-76: 90.8 ms [gauss-76 33.6x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [8.03153141178969e-13 5.09160491688056e-13 4.62206492110789e-12]
rel-76 [1.83868751287866e-03 4.65977239720104e+03 2.82979956402944e+04]
capped_cylinder background=0 radius=100000.0 radius_cap=100000.0 length=0.0 sld=1 sld_solvent=0
! ** capped_cylinder is slow: 2.8 s for 201 points in [1e-05, 1.0]
cpu double: 2768.1 ms gauss-76: 86.2 ms [gauss-76 32.1x faster] ***
! ** gauss-76 is out of tolerance
qvalue [2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [8.03153141178969e-13 5.09160491688056e-13 4.62206492110789e-12]
rel-76 [1.83868751287866e-03 4.65977239720104e+03 2.82979956402944e+04]
core_shell_bicelle_elliptical_belt_rough background=0 radius=90000.0 x_core=1.0 thick_rim=10000.0 thick_face=10000.0 length=180000.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0 sigma=1000.0
cpu double: 711.9 ms gauss-76: 152.7 ms [gauss-76 4.7x faster] ***
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-02]
relerr [6.98454534209187e-11 6.64906671061497e-11 5.34405019005807e-02]
rel-76 [1.16443555404283e-04 8.36590607744081e-03 2.31454001384129e-01]
core_shell_bicelle_elliptical background=0 radius=90000.0 x_core=1.0 thick_rim=10000.0 thick_face=10000.0 length=180000.0 sld_core=0 sld_face=1 sld_rim=1 sld_solvent=0
cpu double: 697.8 ms gauss-76: 160.7 ms [gauss-76 4.3x faster] ***
qvalue [1.00000000000000e-01]
target [8.11642917959339e-05]
actual [5.80430769644655e-05]
relerr [2.84869298060292e-01]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [1.06758449816373e-10 7.39114308251186e-11 4.81114068077031e-02 2.84869298060292e-01]
rel-76 [9.31522529275513e-04 1.55790796139149e-02 2.23895368617276e-01 3.62233351378684e-01]
core_shell_parallelepiped background=0 sld_core=0 sld_a=1 sld_b=1 sld_c=1 sld_solvent=0 length_a=180000.0 length_b=180000.0 length_c=180000.0 thick_rim_a=10000.0 thick_rim_b=10000.0 thick_rim_c=10000.0
cpu double: 430.8 ms gauss-76: 91.1 ms [gauss-76 4.7x faster] ***
qvalue [1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-01]
target [6.98345673054335e+04 3.00783773475261e+02 2.36105089784728e-04]
actual [7.00984945279583e+04 2.61924219691136e+02 4.25855139776527e-05]
relerr [3.77932065320106e-03 1.29194315687782e-01 8.19633223424024e-01]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-01]
relerr [3.77932065320106e-03 1.29194315687782e-01 8.19633223424024e-01]
rel-76 [7.17822024107535e-03 1.66564557817961e-01 9.23838528505826e-01]
elliptical_cylinder background=0 radius_minor=100000.0 axis_ratio=1.0 length=200000 sld=1 sld_solvent=0
cpu double: 421.9 ms gauss-76: 99.6 ms [gauss-76 4.2x faster] ***
qvalue [1.00000000000000e-02]
target [9.35785168289641e-01]
actual [2.47791432535872e-01]
relerr [7.35204787452694e-01]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [7.06622233597165e-11 5.68532462759433e-10 7.35204787452694e-01 3.55485690903060e-02]
rel-76 [3.80598188627008e-03 2.54076483540652e-01 8.77926103428513e-01 7.21416295661587e-01]
flexible_cylinder_elliptical background=0 length=2000000 kuhn_length=200000 radius=200000 axis_ratio=1.0 sld=1 sld_solvent=0
cpu double: 54.6 ms gauss-76: 1.6 ms [gauss-76 35.1x faster] ***
hollow_rectangular_prism_thin_walls background=0 sld=1 sld_solvent=0 length_a=200000 b2a_ratio=1.0 c2a_ratio=1.0
cpu double: 331.8 ms gauss-76: 72.4 ms [gauss-76 4.6x faster] ***
qvalue [1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-01]
target [8.56430787049370e+02 7.05563907715598e+01 1.13816015341818e-01]
actual [8.56144929699486e+02 7.01479891411647e+01 5.53202709676415e-02]
relerr [3.33777526691141e-04 5.78830104444219e-03 5.13950028899704e-01]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-01]
relerr [3.33777526691141e-04 5.78830104444219e-03 5.13950028899704e-01]
rel-76 [9.36321746794216e-04 1.66649391160123e-02 8.48227217539955e-01]
hollow_rectangular_prism background=0 sld=1 sld_solvent=0 length_a=200000 b2a_ratio=1.0 c2a_ratio=1.0 thickness=10000.0
cpu double: 435.5 ms gauss-76: 94.9 ms [gauss-76 4.6x faster] ***
qvalue [1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-01]
target [2.42978141719551e+05 1.05049560360160e+03 9.87931292101115e-04]
actual [2.42955329124302e+05 1.02144394940322e+03 4.06977282690810e-04]
relerr [9.38874381325301e-05 2.76551887497432e-02 5.88051025466298e-01]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-01]
relerr [9.38874381325301e-05 2.76551887497432e-02 5.88051025466298e-01]
rel-76 [2.99875418191674e-04 8.16344362426419e-02 9.37201920567185e-01]
parallelepiped background=0 sld=1 sld_solvent=0 length_a=200000 length_b=200000 length_c=200000
cpu double: 265.9 ms gauss-76: 55.2 ms [gauss-76 4.8x faster] ***
qvalue [1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
target [1.16766869284293e+04 1.81760671001750e+03 2.51648519637758e+00 2.57757351964195e-05]
actual [1.16947856496101e+04 1.81283267798809e+03 1.78496186628208e+00 6.23884113123936e-06]
relerr [1.54998770557014e-03 2.62654841836823e-03 2.90692482971293e-01 7.57956811563381e-01]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [1.54998770557014e-03 2.62654841836823e-03 2.90692482971293e-01 7.57956811563381e-01]
rel-76 [7.44548186547520e-03 2.13338736919987e-01 3.68289786302960e-01 9.21628222984412e-01]
pringle background=0 radius=100000.0 thickness=200000 alpha=0.001 beta=0.001 sld=1 sld_solvent=0
! ** pringle is slow: 10.7 s for 201 points in [1e-05, 1.0]
cpu double: 10683.5 ms gauss-76: 1950.8 ms [gauss-76 5.5x faster] ***
qvalue [5.00000000000000e-04 1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
target [5.31579388633669e+05 5.27862463513970e+04 4.82947216221200e+03 1.68016243213517e+00 6.11226958801765e-05]
actual [4.61165445916087e+05 5.05875699920649e+04 4.71214891387973e+03 5.56221836643101e+00 4.11959261446954e-02]
relerr [1.32461762482117e-01 4.16524475844634e-02 2.42931824413975e-02 2.31052418507090e+00 6.72987388014674e+02]
! ** gauss-76 is out of tolerance
qvalue [5.00000000000000e-04 1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-02 1.00000000000000e-01]
relerr [1.32461762482117e-01 4.16524475844634e-02 2.42931824413975e-02 2.31052418507090e+00 6.72987388014674e+02]
rel-76 [1.16531803286495e-01 1.12688432909316e-01 4.56569948595530e-01 1.33377913520446e+03 1.65262168473276e+04]
rectangular_prism background=0 sld=1 sld_solvent=0 length_a=200000 b2a_ratio=1.0 c2a_ratio=1.0
cpu double: 251.9 ms gauss-76: 56.0 ms [gauss-76 4.5x faster] ***
qvalue [1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-01]
target [1.16766869304387e+04 1.81760671130031e+03 3.55638693177586e-05]
actual [1.16947856493809e+04 1.81283267787616e+03 1.77462986346068e-05]
relerr [1.54998751358641e-03 2.62654918386392e-03 5.01002028883701e-01]
! ** gauss-76 is out of tolerance
qvalue [1.00000000000000e-03 2.00000000000000e-03 1.00000000000000e-01]
relerr [1.54998751358641e-03 2.62654918386392e-03 5.01002028883701e-01]
rel-76 [4.41276355454061e-03 7.50778711594389e-03 8.32131834210397e-01]
triaxial_ellipsoid background=0 sld=1 sld_solvent=0 radius_equat_minor=200000 radius_equat_major=200000 radius_polar=200000
cpu double: 164.0 ms gauss-76: 35.2 ms [gauss-76 4.7x faster] ***