source: sasmodels/Capcyl_Kfun.cpp @ 8faffcd

core_shell_microgelscostrafo411magnetic_modelrelease_v0.94release_v0.95ticket-1257-vesicle-productticket_1156ticket_1265_superballticket_822_more_unit_tests
Last change on this file since 8faffcd was 8faffcd, checked in by HMP1 <helen.park@…>, 10 years ago

Update for Aaron
  • Property mode set to 100644
File size: 4.0 KB
Line 
1real ConvLens_kernel(real len, real rad, real endRad, real x, real tt, real theta)
2{
3    real be = 0;
4        real hDist = -1.0*sqrt(fabs(endRad*endRad-rad*rad));
5        real arg1 = x*cos(theta)*(endRad*tt+hDist+len/2.0);
6        real arg2 = x*endRad*sin(theta)*sqrt(1.0-tt*tt);
7
8        if(arg2 == 0) {be = 0.5;}
9        else {
10                be = NR_BessJ1(arg2)/arg2;
11        }
12        real val = cos(arg1)*(1.0-tt*tt)*be;
13
14        return(val);
15}
16const real Gauss76Z[76]={
17         .999505948362153*(-1.0),               //0
18         .997397786355355*(-1.0),
19         .993608772723527*(-1.0),
20         .988144453359837*(-1.0),
21         .981013938975656*(-1.0),
22         .972229228520377*(-1.0),
23         .961805126758768*(-1.0),
24         .949759207710896*(-1.0),
25         .936111781934811*(-1.0),
26         .92088586125215*(-1.0),
27         .904107119545567*(-1.0),               //10
28         .885803849292083*(-1.0),
29         .866006913771982*(-1.0),
30         .844749694983342*(-1.0),
31         .822068037328975*(-1.0),
32         .7980001871612*(-1.0),
33         .77258672828181*(-1.0),
34         .74587051350361*(-1.0),
35         .717896592387704*(-1.0),
36         .688712135277641*(-1.0),
37         .658366353758143*(-1.0),               //20
38         .626910417672267*(-1.0),
39         .594397368836793*(-1.0),
40         .560882031601237*(-1.0),
41         .526420920401243*(-1.0),
42         .491072144462194*(-1.0),
43         .454895309813726*(-1.0),
44         .417951418780327*(-1.0),
45         .380302767117504*(-1.0),
46         .342012838966962*(-1.0),
47         .303146199807908*(-1.0),               //30
48         .263768387584994*(-1.0),
49         .223945802196474*(-1.0),
50         .183745593528914*(-1.0),
51         .143235548227268*(-1.0),
52         .102483975391227*(-1.0),
53         .0615595913906112*(-1.0),
54         .0205314039939986*(-1.0),
55        -.0205314039939986*(-1.0),
56        -.0615595913906112*(-1.0),
57        -.102483975391227*(-1.0),                       //40
58        -.143235548227268*(-1.0),
59        -.183745593528914*(-1.0),
60        -.223945802196474*(-1.0),
61        -.263768387584994*(-1.0),
62        -.303146199807908*(-1.0),
63        -.342012838966962*(-1.0),
64        -.380302767117504*(-1.0),
65        -.417951418780327*(-1.0),
66        -.454895309813726*(-1.0),
67        -.491072144462194*(-1.0),               //50
68        -.526420920401243*(-1.0),
69        -.560882031601237*(-1.0),
70        -.594397368836793*(-1.0),
71        -.626910417672267*(-1.0),
72        -.658366353758143*(-1.0),
73        -.688712135277641*(-1.0),
74        -.717896592387704*(-1.0),
75        -.74587051350361*(-1.0),
76        -.77258672828181*(-1.0),
77        -.7980001871612*(-1.0), //60
78        -.822068037328975*(-1.0),
79        -.844749694983342*(-1.0),
80        -.866006913771982*(-1.0),
81        -.885803849292083*(-1.0),
82        -.904107119545567*(-1.0),
83        -.92088586125215*(-1.0),
84        -.936111781934811*(-1.0),
85        -.949759207710896*(-1.0),
86        -.961805126758768*(-1.0),
87        -.972229228520377*(-1.0),               //70
88        -.981013938975656*(-1.0),
89        -.988144453359837*(-1.0),
90        -.993608772723527*(-1.0),
91        -.997397786355355*(-1.0),
92        -.999505948362153*(-1.0)                //75
93};
94const real Gauss76Wt[76]={
95        .00126779163408536,             //0
96        .00294910295364247,
97        .00462793522803742,
98        .00629918049732845,
99        .00795984747723973,
100        .00960710541471375,
101        .0112381685696677,
102        .0128502838475101,
103        .0144407317482767,
104        .0160068299122486,
105        .0175459372914742,              //10
106        .0190554584671906,
107        .020532847967908,
108        .0219756145344162,
109        .0233813253070112,
110        .0247476099206597,
111        .026072164497986,
112        .0273527555318275,
113        .028587223650054,
114        .029773487255905,
115        .0309095460374916,              //20
116        .0319934843404216,
117        .0330234743977917,
118        .0339977794120564,
119        .0349147564835508,
120        .0357728593807139,
121        .0365706411473296,
122        .0373067565423816,
123        .0379799643084053,
124        .0385891292645067,
125        .0391332242205184,              //30
126        .0396113317090621,
127        .0400226455325968,
128        .040366472122844,
129        .0406422317102947,
130        .0408494593018285,
131        .040987805464794,
132        .0410570369162294,
133        .0410570369162294,
134        .040987805464794,
135        .0408494593018285,              //40
136        .0406422317102947,
137        .040366472122844,
138        .0400226455325968,
139        .0396113317090621,
140        .0391332242205184,
141        .0385891292645067,
142        .0379799643084053,
143        .0373067565423816,
144        .0365706411473296,
145        .0357728593807139,              //50
146        .0349147564835508,
147        .0339977794120564,
148        .0330234743977917,
149        .0319934843404216,
150        .0309095460374916,
151        .029773487255905,
152        .028587223650054,
153        .0273527555318275,
154        .026072164497986,
155        .0247476099206597,              //60
156        .0233813253070112,
157        .0219756145344162,
158        .020532847967908,
159        .0190554584671906,
160        .0175459372914742,
161        .0160068299122486,
162        .0144407317482767,
163        .0128502838475101,
164        .0112381685696677,
165        .00960710541471375,             //70
166        .00795984747723973,
167        .00629918049732845,
168        .00462793522803742,
169        .00294910295364247,
170        .00126779163408536              //75 (indexed from 0)
171};
Note: See TracBrowser for help on using the repository browser.