sphere.py [71b751d:d277229] in sasmodels

File:

sasmodels/models/sphere.py

-                      r71b751d
+                      rd277229
+             ]
 source = ["lib/sas_3j1x_x.c"]
+source = ["lib/sas_3j1x_x.c","sphere.c"]
 have_Fq = True
+effective_radius_type = ["radius"]
+c_code = """
+static double form_volume(double radius)
+{
+    return M_4PI_3*cube(radius);
+}
+static void Fq(double q, double *f1, double *f2, double sld, double sld_solvent, double radius)
+{
+    const double bes = sas_3j1x_x(q*radius);
+    const double contrast = (sld - sld_solvent);
+    const double form = contrast * form_volume(radius) * bes;
+    *f1 = 1.0e-2*form;
+    *f2 = 1.0e-4*form*form;
+}
+"""
+def ER(radius):
+    """
+    Return equivalent radius (ER)
+    """
+    return radius
+#def ER(radius):
+#    """
+#    Return equivalent radius (ER)
+#    """
+#    return radius
 # VR defaults to 1.0

Note: See TracChangeset for help on using the changeset viewer.