source: sans/Analysis/branches/ajj_23APR07/IGOR_Package_Files/Put in User Procedures/SANS_Models_v3.00/UnifEllipsoid_and_Struct.ipf @ 153

Last change on this file since 153 was 151, checked in by srkline, 15 years ago

(1) - cursors can now be used to select a subrange of USANS data to fit. This is done by th fit wrapper, assigning a subrange of resW to the struct

(2) all of the smeared model functions are now in the latest form of Smear_Model_N() that is NOT a pointwise calculation anymore, since the USANS matrix smearing in inherently not so.

File size: 15.7 KB
Line 
1#pragma rtGlobals=1             // Use modern global access method.
2#pragma IgorVersion = 6.0
3
4// be sure to include all the necessary files...
5
6#include "EffectiveDiameter"
7#include "UniformEllipsoid"
8
9#include "HardSphereStruct"
10#include "HPMSA"
11#include "SquareWellStruct"
12#include "StickyHardSphereStruct"
13
14Proc PlotEllipsoid_HS(num,qmin,qmax)
15        Variable num=128,qmin=0.001,qmax=0.7
16        Prompt num "Enter number of data points for model: "
17        Prompt qmin "Enter minimum q-value (^-1) for model: "
18        Prompt qmax "Enter maximum q-value (^-1) for model: "
19       
20        Make/O/D/n=(num) xwave_EOR_HS,ywave_EOR_HS
21        xwave_EOR_HS =  alog(log(qmin) + x*((log(qmax)-log(qmin))/num))
22        Make/O/D coef_EOR_HS = {0.01,20.,400,3.0e-6,0.01}
23        make/o/t parameters_EOR_HS = {"volume fraction","R(a) rotation axis (A)","R(b) (A)","contrast (A^-2)","incoh. bkg (cm^-1)"}
24        Edit parameters_EOR_HS,coef_EOR_HS
25       
26        Variable/G root:g_EOR_HS
27        g_EOR_HS := Ellipsoid_HS(coef_EOR_HS,ywave_EOR_HS,xwave_EOR_HS)
28        Display ywave_EOR_HS vs xwave_EOR_HS
29        ModifyGraph log=1,marker=29,msize=2,mode=4
30        Label bottom "q (\\S-1\\M)"
31        Label left "Intensity (cm\\S-1\\M)"
32        AutoPositionWindow/M=1/R=$(WinName(0,1)) $WinName(0,2)
33End
34
35// - sets up a dependency to a wrapper, not the actual SmearedModelFunction
36Proc PlotSmearedEllipsoid_HS(str)                                                               
37        String str
38        Prompt str,"Pick the data folder containing the resolution you want",popup,getAList(4)
39       
40        // if any of the resolution waves are missing => abort
41        if(ResolutionWavesMissingDF(str))               //updated to NOT use global strings (in GaussUtils)
42                Abort
43        endif
44       
45        SetDataFolder $("root:"+str)
46       
47        // Setup parameter table for model function
48        Make/O/D smear_coef_EOR_HS = {0.01,20.,400,3.0e-6,0.01}
49        make/o/t smear_parameters_EOR_HS = {"volume fraction","R(a) rotation axis (A)","R(b) (A)","contrast (A^-2)","incoh. bkg (cm^-1)"}
50        Edit smear_parameters_EOR_HS,smear_coef_EOR_HS                                 
51       
52        // output smeared intensity wave, dimensions are identical to experimental QSIG values
53        // make extra copy of experimental q-values for easy plotting
54        Duplicate/O $(str+"_q") smeared_EOR_HS,smeared_qvals                           
55        SetScale d,0,0,"1/cm",smeared_EOR_HS                                                   
56                                       
57        Variable/G gs_EOR_HS=0
58        gs_EOR_HS := fSmearedEllipsoid_HS(smear_coef_EOR_HS,smeared_EOR_HS,smeared_qvals)       //this wrapper fills the STRUCT
59       
60        Display smeared_EOR_HS vs smeared_qvals                                                                 
61        ModifyGraph log=1,marker=29,msize=2,mode=4
62        Label bottom "q (\\S-1\\M)"
63        Label left "Intensity (cm\\S-1\\M)"
64
65        AutoPositionWindow/M=1/R=$(WinName(0,1)) $WinName(0,2)
66       
67        SetDataFolder root:
68End
69       
70
71//AAO function
72Function Ellipsoid_HS(w,yW,xW) : FitFunc
73        Wave w,yW,xW
74       
75        Variable inten,aa,bb
76        aa=w[1]
77        bb=w[2]
78       
79        //setup form factor coefficient wave
80        Make/O/D/N=5 form_EOR_HS
81        form_EOR_HS[0] = 1
82        form_EOR_HS[1] = w[1]
83        form_EOR_HS[2] = w[2]
84        form_EOR_HS[3] = w[3]
85        form_EOR_HS[4] = 0
86       
87        //setup structure factor coefficient wave
88        Make/O/D/N=2 struct_EOR_HS
89        struct_EOR_HS[0] = 0.5*DiamEllip(aa,bb)
90        struct_EOR_HS[1] = w[0]
91       
92        //calculate each and combine
93        Duplicate/O xw, temp_EOR_HS_PQ,temp_EOR_HS_SQ
94        EllipsoidForm(form_EOR_HS,temp_EOR_HS_PQ,xW)
95        HardSphereStruct(struct_EOR_HS,temp_EOR_HS_SQ,xW)
96        yw = temp_EOR_HS_PQ * temp_EOR_HS_SQ
97        yw *= w[0]
98        yw += w[4]
99       
100        //cleanup waves
101//      Killwaves/Z form_EOR_HS,struct_EOR_HS
102       
103        return (0)
104End
105
106Proc PlotEllipsoid_SW(num,qmin,qmax)
107        Variable num=128,qmin=0.001,qmax=0.7
108        Prompt num "Enter number of data points for model: "
109        Prompt qmin "Enter minimum q-value (^-1) for model: "
110        Prompt qmax "Enter maximum q-value (^-1) for model: "
111       
112        Make/O/D/n=(num) xwave_EOR_SW,ywave_EOR_SW
113        xwave_EOR_SW =  alog(log(qmin) + x*((log(qmax)-log(qmin))/num))
114        Make/O/D coef_EOR_SW = {0.01,20.,400,3.0e-6,1.0,1.2,0.01}
115        make/o/t parameters_EOR_SW = {"volume fraction","R(a) rotation axis (A)","R(b) (A)","contrast (A^-2)","well depth (kT)","well width (diam.)","incoh. bkg (cm^-1)"}
116        Edit parameters_EOR_SW,coef_EOR_SW
117       
118        Variable/G root:g_EOR_SW
119        g_EOR_SW := Ellipsoid_SW(coef_EOR_SW,ywave_EOR_SW,xwave_EOR_SW)
120        Display ywave_EOR_SW vs xwave_EOR_SW
121        ModifyGraph log=1,marker=29,msize=2,mode=4
122        Label bottom "q (\\S-1\\M)"
123        Label left "Intensity (cm\\S-1\\M)"
124
125        AutoPositionWindow/M=1/R=$(WinName(0,1)) $WinName(0,2)
126End
127
128// - sets up a dependency to a wrapper, not the actual SmearedModelFunction
129Proc PlotSmearedEllipsoid_SW(str)                                                               
130        String str
131        Prompt str,"Pick the data folder containing the resolution you want",popup,getAList(4)
132       
133        // if any of the resolution waves are missing => abort
134        if(ResolutionWavesMissingDF(str))               //updated to NOT use global strings (in GaussUtils)
135                Abort
136        endif
137       
138        SetDataFolder $("root:"+str)
139       
140        // Setup parameter table for model function
141        Make/O/D smear_coef_EOR_SW = {0.01,20.,400,3.0e-6,1.0,1.2,0.01}
142        make/o/t smear_parameters_EOR_SW = {"volume fraction","R(a) rotation axis (A)","R(b) (A)","contrast (A^-2)","well depth (kT)","well width (diam.)","incoh. bkg (cm^-1)"}
143        Edit smear_parameters_EOR_SW,smear_coef_EOR_SW                                 
144       
145        // output smeared intensity wave, dimensions are identical to experimental QSIG values
146        // make extra copy of experimental q-values for easy plotting
147        Duplicate/O $(str+"_q") smeared_EOR_SW,smeared_qvals                           
148        SetScale d,0,0,"1/cm",smeared_EOR_SW                                                   
149                                       
150        Variable/G gs_EOR_SW=0
151        gs_EOR_SW := fSmearedEllipsoid_SW(smear_coef_EOR_SW,smeared_EOR_SW,smeared_qvals)       //this wrapper fills the STRUCT
152       
153        Display smeared_EOR_SW vs smeared_qvals                                                                 
154        ModifyGraph log=1,marker=29,msize=2,mode=4
155        Label bottom "q (\\S-1\\M)"
156        Label left "Intensity (cm\\S-1\\M)"
157
158        AutoPositionWindow/M=1/R=$(WinName(0,1)) $WinName(0,2)
159       
160        SetDataFolder root:
161End
162       
163
164Function Ellipsoid_SW(w,yw,xw) : FitFunc
165        Wave w,yw,xw
166       
167        Variable inten,aa,bb
168        aa=w[1]
169        bb=w[2]
170       
171        //setup form factor coefficient wave
172        Make/O/D/N=5 form_EOR_SW
173        form_EOR_SW[0] = 1
174        form_EOR_SW[1] = w[1]
175        form_EOR_SW[2] = w[2]
176        form_EOR_SW[3] = w[3]
177        form_EOR_SW[4] = 0
178       
179        //setup structure factor coefficient wave
180        Make/O/D/N=4 struct_EOR_SW
181        struct_EOR_SW[0] = 0.5*DiamEllip(aa,bb)
182        struct_EOR_SW[1] = w[0]
183        struct_EOR_SW[2] = w[4]
184        struct_EOR_SW[3] = w[5]
185       
186        //calculate each and combine
187        Duplicate/O xw, temp_EOR_SW_PQ,temp_EOR_SW_SQ
188        EllipsoidForm(form_EOR_SW,temp_EOR_SW_PQ,xw)
189        SquareWellStruct(struct_EOR_SW,temp_EOR_SW_SQ,xw)
190        yw = temp_EOR_SW_PQ * temp_EOR_SW_SQ
191        yw *= w[0]
192        yw += w[6]
193       
194        //cleanup waves
195//      Killwaves/Z form_EOR_SW,struct_EOR_SW
196       
197        return (0)
198End
199
200Proc PlotEllipsoid_SC(num,qmin,qmax)
201        Variable num=128,qmin=0.001,qmax=0.7
202        Prompt num "Enter number of data points for model: "
203        Prompt qmin "Enter minimum q-value (^-1) for model: "
204        Prompt qmax "Enter maximum q-value (^-1) for model: "
205       
206        if(!DataFolderExists(":HayPenMSA"))
207                NewDataFolder :HayPenMSA
208        endif
209        Make/O/D/N=17 :HayPenMSA:gMSAWave
210       
211        Make/O/D/n=(num) xwave_EOR_SC,ywave_EOR_SC
212        xwave_EOR_SC =  alog(log(qmin) + x*((log(qmax)-log(qmin))/num))
213        Make/O/D coef_EOR_SC = {0.01,20.,400,3.0e-6,20,0,298,78,0.01}
214        make/o/t parameters_EOR_SC = {"volume fraction","R(a) rotation axis (A)","R(b) (A)","contrast (A^-2)","charge","movalent salt(M)","Temperature (K)","dielectric const","incoh. bkg (cm^-1)"}
215        Edit parameters_EOR_SC,coef_EOR_SC
216
217        Variable/G root:g_EOR_SC
218        g_EOR_SC := Ellipsoid_SC(coef_EOR_SC,ywave_EOR_SC,xwave_EOR_SC)
219        Display ywave_EOR_SC vs xwave_EOR_SC
220        ModifyGraph log=1,marker=29,msize=2,mode=4
221        Label bottom "q (\\S-1\\M)"
222        Label left "Intensity (cm\\S-1\\M)"
223
224        AutoPositionWindow/M=1/R=$(WinName(0,1)) $WinName(0,2)
225End
226
227// - sets up a dependency to a wrapper, not the actual SmearedModelFunction
228Proc PlotSmearedEllipsoid_SC(str)                                                               
229        String str
230        Prompt str,"Pick the data folder containing the resolution you want",popup,getAList(4)
231       
232        // if any of the resolution waves are missing => abort
233        if(ResolutionWavesMissingDF(str))               //updated to NOT use global strings (in GaussUtils)
234                Abort
235        endif
236       
237        SetDataFolder $("root:"+str)
238       
239        if(!DataFolderExists(":HayPenMSA"))
240                NewDataFolder :HayPenMSA
241        endif
242        Make/O/D/N=17 :HayPenMSA:gMSAWave
243       
244        // Setup parameter table for model function
245        Make/O/D smear_coef_EOR_SC = {0.01,20.,400,3.0e-6,20,0,298,78,0.01}
246        make/o/t smear_parameters_EOR_SC = {"volume fraction","R(a) rotation axis (A)","R(b) (A)","contrast (A^-2)","charge","movalent salt(M)","Temperature (K)","dielectric const","incoh. bkg (cm^-1)"}
247        Edit smear_parameters_EOR_SC,smear_coef_EOR_SC                                 
248       
249        // output smeared intensity wave, dimensions are identical to experimental QSIG values
250        // make extra copy of experimental q-values for easy plotting
251        Duplicate/O $(str+"_q") smeared_EOR_SC,smeared_qvals                           
252        SetScale d,0,0,"1/cm",smeared_EOR_SC                                                   
253                                       
254        Variable/G gs_EOR_SC=0
255        gs_EOR_SC := fSmearedEllipsoid_SC(smear_coef_EOR_SC,smeared_EOR_SC,smeared_qvals)       //this wrapper fills the STRUCT
256       
257        Display smeared_EOR_SC vs smeared_qvals                                                                 
258        ModifyGraph log=1,marker=29,msize=2,mode=4
259        Label bottom "q (\\S-1\\M)"
260        Label left "Intensity (cm\\S-1\\M)"
261
262        AutoPositionWindow/M=1/R=$(WinName(0,1)) $WinName(0,2)
263       
264        SetDataFolder root:
265End
266
267
268Function Ellipsoid_SC(w,yw,xw) : FitFunc
269        Wave w,yw,xw
270       
271        Variable inten,aa,bb
272        aa=w[1]
273        bb=w[2]
274       
275        //setup form factor coefficient wave
276        Make/O/D/N=5 form_EOR_SC
277        form_EOR_SC[0] = 1
278        form_EOR_SC[1] = w[1]
279        form_EOR_SC[2] = w[2]
280        form_EOR_SC[3] = w[3]
281        form_EOR_SC[4] = 0
282       
283        //setup structure factor coefficient wave
284        Make/O/D/N=6 struct_EOR_SC
285        struct_EOR_SC[0] = DiamEllip(aa,bb)
286        struct_EOR_SC[1] = w[4]
287        struct_EOR_SC[2] = w[0]
288        struct_EOR_SC[3] = w[6]
289        struct_EOR_SC[4] = w[5]
290        struct_EOR_SC[5] = w[7]
291       
292        //calculate each and combine
293        Duplicate/O xw, temp_EOR_SC_PQ,temp_EOR_SC_SQ
294        EllipsoidForm(form_EOR_SC,temp_EOR_SC_PQ,xw)
295        HayterPenfoldMSA(struct_EOR_SC,temp_EOR_SC_SQ,xw)
296        yw = temp_EOR_SC_PQ * temp_EOR_SC_SQ
297        yw *= w[0]
298        yw += w[8]
299       
300        //cleanup waves
301//      Killwaves/Z form_EOR_SC,struct_EOR_SC
302       
303        return (0)
304End
305
306Proc PlotEllipsoid_SHS(num,qmin,qmax)
307        Variable num=128,qmin=0.001,qmax=0.7
308        Prompt num "Enter number of data points for model: "
309        Prompt qmin "Enter minimum q-value (^-1) for model: "
310        Prompt qmax "Enter maximum q-value (^-1) for model: "
311       
312        Make/O/D/n=(num) xwave_EOR_SHS,ywave_EOR_SHS
313        xwave_EOR_SHS =  alog(log(qmin) + x*((log(qmax)-log(qmin))/num))
314        Make/O/D coef_EOR_SHS = {0.01,20.,400,3.0e-6,0.05,0.2,0.01}
315        make/o/t parameters_EOR_SHS = {"volume fraction","R(a) rotation axis (A)","R(b) (A)","contrast (A^-2)","perturbation parameter (0.1)","stickiness, tau","incoh. bkg (cm^-1)"}
316        Edit parameters_EOR_SHS,coef_EOR_SHS
317       
318        Variable/G root:g_EOR_SHS
319        g_EOR_SHS := Ellipsoid_SHS(coef_EOR_SHS,ywave_EOR_SHS,xwave_EOR_SHS)
320        Display ywave_EOR_SHS vs xwave_EOR_SHS
321        ModifyGraph log=1,marker=29,msize=2,mode=4
322        Label bottom "q (\\S-1\\M)"
323        Label left "Intensity (cm\\S-1\\M)"
324
325        AutoPositionWindow/M=1/R=$(WinName(0,1)) $WinName(0,2)
326End
327
328// - sets up a dependency to a wrapper, not the actual SmearedModelFunction
329Proc PlotSmearedEllipsoid_SHS(str)                                                             
330        String str
331        Prompt str,"Pick the data folder containing the resolution you want",popup,getAList(4)
332       
333        // if any of the resolution waves are missing => abort
334        if(ResolutionWavesMissingDF(str))               //updated to NOT use global strings (in GaussUtils)
335                Abort
336        endif
337       
338        SetDataFolder $("root:"+str)
339       
340        // Setup parameter table for model function
341        Make/O/D smear_coef_EOR_SHS = {0.01,20.,400,3.0e-6,0.05,0.2,0.01}
342        make/o/t smear_parameters_EOR_SHS = {"volume fraction","R(a) rotation axis (A)","R(b) (A)","contrast (A^-2)","perturbation parameter (0.1)","stickiness, tau","incoh. bkg (cm^-1)"}
343        Edit smear_parameters_EOR_SHS,smear_coef_EOR_SHS                                       
344       
345        // output smeared intensity wave, dimensions are identical to experimental QSIG values
346        // make extra copy of experimental q-values for easy plotting
347        Duplicate/O $(str+"_q") smeared_EOR_SHS,smeared_qvals                           
348        SetScale d,0,0,"1/cm",smeared_EOR_SHS                                                   
349                                       
350        Variable/G gs_EOR_SHS=0
351        gs_EOR_SHS := fSmearedEllipsoid_SHS(smear_coef_EOR_SHS,smeared_EOR_SHS,smeared_qvals)   //this wrapper fills the STRUCT
352       
353        Display smeared_EOR_SHS vs smeared_qvals                                                                       
354        ModifyGraph log=1,marker=29,msize=2,mode=4
355        Label bottom "q (\\S-1\\M)"
356        Label left "Intensity (cm\\S-1\\M)"
357
358        AutoPositionWindow/M=1/R=$(WinName(0,1)) $WinName(0,2)
359       
360        SetDataFolder root:
361End
362       
363
364Function Ellipsoid_SHS(w,yw,xw) : FitFunc
365        Wave w,yw,xw
366       
367        Variable inten,aa,bb
368        aa=w[1]
369        bb=w[2]
370       
371        //setup form factor coefficient wave
372        Make/O/D/N=5 form_EOR_SHS
373        form_EOR_SHS[0] = 1
374        form_EOR_SHS[1] = w[1]
375        form_EOR_SHS[2] = w[2]
376        form_EOR_SHS[3] = w[3]
377        form_EOR_SHS[4] = 0
378       
379        //setup structure factor coefficient wave
380        Make/O/D/N=4 struct_EOR_SHS
381        struct_EOR_SHS[0] = 0.5*DiamEllip(aa,bb)
382        struct_EOR_SHS[1] = w[0]
383        struct_EOR_SHS[2] = w[4]
384        struct_EOR_SHS[3] = w[5]
385       
386        //calculate each and combine
387        Duplicate/O xw, temp_EOR_SHS_PQ,temp_EOR_SHS_SQ
388        EllipsoidForm(form_EOR_SHS,temp_EOR_SHS_PQ,xw)
389        StickyHS_Struct(struct_EOR_SHS,temp_EOR_SHS_SQ,xw)
390        yw = temp_EOR_SHS_PQ *temp_EOR_SHS_SQ
391        yw *= w[0]
392        yw += w[6]
393       
394        //cleanup waves
395//      Killwaves/Z form_EOR_SHS,struct_EOR_SHS
396       
397        return (0)
398End
399
400
401
402// this is all there is to the smeared calculation!
403Function SmearedEllipsoid_HS(s) :FitFunc
404        Struct ResSmearAAOStruct &s
405
406//      the name of your unsmeared model is the first argument
407        Smear_Model_20(Ellipsoid_HS,s.coefW,s.xW,s.yW,s.resW)
408
409        return(0)
410End
411
412// this is all there is to the smeared calculation!
413Function SmearedEllipsoid_SW(s) :FitFunc
414        Struct ResSmearAAOStruct &s
415
416//      the name of your unsmeared model is the first argument
417        Smear_Model_20(Ellipsoid_SW,s.coefW,s.xW,s.yW,s.resW)
418
419        return(0)
420End
421
422// this is all there is to the smeared calculation!
423Function SmearedEllipsoid_SC(s) :FitFunc
424        Struct ResSmearAAOStruct &s
425
426//      the name of your unsmeared model is the first argument
427        Smear_Model_20(Ellipsoid_SC,s.coefW,s.xW,s.yW,s.resW)
428
429        return(0)
430End
431
432// this is all there is to the smeared calculation!
433Function SmearedEllipsoid_SHS(s) :FitFunc
434        Struct ResSmearAAOStruct &s
435
436//      the name of your unsmeared model is the first argument
437        Smear_Model_20(Ellipsoid_SHS,s.coefW,s.xW,s.yW,s.resW)
438
439        return(0)
440End
441
442//wrapper to calculate the smeared model as an AAO-Struct
443// fills the struct and calls the ususal function with the STRUCT parameter
444//
445// used only for the dependency, not for fitting
446//
447Function fSmearedEllipsoid_HS(coefW,yW,xW)
448        Wave coefW,yW,xW
449       
450        String str = getWavesDataFolder(yW,0)
451        String DF="root:"+str+":"
452       
453        WAVE resW = $(DF+str+"_res")
454       
455        STRUCT ResSmearAAOStruct fs
456        WAVE fs.coefW = coefW   
457        WAVE fs.yW = yW
458        WAVE fs.xW = xW
459        WAVE fs.resW = resW
460       
461        Variable err
462        err = SmearedEllipsoid_HS(fs)
463       
464        return (0)
465End
466
467//wrapper to calculate the smeared model as an AAO-Struct
468// fills the struct and calls the ususal function with the STRUCT parameter
469//
470// used only for the dependency, not for fitting
471//
472Function fSmearedEllipsoid_SW(coefW,yW,xW)
473        Wave coefW,yW,xW
474       
475        String str = getWavesDataFolder(yW,0)
476        String DF="root:"+str+":"
477       
478        WAVE resW = $(DF+str+"_res")
479       
480        STRUCT ResSmearAAOStruct fs
481        WAVE fs.coefW = coefW   
482        WAVE fs.yW = yW
483        WAVE fs.xW = xW
484        WAVE fs.resW = resW
485       
486        Variable err
487        err = SmearedEllipsoid_SW(fs)
488       
489        return (0)
490End
491
492//wrapper to calculate the smeared model as an AAO-Struct
493// fills the struct and calls the ususal function with the STRUCT parameter
494//
495// used only for the dependency, not for fitting
496//
497Function fSmearedEllipsoid_SC(coefW,yW,xW)
498        Wave coefW,yW,xW
499       
500        String str = getWavesDataFolder(yW,0)
501        String DF="root:"+str+":"
502       
503        WAVE resW = $(DF+str+"_res")
504       
505        STRUCT ResSmearAAOStruct fs
506        WAVE fs.coefW = coefW   
507        WAVE fs.yW = yW
508        WAVE fs.xW = xW
509        WAVE fs.resW = resW
510       
511        Variable err
512        err = SmearedEllipsoid_SC(fs)
513       
514        return (0)
515End
516
517//wrapper to calculate the smeared model as an AAO-Struct
518// fills the struct and calls the ususal function with the STRUCT parameter
519//
520// used only for the dependency, not for fitting
521//
522Function fSmearedEllipsoid_SHS(coefW,yW,xW)
523        Wave coefW,yW,xW
524       
525        String str = getWavesDataFolder(yW,0)
526        String DF="root:"+str+":"
527       
528        WAVE resW = $(DF+str+"_res")
529       
530        STRUCT ResSmearAAOStruct fs
531        WAVE fs.coefW = coefW   
532        WAVE fs.yW = yW
533        WAVE fs.xW = xW
534        WAVE fs.resW = resW
535       
536        Variable err
537        err = SmearedEllipsoid_SHS(fs)
538       
539        return (0)
540End
Note: See TracBrowser for help on using the repository browser.