Changeset 66dbbfb in sasmodels for explore/asymint.py
 Timestamp:
 Mar 18, 2019 5:08:46 PM (5 months ago)
 Branches:
 master, core_shell_microgels, magnetic_model, ticket1257vesicleproduct, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests
 Children:
 edc783e
 Parents:
 242b361
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

explore/asymint.py
r242b361 r66dbbfb 461 461 fn = lambda azimuthal, polar: kernel_2d(q=q, theta=polar, phi=azimuthal) 462 462 Iq = integrate_spherical(fn, rule=rule_obj)/(4*pi) 463 print("%s degree=%d points=%s => % g"463 print("%s degree=%d points=%s => %.15g" 464 464 % (rule, rule_obj.degree, len(rule_obj.points), Iq)) 465 465 … … 525 525 quadpy_method(Q, "Stroud:U3 141") 526 526 527 print("gauss20 ", *gauss_quad_2d(Q, n=20))528 print("gauss76 ", *gauss_quad_2d(Q, n=76))529 print("gauss150 ", *gauss_quad_2d(Q, n=150))530 print("gauss500 ", *gauss_quad_2d(Q, n=500))531 print("gauss1025 ", *gauss_quad_2d(Q, n=1025))532 print("gauss2049 ", *gauss_quad_2d(Q, n=2049))533 print("gauss20 usub ", *gauss_quad_usub(Q, n=20))534 print("gauss76 usub ", *gauss_quad_usub(Q, n=76))535 print("gauss150 usub ", *gauss_quad_usub(Q, n=150))527 print("gauss20 points=%d => %.15g" % gauss_quad_2d(Q, n=20)) 528 print("gauss76 points=%d => %.15g" % gauss_quad_2d(Q, n=76)) 529 print("gauss150 points=%d => %.15g" % gauss_quad_2d(Q, n=150)) 530 print("gauss500 points=%d => %.15g" % gauss_quad_2d(Q, n=500)) 531 print("gauss1025 points=%d => %.15g" % gauss_quad_2d(Q, n=1025)) 532 print("gauss2049 points=%d => %.15g" % gauss_quad_2d(Q, n=2049)) 533 print("gauss20 usub points=%d => %.15g" % gauss_quad_usub(Q, n=20)) 534 print("gauss76 usub points=%d => %.15g" % gauss_quad_usub(Q, n=76)) 535 print("gauss150 usub points=%d => %.15g" % gauss_quad_usub(Q, n=150)) 536 536 537 537 #gridded_2d(Q, n=2**8+1)
Note: See TracChangeset
for help on using the changeset viewer.