Changeset 706f466 in sasmodels


Ignore:
Timestamp:
Sep 28, 2017 2:40:23 PM (7 years 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:
6ab64c9
Parents:
62d7601
Message:

fix doc strings for python-based sas special functions

Files:
2 edited

Legend:

Unmodified
Added
Removed
  • doc/genapi.py

    r2e66ef5 r706f466  
    5959    #('alignment', 'GPU data alignment [unused]'), 
    6060    ('bumps_model', 'Bumps interface'), 
     61    ('compare', 'Compare models on different compute engines'), 
    6162    ('compare_many', 'Batch compare models on different compute engines'), 
    62     ('compare', 'Compare models on different compute engines'), 
     63    ('conversion_table', 'Model conversion table'), 
    6364    ('convert', 'Sasview to sasmodel converter'), 
    6465    ('core', 'Model access'), 
     
    8283    ('sasview_model', 'Sasview interface'), 
    8384    ('sesans', 'SESANS calculation routines'), 
     85    ('special', 'Special functions library'), 
    8486    ('weights', 'Distribution functions'), 
    8587] 
  • sasmodels/special.py

    r4f611f1 r706f466  
    1 """ 
     1r""" 
    22Special Functions 
    33................. 
     
    1010    M_PI, M_PI_2, M_PI_4, M_SQRT1_2, M_E: 
    1111        $\pi$, $\pi/2$, $\pi/4$, $1/\sqrt{2}$ and Euler's constant $e$ 
     12 
    1213    exp, log, pow(x,y), expm1, sqrt: 
    1314        Power functions $e^x$, $\ln x$, $x^y$, $e^x - 1$, $\sqrt{x}$. 
    1415        The function expm1(x) is accurate across all $x$, including $x$ 
    1516        very close to zero. 
     17 
    1618    sin, cos, tan, asin, acos, atan: 
    1719        Trigonometry functions and inverses, operating on radians. 
     20 
    1821    sinh, cosh, tanh, asinh, acosh, atanh: 
    1922        Hyperbolic trigonometry functions. 
     23 
    2024    atan2(y,x): 
    2125        Angle from the $x$\ -axis to the point $(x,y)$, which is equal to 
     
    2428        atan(y/x) would return a value in quadrant I. Similarly for 
    2529        quadrants II and IV when $x$ and $y$ have opposite sign. 
     30 
    2631    fmin(x,y), fmax(x,y), trunc, rint: 
    2732        Floating point functions.  rint(x) returns the nearest integer. 
     33 
    2834    NAN: 
    2935        NaN, Not a Number, $0/0$.  Use isnan(x) to test for NaN.  Note that 
    3036        you cannot use :code:`x == NAN` to test for NaN values since that 
    3137        will always return false.  NAN does not equal NAN! 
     38 
    3239    INFINITY: 
    3340        $\infty, 1/0$.  Use isinf(x) to test for infinity, or isfinite(x) 
    3441        to test for finite and not NaN. 
     42 
    3543    erf, erfc, tgamma, lgamma:  **do not use** 
    3644        Special functions that should be part of the standard, but are missing 
     
    4250    M_PI_180, M_4PI_3: 
    4351        $\frac{\pi}{180}$, $\frac{4\pi}{3}$ 
     52 
    4453    SINCOS(x, s, c): 
    4554        Macro which sets s=sin(x) and c=cos(x). The variables *c* and *s* 
    4655        must be declared first. 
     56 
    4757    square(x): 
    4858        $x^2$ 
     59 
    4960    cube(x): 
    5061        $x^3$ 
     62 
    5163    sas_sinx_x(x): 
    5264        $\sin(x)/x$, with limit $\sin(0)/0 = 1$. 
     65 
    5366    powr(x, y): 
    5467        $x^y$ for $x \ge 0$; this is faster than general $x^y$ on some GPUs. 
     68 
    5569    pown(x, n): 
    5670        $x^n$ for $n$ integer; this is faster than general $x^n$ on some GPUs. 
     71 
    5772    FLOAT_SIZE: 
    5873        The number of bytes in a floating point value.  Even though all 
     
    6782            ... code for single precision ... 
    6883            #endif 
     84 
    6985    SAS_DOUBLE: 
    7086        A replacement for :code:`double` so that the declared variable will 
     
    88104        sorted from highest to lowest. 
    89105 
    90         :code:`source = ["lib/polevl.c", ...]` (`link to code <https://github.com/SasView/sasmodels/tree/master/sasmodels/models/lib/polevl.c>`_) 
    91  
    92106    p1evl(x, c, n): 
    93107        Evaluation of normalized polynomial $p(x) = x^n + \sum_{i=0}^{n-1} c_i x^i$ 
     
    97111        sorted from highest to lowest. 
    98112 
    99         :code:`source = ["lib/polevl.c", ...]` 
    100         (`link to code <https://github.com/SasView/sasmodels/tree/master/sasmodels/models/lib/polevl.c>`_) 
    101  
    102113    sas_gamma(x): 
    103114        Gamma function $\text{sas_gamma}(x) = \Gamma(x)$. 
     
    105116        The standard math function, tgamma(x) is unstable for $x < 1$ 
    106117        on some platforms. 
    107  
    108         :code:`source = ["lib/sasgamma.c", ...]` 
    109         (`link to code <https://github.com/SasView/sasmodels/tree/master/sasmodels/models/lib/sas_gamma.c>`_) 
    110118 
    111119    sas_erf(x), sas_erfc(x): 
     
    118126        on some platforms. 
    119127 
    120         :code:`source = ["lib/polevl.c", "lib/sas_erf.c", ...]` 
    121         (`link to error functions' code <https://github.com/SasView/sasmodels/tree/master/sasmodels/models/lib/sas_erf.c>`_) 
    122  
    123128    sas_J0(x): 
    124129        Bessel function of the first kind $\text{sas_J0}(x)=J_0(x)$ where 
     
    127132        The standard math function j0(x) is not available on all platforms. 
    128133 
    129         :code:`source = ["lib/polevl.c", "lib/sas_J0.c", ...]` 
    130         (`link to Bessel function's code <https://github.com/SasView/sasmodels/tree/master/sasmodels/models/lib/sas_J0.c>`_) 
    131  
    132134    sas_J1(x): 
    133135        Bessel function of the first kind  $\text{sas_J1}(x)=J_1(x)$ where 
     
    135137 
    136138        The standard math function j1(x) is not available on all platforms. 
    137  
    138         :code:`source = ["lib/polevl.c", "lib/sas_J1.c", ...]` 
    139         (`link to Bessel function's code <https://github.com/SasView/sasmodels/tree/master/sasmodels/models/lib/sas_J1.c>`_) 
    140139 
    141140    sas_JN(n, x): 
     
    147146        The standard math function jn(n, x) is not available on all platforms. 
    148147 
    149         :code:`source = ["lib/polevl.c", "lib/sas_J0.c", "lib/sas_J1.c", "lib/sas_JN.c", ...]` 
    150         (`link to Bessel function's code <https://github.com/SasView/sasmodels/tree/master/sasmodels/models/lib/sas_JN.c>`_) 
    151  
    152148    sas_Si(x): 
    153149        Sine integral $\text{Si}(x) = \int_0^x \tfrac{\sin t}{t}\,dt$. 
     
    170166           - \frac{x^3}{3\times 3!} + \frac{x^5}{5 \times 5!} - \frac{x^7}{7 \times 7!} 
    171167           + \frac{x^9}{9\times 9!} - \frac{x^{11}}{11\times 11!} 
    172  
    173         :code:`source = ["lib/Si.c", ...]` 
    174         (`link to code <https://github.com/SasView/sasmodels/tree/master/sasmodels/models/lib/Si.c>`_) 
    175168 
    176169    sas_3j1x_x(x): 
     
    182175        This function uses a Taylor series for small $x$ for numerical accuracy. 
    183176 
    184         :code:`source = ["lib/sas_3j1x_x.c", ...]` 
    185         (`link to code <https://github.com/SasView/sasmodels/tree/master/sasmodels/models/lib/sas_3j1x_x.c>`_) 
    186  
    187177 
    188178    sas_2J1x_x(x): 
     
    191181        and first order. 
    192182 
    193         :code:`source = ["lib/polevl.c", "lib/sas_J1.c", ...]` 
    194         (`link to Bessel form's code <https://github.com/SasView/sasmodels/tree/master/sasmodels/models/lib/sas_J1.c>`_) 
    195  
    196183 
    197184    Gauss76Z[i], Gauss76Wt[i]: 
     
    201188        Similar arrays are available in :code:`gauss20.c` for 20-point 
    202189        quadrature and in :code:`gauss150.c` for 150-point quadrature. 
    203  
    204         :code:`source = ["lib/gauss76.c", ...]` 
    205         (`link to code <https://github.com/SasView/sasmodels/tree/master/sasmodels/models/lib/gauss76.c>`_) 
    206  
    207190 
    208191""" 
     
    290273            retvalue = (sin(x) - x*cos(x))/x**2 
    291274        retvalue[x == 0.] = 0. 
    292          
     275 
    293276def sas_3j1x_x(x): 
    294277    if np.isscalar(x): 
     
    313296 
    314297Gauss20Wt = np.array([ 
    315         .0176140071391521, 
    316         .0406014298003869, 
    317         .0626720483341091, 
    318         .0832767415767047, 
    319         .10193011981724, 
    320         .118194531961518, 
    321         .131688638449177, 
    322         .142096109318382, 
    323         .149172986472604, 
    324         .152753387130726, 
    325         .152753387130726, 
    326         .149172986472604, 
    327         .142096109318382, 
    328         .131688638449177, 
    329         .118194531961518, 
    330         .10193011981724, 
    331         .0832767415767047, 
    332         .0626720483341091, 
    333         .0406014298003869, 
    334         .0176140071391521 
     298    .0176140071391521, 
     299    .0406014298003869, 
     300    .0626720483341091, 
     301    .0832767415767047, 
     302    .10193011981724, 
     303    .118194531961518, 
     304    .131688638449177, 
     305    .142096109318382, 
     306    .149172986472604, 
     307    .152753387130726, 
     308    .152753387130726, 
     309    .149172986472604, 
     310    .142096109318382, 
     311    .131688638449177, 
     312    .118194531961518, 
     313    .10193011981724, 
     314    .0832767415767047, 
     315    .0626720483341091, 
     316    .0406014298003869, 
     317    .0176140071391521 
    335318]) 
    336319 
    337320Gauss20Z = np.array([ 
    338         -.993128599185095, 
    339         -.963971927277914, 
    340         -.912234428251326, 
    341         -.839116971822219, 
    342         -.746331906460151, 
    343         -.636053680726515, 
    344         -.510867001950827, 
    345         -.37370608871542, 
    346         -.227785851141645, 
    347         -.076526521133497, 
    348         .0765265211334973, 
    349         .227785851141645, 
    350         .37370608871542, 
    351         .510867001950827, 
    352         .636053680726515, 
    353         .746331906460151, 
    354         .839116971822219, 
    355         .912234428251326, 
    356         .963971927277914, 
    357         .993128599185095 
     321    -.993128599185095, 
     322    -.963971927277914, 
     323    -.912234428251326, 
     324    -.839116971822219, 
     325    -.746331906460151, 
     326    -.636053680726515, 
     327    -.510867001950827, 
     328    -.37370608871542, 
     329    -.227785851141645, 
     330    -.076526521133497, 
     331    .0765265211334973, 
     332    .227785851141645, 
     333    .37370608871542, 
     334    .510867001950827, 
     335    .636053680726515, 
     336    .746331906460151, 
     337    .839116971822219, 
     338    .912234428251326, 
     339    .963971927277914, 
     340    .993128599185095 
    358341]) 
    359342 
    360343Gauss76Wt = np.array([ 
    361         .00126779163408536,             #0 
    362         .00294910295364247, 
    363         .00462793522803742, 
    364         .00629918049732845, 
    365         .00795984747723973, 
    366         .00960710541471375, 
    367         .0112381685696677, 
    368         .0128502838475101, 
    369         .0144407317482767, 
    370         .0160068299122486, 
    371         .0175459372914742,              #10 
    372         .0190554584671906, 
    373         .020532847967908, 
    374         .0219756145344162, 
    375         .0233813253070112, 
    376         .0247476099206597, 
    377         .026072164497986, 
    378         .0273527555318275, 
    379         .028587223650054, 
    380         .029773487255905, 
    381         .0309095460374916,              #20 
    382         .0319934843404216, 
    383         .0330234743977917, 
    384         .0339977794120564, 
    385         .0349147564835508, 
    386         .0357728593807139, 
    387         .0365706411473296, 
    388         .0373067565423816, 
    389         .0379799643084053, 
    390         .0385891292645067, 
    391         .0391332242205184,              #30 
    392         .0396113317090621, 
    393         .0400226455325968, 
    394         .040366472122844, 
    395         .0406422317102947, 
    396         .0408494593018285, 
    397         .040987805464794, 
    398         .0410570369162294, 
    399         .0410570369162294, 
    400         .040987805464794, 
    401         .0408494593018285,              #40 
    402         .0406422317102947, 
    403         .040366472122844, 
    404         .0400226455325968, 
    405         .0396113317090621, 
    406         .0391332242205184, 
    407         .0385891292645067, 
    408         .0379799643084053, 
    409         .0373067565423816, 
    410         .0365706411473296, 
    411         .0357728593807139,              #50 
    412         .0349147564835508, 
    413         .0339977794120564, 
    414         .0330234743977917, 
    415         .0319934843404216, 
    416         .0309095460374916, 
    417         .029773487255905, 
    418         .028587223650054, 
    419         .0273527555318275, 
    420         .026072164497986, 
    421         .0247476099206597,              #60 
    422         .0233813253070112, 
    423         .0219756145344162, 
    424         .020532847967908, 
    425         .0190554584671906, 
    426         .0175459372914742, 
    427         .0160068299122486, 
    428         .0144407317482767, 
    429         .0128502838475101, 
    430         .0112381685696677, 
    431         .00960710541471375,             #70 
    432         .00795984747723973, 
    433         .00629918049732845, 
    434         .00462793522803742, 
    435         .00294910295364247, 
    436         .00126779163408536              #75 (indexed from 0) 
     344    .00126779163408536,         #0 
     345    .00294910295364247, 
     346    .00462793522803742, 
     347    .00629918049732845, 
     348    .00795984747723973, 
     349    .00960710541471375, 
     350    .0112381685696677, 
     351    .0128502838475101, 
     352    .0144407317482767, 
     353    .0160068299122486, 
     354    .0175459372914742,          #10 
     355    .0190554584671906, 
     356    .020532847967908, 
     357    .0219756145344162, 
     358    .0233813253070112, 
     359    .0247476099206597, 
     360    .026072164497986, 
     361    .0273527555318275, 
     362    .028587223650054, 
     363    .029773487255905, 
     364    .0309095460374916,          #20 
     365    .0319934843404216, 
     366    .0330234743977917, 
     367    .0339977794120564, 
     368    .0349147564835508, 
     369    .0357728593807139, 
     370    .0365706411473296, 
     371    .0373067565423816, 
     372    .0379799643084053, 
     373    .0385891292645067, 
     374    .0391332242205184,          #30 
     375    .0396113317090621, 
     376    .0400226455325968, 
     377    .040366472122844, 
     378    .0406422317102947, 
     379    .0408494593018285, 
     380    .040987805464794, 
     381    .0410570369162294, 
     382    .0410570369162294, 
     383    .040987805464794, 
     384    .0408494593018285,          #40 
     385    .0406422317102947, 
     386    .040366472122844, 
     387    .0400226455325968, 
     388    .0396113317090621, 
     389    .0391332242205184, 
     390    .0385891292645067, 
     391    .0379799643084053, 
     392    .0373067565423816, 
     393    .0365706411473296, 
     394    .0357728593807139,          #50 
     395    .0349147564835508, 
     396    .0339977794120564, 
     397    .0330234743977917, 
     398    .0319934843404216, 
     399    .0309095460374916, 
     400    .029773487255905, 
     401    .028587223650054, 
     402    .0273527555318275, 
     403    .026072164497986, 
     404    .0247476099206597,          #60 
     405    .0233813253070112, 
     406    .0219756145344162, 
     407    .020532847967908, 
     408    .0190554584671906, 
     409    .0175459372914742, 
     410    .0160068299122486, 
     411    .0144407317482767, 
     412    .0128502838475101, 
     413    .0112381685696677, 
     414    .00960710541471375,         #70 
     415    .00795984747723973, 
     416    .00629918049732845, 
     417    .00462793522803742, 
     418    .00294910295364247, 
     419    .00126779163408536          #75 (indexed from 0) 
    437420]) 
    438421 
    439422Gauss76Z = np.array([ 
    440         -.999505948362153,              #0 
    441         -.997397786355355, 
    442         -.993608772723527, 
    443         -.988144453359837, 
    444         -.981013938975656, 
    445         -.972229228520377, 
    446         -.961805126758768, 
    447         -.949759207710896, 
    448         -.936111781934811, 
    449         -.92088586125215, 
    450         -.904107119545567,              #10 
    451         -.885803849292083, 
    452         -.866006913771982, 
    453         -.844749694983342, 
    454         -.822068037328975, 
    455         -.7980001871612, 
    456         -.77258672828181, 
    457         -.74587051350361, 
    458         -.717896592387704, 
    459         -.688712135277641, 
    460         -.658366353758143,              #20 
    461         -.626910417672267, 
    462         -.594397368836793, 
    463         -.560882031601237, 
    464         -.526420920401243, 
    465         -.491072144462194, 
    466         -.454895309813726, 
    467         -.417951418780327, 
    468         -.380302767117504, 
    469         -.342012838966962, 
    470         -.303146199807908,              #30 
    471         -.263768387584994, 
    472         -.223945802196474, 
    473         -.183745593528914, 
    474         -.143235548227268, 
    475         -.102483975391227, 
    476         -.0615595913906112, 
    477         -.0205314039939986, 
    478         .0205314039939986, 
    479         .0615595913906112, 
    480         .102483975391227,                       #40 
    481         .143235548227268, 
    482         .183745593528914, 
    483         .223945802196474, 
    484         .263768387584994, 
    485         .303146199807908, 
    486         .342012838966962, 
    487         .380302767117504, 
    488         .417951418780327, 
    489         .454895309813726, 
    490         .491072144462194,               #50 
    491         .526420920401243, 
    492         .560882031601237, 
    493         .594397368836793, 
    494         .626910417672267, 
    495         .658366353758143, 
    496         .688712135277641, 
    497         .717896592387704, 
    498         .74587051350361, 
    499         .77258672828181, 
    500         .7980001871612, #60 
    501         .822068037328975, 
    502         .844749694983342, 
    503         .866006913771982, 
    504         .885803849292083, 
    505         .904107119545567, 
    506         .92088586125215, 
    507         .936111781934811, 
    508         .949759207710896, 
    509         .961805126758768, 
    510         .972229228520377,               #70 
    511         .981013938975656, 
    512         .988144453359837, 
    513         .993608772723527, 
    514         .997397786355355, 
    515         .999505948362153                #75 
     423    -.999505948362153,          #0 
     424    -.997397786355355, 
     425    -.993608772723527, 
     426    -.988144453359837, 
     427    -.981013938975656, 
     428    -.972229228520377, 
     429    -.961805126758768, 
     430    -.949759207710896, 
     431    -.936111781934811, 
     432    -.92088586125215, 
     433    -.904107119545567,          #10 
     434    -.885803849292083, 
     435    -.866006913771982, 
     436    -.844749694983342, 
     437    -.822068037328975, 
     438    -.7980001871612, 
     439    -.77258672828181, 
     440    -.74587051350361, 
     441    -.717896592387704, 
     442    -.688712135277641, 
     443    -.658366353758143,          #20 
     444    -.626910417672267, 
     445    -.594397368836793, 
     446    -.560882031601237, 
     447    -.526420920401243, 
     448    -.491072144462194, 
     449    -.454895309813726, 
     450    -.417951418780327, 
     451    -.380302767117504, 
     452    -.342012838966962, 
     453    -.303146199807908,          #30 
     454    -.263768387584994, 
     455    -.223945802196474, 
     456    -.183745593528914, 
     457    -.143235548227268, 
     458    -.102483975391227, 
     459    -.0615595913906112, 
     460    -.0205314039939986, 
     461    .0205314039939986, 
     462    .0615595913906112, 
     463    .102483975391227,                   #40 
     464    .143235548227268, 
     465    .183745593528914, 
     466    .223945802196474, 
     467    .263768387584994, 
     468    .303146199807908, 
     469    .342012838966962, 
     470    .380302767117504, 
     471    .417951418780327, 
     472    .454895309813726, 
     473    .491072144462194,           #50 
     474    .526420920401243, 
     475    .560882031601237, 
     476    .594397368836793, 
     477    .626910417672267, 
     478    .658366353758143, 
     479    .688712135277641, 
     480    .717896592387704, 
     481    .74587051350361, 
     482    .77258672828181, 
     483    .7980001871612,     #60 
     484    .822068037328975, 
     485    .844749694983342, 
     486    .866006913771982, 
     487    .885803849292083, 
     488    .904107119545567, 
     489    .92088586125215, 
     490    .936111781934811, 
     491    .949759207710896, 
     492    .961805126758768, 
     493    .972229228520377,           #70 
     494    .981013938975656, 
     495    .988144453359837, 
     496    .993608772723527, 
     497    .997397786355355, 
     498    .999505948362153            #75 
    516499]) 
    517500 
    518501Gauss150Z = np.array([ 
    519         -0.9998723404457334, 
    520         -0.9993274305065947, 
    521         -0.9983473449340834, 
    522         -0.9969322929775997, 
    523         -0.9950828645255290, 
    524         -0.9927998590434373, 
    525         -0.9900842691660192, 
    526         -0.9869372772712794, 
    527         -0.9833602541697529, 
    528         -0.9793547582425894, 
    529         -0.9749225346595943, 
    530         -0.9700655145738374, 
    531         -0.9647858142586956, 
    532         -0.9590857341746905, 
    533         -0.9529677579610971, 
    534         -0.9464345513503147, 
    535         -0.9394889610042837, 
    536         -0.9321340132728527, 
    537         -0.9243729128743136, 
    538         -0.9162090414984952, 
    539         -0.9076459563329236, 
    540         -0.8986873885126239, 
    541         -0.8893372414942055, 
    542         -0.8795995893549102, 
    543         -0.8694786750173527, 
    544         -0.8589789084007133, 
    545         -0.8481048644991847, 
    546         -0.8368612813885015, 
    547         -0.8252530581614230, 
    548         -0.8132852527930605, 
    549         -0.8009630799369827, 
    550         -0.7882919086530552, 
    551         -0.7752772600680049, 
    552         -0.7619248049697269, 
    553         -0.7482403613363824, 
    554         -0.7342298918013638, 
    555         -0.7198995010552305, 
    556         -0.7052554331857488, 
    557         -0.6903040689571928, 
    558         -0.6750519230300931, 
    559         -0.6595056411226444, 
    560         -0.6436719971150083, 
    561         -0.6275578900977726, 
    562         -0.6111703413658551, 
    563         -0.5945164913591590, 
    564         -0.5776035965513142, 
    565         -0.5604390262878617, 
    566         -0.5430302595752546, 
    567         -0.5253848818220803, 
    568         -0.5075105815339176, 
    569         -0.4894151469632753, 
    570         -0.4711064627160663, 
    571         -0.4525925063160997, 
    572         -0.4338813447290861, 
    573         -0.4149811308476706, 
    574         -0.3959000999390257, 
    575         -0.3766465660565522, 
    576         -0.3572289184172501, 
    577         -0.3376556177463400, 
    578         -0.3179351925907259, 
    579         -0.2980762356029071, 
    580         -0.2780873997969574, 
    581         -0.2579773947782034, 
    582         -0.2377549829482451, 
    583         -0.2174289756869712, 
    584         -0.1970082295132342, 
    585         -0.1765016422258567, 
    586         -0.1559181490266516, 
    587         -0.1352667186271445, 
    588         -0.1145563493406956, 
    589         -0.0937960651617229, 
    590         -0.0729949118337358, 
    591         -0.0521619529078925, 
    592         -0.0313062657937972, 
    593         -0.0104369378042598, 
    594         0.0104369378042598, 
    595         0.0313062657937972, 
    596         0.0521619529078925, 
    597         0.0729949118337358, 
    598         0.0937960651617229, 
    599         0.1145563493406956, 
    600         0.1352667186271445, 
    601         0.1559181490266516, 
    602         0.1765016422258567, 
    603         0.1970082295132342, 
    604         0.2174289756869712, 
    605         0.2377549829482451, 
    606         0.2579773947782034, 
    607         0.2780873997969574, 
    608         0.2980762356029071, 
    609         0.3179351925907259, 
    610         0.3376556177463400, 
    611         0.3572289184172501, 
    612         0.3766465660565522, 
    613         0.3959000999390257, 
    614         0.4149811308476706, 
    615         0.4338813447290861, 
    616         0.4525925063160997, 
    617         0.4711064627160663, 
    618         0.4894151469632753, 
    619         0.5075105815339176, 
    620         0.5253848818220803, 
    621         0.5430302595752546, 
    622         0.5604390262878617, 
    623         0.5776035965513142, 
    624         0.5945164913591590, 
    625         0.6111703413658551, 
    626         0.6275578900977726, 
    627         0.6436719971150083, 
    628         0.6595056411226444, 
    629         0.6750519230300931, 
    630         0.6903040689571928, 
    631         0.7052554331857488, 
    632         0.7198995010552305, 
    633         0.7342298918013638, 
    634         0.7482403613363824, 
    635         0.7619248049697269, 
    636         0.7752772600680049, 
    637         0.7882919086530552, 
    638         0.8009630799369827, 
    639         0.8132852527930605, 
    640         0.8252530581614230, 
    641         0.8368612813885015, 
    642         0.8481048644991847, 
    643         0.8589789084007133, 
    644         0.8694786750173527, 
    645         0.8795995893549102, 
    646         0.8893372414942055, 
    647         0.8986873885126239, 
    648         0.9076459563329236, 
    649         0.9162090414984952, 
    650         0.9243729128743136, 
    651         0.9321340132728527, 
    652         0.9394889610042837, 
    653         0.9464345513503147, 
    654         0.9529677579610971, 
    655         0.9590857341746905, 
    656         0.9647858142586956, 
    657         0.9700655145738374, 
    658         0.9749225346595943, 
    659         0.9793547582425894, 
    660         0.9833602541697529, 
    661         0.9869372772712794, 
    662         0.9900842691660192, 
    663         0.9927998590434373, 
    664         0.9950828645255290, 
    665         0.9969322929775997, 
    666         0.9983473449340834, 
    667         0.9993274305065947, 
    668         0.9998723404457334 
     502      -0.9998723404457334, 
     503      -0.9993274305065947, 
     504      -0.9983473449340834, 
     505      -0.9969322929775997, 
     506      -0.9950828645255290, 
     507      -0.9927998590434373, 
     508      -0.9900842691660192, 
     509      -0.9869372772712794, 
     510      -0.9833602541697529, 
     511      -0.9793547582425894, 
     512      -0.9749225346595943, 
     513      -0.9700655145738374, 
     514      -0.9647858142586956, 
     515      -0.9590857341746905, 
     516      -0.9529677579610971, 
     517      -0.9464345513503147, 
     518      -0.9394889610042837, 
     519      -0.9321340132728527, 
     520      -0.9243729128743136, 
     521      -0.9162090414984952, 
     522      -0.9076459563329236, 
     523      -0.8986873885126239, 
     524      -0.8893372414942055, 
     525      -0.8795995893549102, 
     526      -0.8694786750173527, 
     527      -0.8589789084007133, 
     528      -0.8481048644991847, 
     529      -0.8368612813885015, 
     530      -0.8252530581614230, 
     531      -0.8132852527930605, 
     532      -0.8009630799369827, 
     533      -0.7882919086530552, 
     534      -0.7752772600680049, 
     535      -0.7619248049697269, 
     536      -0.7482403613363824, 
     537      -0.7342298918013638, 
     538      -0.7198995010552305, 
     539      -0.7052554331857488, 
     540      -0.6903040689571928, 
     541      -0.6750519230300931, 
     542      -0.6595056411226444, 
     543      -0.6436719971150083, 
     544      -0.6275578900977726, 
     545      -0.6111703413658551, 
     546      -0.5945164913591590, 
     547      -0.5776035965513142, 
     548      -0.5604390262878617, 
     549      -0.5430302595752546, 
     550      -0.5253848818220803, 
     551      -0.5075105815339176, 
     552      -0.4894151469632753, 
     553      -0.4711064627160663, 
     554      -0.4525925063160997, 
     555      -0.4338813447290861, 
     556      -0.4149811308476706, 
     557      -0.3959000999390257, 
     558      -0.3766465660565522, 
     559      -0.3572289184172501, 
     560      -0.3376556177463400, 
     561      -0.3179351925907259, 
     562      -0.2980762356029071, 
     563      -0.2780873997969574, 
     564      -0.2579773947782034, 
     565      -0.2377549829482451, 
     566      -0.2174289756869712, 
     567      -0.1970082295132342, 
     568      -0.1765016422258567, 
     569      -0.1559181490266516, 
     570      -0.1352667186271445, 
     571      -0.1145563493406956, 
     572      -0.0937960651617229, 
     573      -0.0729949118337358, 
     574      -0.0521619529078925, 
     575      -0.0313062657937972, 
     576      -0.0104369378042598, 
     577      0.0104369378042598, 
     578      0.0313062657937972, 
     579      0.0521619529078925, 
     580      0.0729949118337358, 
     581      0.0937960651617229, 
     582      0.1145563493406956, 
     583      0.1352667186271445, 
     584      0.1559181490266516, 
     585      0.1765016422258567, 
     586      0.1970082295132342, 
     587      0.2174289756869712, 
     588      0.2377549829482451, 
     589      0.2579773947782034, 
     590      0.2780873997969574, 
     591      0.2980762356029071, 
     592      0.3179351925907259, 
     593      0.3376556177463400, 
     594      0.3572289184172501, 
     595      0.3766465660565522, 
     596      0.3959000999390257, 
     597      0.4149811308476706, 
     598      0.4338813447290861, 
     599      0.4525925063160997, 
     600      0.4711064627160663, 
     601      0.4894151469632753, 
     602      0.5075105815339176, 
     603      0.5253848818220803, 
     604      0.5430302595752546, 
     605      0.5604390262878617, 
     606      0.5776035965513142, 
     607      0.5945164913591590, 
     608      0.6111703413658551, 
     609      0.6275578900977726, 
     610      0.6436719971150083, 
     611      0.6595056411226444, 
     612      0.6750519230300931, 
     613      0.6903040689571928, 
     614      0.7052554331857488, 
     615      0.7198995010552305, 
     616      0.7342298918013638, 
     617      0.7482403613363824, 
     618      0.7619248049697269, 
     619      0.7752772600680049, 
     620      0.7882919086530552, 
     621      0.8009630799369827, 
     622      0.8132852527930605, 
     623      0.8252530581614230, 
     624      0.8368612813885015, 
     625      0.8481048644991847, 
     626      0.8589789084007133, 
     627      0.8694786750173527, 
     628      0.8795995893549102, 
     629      0.8893372414942055, 
     630      0.8986873885126239, 
     631      0.9076459563329236, 
     632      0.9162090414984952, 
     633      0.9243729128743136, 
     634      0.9321340132728527, 
     635      0.9394889610042837, 
     636      0.9464345513503147, 
     637      0.9529677579610971, 
     638      0.9590857341746905, 
     639      0.9647858142586956, 
     640      0.9700655145738374, 
     641      0.9749225346595943, 
     642      0.9793547582425894, 
     643      0.9833602541697529, 
     644      0.9869372772712794, 
     645      0.9900842691660192, 
     646      0.9927998590434373, 
     647      0.9950828645255290, 
     648      0.9969322929775997, 
     649      0.9983473449340834, 
     650      0.9993274305065947, 
     651      0.9998723404457334 
    669652]) 
    670653 
    671654Gauss150Wt = np.array([ 
    672         0.0003276086705538, 
    673         0.0007624720924706, 
    674         0.0011976474864367, 
    675         0.0016323569986067, 
    676         0.0020663664924131, 
    677         0.0024994789888943, 
    678         0.0029315036836558, 
    679         0.0033622516236779, 
    680         0.0037915348363451, 
    681         0.0042191661429919, 
    682         0.0046449591497966, 
    683         0.0050687282939456, 
    684         0.0054902889094487, 
    685         0.0059094573005900, 
    686         0.0063260508184704, 
    687         0.0067398879387430, 
    688         0.0071507883396855, 
    689         0.0075585729801782, 
    690         0.0079630641773633, 
    691         0.0083640856838475, 
    692         0.0087614627643580, 
    693         0.0091550222717888, 
    694         0.0095445927225849, 
    695         0.0099300043714212, 
    696         0.0103110892851360, 
    697         0.0106876814158841, 
    698         0.0110596166734735, 
    699         0.0114267329968529, 
    700         0.0117888704247183, 
    701         0.0121458711652067, 
    702         0.0124975796646449, 
    703         0.0128438426753249, 
    704         0.0131845093222756, 
    705         0.0135194311690004, 
    706         0.0138484622795371, 
    707         0.0141714592928592, 
    708         0.0144882814685445, 
    709         0.0147987907597169, 
    710         0.0151028518701744, 
    711         0.0154003323133401, 
    712         0.0156911024699895, 
    713         0.0159750356447283, 
    714         0.0162520081211971, 
    715         0.0165218992159766, 
    716         0.0167845913311726, 
    717         0.0170399700056559, 
    718         0.0172879239649355, 
    719         0.0175283451696437, 
    720         0.0177611288626114, 
    721         0.0179861736145128, 
    722         0.0182033813680609, 
    723         0.0184126574807331, 
    724         0.0186139107660094, 
    725         0.0188070535331042, 
    726         0.0189920016251754, 
    727         0.0191686744559934, 
    728         0.0193369950450545, 
    729         0.0194968900511231, 
    730         0.0196482898041878, 
    731         0.0197911283358190, 
    732         0.0199253434079123, 
    733         0.0200508765398072, 
    734         0.0201676730337687, 
    735         0.0202756819988200, 
    736         0.0203748563729175, 
    737         0.0204651529434560, 
    738         0.0205465323660984, 
    739         0.0206189591819181, 
    740         0.0206824018328499, 
    741         0.0207368326754401, 
    742         0.0207822279928917, 
    743         0.0208185680053983, 
    744         0.0208458368787627, 
    745         0.0208640227312962, 
    746         0.0208731176389954, 
    747         0.0208731176389954, 
    748         0.0208640227312962, 
    749         0.0208458368787627, 
    750         0.0208185680053983, 
    751         0.0207822279928917, 
    752         0.0207368326754401, 
    753         0.0206824018328499, 
    754         0.0206189591819181, 
    755         0.0205465323660984, 
    756         0.0204651529434560, 
    757         0.0203748563729175, 
    758         0.0202756819988200, 
    759         0.0201676730337687, 
    760         0.0200508765398072, 
    761         0.0199253434079123, 
    762         0.0197911283358190, 
    763         0.0196482898041878, 
    764         0.0194968900511231, 
    765         0.0193369950450545, 
    766         0.0191686744559934, 
    767         0.0189920016251754, 
    768         0.0188070535331042, 
    769         0.0186139107660094, 
    770         0.0184126574807331, 
    771         0.0182033813680609, 
    772         0.0179861736145128, 
    773         0.0177611288626114, 
    774         0.0175283451696437, 
    775         0.0172879239649355, 
    776         0.0170399700056559, 
    777         0.0167845913311726, 
    778         0.0165218992159766, 
    779         0.0162520081211971, 
    780         0.0159750356447283, 
    781         0.0156911024699895, 
    782         0.0154003323133401, 
    783         0.0151028518701744, 
    784         0.0147987907597169, 
    785         0.0144882814685445, 
    786         0.0141714592928592, 
    787         0.0138484622795371, 
    788         0.0135194311690004, 
    789         0.0131845093222756, 
    790         0.0128438426753249, 
    791         0.0124975796646449, 
    792         0.0121458711652067, 
    793         0.0117888704247183, 
    794         0.0114267329968529, 
    795         0.0110596166734735, 
    796         0.0106876814158841, 
    797         0.0103110892851360, 
    798         0.0099300043714212, 
    799         0.0095445927225849, 
    800         0.0091550222717888, 
    801         0.0087614627643580, 
    802         0.0083640856838475, 
    803         0.0079630641773633, 
    804         0.0075585729801782, 
    805         0.0071507883396855, 
    806         0.0067398879387430, 
    807         0.0063260508184704, 
    808         0.0059094573005900, 
    809         0.0054902889094487, 
    810         0.0050687282939456, 
    811         0.0046449591497966, 
    812         0.0042191661429919, 
    813         0.0037915348363451, 
    814         0.0033622516236779, 
    815         0.0029315036836558, 
    816         0.0024994789888943, 
    817         0.0020663664924131, 
    818         0.0016323569986067, 
    819         0.0011976474864367, 
    820         0.0007624720924706, 
    821         0.0003276086705538 
     655      0.0003276086705538, 
     656      0.0007624720924706, 
     657      0.0011976474864367, 
     658      0.0016323569986067, 
     659      0.0020663664924131, 
     660      0.0024994789888943, 
     661      0.0029315036836558, 
     662      0.0033622516236779, 
     663      0.0037915348363451, 
     664      0.0042191661429919, 
     665      0.0046449591497966, 
     666      0.0050687282939456, 
     667      0.0054902889094487, 
     668      0.0059094573005900, 
     669      0.0063260508184704, 
     670      0.0067398879387430, 
     671      0.0071507883396855, 
     672      0.0075585729801782, 
     673      0.0079630641773633, 
     674      0.0083640856838475, 
     675      0.0087614627643580, 
     676      0.0091550222717888, 
     677      0.0095445927225849, 
     678      0.0099300043714212, 
     679      0.0103110892851360, 
     680      0.0106876814158841, 
     681      0.0110596166734735, 
     682      0.0114267329968529, 
     683      0.0117888704247183, 
     684      0.0121458711652067, 
     685      0.0124975796646449, 
     686      0.0128438426753249, 
     687      0.0131845093222756, 
     688      0.0135194311690004, 
     689      0.0138484622795371, 
     690      0.0141714592928592, 
     691      0.0144882814685445, 
     692      0.0147987907597169, 
     693      0.0151028518701744, 
     694      0.0154003323133401, 
     695      0.0156911024699895, 
     696      0.0159750356447283, 
     697      0.0162520081211971, 
     698      0.0165218992159766, 
     699      0.0167845913311726, 
     700      0.0170399700056559, 
     701      0.0172879239649355, 
     702      0.0175283451696437, 
     703      0.0177611288626114, 
     704      0.0179861736145128, 
     705      0.0182033813680609, 
     706      0.0184126574807331, 
     707      0.0186139107660094, 
     708      0.0188070535331042, 
     709      0.0189920016251754, 
     710      0.0191686744559934, 
     711      0.0193369950450545, 
     712      0.0194968900511231, 
     713      0.0196482898041878, 
     714      0.0197911283358190, 
     715      0.0199253434079123, 
     716      0.0200508765398072, 
     717      0.0201676730337687, 
     718      0.0202756819988200, 
     719      0.0203748563729175, 
     720      0.0204651529434560, 
     721      0.0205465323660984, 
     722      0.0206189591819181, 
     723      0.0206824018328499, 
     724      0.0207368326754401, 
     725      0.0207822279928917, 
     726      0.0208185680053983, 
     727      0.0208458368787627, 
     728      0.0208640227312962, 
     729      0.0208731176389954, 
     730      0.0208731176389954, 
     731      0.0208640227312962, 
     732      0.0208458368787627, 
     733      0.0208185680053983, 
     734      0.0207822279928917, 
     735      0.0207368326754401, 
     736      0.0206824018328499, 
     737      0.0206189591819181, 
     738      0.0205465323660984, 
     739      0.0204651529434560, 
     740      0.0203748563729175, 
     741      0.0202756819988200, 
     742      0.0201676730337687, 
     743      0.0200508765398072, 
     744      0.0199253434079123, 
     745      0.0197911283358190, 
     746      0.0196482898041878, 
     747      0.0194968900511231, 
     748      0.0193369950450545, 
     749      0.0191686744559934, 
     750      0.0189920016251754, 
     751      0.0188070535331042, 
     752      0.0186139107660094, 
     753      0.0184126574807331, 
     754      0.0182033813680609, 
     755      0.0179861736145128, 
     756      0.0177611288626114, 
     757      0.0175283451696437, 
     758      0.0172879239649355, 
     759      0.0170399700056559, 
     760      0.0167845913311726, 
     761      0.0165218992159766, 
     762      0.0162520081211971, 
     763      0.0159750356447283, 
     764      0.0156911024699895, 
     765      0.0154003323133401, 
     766      0.0151028518701744, 
     767      0.0147987907597169, 
     768      0.0144882814685445, 
     769      0.0141714592928592, 
     770      0.0138484622795371, 
     771      0.0135194311690004, 
     772      0.0131845093222756, 
     773      0.0128438426753249, 
     774      0.0124975796646449, 
     775      0.0121458711652067, 
     776      0.0117888704247183, 
     777      0.0114267329968529, 
     778      0.0110596166734735, 
     779      0.0106876814158841, 
     780      0.0103110892851360, 
     781      0.0099300043714212, 
     782      0.0095445927225849, 
     783      0.0091550222717888, 
     784      0.0087614627643580, 
     785      0.0083640856838475, 
     786      0.0079630641773633, 
     787      0.0075585729801782, 
     788      0.0071507883396855, 
     789      0.0067398879387430, 
     790      0.0063260508184704, 
     791      0.0059094573005900, 
     792      0.0054902889094487, 
     793      0.0050687282939456, 
     794      0.0046449591497966, 
     795      0.0042191661429919, 
     796      0.0037915348363451, 
     797      0.0033622516236779, 
     798      0.0029315036836558, 
     799      0.0024994789888943, 
     800      0.0020663664924131, 
     801      0.0016323569986067, 
     802      0.0011976474864367, 
     803      0.0007624720924706, 
     804      0.0003276086705538 
    822805]) 
Note: See TracChangeset for help on using the changeset viewer.