Changeset 821


Ignore:
Timestamp:
Jul 19, 2011 7:24:27 PM (11 years ago)
Author:
ajj
Message:
  • Turned off file-check for me as well as Steve. DoCheck?(0) is getting boring!
Location:
sans
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • sans/Dev/trunk/NCNR_User_Procedures/Analysis/Models/NewModels_2008/PolyCore_and_NShells_v40.ipf

    r819 r821  
    289289        Wave cw,yw,xw 
    290290         
    291 //#if exists("PolyOneShellX") 
    292 //      yw = PolyOneShellX(cw,xw) 
    293 //#else 
     291#if exists("PolyOneShellX") 
     292        yw = PolyOneShellX(cw,xw) 
     293#else 
    294294        yw = fPolyOneShell(cw,xw) 
    295 //#endif 
     295#endif 
    296296        return(0) 
    297297End 
     
    300300        Wave cw,yw,xw 
    301301         
    302 //#if exists("PolyTwoShellX") 
     302#if exists("PolyTwoShellX") 
    303303        yw = PolyTwoShellX(cw,xw) 
    304 //#else 
     304#else 
    305305        yw = fPolyTwoShell(cw,xw) 
    306 //#endif 
     306#endif 
    307307        return(0) 
    308308End 
     
    311311        Wave cw,yw,xw 
    312312         
    313 //#if exists("PolyThreeShellX") 
    314 //      yw =PolyThreeShellX(cw,xw) 
    315 //#else 
     313#if exists("PolyThreeShellX") 
     314        yw =PolyThreeShellX(cw,xw) 
     315#else 
    316316        yw = fPolyThreeShell(cw,xw) 
    317 //#endif 
     317#endif 
    318318        return(0) 
    319319End 
     
    322322        Wave cw,yw,xw 
    323323         
    324 //#if exists("PolyFourShellX") 
    325 //      yw = PolyFourShellX(cw,xw) 
    326 //#else 
     324#if exists("PolyFourShellX") 
     325        yw = PolyFourShellX(cw,xw) 
     326#else 
    327327        yw = fPolyFourShell(cw,xw) 
    328 //#endif 
     328#endif 
    329329        return(0) 
    330330End 
     
    358358        NVAR/Z gNord=gNord 
    359359        if(! NVAR_Exists(gNord) ) 
    360                 nord=120                //use 76 pts as default 
     360                nord=76         //use 76 pts as default 
    361361        else 
    362362                if( (gNord == 20) || (gNord ==76) ) 
     
    465465//      NVAR/Z gNord=gNord 
    466466//      if(! NVAR_Exists(gNord) ) 
    467                 nord=120                //use 76 pts as default 
     467                nord=76         //use 76 pts as default 
    468468//      else 
    469469//              if( (gNord == 20) || (gNord ==76) ) 
     
    575575//      NVAR/Z gNord=gNord 
    576576//      if(! NVAR_Exists(gNord) ) 
    577                 nord=120                //use 76 pts as default 
     577                nord=150                //use 150 pts as default - same as XOP 
    578578//      else 
    579579//              if( (gNord == 20) || (gNord ==76) ) 
     
    690690        NVAR/Z gNord=gNord 
    691691//      if(! NVAR_Exists(gNord) ) 
    692                 nord=120                //use 76 pts as default 
     692                nord=150                //use 150 pts as default - same as XOP 
    693693//      else 
    694694//              if( (gNord == 20) || (gNord ==76) ) 
  • sans/Dev/trunk/NCNR_User_Procedures/Analysis/Packages/ModelPicker/SANSModelPicker_v40.ipf

    r775 r821  
    8181                // always turn off file checking for me 
    8282                checkForFiles = !(stringmatch(ParseFilePath(0,SpecialDirPath("Desktop",0,0,0),":",1,1),"s*ine"))                //zero for me 
     83                checkForFiles = !(stringmatch(ParseFilePath(0,SpecialDirPath("Desktop",0,0,0),":",1,1),"ajj"))          //and zero for me 
     84 
    8385                 
    8486                // turn off file checking if the proper alias to the NCNR procedures is there 
  • sans/XOP_Dev/SANSAnalysis/lib/libSphere.c

    r678 r821  
    17211721} 
    17221722 
     1723// Use 150 point integral. This is probably excessive, but we need at least 100 points  
     1724// and this avoids yet another set of abcissae and weights 
    17231725double 
    17241726PolyThreeShell(double dp[], double x) 
     
    17271729        double va,vb,summ,yyy,zi; 
    17281730        double answer,zp1,zp2,zp3,vpoly,range,temp_3sf[11],pi; 
    1729         int nord=76,ii; 
     1731        int nord=150,ii; 
    17301732        double thick1,thick2,thick3; 
    17311733        double rhoshel1,rhoshel2,rhoshel3; 
     
    17741776        for(ii=0;ii<nord;ii+=1) { 
    17751777                // calculate Gauss points on integration interval (r-value for evaluation) 
    1776                 zi = ( Gauss76Z[ii]*(vb-va) + vb + va )/2.0; 
     1778                zi = ( Gauss150Z[ii]*(vb-va) + vb + va )/2.0; 
    17771779                temp_3sf[1] = zi; 
    1778                 yyy = Gauss76Wt[ii] * SchulzPoint(zi,rcore,zz) * ThreeShell(temp_3sf,x); 
     1780                yyy = Gauss150Wt[ii] * SchulzPoint(zi,rcore,zz) * ThreeShell(temp_3sf,x); 
    17791781                //un-normalize by volume 
    17801782                yyy *= 4.0*pi/3.0*pow((zi+thick1+thick2+thick3),3); 
     
    17981800} 
    17991801 
     1802 
     1803// Use 150 point integral. This is probably excessive, but we need at least 100 points  
     1804// and this avoids yet another set of abcissae and weights 
    18001805double 
    18011806PolyFourShell(double dp[], double x) 
     
    18041809        double va,vb,summ,yyy,zi; 
    18051810        double answer,zp1,zp2,zp3,vpoly,range,temp_4sf[13],pi; 
    1806         int nord=76,ii; 
     1811        int nord=150,ii; 
    18071812        double thick1,thick2,thick3,thick4; 
    18081813        double rhoshel1,rhoshel2,rhoshel3,rhoshel4; 
     
    18551860        for(ii=0;ii<nord;ii+=1) { 
    18561861                // calculate Gauss points on integration interval (r-value for evaluation) 
    1857                 zi = ( Gauss76Z[ii]*(vb-va) + vb + va )/2.0; 
     1862                zi = ( Gauss150Z[ii]*(vb-va) + vb + va )/2.0; 
    18581863                temp_4sf[1] = zi; 
    1859                 yyy = Gauss76Wt[ii] * SchulzPoint(zi,rcore,zz) * FourShell(temp_4sf,x); 
     1864                yyy = Gauss150Wt[ii] * SchulzPoint(zi,rcore,zz) * FourShell(temp_4sf,x); 
    18601865                //un-normalize by volume 
    18611866                yyy *= 4.0*pi/3.0*pow((zi+thick1+thick2+thick3+thick4),3); 
Note: See TracChangeset for help on using the changeset viewer.