Changeset 0ba3b08 in sasview for sansmodels/src/c_models/fuzzysphere.cpp

Timestamp:

Jan 5, 2012 10:16:29 AM (13 years ago)

Author:

Mathieu Doucet <doucetm@…>

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:

011e0e4

Parents:

bbbed8c

Message:

refactored bunch of models

File:

• sansmodels/src/c_models/fuzzysphere.cpp (modified) (1 diff)

sansmodels/src/c_models/fuzzysphere.cpp

@@ -15 +15 @@
 
 #include <math.h>
-#include "models.hh"
 #include "parameters.hh"
 #include <stdio.h>
+#include <stdlib.h>
 using namespace std;
+#include "fuzzysphere.h"
 
 extern "C" {
-        #include "fuzzysphere.h"
+#include "libSphere.h"
+}
+
+// scattering from a uniform sphere w/ fuzzy surface
+// Modified from FuzzySpheres in libigor/libSphere.c without polydispersion: JHC
+static double fuzzysphere_kernel(double dp[], double q){
+  double pi,x,xr;
+  double radius,sldSph,sldSolv,scale,bkg,delrho,fuzziness,f2,bes,vol,f;   //my local names
+
+  pi = 4.0*atan(1.0);
+  x= q;
+  scale = dp[0];
+  radius = dp[1];
+  fuzziness = dp[2];
+  sldSph = dp[3];
+  sldSolv = dp[4];
+  bkg = dp[5];
+  delrho=sldSph-sldSolv;
+
+  xr = x*radius;
+  //handle xr==0 separately
+  if(xr == 0.0){
+    bes = 1.0;
+  }else{
+    bes = 3.0*(sin(xr)-xr*cos(xr))/(xr*xr*xr);
+  }
+  vol = 4.0*pi/3.0*radius*radius*radius;
+  f = vol*bes*delrho;   // [=] A
+  f *= exp(-0.5*fuzziness*fuzziness*x*x);
+  // normalize to single particle volume, convert to 1/cm
+  f2 = f * f / vol * 1.0e8;   // [=] 1/cm
+
+  f2 *= scale;
+  f2 += bkg;
+
+    return(f2); //scale, and add in the background
 }