Changeset 5f0dcab in sasview for sansmodels/src/sans/models/libigor/libSphere.h
- Timestamp:
- Dec 7, 2011 3:39:05 PM (12 years ago)
- Branches:
- master, ESS_GUI, ESS_GUI_Docs, ESS_GUI_batch_fitting, ESS_GUI_bumps_abstraction, ESS_GUI_iss1116, ESS_GUI_iss879, ESS_GUI_iss959, ESS_GUI_opencl, ESS_GUI_ordering, ESS_GUI_sync_sascalc, costrafo411, magnetic_scatt, release-4.1.1, release-4.1.2, release-4.2.2, release_4.0.1, ticket-1009, ticket-1094-headless, ticket-1242-2d-resolution, ticket-1243, ticket-1249, ticket885, unittest-saveload
- Children:
- feadd6f
- Parents:
- 2605da22
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sansmodels/src/sans/models/libigor/libSphere.h
r6a2c931 r5f0dcab 49 49 double SC_Integrand(double w[], double qq, double xx, double yy); 50 50 double SCeval(double Theta, double Phi, double temp3, double temp4, double temp5); 51 52 53 54 static double55 gammln(double xx) {56 double x,y,tmp,ser;57 static double cof[6]={76.18009172947146,-86.50532032941677,58 24.01409824083091,-1.231739572450155,59 0.1208650973866179e-2,-0.5395239384953e-5};60 int j;61 62 y=x=xx;63 tmp=x+5.5;64 tmp -= (x+0.5)*log(tmp);65 ser=1.000000000190015;66 for (j=0;j<=5;j++) ser += cof[j]/++y;67 return -tmp+log(2.5066282746310005*ser/x);68 }69 70 static double71 LogNormal_distr(double sig, double mu, double pt)72 {73 double retval,pi;74 75 pi = 4.0*atan(1.0);76 retval = (1.0/ (sig*pt*sqrt(2.0*pi)) )*exp( -0.5*(log(pt) - mu)*(log(pt) - mu)/sig/sig );77 return(retval);78 }79 80 static double81 Gauss_distr(double sig, double avg, double pt)82 {83 double retval,Pi;84 85 Pi = 4.0*atan(1.0);86 retval = (1.0/ (sig*sqrt(2.0*Pi)) )*exp(-(avg-pt)*(avg-pt)/sig/sig/2.0);87 return(retval);88 }89 90 static double SchulzPoint(double x, double avg, double zz) {91 double dr;92 dr = zz*log(x) - gammln(zz+1.0)+(zz+1.0)*log((zz+1.0)/avg)-(x/avg*(zz+1.0));93 return (exp(dr));94 };95
Note: See TracChangeset
for help on using the changeset viewer.