Changeset becded3 in sasmodels for sasmodels/models/capped_cylinder.c

Timestamp:

Oct 17, 2017 6:21:27 PM (7 years ago)

Author:

Paul Kienzle <pkienzle@…>

Branches:

master, core_shell_microgels, magnetic_model, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests

Children:

8698a0d

Parents:

2a602c7

Message:

update oriented models to new interface (which will be in the next commit)

File:

• sasmodels/models/capped_cylinder.c (modified) (5 diffs)

sasmodels/models/capped_cylinder.c

@@ -1 @@
-double form_volume(double radius, double radius_cap, double length);
-double Iq(double q, double sld, double solvent_sld,
-    double radius, double radius_cap, double length);
-double Iqxy(double qx, double qy, double sld, double solvent_sld,
-    double radius, double radius_cap, double length, double theta, double phi);
-
 #define INVALID(v) (v.radius_cap < v.radius)
 
@@ -17 +11 @@
 static double
 _cap_kernel(double qab, double qc, double h, double radius_cap,
-            double half_length)
+    double half_length)
 {
     // translate a point in [-1,1] to a point in [lower,upper]
@@ -60 +54 @@
 }
 
-double form_volume(double radius, double radius_cap, double length)
+static double
+form_volume(double radius, double radius_cap, double length)
 {
     // cap radius should never be less than radius when this is called
@@ -89 +84 @@
 }
 
-double Iq(double q, double sld, double solvent_sld,
-          double radius, double radius_cap, double length)
+static double
+Iq(double q, double sld, double solvent_sld,
+    double radius, double radius_cap, double length)
 {
     const double h = sqrt(radius_cap*radius_cap - radius*radius);
@@ -118 +114 @@
 
 
-double Iqxy(double qx, double qy,
+static double
+Iqxy(double qab, double qc,
     double sld, double solvent_sld, double radius,
-    double radius_cap, double length,
-    double theta, double phi)
+    double radius_cap, double length)
 {
-    double q, sin_alpha, cos_alpha;
-    ORIENT_SYMMETRIC(qx, qy, theta, phi, q, sin_alpha, cos_alpha);
-    const double qab = q*sin_alpha;
-    const double qc = q*cos_alpha;
-
     const double h = sqrt(radius_cap*radius_cap - radius*radius);
     const double Aq = _fq(qab, qc, h, radius_cap, radius, 0.5*length);