source: sasmodels/Kernel/Capcyl_Kfun.cpp @ 099e053

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

further organizing

  • Property mode set to 100644
File size: 3.7 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}
16real Gauss76Z(int i)
17{
18real array[76] = {.999505948362153*(-1.0),.997397786355355*(-1.0),.993608772723527*(-1.0),.988144453359837*(-1.0),
19.981013938975656*(-1.0),.972229228520377*(-1.0),.961805126758768*(-1.0),.949759207710896*(-1.0),
20.936111781934811*(-1.0),.92088586125215*(-1.0),.904107119545567*(-1.0),.885803849292083*(-1.0),
21.866006913771982*(-1.0),.844749694983342*(-1.0),.822068037328975*(-1.0),.7980001871612*(-1.0),
22.77258672828181*(-1.0),.74587051350361*(-1.0),.717896592387704*(-1.0),.688712135277641*(-1.0),
23.658366353758143*(-1.0),.626910417672267*(-1.0),.594397368836793*(-1.0),.560882031601237*(-1.0),
24.526420920401243*(-1.0),.491072144462194*(-1.0),.454895309813726*(-1.0),.417951418780327*(-1.0),
25.380302767117504*(-1.0),.342012838966962*(-1.0),.303146199807908*(-1.0),.263768387584994*(-1.0),
26.223945802196474*(-1.0),.183745593528914*(-1.0),.143235548227268*(-1.0),.102483975391227*(-1.0),
27.0615595913906112*(-1.0),.0205314039939986*(-1.0),-.0205314039939986*(-1.0),-.0615595913906112*(-1.0),
28-.102483975391227*(-1.0),-.143235548227268*(-1.0),-.183745593528914*(-1.0),-.223945802196474*(-1.0),
29-.263768387584994*(-1.0),-.303146199807908*(-1.0),-.342012838966962*(-1.0),-.380302767117504*(-1.0),
30-.417951418780327*(-1.0),-.454895309813726*(-1.0),-.491072144462194*(-1.0),-.526420920401243*(-1.0),
31-.560882031601237*(-1.0),-.594397368836793*(-1.0),-.626910417672267*(-1.0),-.658366353758143*(-1.0),
32-.688712135277641*(-1.0),-.717896592387704*(-1.0),-.74587051350361*(-1.0),-.77258672828181*(-1.0),
33-.7980001871612*(-1.0), -.822068037328975*(-1.0),-.844749694983342*(-1.0),-.866006913771982*(-1.0),
34-.885803849292083*(-1.0),-.904107119545567*(-1.0),-.92088586125215*(-1.0),-.936111781934811*(-1.0),
35-.949759207710896*(-1.0),-.961805126758768*(-1.0),-.972229228520377*(-1.0),-.981013938975656*(-1.0),
36-.988144453359837*(-1.0),-.993608772723527*(-1.0),-.997397786355355*(-1.0),-.999505948362153*(-1.0) };
37
38return array[i];
39}
40real Gauss76Wt(int i)
41{
42real array[76] = {.00126779163408536,.00294910295364247,.00462793522803742,.00629918049732845,
43.00795984747723973,.00960710541471375,.0112381685696677,.0128502838475101,
44.0144407317482767,.0160068299122486,.0175459372914742,.0190554584671906,
45.020532847967908,.0219756145344162,.0233813253070112,.0247476099206597,
46.026072164497986,.0273527555318275,.028587223650054,.029773487255905,
47.0309095460374916,.0319934843404216,.0330234743977917,.0339977794120564,
48.0349147564835508,.0357728593807139,.0365706411473296,.0373067565423816,
49.0379799643084053,.0385891292645067,.0391332242205184,.0396113317090621,
50.0400226455325968,.040366472122844,.0406422317102947,.0408494593018285,
51.040987805464794,.0410570369162294,.0410570369162294,.040987805464794,
52.0408494593018285,.0406422317102947,.040366472122844,.0400226455325968,
53.0396113317090621,.0391332242205184,.0385891292645067,.0379799643084053,
54.0373067565423816,.0365706411473296,.0357728593807139,.0349147564835508,
55.0339977794120564,.0330234743977917,.0319934843404216,.0309095460374916,
56.029773487255905,.028587223650054,.0273527555318275,.026072164497986,
57.0247476099206597,.0233813253070112,.0219756145344162,.020532847967908,
58.0190554584671906,.0175459372914742,.0160068299122486,.0144407317482767,
59.0128502838475101,.0112381685696677,.00960710541471375,.00795984747723973,
60.00629918049732845,.00462793522803742,.00294910295364247,.00126779163408536 };
61
62return array[i];
63}
Note: See TracBrowser for help on using the repository browser.