Changeset 994d77f in sasmodels for sasmodels/models/lib


Ignore:
Timestamp:
Oct 30, 2014 12:33:53 PM (10 years ago)
Author:
Paul Kienzle <pkienzle@…>
Branches:
master, core_shell_microgels, costrafo411, magnetic_model, release_v0.94, release_v0.95, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests
Children:
ef2861b
Parents:
d087487b
Message:

Convert double to float rather than using real

Location:
sasmodels/models/lib
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • sasmodels/models/lib/J1.c

    r14de349 r994d77f  
    1 real J1(real x); 
    2 real J1(real x) 
     1double J1(double x); 
     2double J1(double x) 
    33{ 
    4   const real ax = fabs(x); 
    5   if (ax < REAL(8.0)) { 
    6     const real y = x*x; 
    7     const real ans1 = x*(REAL(72362614232.0) 
    8               + y*(REAL(-7895059235.0) 
    9               + y*(REAL(242396853.1) 
    10               + y*(REAL(-2972611.439) 
    11               + y*(REAL(15704.48260) 
    12               + y*(REAL(-30.16036606))))))); 
    13     const real ans2 = REAL(144725228442.0) 
    14               + y*(REAL(2300535178.0) 
    15               + y*(REAL(18583304.74) 
    16               + y*(REAL(99447.43394) 
    17               + y*(REAL(376.9991397) 
     4  const double ax = fabs(x); 
     5  if (ax < 8.0) { 
     6    const double y = x*x; 
     7    const double ans1 = x*(72362614232.0 
     8              + y*(-7895059235.0 
     9              + y*(242396853.1 
     10              + y*(-2972611.439 
     11              + y*(15704.48260 
     12              + y*(-30.16036606)))))); 
     13    const double ans2 = 144725228442.0 
     14              + y*(2300535178.0 
     15              + y*(18583304.74 
     16              + y*(99447.43394 
     17              + y*(376.9991397 
    1818              + y)))); 
    1919    return ans1/ans2; 
    2020  } else { 
    21     const real y = REAL(64.0)/(ax*ax); 
    22     const real xx = ax - REAL(2.356194491); 
    23     const real ans1 = REAL(1.0) 
    24               + y*(REAL(0.183105e-2) 
    25               + y*(REAL(-0.3516396496e-4) 
    26               + y*(REAL(0.2457520174e-5) 
    27               + y*REAL(-0.240337019e-6)))); 
    28     const real ans2 = REAL(0.04687499995) 
    29               + y*(REAL(-0.2002690873e-3) 
    30               + y*(REAL(0.8449199096e-5) 
    31               + y*(REAL(-0.88228987e-6) 
    32               + y*REAL(0.105787412e-6)))); 
    33     real sn,cn; 
     21    const double y = 64.0/(ax*ax); 
     22    const double xx = ax - 2.356194491; 
     23    const double ans1 = 1.0 
     24              + y*(0.183105e-2 
     25              + y*(-0.3516396496e-4 
     26              + y*(0.2457520174e-5 
     27              + y*-0.240337019e-6))); 
     28    const double ans2 = 0.04687499995 
     29              + y*(-0.2002690873e-3 
     30              + y*(0.8449199096e-5 
     31              + y*(-0.88228987e-6 
     32              + y*0.105787412e-6))); 
     33    double sn,cn; 
    3434    SINCOS(xx, sn, cn); 
    35     const real ans = sqrt(REAL(0.636619772)/ax) * (cn*ans1 - (REAL(8.0)/ax)*sn*ans2); 
    36     return (x < REAL(0.0)) ? -ans : ans; 
     35    const double ans = sqrt(0.636619772/ax) * (cn*ans1 - (8.0/ax)*sn*ans2); 
     36    return (x < 0.0) ? -ans : ans; 
    3737  } 
    3838} 

  • sasmodels/models/lib/gauss150.c

    r14de349 r994d77f  
    99 
    1010// Gaussians 
    11 constant real Gauss150Z[150]={ 
    12         REAL(-0.9998723404457334), 
    13         REAL(-0.9993274305065947), 
    14         REAL(-0.9983473449340834), 
    15         REAL(-0.9969322929775997), 
    16         REAL(-0.9950828645255290), 
    17         REAL(-0.9927998590434373), 
    18         REAL(-0.9900842691660192), 
    19         REAL(-0.9869372772712794), 
    20         REAL(-0.9833602541697529), 
    21         REAL(-0.9793547582425894), 
    22         REAL(-0.9749225346595943), 
    23         REAL(-0.9700655145738374), 
    24         REAL(-0.9647858142586956), 
    25         REAL(-0.9590857341746905), 
    26         REAL(-0.9529677579610971), 
    27         REAL(-0.9464345513503147), 
    28         REAL(-0.9394889610042837), 
    29         REAL(-0.9321340132728527), 
    30         REAL(-0.9243729128743136), 
    31         REAL(-0.9162090414984952), 
    32         REAL(-0.9076459563329236), 
    33         REAL(-0.8986873885126239), 
    34         REAL(-0.8893372414942055), 
    35         REAL(-0.8795995893549102), 
    36         REAL(-0.8694786750173527), 
    37         REAL(-0.8589789084007133), 
    38         REAL(-0.8481048644991847), 
    39         REAL(-0.8368612813885015), 
    40         REAL(-0.8252530581614230), 
    41         REAL(-0.8132852527930605), 
    42         REAL(-0.8009630799369827), 
    43         REAL(-0.7882919086530552), 
    44         REAL(-0.7752772600680049), 
    45         REAL(-0.7619248049697269), 
    46         REAL(-0.7482403613363824), 
    47         REAL(-0.7342298918013638), 
    48         REAL(-0.7198995010552305), 
    49         REAL(-0.7052554331857488), 
    50         REAL(-0.6903040689571928), 
    51         REAL(-0.6750519230300931), 
    52         REAL(-0.6595056411226444), 
    53         REAL(-0.6436719971150083), 
    54         REAL(-0.6275578900977726), 
    55         REAL(-0.6111703413658551), 
    56         REAL(-0.5945164913591590), 
    57         REAL(-0.5776035965513142), 
    58         REAL(-0.5604390262878617), 
    59         REAL(-0.5430302595752546), 
    60         REAL(-0.5253848818220803), 
    61         REAL(-0.5075105815339176), 
    62         REAL(-0.4894151469632753), 
    63         REAL(-0.4711064627160663), 
    64         REAL(-0.4525925063160997), 
    65         REAL(-0.4338813447290861), 
    66         REAL(-0.4149811308476706), 
    67         REAL(-0.3959000999390257), 
    68         REAL(-0.3766465660565522), 
    69         REAL(-0.3572289184172501), 
    70         REAL(-0.3376556177463400), 
    71         REAL(-0.3179351925907259), 
    72         REAL(-0.2980762356029071), 
    73         REAL(-0.2780873997969574), 
    74         REAL(-0.2579773947782034), 
    75         REAL(-0.2377549829482451), 
    76         REAL(-0.2174289756869712), 
    77         REAL(-0.1970082295132342), 
    78         REAL(-0.1765016422258567), 
    79         REAL(-0.1559181490266516), 
    80         REAL(-0.1352667186271445), 
    81         REAL(-0.1145563493406956), 
    82         REAL(-0.0937960651617229), 
    83         REAL(-0.0729949118337358), 
    84         REAL(-0.0521619529078925), 
    85         REAL(-0.0313062657937972), 
    86         REAL(-0.0104369378042598), 
    87         REAL(0.0104369378042598), 
    88         REAL(0.0313062657937972), 
    89         REAL(0.0521619529078925), 
    90         REAL(0.0729949118337358), 
    91         REAL(0.0937960651617229), 
    92         REAL(0.1145563493406956), 
    93         REAL(0.1352667186271445), 
    94         REAL(0.1559181490266516), 
    95         REAL(0.1765016422258567), 
    96         REAL(0.1970082295132342), 
    97         REAL(0.2174289756869712), 
    98         REAL(0.2377549829482451), 
    99         REAL(0.2579773947782034), 
    100         REAL(0.2780873997969574), 
    101         REAL(0.2980762356029071), 
    102         REAL(0.3179351925907259), 
    103         REAL(0.3376556177463400), 
    104         REAL(0.3572289184172501), 
    105         REAL(0.3766465660565522), 
    106         REAL(0.3959000999390257), 
    107         REAL(0.4149811308476706), 
    108         REAL(0.4338813447290861), 
    109         REAL(0.4525925063160997), 
    110         REAL(0.4711064627160663), 
    111         REAL(0.4894151469632753), 
    112         REAL(0.5075105815339176), 
    113         REAL(0.5253848818220803), 
    114         REAL(0.5430302595752546), 
    115         REAL(0.5604390262878617), 
    116         REAL(0.5776035965513142), 
    117         REAL(0.5945164913591590), 
    118         REAL(0.6111703413658551), 
    119         REAL(0.6275578900977726), 
    120         REAL(0.6436719971150083), 
    121         REAL(0.6595056411226444), 
    122         REAL(0.6750519230300931), 
    123         REAL(0.6903040689571928), 
    124         REAL(0.7052554331857488), 
    125         REAL(0.7198995010552305), 
    126         REAL(0.7342298918013638), 
    127         REAL(0.7482403613363824), 
    128         REAL(0.7619248049697269), 
    129         REAL(0.7752772600680049), 
    130         REAL(0.7882919086530552), 
    131         REAL(0.8009630799369827), 
    132         REAL(0.8132852527930605), 
    133         REAL(0.8252530581614230), 
    134         REAL(0.8368612813885015), 
    135         REAL(0.8481048644991847), 
    136         REAL(0.8589789084007133), 
    137         REAL(0.8694786750173527), 
    138         REAL(0.8795995893549102), 
    139         REAL(0.8893372414942055), 
    140         REAL(0.8986873885126239), 
    141         REAL(0.9076459563329236), 
    142         REAL(0.9162090414984952), 
    143         REAL(0.9243729128743136), 
    144         REAL(0.9321340132728527), 
    145         REAL(0.9394889610042837), 
    146         REAL(0.9464345513503147), 
    147         REAL(0.9529677579610971), 
    148         REAL(0.9590857341746905), 
    149         REAL(0.9647858142586956), 
    150         REAL(0.9700655145738374), 
    151         REAL(0.9749225346595943), 
    152         REAL(0.9793547582425894), 
    153         REAL(0.9833602541697529), 
    154         REAL(0.9869372772712794), 
    155         REAL(0.9900842691660192), 
    156         REAL(0.9927998590434373), 
    157         REAL(0.9950828645255290), 
    158         REAL(0.9969322929775997), 
    159         REAL(0.9983473449340834), 
    160         REAL(0.9993274305065947), 
    161         REAL(0.9998723404457334) 
     11constant double Gauss150Z[150]={ 
     12        -0.9998723404457334, 
     13        -0.9993274305065947, 
     14        -0.9983473449340834, 
     15        -0.9969322929775997, 
     16        -0.9950828645255290, 
     17        -0.9927998590434373, 
     18        -0.9900842691660192, 
     19        -0.9869372772712794, 
     20        -0.9833602541697529, 
     21        -0.9793547582425894, 
     22        -0.9749225346595943, 
     23        -0.9700655145738374, 
     24        -0.9647858142586956, 
     25        -0.9590857341746905, 
     26        -0.9529677579610971, 
     27        -0.9464345513503147, 
     28        -0.9394889610042837, 
     29        -0.9321340132728527, 
     30        -0.9243729128743136, 
     31        -0.9162090414984952, 
     32        -0.9076459563329236, 
     33        -0.8986873885126239, 
     34        -0.8893372414942055, 
     35        -0.8795995893549102, 
     36        -0.8694786750173527, 
     37        -0.8589789084007133, 
     38        -0.8481048644991847, 
     39        -0.8368612813885015, 
     40        -0.8252530581614230, 
     41        -0.8132852527930605, 
     42        -0.8009630799369827, 
     43        -0.7882919086530552, 
     44        -0.7752772600680049, 
     45        -0.7619248049697269, 
     46        -0.7482403613363824, 
     47        -0.7342298918013638, 
     48        -0.7198995010552305, 
     49        -0.7052554331857488, 
     50        -0.6903040689571928, 
     51        -0.6750519230300931, 
     52        -0.6595056411226444, 
     53        -0.6436719971150083, 
     54        -0.6275578900977726, 
     55        -0.6111703413658551, 
     56        -0.5945164913591590, 
     57        -0.5776035965513142, 
     58        -0.5604390262878617, 
     59        -0.5430302595752546, 
     60        -0.5253848818220803, 
     61        -0.5075105815339176, 
     62        -0.4894151469632753, 
     63        -0.4711064627160663, 
     64        -0.4525925063160997, 
     65        -0.4338813447290861, 
     66        -0.4149811308476706, 
     67        -0.3959000999390257, 
     68        -0.3766465660565522, 
     69        -0.3572289184172501, 
     70        -0.3376556177463400, 
     71        -0.3179351925907259, 
     72        -0.2980762356029071, 
     73        -0.2780873997969574, 
     74        -0.2579773947782034, 
     75        -0.2377549829482451, 
     76        -0.2174289756869712, 
     77        -0.1970082295132342, 
     78        -0.1765016422258567, 
     79        -0.1559181490266516, 
     80        -0.1352667186271445, 
     81        -0.1145563493406956, 
     82        -0.0937960651617229, 
     83        -0.0729949118337358, 
     84        -0.0521619529078925, 
     85        -0.0313062657937972, 
     86        -0.0104369378042598, 
     87        0.0104369378042598, 
     88        0.0313062657937972, 
     89        0.0521619529078925, 
     90        0.0729949118337358, 
     91        0.0937960651617229, 
     92        0.1145563493406956, 
     93        0.1352667186271445, 
     94        0.1559181490266516, 
     95        0.1765016422258567, 
     96        0.1970082295132342, 
     97        0.2174289756869712, 
     98        0.2377549829482451, 
     99        0.2579773947782034, 
     100        0.2780873997969574, 
     101        0.2980762356029071, 
     102        0.3179351925907259, 
     103        0.3376556177463400, 
     104        0.3572289184172501, 
     105        0.3766465660565522, 
     106        0.3959000999390257, 
     107        0.4149811308476706, 
     108        0.4338813447290861, 
     109        0.4525925063160997, 
     110        0.4711064627160663, 
     111        0.4894151469632753, 
     112        0.5075105815339176, 
     113        0.5253848818220803, 
     114        0.5430302595752546, 
     115        0.5604390262878617, 
     116        0.5776035965513142, 
     117        0.5945164913591590, 
     118        0.6111703413658551, 
     119        0.6275578900977726, 
     120        0.6436719971150083, 
     121        0.6595056411226444, 
     122        0.6750519230300931, 
     123        0.6903040689571928, 
     124        0.7052554331857488, 
     125        0.7198995010552305, 
     126        0.7342298918013638, 
     127        0.7482403613363824, 
     128        0.7619248049697269, 
     129        0.7752772600680049, 
     130        0.7882919086530552, 
     131        0.8009630799369827, 
     132        0.8132852527930605, 
     133        0.8252530581614230, 
     134        0.8368612813885015, 
     135        0.8481048644991847, 
     136        0.8589789084007133, 
     137        0.8694786750173527, 
     138        0.8795995893549102, 
     139        0.8893372414942055, 
     140        0.8986873885126239, 
     141        0.9076459563329236, 
     142        0.9162090414984952, 
     143        0.9243729128743136, 
     144        0.9321340132728527, 
     145        0.9394889610042837, 
     146        0.9464345513503147, 
     147        0.9529677579610971, 
     148        0.9590857341746905, 
     149        0.9647858142586956, 
     150        0.9700655145738374, 
     151        0.9749225346595943, 
     152        0.9793547582425894, 
     153        0.9833602541697529, 
     154        0.9869372772712794, 
     155        0.9900842691660192, 
     156        0.9927998590434373, 
     157        0.9950828645255290, 
     158        0.9969322929775997, 
     159        0.9983473449340834, 
     160        0.9993274305065947, 
     161        0.9998723404457334 
    162162}; 
    163163 
    164 constant real Gauss150Wt[150]={ 
    165         REAL(0.0003276086705538), 
    166         REAL(0.0007624720924706), 
    167         REAL(0.0011976474864367), 
    168         REAL(0.0016323569986067), 
    169         REAL(0.0020663664924131), 
    170         REAL(0.0024994789888943), 
    171         REAL(0.0029315036836558), 
    172         REAL(0.0033622516236779), 
    173         REAL(0.0037915348363451), 
    174         REAL(0.0042191661429919), 
    175         REAL(0.0046449591497966), 
    176         REAL(0.0050687282939456), 
    177         REAL(0.0054902889094487), 
    178         REAL(0.0059094573005900), 
    179         REAL(0.0063260508184704), 
    180         REAL(0.0067398879387430), 
    181         REAL(0.0071507883396855), 
    182         REAL(0.0075585729801782), 
    183         REAL(0.0079630641773633), 
    184         REAL(0.0083640856838475), 
    185         REAL(0.0087614627643580), 
    186         REAL(0.0091550222717888), 
    187         REAL(0.0095445927225849), 
    188         REAL(0.0099300043714212), 
    189         REAL(0.0103110892851360), 
    190         REAL(0.0106876814158841), 
    191         REAL(0.0110596166734735), 
    192         REAL(0.0114267329968529), 
    193         REAL(0.0117888704247183), 
    194         REAL(0.0121458711652067), 
    195         REAL(0.0124975796646449), 
    196         REAL(0.0128438426753249), 
    197         REAL(0.0131845093222756), 
    198         REAL(0.0135194311690004), 
    199         REAL(0.0138484622795371), 
    200         REAL(0.0141714592928592), 
    201         REAL(0.0144882814685445), 
    202         REAL(0.0147987907597169), 
    203         REAL(0.0151028518701744), 
    204         REAL(0.0154003323133401), 
    205         REAL(0.0156911024699895), 
    206         REAL(0.0159750356447283), 
    207         REAL(0.0162520081211971), 
    208         REAL(0.0165218992159766), 
    209         REAL(0.0167845913311726), 
    210         REAL(0.0170399700056559), 
    211         REAL(0.0172879239649355), 
    212         REAL(0.0175283451696437), 
    213         REAL(0.0177611288626114), 
    214         REAL(0.0179861736145128), 
    215         REAL(0.0182033813680609), 
    216         REAL(0.0184126574807331), 
    217         REAL(0.0186139107660094), 
    218         REAL(0.0188070535331042), 
    219         REAL(0.0189920016251754), 
    220         REAL(0.0191686744559934), 
    221         REAL(0.0193369950450545), 
    222         REAL(0.0194968900511231), 
    223         REAL(0.0196482898041878), 
    224         REAL(0.0197911283358190), 
    225         REAL(0.0199253434079123), 
    226         REAL(0.0200508765398072), 
    227         REAL(0.0201676730337687), 
    228         REAL(0.0202756819988200), 
    229         REAL(0.0203748563729175), 
    230         REAL(0.0204651529434560), 
    231         REAL(0.0205465323660984), 
    232         REAL(0.0206189591819181), 
    233         REAL(0.0206824018328499), 
    234         REAL(0.0207368326754401), 
    235         REAL(0.0207822279928917), 
    236         REAL(0.0208185680053983), 
    237         REAL(0.0208458368787627), 
    238         REAL(0.0208640227312962), 
    239         REAL(0.0208731176389954), 
    240         REAL(0.0208731176389954), 
    241         REAL(0.0208640227312962), 
    242         REAL(0.0208458368787627), 
    243         REAL(0.0208185680053983), 
    244         REAL(0.0207822279928917), 
    245         REAL(0.0207368326754401), 
    246         REAL(0.0206824018328499), 
    247         REAL(0.0206189591819181), 
    248         REAL(0.0205465323660984), 
    249         REAL(0.0204651529434560), 
    250         REAL(0.0203748563729175), 
    251         REAL(0.0202756819988200), 
    252         REAL(0.0201676730337687), 
    253         REAL(0.0200508765398072), 
    254         REAL(0.0199253434079123), 
    255         REAL(0.0197911283358190), 
    256         REAL(0.0196482898041878), 
    257         REAL(0.0194968900511231), 
    258         REAL(0.0193369950450545), 
    259         REAL(0.0191686744559934), 
    260         REAL(0.0189920016251754), 
    261         REAL(0.0188070535331042), 
    262         REAL(0.0186139107660094), 
    263         REAL(0.0184126574807331), 
    264         REAL(0.0182033813680609), 
    265         REAL(0.0179861736145128), 
    266         REAL(0.0177611288626114), 
    267         REAL(0.0175283451696437), 
    268         REAL(0.0172879239649355), 
    269         REAL(0.0170399700056559), 
    270         REAL(0.0167845913311726), 
    271         REAL(0.0165218992159766), 
    272         REAL(0.0162520081211971), 
    273         REAL(0.0159750356447283), 
    274         REAL(0.0156911024699895), 
    275         REAL(0.0154003323133401), 
    276         REAL(0.0151028518701744), 
    277         REAL(0.0147987907597169), 
    278         REAL(0.0144882814685445), 
    279         REAL(0.0141714592928592), 
    280         REAL(0.0138484622795371), 
    281         REAL(0.0135194311690004), 
    282         REAL(0.0131845093222756), 
    283         REAL(0.0128438426753249), 
    284         REAL(0.0124975796646449), 
    285         REAL(0.0121458711652067), 
    286         REAL(0.0117888704247183), 
    287         REAL(0.0114267329968529), 
    288         REAL(0.0110596166734735), 
    289         REAL(0.0106876814158841), 
    290         REAL(0.0103110892851360), 
    291         REAL(0.0099300043714212), 
    292         REAL(0.0095445927225849), 
    293         REAL(0.0091550222717888), 
    294         REAL(0.0087614627643580), 
    295         REAL(0.0083640856838475), 
    296         REAL(0.0079630641773633), 
    297         REAL(0.0075585729801782), 
    298         REAL(0.0071507883396855), 
    299         REAL(0.0067398879387430), 
    300         REAL(0.0063260508184704), 
    301         REAL(0.0059094573005900), 
    302         REAL(0.0054902889094487), 
    303         REAL(0.0050687282939456), 
    304         REAL(0.0046449591497966), 
    305         REAL(0.0042191661429919), 
    306         REAL(0.0037915348363451), 
    307         REAL(0.0033622516236779), 
    308         REAL(0.0029315036836558), 
    309         REAL(0.0024994789888943), 
    310         REAL(0.0020663664924131), 
    311         REAL(0.0016323569986067), 
    312         REAL(0.0011976474864367), 
    313         REAL(0.0007624720924706), 
    314         REAL(0.0003276086705538) 
     164constant double Gauss150Wt[150]={ 
     165        0.0003276086705538, 
     166        0.0007624720924706, 
     167        0.0011976474864367, 
     168        0.0016323569986067, 
     169        0.0020663664924131, 
     170        0.0024994789888943, 
     171        0.0029315036836558, 
     172        0.0033622516236779, 
     173        0.0037915348363451, 
     174        0.0042191661429919, 
     175        0.0046449591497966, 
     176        0.0050687282939456, 
     177        0.0054902889094487, 
     178        0.0059094573005900, 
     179        0.0063260508184704, 
     180        0.0067398879387430, 
     181        0.0071507883396855, 
     182        0.0075585729801782, 
     183        0.0079630641773633, 
     184        0.0083640856838475, 
     185        0.0087614627643580, 
     186        0.0091550222717888, 
     187        0.0095445927225849, 
     188        0.0099300043714212, 
     189        0.0103110892851360, 
     190        0.0106876814158841, 
     191        0.0110596166734735, 
     192        0.0114267329968529, 
     193        0.0117888704247183, 
     194        0.0121458711652067, 
     195        0.0124975796646449, 
     196        0.0128438426753249, 
     197        0.0131845093222756, 
     198        0.0135194311690004, 
     199        0.0138484622795371, 
     200        0.0141714592928592, 
     201        0.0144882814685445, 
     202        0.0147987907597169, 
     203        0.0151028518701744, 
     204        0.0154003323133401, 
     205        0.0156911024699895, 
     206        0.0159750356447283, 
     207        0.0162520081211971, 
     208        0.0165218992159766, 
     209        0.0167845913311726, 
     210        0.0170399700056559, 
     211        0.0172879239649355, 
     212        0.0175283451696437, 
     213        0.0177611288626114, 
     214        0.0179861736145128, 
     215        0.0182033813680609, 
     216        0.0184126574807331, 
     217        0.0186139107660094, 
     218        0.0188070535331042, 
     219        0.0189920016251754, 
     220        0.0191686744559934, 
     221        0.0193369950450545, 
     222        0.0194968900511231, 
     223        0.0196482898041878, 
     224        0.0197911283358190, 
     225        0.0199253434079123, 
     226        0.0200508765398072, 
     227        0.0201676730337687, 
     228        0.0202756819988200, 
     229        0.0203748563729175, 
     230        0.0204651529434560, 
     231        0.0205465323660984, 
     232        0.0206189591819181, 
     233        0.0206824018328499, 
     234        0.0207368326754401, 
     235        0.0207822279928917, 
     236        0.0208185680053983, 
     237        0.0208458368787627, 
     238        0.0208640227312962, 
     239        0.0208731176389954, 
     240        0.0208731176389954, 
     241        0.0208640227312962, 
     242        0.0208458368787627, 
     243        0.0208185680053983, 
     244        0.0207822279928917, 
     245        0.0207368326754401, 
     246        0.0206824018328499, 
     247        0.0206189591819181, 
     248        0.0205465323660984, 
     249        0.0204651529434560, 
     250        0.0203748563729175, 
     251        0.0202756819988200, 
     252        0.0201676730337687, 
     253        0.0200508765398072, 
     254        0.0199253434079123, 
     255        0.0197911283358190, 
     256        0.0196482898041878, 
     257        0.0194968900511231, 
     258        0.0193369950450545, 
     259        0.0191686744559934, 
     260        0.0189920016251754, 
     261        0.0188070535331042, 
     262        0.0186139107660094, 
     263        0.0184126574807331, 
     264        0.0182033813680609, 
     265        0.0179861736145128, 
     266        0.0177611288626114, 
     267        0.0175283451696437, 
     268        0.0172879239649355, 
     269        0.0170399700056559, 
     270        0.0167845913311726, 
     271        0.0165218992159766, 
     272        0.0162520081211971, 
     273        0.0159750356447283, 
     274        0.0156911024699895, 
     275        0.0154003323133401, 
     276        0.0151028518701744, 
     277        0.0147987907597169, 
     278        0.0144882814685445, 
     279        0.0141714592928592, 
     280        0.0138484622795371, 
     281        0.0135194311690004, 
     282        0.0131845093222756, 
     283        0.0128438426753249, 
     284        0.0124975796646449, 
     285        0.0121458711652067, 
     286        0.0117888704247183, 
     287        0.0114267329968529, 
     288        0.0110596166734735, 
     289        0.0106876814158841, 
     290        0.0103110892851360, 
     291        0.0099300043714212, 
     292        0.0095445927225849, 
     293        0.0091550222717888, 
     294        0.0087614627643580, 
     295        0.0083640856838475, 
     296        0.0079630641773633, 
     297        0.0075585729801782, 
     298        0.0071507883396855, 
     299        0.0067398879387430, 
     300        0.0063260508184704, 
     301        0.0059094573005900, 
     302        0.0054902889094487, 
     303        0.0050687282939456, 
     304        0.0046449591497966, 
     305        0.0042191661429919, 
     306        0.0037915348363451, 
     307        0.0033622516236779, 
     308        0.0029315036836558, 
     309        0.0024994789888943, 
     310        0.0020663664924131, 
     311        0.0016323569986067, 
     312        0.0011976474864367, 
     313        0.0007624720924706, 
     314        0.0003276086705538 
    315315}; 

  • sasmodels/models/lib/gauss20.c

    r14de349 r994d77f  
    99 
    1010// Gaussians 
    11 constant real Gauss20Wt[20]={ 
    12         REAL(.0176140071391521), 
    13         REAL(.0406014298003869), 
    14         REAL(.0626720483341091), 
    15         REAL(.0832767415767047), 
    16         REAL(.10193011981724), 
    17         REAL(.118194531961518), 
    18         REAL(.131688638449177), 
    19         REAL(.142096109318382), 
    20         REAL(.149172986472604), 
    21         REAL(.152753387130726), 
    22         REAL(.152753387130726), 
    23         REAL(.149172986472604), 
    24         REAL(.142096109318382), 
    25         REAL(.131688638449177), 
    26         REAL(.118194531961518), 
    27         REAL(.10193011981724), 
    28         REAL(.0832767415767047), 
    29         REAL(.0626720483341091), 
    30         REAL(.0406014298003869), 
    31         REAL(.0176140071391521) 
     11constant double Gauss20Wt[20]={ 
     12        .0176140071391521, 
     13        .0406014298003869, 
     14        .0626720483341091, 
     15        .0832767415767047, 
     16        .10193011981724, 
     17        .118194531961518, 
     18        .131688638449177, 
     19        .142096109318382, 
     20        .149172986472604, 
     21        .152753387130726, 
     22        .152753387130726, 
     23        .149172986472604, 
     24        .142096109318382, 
     25        .131688638449177, 
     26        .118194531961518, 
     27        .10193011981724, 
     28        .0832767415767047, 
     29        .0626720483341091, 
     30        .0406014298003869, 
     31        .0176140071391521 
    3232}; 
    3333 
    34 constant real Gauss20Z[20]={ 
    35         REAL(-.993128599185095), 
    36         REAL(-.963971927277914), 
    37         REAL(-.912234428251326), 
    38         REAL(-.839116971822219), 
    39         REAL(-.746331906460151), 
    40         REAL(-.636053680726515), 
    41         REAL(-.510867001950827), 
    42         REAL(-.37370608871542), 
    43         REAL(-.227785851141645), 
    44         REAL(-.076526521133497), 
    45         REAL(.0765265211334973), 
    46         REAL(.227785851141645), 
    47         REAL(.37370608871542), 
    48         REAL(.510867001950827), 
    49         REAL(.636053680726515), 
    50         REAL(.746331906460151), 
    51         REAL(.839116971822219), 
    52         REAL(.912234428251326), 
    53         REAL(.963971927277914), 
    54         REAL(.993128599185095) 
     34constant double Gauss20Z[20]={ 
     35        -.993128599185095, 
     36        -.963971927277914, 
     37        -.912234428251326, 
     38        -.839116971822219, 
     39        -.746331906460151, 
     40        -.636053680726515, 
     41        -.510867001950827, 
     42        -.37370608871542, 
     43        -.227785851141645, 
     44        -.076526521133497, 
     45        .0765265211334973, 
     46        .227785851141645, 
     47        .37370608871542, 
     48        .510867001950827, 
     49        .636053680726515, 
     50        .746331906460151, 
     51        .839116971822219, 
     52        .912234428251326, 
     53        .963971927277914, 
     54        .993128599185095 
    5555}; 

  • sasmodels/models/lib/gauss76.c

    r14de349 r994d77f  
    99 
    1010// Gaussians 
    11 constant real Gauss76Wt[76]={ 
    12         REAL(.00126779163408536),               //0 
    13         REAL(.00294910295364247), 
    14         REAL(.00462793522803742), 
    15         REAL(.00629918049732845), 
    16         REAL(.00795984747723973), 
    17         REAL(.00960710541471375), 
    18         REAL(.0112381685696677), 
    19         REAL(.0128502838475101), 
    20         REAL(.0144407317482767), 
    21         REAL(.0160068299122486), 
    22         REAL(.0175459372914742),                //10 
    23         REAL(.0190554584671906), 
    24         REAL(.020532847967908), 
    25         REAL(.0219756145344162), 
    26         REAL(.0233813253070112), 
    27         REAL(.0247476099206597), 
    28         REAL(.026072164497986), 
    29         REAL(.0273527555318275), 
    30         REAL(.028587223650054), 
    31         REAL(.029773487255905), 
    32         REAL(.0309095460374916),                //20 
    33         REAL(.0319934843404216), 
    34         REAL(.0330234743977917), 
    35         REAL(.0339977794120564), 
    36         REAL(.0349147564835508), 
    37         REAL(.0357728593807139), 
    38         REAL(.0365706411473296), 
    39         REAL(.0373067565423816), 
    40         REAL(.0379799643084053), 
    41         REAL(.0385891292645067), 
    42         REAL(.0391332242205184),                //30 
    43         REAL(.0396113317090621), 
    44         REAL(.0400226455325968), 
    45         REAL(.040366472122844), 
    46         REAL(.0406422317102947), 
    47         REAL(.0408494593018285), 
    48         REAL(.040987805464794), 
    49         REAL(.0410570369162294), 
    50         REAL(.0410570369162294), 
    51         REAL(.040987805464794), 
    52         REAL(.0408494593018285),                //40 
    53         REAL(.0406422317102947), 
    54         REAL(.040366472122844), 
    55         REAL(.0400226455325968), 
    56         REAL(.0396113317090621), 
    57         REAL(.0391332242205184), 
    58         REAL(.0385891292645067), 
    59         REAL(.0379799643084053), 
    60         REAL(.0373067565423816), 
    61         REAL(.0365706411473296), 
    62         REAL(.0357728593807139),                //50 
    63         REAL(.0349147564835508), 
    64         REAL(.0339977794120564), 
    65         REAL(.0330234743977917), 
    66         REAL(.0319934843404216), 
    67         REAL(.0309095460374916), 
    68         REAL(.029773487255905), 
    69         REAL(.028587223650054), 
    70         REAL(.0273527555318275), 
    71         REAL(.026072164497986), 
    72         REAL(.0247476099206597),                //60 
    73         REAL(.0233813253070112), 
    74         REAL(.0219756145344162), 
    75         REAL(.020532847967908), 
    76         REAL(.0190554584671906), 
    77         REAL(.0175459372914742), 
    78         REAL(.0160068299122486), 
    79         REAL(.0144407317482767), 
    80         REAL(.0128502838475101), 
    81         REAL(.0112381685696677), 
    82         REAL(.00960710541471375),               //70 
    83         REAL(.00795984747723973), 
    84         REAL(.00629918049732845), 
    85         REAL(.00462793522803742), 
    86         REAL(.00294910295364247), 
    87         REAL(.00126779163408536)                //75 (indexed from 0) 
     11constant double Gauss76Wt[76]={ 
     12        .00126779163408536,             //0 
     13        .00294910295364247, 
     14        .00462793522803742, 
     15        .00629918049732845, 
     16        .00795984747723973, 
     17        .00960710541471375, 
     18        .0112381685696677, 
     19        .0128502838475101, 
     20        .0144407317482767, 
     21        .0160068299122486, 
     22        .0175459372914742,              //10 
     23        .0190554584671906, 
     24        .020532847967908, 
     25        .0219756145344162, 
     26        .0233813253070112, 
     27        .0247476099206597, 
     28        .026072164497986, 
     29        .0273527555318275, 
     30        .028587223650054, 
     31        .029773487255905, 
     32        .0309095460374916,              //20 
     33        .0319934843404216, 
     34        .0330234743977917, 
     35        .0339977794120564, 
     36        .0349147564835508, 
     37        .0357728593807139, 
     38        .0365706411473296, 
     39        .0373067565423816, 
     40        .0379799643084053, 
     41        .0385891292645067, 
     42        .0391332242205184,              //30 
     43        .0396113317090621, 
     44        .0400226455325968, 
     45        .040366472122844, 
     46        .0406422317102947, 
     47        .0408494593018285, 
     48        .040987805464794, 
     49        .0410570369162294, 
     50        .0410570369162294, 
     51        .040987805464794, 
     52        .0408494593018285,              //40 
     53        .0406422317102947, 
     54        .040366472122844, 
     55        .0400226455325968, 
     56        .0396113317090621, 
     57        .0391332242205184, 
     58        .0385891292645067, 
     59        .0379799643084053, 
     60        .0373067565423816, 
     61        .0365706411473296, 
     62        .0357728593807139,              //50 
     63        .0349147564835508, 
     64        .0339977794120564, 
     65        .0330234743977917, 
     66        .0319934843404216, 
     67        .0309095460374916, 
     68        .029773487255905, 
     69        .028587223650054, 
     70        .0273527555318275, 
     71        .026072164497986, 
     72        .0247476099206597,              //60 
     73        .0233813253070112, 
     74        .0219756145344162, 
     75        .020532847967908, 
     76        .0190554584671906, 
     77        .0175459372914742, 
     78        .0160068299122486, 
     79        .0144407317482767, 
     80        .0128502838475101, 
     81        .0112381685696677, 
     82        .00960710541471375,             //70 
     83        .00795984747723973, 
     84        .00629918049732845, 
     85        .00462793522803742, 
     86        .00294910295364247, 
     87        .00126779163408536              //75 (indexed from 0) 
    8888}; 
    8989 
    90 constant real Gauss76Z[76]={ 
    91         REAL(-.999505948362153),                //0 
    92         REAL(-.997397786355355), 
    93         REAL(-.993608772723527), 
    94         REAL(-.988144453359837), 
    95         REAL(-.981013938975656), 
    96         REAL(-.972229228520377), 
    97         REAL(-.961805126758768), 
    98         REAL(-.949759207710896), 
    99         REAL(-.936111781934811), 
    100         REAL(-.92088586125215), 
    101         REAL(-.904107119545567),                //10 
    102         REAL(-.885803849292083), 
    103         REAL(-.866006913771982), 
    104         REAL(-.844749694983342), 
    105         REAL(-.822068037328975), 
    106         REAL(-.7980001871612), 
    107         REAL(-.77258672828181), 
    108         REAL(-.74587051350361), 
    109         REAL(-.717896592387704), 
    110         REAL(-.688712135277641), 
    111         REAL(-.658366353758143),                //20 
    112         REAL(-.626910417672267), 
    113         REAL(-.594397368836793), 
    114         REAL(-.560882031601237), 
    115         REAL(-.526420920401243), 
    116         REAL(-.491072144462194), 
    117         REAL(-.454895309813726), 
    118         REAL(-.417951418780327), 
    119         REAL(-.380302767117504), 
    120         REAL(-.342012838966962), 
    121         REAL(-.303146199807908),                //30 
    122         REAL(-.263768387584994), 
    123         REAL(-.223945802196474), 
    124         REAL(-.183745593528914), 
    125         REAL(-.143235548227268), 
    126         REAL(-.102483975391227), 
    127         REAL(-.0615595913906112), 
    128         REAL(-.0205314039939986), 
    129         REAL(.0205314039939986), 
    130         REAL(.0615595913906112), 
    131         REAL(.102483975391227),                 //40 
    132         REAL(.143235548227268), 
    133         REAL(.183745593528914), 
    134         REAL(.223945802196474), 
    135         REAL(.263768387584994), 
    136         REAL(.303146199807908), 
    137         REAL(.342012838966962), 
    138         REAL(.380302767117504), 
    139         REAL(.417951418780327), 
    140         REAL(.454895309813726), 
    141         REAL(.491072144462194),         //50 
    142         REAL(.526420920401243), 
    143         REAL(.560882031601237), 
    144         REAL(.594397368836793), 
    145         REAL(.626910417672267), 
    146         REAL(.658366353758143), 
    147         REAL(.688712135277641), 
    148         REAL(.717896592387704), 
    149         REAL(.74587051350361), 
    150         REAL(.77258672828181), 
    151         REAL(.7980001871612),   //60 
    152         REAL(.822068037328975), 
    153         REAL(.844749694983342), 
    154         REAL(.866006913771982), 
    155         REAL(.885803849292083), 
    156         REAL(.904107119545567), 
    157         REAL(.92088586125215), 
    158         REAL(.936111781934811), 
    159         REAL(.949759207710896), 
    160         REAL(.961805126758768), 
    161         REAL(.972229228520377),         //70 
    162         REAL(.981013938975656), 
    163         REAL(.988144453359837), 
    164         REAL(.993608772723527), 
    165         REAL(.997397786355355), 
    166         REAL(.999505948362153)          //75 
     90constant double Gauss76Z[76]={ 
     91        -.999505948362153,              //0 
     92        -.997397786355355, 
     93        -.993608772723527, 
     94        -.988144453359837, 
     95        -.981013938975656, 
     96        -.972229228520377, 
     97        -.961805126758768, 
     98        -.949759207710896, 
     99        -.936111781934811, 
     100        -.92088586125215, 
     101        -.904107119545567,              //10 
     102        -.885803849292083, 
     103        -.866006913771982, 
     104        -.844749694983342, 
     105        -.822068037328975, 
     106        -.7980001871612, 
     107        -.77258672828181, 
     108        -.74587051350361, 
     109        -.717896592387704, 
     110        -.688712135277641, 
     111        -.658366353758143,              //20 
     112        -.626910417672267, 
     113        -.594397368836793, 
     114        -.560882031601237, 
     115        -.526420920401243, 
     116        -.491072144462194, 
     117        -.454895309813726, 
     118        -.417951418780327, 
     119        -.380302767117504, 
     120        -.342012838966962, 
     121        -.303146199807908,              //30 
     122        -.263768387584994, 
     123        -.223945802196474, 
     124        -.183745593528914, 
     125        -.143235548227268, 
     126        -.102483975391227, 
     127        -.0615595913906112, 
     128        -.0205314039939986, 
     129        .0205314039939986, 
     130        .0615595913906112, 
     131        .102483975391227,                       //40 
     132        .143235548227268, 
     133        .183745593528914, 
     134        .223945802196474, 
     135        .263768387584994, 
     136        .303146199807908, 
     137        .342012838966962, 
     138        .380302767117504, 
     139        .417951418780327, 
     140        .454895309813726, 
     141        .491072144462194,               //50 
     142        .526420920401243, 
     143        .560882031601237, 
     144        .594397368836793, 
     145        .626910417672267, 
     146        .658366353758143, 
     147        .688712135277641, 
     148        .717896592387704, 
     149        .74587051350361, 
     150        .77258672828181, 
     151        .7980001871612, //60 
     152        .822068037328975, 
     153        .844749694983342, 
     154        .866006913771982, 
     155        .885803849292083, 
     156        .904107119545567, 
     157        .92088586125215, 
     158        .936111781934811, 
     159        .949759207710896, 
     160        .961805126758768, 
     161        .972229228520377,               //70 
     162        .981013938975656, 
     163        .988144453359837, 
     164        .993608772723527, 
     165        .997397786355355, 
     166        .999505948362153                //75 
    167167}; 
Note: See TracChangeset for help on using the changeset viewer.