Changeset 34c2649 in sasview for sansmodels/src/sans/models/libigor/libSphere.c
- Timestamp:
- Nov 21, 2011 4:52:25 PM (13 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:
- 1810613
- Parents:
- afea8fe
- File:
-
- 1 edited
-
sansmodels/src/sans/models/libigor/libSphere.c (modified) (13 diffs)
Legend:
- Unmodified
- Added
- Removed
-
sansmodels/src/sans/models/libigor/libSphere.c
r3a084629 r34c2649 653 653 } 654 654 655 /* 655 656 static double 656 657 LogNormal_distr(double sig, double mu, double pt) … … 672 673 return(retval); 673 674 } 675 */ 674 676 675 677 // scattering from a core shell sphere with a (Schulz) polydisperse core and constant ratio (shell thickness)/(core radius) … … 1220 1222 } 1221 1223 1224 /* 1222 1225 static double 1223 1226 SchulzPoint(double x, double avg, double zz) { … … 1245 1248 return -tmp+log(2.5066282746310005*ser/x); 1246 1249 } 1250 */ 1247 1251 1248 1252 double … … 1262 1266 // variables are: 1263 1267 //[0] scale factor 1264 //[1] radius of core [ ]1265 //[2] SLD of the core [ -2]1266 //[3] thickness of the shell [ ]1268 //[1] radius of core [ï¿œ] 1269 //[2] SLD of the core [ï¿œ-2] 1270 //[3] thickness of the shell [ï¿œ] 1267 1271 //[4] SLD of the shell 1268 1272 //[5] SLD of the solvent … … 1305 1309 f += vol*bes*contr; 1306 1310 1307 // normalize to particle volume and rescale from [ -1] to [cm-1]1311 // normalize to particle volume and rescale from [ï¿œ-1] to [cm-1] 1308 1312 f2 = f*f/vol*1.0e8; 1309 1313 … … 1321 1325 // variables are: 1322 1326 //[0] scale factor 1323 //[1] radius of core [ ]1324 //[2] SLD of the core [ -2]1325 //[3] thickness of shell 1 [ ]1327 //[1] radius of core [ï¿œ] 1328 //[2] SLD of the core [ï¿œ-2] 1329 //[3] thickness of shell 1 [ï¿œ] 1326 1330 //[4] SLD of shell 1 1327 //[5] thickness of shell 2 [ ]1331 //[5] thickness of shell 2 [ï¿œ] 1328 1332 //[6] SLD of shell 2 1329 1333 //[7] SLD of the solvent … … 1379 1383 1380 1384 1381 // normalize to particle volume and rescale from [ -1] to [cm-1]1385 // normalize to particle volume and rescale from [ï¿œ-1] to [cm-1] 1382 1386 f2 = f*f/vol*1.0e8; 1383 1387 … … 1395 1399 // variables are: 1396 1400 //[0] scale factor 1397 //[1] radius of core [ ]1398 //[2] SLD of the core [ -2]1399 //[3] thickness of shell 1 [ ]1401 //[1] radius of core [ï¿œ] 1402 //[2] SLD of the core [ï¿œ-2] 1403 //[3] thickness of shell 1 [ï¿œ] 1400 1404 //[4] SLD of shell 1 1401 //[5] thickness of shell 2 [ ]1405 //[5] thickness of shell 2 [ï¿œ] 1402 1406 //[6] SLD of shell 2 1403 1407 //[7] thickness of shell 3 … … 1467 1471 f += vol*bes*contr; 1468 1472 1469 // normalize to particle volume and rescale from [ -1] to [cm-1]1473 // normalize to particle volume and rescale from [ï¿œ-1] to [cm-1] 1470 1474 f2 = f*f/vol*1.0e8; 1471 1475 … … 1483 1487 // variables are: 1484 1488 //[0] scale factor 1485 //[1] radius of core [ ]1486 //[2] SLD of the core [ -2]1487 //[3] thickness of shell 1 [ ]1489 //[1] radius of core [ï¿œ] 1490 //[2] SLD of the core [ï¿œ-2] 1491 //[3] thickness of shell 1 [ï¿œ] 1488 1492 //[4] SLD of shell 1 1489 //[5] thickness of shell 2 [ ]1493 //[5] thickness of shell 2 [ï¿œ] 1490 1494 //[6] SLD of shell 2 1491 1495 //[7] thickness of shell 3 … … 1570 1574 1571 1575 1572 // normalize to particle volume and rescale from [ -1] to [cm-1]1576 // normalize to particle volume and rescale from [ï¿œ-1] to [cm-1] 1573 1577 f2 = f*f/vol*1.0e8; 1574 1578 … … 1997 2001 bes = 3.0*(sin(x*radius)-x*radius*cos(x*radius))/(x*x*x)/(radius*radius*radius); 1998 2002 vol = 4.0*pi/3.0*radius*radius*radius; 1999 f = vol*bes*delrho ; // [=] 2003 f = vol*bes*delrho ; // [=] ï¿œ 2000 2004 // normalize to single particle volume, convert to 1/cm 2001 2005 f2 = f * f / vol * 1.0e8; // [=] 1/cm
Note: See TracChangeset
for help on using the changeset viewer.
