Changeset 62


Ignore:
Timestamp:
Mar 26, 2007 1:49:49 PM (15 years ago)
Author:
srkline
Message:

scale factor added to final Beaucage model

File:
1 edited

Legend:

Unmodified
Added
Removed
  • sans/Analysis/trunk/Put in User Procedures/SANS_Models_v3.00/NewModels_2006/Beaucage.ipf

    r61 r62  
    198198         
    199199        Variable ans,erf1,prec=1e-15 
    200         Variable G1,Rg1,B1,Pow1,bkg 
    201          
    202         G1 = w[0] 
    203         Rg1 = w[1] 
    204         B1 = w[2] 
    205         Pow1 = w[3] 
    206         bkg = w[4] 
     200        Variable G1,Rg1,B1,Pow1,bkg,scale 
     201         
     202        scale = w[0] 
     203        G1 = w[1] 
     204        Rg1 = w[2] 
     205        B1 = w[3] 
     206        Pow1 = w[4] 
     207        bkg = w[5] 
    207208         
    208209        erf1 = erf( (x*Rg1/sqrt(6)) ,prec) 
     
    210211        ans = G1*exp(-x*x*Rg1*Rg1/3) 
    211212        ans += B1*(erf1^3/x)^Pow1 
     213        ans *= scale 
    212214        ans += bkg 
    213215         
     
    220222         
    221223        Variable ans,G1,Rg1,B1,G2,Rg2,B2,Pow1,Pow2,bkg 
    222         Variable erf1,erf2,prec=1e-15 
     224        Variable erf1,erf2,prec=1e-15,scale 
    223225         
    224226        //Rsub = Rs 
    225         G1 = w[0]       //equivalent to I(0) 
    226         Rg1 = w[1] 
    227         B1 = w[2] 
    228         Pow1 = w[3] 
    229         G2 = w[4] 
    230         Rg2 = w[5] 
    231         B2 = w[6] 
    232         Pow2 = w[7] 
    233         bkg = w[8] 
     227        scale = w[0] 
     228        G1 = w[1]       //equivalent to I(0) 
     229        Rg1 = w[2] 
     230        B1 = w[3] 
     231        Pow1 = w[4] 
     232        G2 = w[5] 
     233        Rg2 = w[6] 
     234        B2 = w[7] 
     235        Pow2 = w[8] 
     236        bkg = w[9] 
    234237         
    235238        erf1 = erf( (x*Rg1/sqrt(6)) ,prec) 
     
    241244        ans += G2*exp(-x*x*Rg2*Rg2/3) 
    242245        ans += B2*(erf2^3/x)^Pow2 
    243          
     246        ans *= scale 
    244247        ans += bkg 
    245248         
     
    254257        Variable ans,G1,Rg1,B1,G2,Rg2,B2,Pow1,Pow2,bkg 
    255258        Variable G3,Rg3,B3,Pow3,erf3 
    256         Variable erf1,erf2,prec=1e-15 
     259        Variable erf1,erf2,prec=1e-15,scale 
    257260         
    258261        //Rsub = Rs 
    259         G1 = w[0]       //equivalent to I(0) 
    260         Rg1 = w[1] 
    261         B1 = w[2] 
    262         Pow1 = w[3] 
    263         G2 = w[4] 
    264         Rg2 = w[5] 
    265         B2 = w[6] 
    266         Pow2 = w[7] 
    267         G3 = w[8] 
    268         Rg3 = w[9] 
    269         B3 = w[10] 
    270         Pow3 = w[11] 
    271         bkg = w[12] 
     262        scale = w[0] 
     263        G1 = w[1]       //equivalent to I(0) 
     264        Rg1 = w[2] 
     265        B1 = w[3] 
     266        Pow1 = w[4] 
     267        G2 = w[5] 
     268        Rg2 = w[6] 
     269        B2 = w[7] 
     270        Pow2 = w[8] 
     271        G3 = w[9] 
     272        Rg3 = w[10] 
     273        B3 = w[11] 
     274        Pow3 = w[12] 
     275        bkg = w[13] 
    272276         
    273277        erf1 = erf( (x*Rg1/sqrt(6)) ,prec) 
     
    279283        ans += G2*exp(-x*x*Rg2*Rg2/3) + B2*exp(-x*x*Rg3*Rg3/3)*(erf2^3/x)^Pow2 
    280284        ans += G3*exp(-x*x*Rg3*Rg3/3) + B3*(erf3^3/x)^Pow3 
    281          
     285        ans *= scale 
    282286        ans += bkg 
    283287         
     
    292296        Variable G3,Rg3,B3,Pow3,erf3 
    293297        Variable G4,Rg4,B4,Pow4,erf4 
    294         Variable erf1,erf2,prec=1e-15 
     298        Variable erf1,erf2,prec=1e-15,scale 
    295299         
    296300        //Rsub = Rs 
    297         G1 = w[0]       //equivalent to I(0) 
    298         Rg1 = w[1] 
    299         B1 = w[2] 
    300         Pow1 = w[3] 
    301         G2 = w[4] 
    302         Rg2 = w[5] 
    303         B2 = w[6] 
    304         Pow2 = w[7] 
    305         G3 = w[8] 
    306         Rg3 = w[9] 
    307         B3 = w[10] 
    308         Pow3 = w[11] 
    309         G4 = w[12] 
    310         Rg4 = w[13] 
    311         B4 = w[14] 
    312         Pow4 = w[15] 
    313         bkg = w[16] 
     301        scale = w[0] 
     302        G1 = w[1]       //equivalent to I(0) 
     303        Rg1 = w[2] 
     304        B1 = w[3] 
     305        Pow1 = w[4] 
     306        G2 = w[5] 
     307        Rg2 = w[6] 
     308        B2 = w[7] 
     309        Pow2 = w[8] 
     310        G3 = w[9] 
     311        Rg3 = w[10] 
     312        B3 = w[11] 
     313        Pow3 = w[12] 
     314        G4 = w[13] 
     315        Rg4 = w[14] 
     316        B4 = w[15] 
     317        Pow4 = w[16] 
     318        bkg = w[17] 
    314319         
    315320        erf1 = erf( (x*Rg1/sqrt(6)) ,prec) 
     
    322327        ans += G3*exp(-x*x*Rg3*Rg3/3) + B3*exp(-x*x*Rg4*Rg4/3)*(erf3^3/x)^Pow3 
    323328        ans += G4*exp(-x*x*Rg4*Rg4/3) + B4*(erf4^3/x)^Pow4 
    324          
     329        ans *= scale 
    325330        ans += bkg 
    326331         
Note: See TracChangeset for help on using the changeset viewer.