Changeset 95ce773 in sasmodels for sasmodels/models/lib/polevl.c


Ignore:
Timestamp:
Mar 18, 2016 12:45:38 PM (8 years ago)
Author:
wojciech
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:
9f7a852
Parents:
a629d8e
Message:

Bessel functions clean-up

File:
1 edited

Legend:

Unmodified
Added
Removed

  • sasmodels/models/lib/polevl.c

    rfad5dc1 r95ce773  
    5151*/ 
    5252 
     53constant double RPJ1[8] = { 
     54    -8.99971225705559398224E8, 
     55    4.52228297998194034323E11, 
     56    -7.27494245221818276015E13, 
     57    3.68295732863852883286E15, 
     58    0.0, 
     59    0.0, 
     60    0.0, 
     61    0.0 }; 
     62 
     63constant double RQJ1[8] = { 
     64    6.20836478118054335476E2, 
     65    2.56987256757748830383E5, 
     66    8.35146791431949253037E7, 
     67    2.21511595479792499675E10, 
     68    4.74914122079991414898E12, 
     69    7.84369607876235854894E14, 
     70    8.95222336184627338078E16, 
     71    5.32278620332680085395E18 
     72    }; 
     73 
     74constant double PPJ1[8] = { 
     75    7.62125616208173112003E-4, 
     76    7.31397056940917570436E-2, 
     77    1.12719608129684925192E0, 
     78    5.11207951146807644818E0, 
     79    8.42404590141772420927E0, 
     80    5.21451598682361504063E0, 
     81    1.00000000000000000254E0, 
     82    0.0} ; 
     83 
     84 
     85constant double PQJ1[8] = { 
     86    5.71323128072548699714E-4, 
     87    6.88455908754495404082E-2, 
     88    1.10514232634061696926E0, 
     89    5.07386386128601488557E0, 
     90    8.39985554327604159757E0, 
     91    5.20982848682361821619E0, 
     92    9.99999999999999997461E-1, 
     93    0.0 }; 
     94 
     95constant double QPJ1[8] = { 
     96    5.10862594750176621635E-2, 
     97    4.98213872951233449420E0, 
     98    7.58238284132545283818E1, 
     99    3.66779609360150777800E2, 
     100    7.10856304998926107277E2, 
     101    5.97489612400613639965E2, 
     102    2.11688757100572135698E2, 
     103    2.52070205858023719784E1 }; 
     104 
     105constant double QQJ1[8] = { 
     106    7.42373277035675149943E1, 
     107    1.05644886038262816351E3, 
     108    4.98641058337653607651E3, 
     109    9.56231892404756170795E3, 
     110    7.99704160447350683650E3, 
     111    2.82619278517639096600E3, 
     112    3.36093607810698293419E2, 
     113    0.0 }; 
     114 
     115constant double JPJ1[8] = { 
     116    -4.878788132172128E-009, 
     117    6.009061827883699E-007, 
     118    -4.541343896997497E-005, 
     119    1.937383947804541E-003, 
     120    -3.405537384615824E-002, 
     121    0.0, 
     122    0.0, 
     123    0.0 
     124    }; 
     125 
     126constant double MO1J1[8] = { 
     127    6.913942741265801E-002, 
     128    -2.284801500053359E-001, 
     129    3.138238455499697E-001, 
     130    -2.102302420403875E-001, 
     131    5.435364690523026E-003, 
     132    1.493389585089498E-001, 
     133    4.976029650847191E-006, 
     134    7.978845453073848E-001 
     135    }; 
     136 
     137constant double PH1J1[8] = { 
     138    -4.497014141919556E+001, 
     139    5.073465654089319E+001, 
     140    -2.485774108720340E+001, 
     141    7.222973196770240E+000, 
     142    -1.544842782180211E+000, 
     143    3.503787691653334E-001, 
     144    -1.637986776941202E-001, 
     145    3.749989509080821E-001 
     146    }; 
     147 
     148 constant double PPJ0[8] = { 
     149        7.96936729297347051624E-4, 
     150        8.28352392107440799803E-2, 
     151        1.23953371646414299388E0, 
     152        5.44725003058768775090E0, 
     153        8.74716500199817011941E0, 
     154        5.30324038235394892183E0, 
     155        9.99999999999999997821E-1, 
     156        0.0 
     157    }; 
     158 
     159 constant double PQJ0[8] = { 
     160        9.24408810558863637013E-4, 
     161        8.56288474354474431428E-2, 
     162        1.25352743901058953537E0, 
     163        5.47097740330417105182E0, 
     164        8.76190883237069594232E0, 
     165        5.30605288235394617618E0, 
     166        1.00000000000000000218E0, 
     167        0.0 
     168    }; 
     169 
     170 constant double QPJ0[8] = { 
     171        -1.13663838898469149931E-2, 
     172        -1.28252718670509318512E0, 
     173        -1.95539544257735972385E1, 
     174        -9.32060152123768231369E1, 
     175        -1.77681167980488050595E2, 
     176        -1.47077505154951170175E2, 
     177        -5.14105326766599330220E1, 
     178        -6.05014350600728481186E0, 
     179    }; 
     180 
     181 constant double QQJ0[8] = { 
     182        /*  1.00000000000000000000E0,*/ 
     183        6.43178256118178023184E1, 
     184        8.56430025976980587198E2, 
     185        3.88240183605401609683E3, 
     186        7.24046774195652478189E3, 
     187        5.93072701187316984827E3, 
     188        2.06209331660327847417E3, 
     189        2.42005740240291393179E2, 
     190    }; 
     191 
     192 constant double YPJ0[8] = { 
     193        1.55924367855235737965E4, 
     194        -1.46639295903971606143E7, 
     195        5.43526477051876500413E9, 
     196        -9.82136065717911466409E11, 
     197        8.75906394395366999549E13, 
     198        -3.46628303384729719441E15, 
     199        4.42733268572569800351E16, 
     200        -1.84950800436986690637E16, 
     201 }; 
     202 
     203 
     204 constant double YQJ0[7] = { 
     205        /* 1.00000000000000000000E0,*/ 
     206        1.04128353664259848412E3, 
     207        6.26107330137134956842E5, 
     208        2.68919633393814121987E8, 
     209        8.64002487103935000337E10, 
     210        2.02979612750105546709E13, 
     211        3.17157752842975028269E15, 
     212        2.50596256172653059228E17, 
     213  }; 
     214 
     215 constant double RPJ0[8] = { 
     216        -4.79443220978201773821E9, 
     217        1.95617491946556577543E12, 
     218        -2.49248344360967716204E14, 
     219        9.70862251047306323952E15, 
     220        0.0, 
     221        0.0, 
     222        0.0, 
     223        0.0 
     224  }; 
     225 
     226 constant double RQJ0[8] = { 
     227        /* 1.00000000000000000000E0,*/ 
     228        4.99563147152651017219E2, 
     229        1.73785401676374683123E5, 
     230        4.84409658339962045305E7, 
     231        1.11855537045356834862E10, 
     232        2.11277520115489217587E12, 
     233        3.10518229857422583814E14, 
     234        3.18121955943204943306E16, 
     235        1.71086294081043136091E18, 
     236  }; 
     237 
     238 constant double MOJ0[8] = { 
     239        -6.838999669318810E-002, 
     240        1.864949361379502E-001, 
     241        -2.145007480346739E-001, 
     242        1.197549369473540E-001, 
     243        -3.560281861530129E-003, 
     244        -4.969382655296620E-002, 
     245        -3.355424622293709E-006, 
     246        7.978845717621440E-001 
     247  }; 
     248 
     249 constant double PHJ0[8] = { 
     250        3.242077816988247E+001, 
     251        -3.630592630518434E+001, 
     252        1.756221482109099E+001, 
     253        -4.974978466280903E+000, 
     254        1.001973420681837E+000, 
     255        -1.939906941791308E-001, 
     256        6.490598792654666E-002, 
     257        -1.249992184872738E-001 
     258  }; 
     259 
     260 constant double JPJ0[8] = { 
     261        -6.068350350393235E-008, 
     262        6.388945720783375E-006, 
     263        -3.969646342510940E-004, 
     264        1.332913422519003E-002, 
     265        -1.729150680240724E-001, 
     266        0.0, 
     267        0.0, 
     268        0.0 
     269 }; 
    53270 
    54271double polevl( double x, constant double *coef, int N ) { 
Note: See TracChangeset for help on using the changeset viewer.