Changeset 74e9b5f in sasmodels for sasmodels/models/lib/sas_J1.c
- Timestamp:
- Oct 12, 2018 8:52:48 PM (6 years ago)
- Branches:
- master, core_shell_microgels, magnetic_model, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests
- Children:
- 4de14584
- Parents:
- b0de252
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sasmodels/models/lib/sas_J1.c
r0db7dbd r74e9b5f 43 43 //Cephes double pression function 44 44 45 constant _vardouble RPJ1[8] = {45 constant double RPJ1[8] = { 46 46 -8.99971225705559398224E8, 47 47 4.52228297998194034323E11, … … 53 53 0.0 }; 54 54 55 constant _vardouble RQJ1[8] = {55 constant double RQJ1[8] = { 56 56 6.20836478118054335476E2, 57 57 2.56987256757748830383E5, … … 64 64 }; 65 65 66 constant _vardouble PPJ1[8] = {66 constant double PPJ1[8] = { 67 67 7.62125616208173112003E-4, 68 68 7.31397056940917570436E-2, … … 75 75 76 76 77 constant _vardouble PQJ1[8] = {77 constant double PQJ1[8] = { 78 78 5.71323128072548699714E-4, 79 79 6.88455908754495404082E-2, … … 85 85 0.0 }; 86 86 87 constant _vardouble QPJ1[8] = {87 constant double QPJ1[8] = { 88 88 5.10862594750176621635E-2, 89 89 4.98213872951233449420E0, … … 95 95 2.52070205858023719784E1 }; 96 96 97 constant _vardouble QQJ1[8] = {97 constant double QQJ1[8] = { 98 98 7.42373277035675149943E1, 99 99 1.05644886038262816351E3, … … 105 105 0.0 }; 106 106 107 __device__static107 static 108 108 double cephes_j1(double x) 109 109 { … … 155 155 #else 156 156 //Single precission version of cephes 157 constant _varfloat JPJ1[8] = {157 constant float JPJ1[8] = { 158 158 -4.878788132172128E-009, 159 159 6.009061827883699E-007, … … 166 166 }; 167 167 168 constant _varfloat MO1J1[8] = {168 constant float MO1J1[8] = { 169 169 6.913942741265801E-002, 170 170 -2.284801500053359E-001, … … 177 177 }; 178 178 179 constant _varfloat PH1J1[8] = {179 constant float PH1J1[8] = { 180 180 -4.497014141919556E+001, 181 181 5.073465654089319E+001, … … 188 188 }; 189 189 190 __device__static190 static 191 191 float cephes_j1f(float xx) 192 192 { … … 239 239 240 240 //Finally J1c function that equals 2*J1(x)/x 241 __device__static241 static 242 242 double sas_2J1x_x(double x) 243 243 {
Note: See TracChangeset
for help on using the changeset viewer.