Changeset 66dbbfb in sasmodels


Ignore:
Timestamp:
Mar 18, 2019 5:08:46 PM (4 months ago)
Author:
Paul Kienzle <pkienzle@…>
Branches:
master, core_shell_microgels, magnetic_model, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests
Children:
edc783e
Parents:
242b361
Message:

prettier output for explore/asymint.py

File:
1 edited

Legend:

Unmodified
Added
Removed
  • explore/asymint.py

    r242b361 r66dbbfb  
    461461    fn = lambda azimuthal, polar: kernel_2d(q=q, theta=polar, phi=azimuthal) 
    462462    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" 
    464464          % (rule, rule_obj.degree, len(rule_obj.points), Iq)) 
    465465 
     
    525525    quadpy_method(Q, "Stroud:U3 14-1") 
    526526 
    527     print("gauss-20", *gauss_quad_2d(Q, n=20)) 
    528     print("gauss-76", *gauss_quad_2d(Q, n=76)) 
    529     print("gauss-150", *gauss_quad_2d(Q, n=150)) 
    530     print("gauss-500", *gauss_quad_2d(Q, n=500)) 
    531     print("gauss-1025", *gauss_quad_2d(Q, n=1025)) 
    532     print("gauss-2049", *gauss_quad_2d(Q, n=2049)) 
    533     print("gauss-20 usub", *gauss_quad_usub(Q, n=20)) 
    534     print("gauss-76 usub", *gauss_quad_usub(Q, n=76)) 
    535     print("gauss-150 usub", *gauss_quad_usub(Q, n=150)) 
     527    print("gauss-20 points=%d => %.15g" % gauss_quad_2d(Q, n=20)) 
     528    print("gauss-76 points=%d => %.15g" % gauss_quad_2d(Q, n=76)) 
     529    print("gauss-150 points=%d => %.15g" % gauss_quad_2d(Q, n=150)) 
     530    print("gauss-500 points=%d => %.15g" % gauss_quad_2d(Q, n=500)) 
     531    print("gauss-1025 points=%d => %.15g" % gauss_quad_2d(Q, n=1025)) 
     532    print("gauss-2049 points=%d => %.15g" % gauss_quad_2d(Q, n=2049)) 
     533    print("gauss-20 usub points=%d => %.15g" % gauss_quad_usub(Q, n=20)) 
     534    print("gauss-76 usub points=%d => %.15g" % gauss_quad_usub(Q, n=76)) 
     535    print("gauss-150 usub points=%d => %.15g" % gauss_quad_usub(Q, n=150)) 
    536536 
    537537    #gridded_2d(Q, n=2**8+1) 
Note: See TracChangeset for help on using the changeset viewer.