source: sans/XOP_Dev/SANSAnalysis/lib/GaussWeights.h @ 453

Last change on this file since 453 was 453, checked in by srkline, 14 years ago

Additions to the library of the 2008 model functions. Direct proting of the Igor code, duplicated by these XOPs

File size: 11.5 KB
Line 
1/*
2 *  GaussWeights.h
3 *  SANSAnalysis
4 *
5 *  Created by Andrew Jackson on 4/23/07.
6 *  Copyright 2007 __MyCompanyName__. All rights reserved.
7 *
8 */
9
10// Gaussians
11static 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
32};
33
34static 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
55};
56
57static double Gauss76Wt[76]={
58        .00126779163408536,             //0
59        .00294910295364247,
60        .00462793522803742,
61        .00629918049732845,
62        .00795984747723973,
63        .00960710541471375,
64        .0112381685696677,
65        .0128502838475101,
66        .0144407317482767,
67        .0160068299122486,
68        .0175459372914742,              //10
69        .0190554584671906,
70        .020532847967908,
71        .0219756145344162,
72        .0233813253070112,
73        .0247476099206597,
74        .026072164497986,
75        .0273527555318275,
76        .028587223650054,
77        .029773487255905,
78        .0309095460374916,              //20
79        .0319934843404216,
80        .0330234743977917,
81        .0339977794120564,
82        .0349147564835508,
83        .0357728593807139,
84        .0365706411473296,
85        .0373067565423816,
86        .0379799643084053,
87        .0385891292645067,
88        .0391332242205184,              //30
89        .0396113317090621,
90        .0400226455325968,
91        .040366472122844,
92        .0406422317102947,
93        .0408494593018285,
94        .040987805464794,
95        .0410570369162294,
96        .0410570369162294,
97        .040987805464794,
98        .0408494593018285,              //40
99        .0406422317102947,
100        .040366472122844,
101        .0400226455325968,
102        .0396113317090621,
103        .0391332242205184,
104        .0385891292645067,
105        .0379799643084053,
106        .0373067565423816,
107        .0365706411473296,
108        .0357728593807139,              //50
109        .0349147564835508,
110        .0339977794120564,
111        .0330234743977917,
112        .0319934843404216,
113        .0309095460374916,
114        .029773487255905,
115        .028587223650054,
116        .0273527555318275,
117        .026072164497986,
118        .0247476099206597,              //60
119        .0233813253070112,
120        .0219756145344162,
121        .020532847967908,
122        .0190554584671906,
123        .0175459372914742,
124        .0160068299122486,
125        .0144407317482767,
126        .0128502838475101,
127        .0112381685696677,
128        .00960710541471375,             //70
129        .00795984747723973,
130        .00629918049732845,
131        .00462793522803742,
132        .00294910295364247,
133        .00126779163408536              //75 (indexed from 0)
134};
135static double Gauss76Z[76]={
136         .999505948362153*(-1.0),               //0
137         .997397786355355*(-1.0),
138         .993608772723527*(-1.0),
139         .988144453359837*(-1.0),
140         .981013938975656*(-1.0),
141         .972229228520377*(-1.0),
142         .961805126758768*(-1.0),
143         .949759207710896*(-1.0),
144         .936111781934811*(-1.0),
145         .92088586125215*(-1.0),
146         .904107119545567*(-1.0),               //10
147         .885803849292083*(-1.0),
148         .866006913771982*(-1.0),
149         .844749694983342*(-1.0),
150         .822068037328975*(-1.0),
151         .7980001871612*(-1.0),
152         .77258672828181*(-1.0),
153         .74587051350361*(-1.0),
154         .717896592387704*(-1.0),
155         .688712135277641*(-1.0),
156         .658366353758143*(-1.0),               //20
157         .626910417672267*(-1.0),
158         .594397368836793*(-1.0),
159         .560882031601237*(-1.0),
160         .526420920401243*(-1.0),
161         .491072144462194*(-1.0),
162         .454895309813726*(-1.0),
163         .417951418780327*(-1.0),
164         .380302767117504*(-1.0),
165         .342012838966962*(-1.0),
166         .303146199807908*(-1.0),               //30
167         .263768387584994*(-1.0),
168         .223945802196474*(-1.0),
169         .183745593528914*(-1.0),
170         .143235548227268*(-1.0),
171         .102483975391227*(-1.0),
172         .0615595913906112*(-1.0),
173         .0205314039939986*(-1.0),
174        -.0205314039939986*(-1.0),
175        -.0615595913906112*(-1.0),
176        -.102483975391227*(-1.0),                       //40
177        -.143235548227268*(-1.0),
178        -.183745593528914*(-1.0),
179        -.223945802196474*(-1.0),
180        -.263768387584994*(-1.0),
181        -.303146199807908*(-1.0),
182        -.342012838966962*(-1.0),
183        -.380302767117504*(-1.0),
184        -.417951418780327*(-1.0),
185        -.454895309813726*(-1.0),
186        -.491072144462194*(-1.0),               //50
187        -.526420920401243*(-1.0),
188        -.560882031601237*(-1.0),
189        -.594397368836793*(-1.0),
190        -.626910417672267*(-1.0),
191        -.658366353758143*(-1.0),
192        -.688712135277641*(-1.0),
193        -.717896592387704*(-1.0),
194        -.74587051350361*(-1.0),
195        -.77258672828181*(-1.0),
196        -.7980001871612*(-1.0), //60
197        -.822068037328975*(-1.0),
198        -.844749694983342*(-1.0),
199        -.866006913771982*(-1.0),
200        -.885803849292083*(-1.0),
201        -.904107119545567*(-1.0),
202        -.92088586125215*(-1.0),
203        -.936111781934811*(-1.0),
204        -.949759207710896*(-1.0),
205        -.961805126758768*(-1.0),
206        -.972229228520377*(-1.0),               //70
207        -.981013938975656*(-1.0),
208        -.988144453359837*(-1.0),
209        -.993608772723527*(-1.0),
210        -.997397786355355*(-1.0),
211        -.999505948362153*(-1.0)                //75
212};
213
214static double Gauss150Z[150]={
215        -0.9998723404457334,
216        -0.9993274305065947,
217        -0.9983473449340834,
218        -0.9969322929775997,
219        -0.9950828645255290,
220        -0.9927998590434373,
221        -0.9900842691660192,
222        -0.9869372772712794,
223        -0.9833602541697529,
224        -0.9793547582425894,
225        -0.9749225346595943,
226        -0.9700655145738374,
227        -0.9647858142586956,
228        -0.9590857341746905,
229        -0.9529677579610971,
230        -0.9464345513503147,
231        -0.9394889610042837,
232        -0.9321340132728527,
233        -0.9243729128743136,
234        -0.9162090414984952,
235        -0.9076459563329236,
236        -0.8986873885126239,
237        -0.8893372414942055,
238        -0.8795995893549102,
239        -0.8694786750173527,
240        -0.8589789084007133,
241        -0.8481048644991847,
242        -0.8368612813885015,
243        -0.8252530581614230,
244        -0.8132852527930605,
245        -0.8009630799369827,
246        -0.7882919086530552,
247        -0.7752772600680049,
248        -0.7619248049697269,
249        -0.7482403613363824,
250        -0.7342298918013638,
251        -0.7198995010552305,
252        -0.7052554331857488,
253        -0.6903040689571928,
254        -0.6750519230300931,
255        -0.6595056411226444,
256        -0.6436719971150083,
257        -0.6275578900977726,
258        -0.6111703413658551,
259        -0.5945164913591590,
260        -0.5776035965513142,
261        -0.5604390262878617,
262        -0.5430302595752546,
263        -0.5253848818220803,
264        -0.5075105815339176,
265        -0.4894151469632753,
266        -0.4711064627160663,
267        -0.4525925063160997,
268        -0.4338813447290861,
269        -0.4149811308476706,
270        -0.3959000999390257,
271        -0.3766465660565522,
272        -0.3572289184172501,
273        -0.3376556177463400,
274        -0.3179351925907259,
275        -0.2980762356029071,
276        -0.2780873997969574,
277        -0.2579773947782034,
278        -0.2377549829482451,
279        -0.2174289756869712,
280        -0.1970082295132342,
281        -0.1765016422258567,
282        -0.1559181490266516,
283        -0.1352667186271445,
284        -0.1145563493406956,
285        -0.0937960651617229,
286        -0.0729949118337358,
287        -0.0521619529078925,
288        -0.0313062657937972,
289        -0.0104369378042598,
290        0.0104369378042598,
291        0.0313062657937972,
292        0.0521619529078925,
293        0.0729949118337358,
294        0.0937960651617229,
295        0.1145563493406956,
296        0.1352667186271445,
297        0.1559181490266516,
298        0.1765016422258567,
299        0.1970082295132342,
300        0.2174289756869712,
301        0.2377549829482451,
302        0.2579773947782034,
303        0.2780873997969574,
304        0.2980762356029071,
305        0.3179351925907259,
306        0.3376556177463400,
307        0.3572289184172501,
308        0.3766465660565522,
309        0.3959000999390257,
310        0.4149811308476706,
311        0.4338813447290861,
312        0.4525925063160997,
313        0.4711064627160663,
314        0.4894151469632753,
315        0.5075105815339176,
316        0.5253848818220803,
317        0.5430302595752546,
318        0.5604390262878617,
319        0.5776035965513142,
320        0.5945164913591590,
321        0.6111703413658551,
322        0.6275578900977726,
323        0.6436719971150083,
324        0.6595056411226444,
325        0.6750519230300931,
326        0.6903040689571928,
327        0.7052554331857488,
328        0.7198995010552305,
329        0.7342298918013638,
330        0.7482403613363824,
331        0.7619248049697269,
332        0.7752772600680049,
333        0.7882919086530552,
334        0.8009630799369827,
335        0.8132852527930605,
336        0.8252530581614230,
337        0.8368612813885015,
338        0.8481048644991847,
339        0.8589789084007133,
340        0.8694786750173527,
341        0.8795995893549102,
342        0.8893372414942055,
343        0.8986873885126239,
344        0.9076459563329236,
345        0.9162090414984952,
346        0.9243729128743136,
347        0.9321340132728527,
348        0.9394889610042837,
349        0.9464345513503147,
350        0.9529677579610971,
351        0.9590857341746905,
352        0.9647858142586956,
353        0.9700655145738374,
354        0.9749225346595943,
355        0.9793547582425894,
356        0.9833602541697529,
357        0.9869372772712794,
358        0.9900842691660192,
359        0.9927998590434373,
360        0.9950828645255290,
361        0.9969322929775997,
362        0.9983473449340834,
363        0.9993274305065947,
364        0.9998723404457334
365};
366
367static double Gauss150Wt[150]={
368        0.0003276086705538,
369        0.0007624720924706,
370        0.0011976474864367,
371        0.0016323569986067,
372        0.0020663664924131,
373        0.0024994789888943,
374        0.0029315036836558,
375        0.0033622516236779,
376        0.0037915348363451,
377        0.0042191661429919,
378        0.0046449591497966,
379        0.0050687282939456,
380        0.0054902889094487,
381        0.0059094573005900,
382        0.0063260508184704,
383        0.0067398879387430,
384        0.0071507883396855,
385        0.0075585729801782,
386        0.0079630641773633,
387        0.0083640856838475,
388        0.0087614627643580,
389        0.0091550222717888,
390        0.0095445927225849,
391        0.0099300043714212,
392        0.0103110892851360,
393        0.0106876814158841,
394        0.0110596166734735,
395        0.0114267329968529,
396        0.0117888704247183,
397        0.0121458711652067,
398        0.0124975796646449,
399        0.0128438426753249,
400        0.0131845093222756,
401        0.0135194311690004,
402        0.0138484622795371,
403        0.0141714592928592,
404        0.0144882814685445,
405        0.0147987907597169,
406        0.0151028518701744,
407        0.0154003323133401,
408        0.0156911024699895,
409        0.0159750356447283,
410        0.0162520081211971,
411        0.0165218992159766,
412        0.0167845913311726,
413        0.0170399700056559,
414        0.0172879239649355,
415        0.0175283451696437,
416        0.0177611288626114,
417        0.0179861736145128,
418        0.0182033813680609,
419        0.0184126574807331,
420        0.0186139107660094,
421        0.0188070535331042,
422        0.0189920016251754,
423        0.0191686744559934,
424        0.0193369950450545,
425        0.0194968900511231,
426        0.0196482898041878,
427        0.0197911283358190,
428        0.0199253434079123,
429        0.0200508765398072,
430        0.0201676730337687,
431        0.0202756819988200,
432        0.0203748563729175,
433        0.0204651529434560,
434        0.0205465323660984,
435        0.0206189591819181,
436        0.0206824018328499,
437        0.0207368326754401,
438        0.0207822279928917,
439        0.0208185680053983,
440        0.0208458368787627,
441        0.0208640227312962,
442        0.0208731176389954,
443        0.0208731176389954,
444        0.0208640227312962,
445        0.0208458368787627,
446        0.0208185680053983,
447        0.0207822279928917,
448        0.0207368326754401,
449        0.0206824018328499,
450        0.0206189591819181,
451        0.0205465323660984,
452        0.0204651529434560,
453        0.0203748563729175,
454        0.0202756819988200,
455        0.0201676730337687,
456        0.0200508765398072,
457        0.0199253434079123,
458        0.0197911283358190,
459        0.0196482898041878,
460        0.0194968900511231,
461        0.0193369950450545,
462        0.0191686744559934,
463        0.0189920016251754,
464        0.0188070535331042,
465        0.0186139107660094,
466        0.0184126574807331,
467        0.0182033813680609,
468        0.0179861736145128,
469        0.0177611288626114,
470        0.0175283451696437,
471        0.0172879239649355,
472        0.0170399700056559,
473        0.0167845913311726,
474        0.0165218992159766,
475        0.0162520081211971,
476        0.0159750356447283,
477        0.0156911024699895,
478        0.0154003323133401,
479        0.0151028518701744,
480        0.0147987907597169,
481        0.0144882814685445,
482        0.0141714592928592,
483        0.0138484622795371,
484        0.0135194311690004,
485        0.0131845093222756,
486        0.0128438426753249,
487        0.0124975796646449,
488        0.0121458711652067,
489        0.0117888704247183,
490        0.0114267329968529,
491        0.0110596166734735,
492        0.0106876814158841,
493        0.0103110892851360,
494        0.0099300043714212,
495        0.0095445927225849,
496        0.0091550222717888,
497        0.0087614627643580,
498        0.0083640856838475,
499        0.0079630641773633,
500        0.0075585729801782,
501        0.0071507883396855,
502        0.0067398879387430,
503        0.0063260508184704,
504        0.0059094573005900,
505        0.0054902889094487,
506        0.0050687282939456,
507        0.0046449591497966,
508        0.0042191661429919,
509        0.0037915348363451,
510        0.0033622516236779,
511        0.0029315036836558,
512        0.0024994789888943,
513        0.0020663664924131,
514        0.0016323569986067,
515        0.0011976474864367,
516        0.0007624720924706,
517        0.0003276086705538
518};
519
Note: See TracBrowser for help on using the repository browser.