Changeset 99658f6 in sasmodels for sasmodels/models/capped_cylinder.c

Timestamp:

Nov 6, 2018 2:10:43 PM (5 years ago)

Author:

grethevj

Branches:

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

Children:

cf3d0ce

Parents:

5024a56

Message:

updated ER functions including cylinder excluded volume, to match 4.x

File:

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

sasmodels/models/capped_cylinder.c

@@ -85 +85 @@
 
 static double
+radius_from_excluded_volume(double radius, double radius_cap, double length)
+{
+    const double hc = radius_cap - sqrt(radius_cap*radius_cap - radius*radius);
+    const double length_tot = length + 2.0*hc;
+    return 0.5*cbrt(0.75*radius*(2.0*radius*length_tot + (radius + length_tot)*(M_PI*radius + length_tot)));
+}
+
+static double
 radius_from_volume(double radius, double radius_cap, double length)
 {
@@ -103 +111 @@
     switch (mode) {
     default:
-    case 1: // equivalent sphere
+    case 1: // equivalent cylinder excluded volume
+        return radius_from_excluded_volume(radius, radius_cap, length);
+    case 2: // equivalent volume sphere
         return radius_from_volume(radius, radius_cap, length);
-    case 2: // radius
+    case 3: // radius
         return radius;
-    case 3: // half length
+    case 4: // half length
         return 0.5*length;
-    case 4: // half total length
+    case 5: // half total length
         return radius_from_totallength(radius, radius_cap,length);
     }