Changeset 2cc633b in sasview for sansmodels/src/sans/models/c_models

Timestamp:

Aug 3, 2009 7:44:15 PM (15 years ago)

Author:

Jae Cho <jhjcho@…>

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:

42f193a

Parents:

0f3fefe

Message:

FIXED A BUG

Location:

sansmodels/src/sans/models/c_models

Files:

• CMultiShellModel.cpp (modified) (4 diffs)
• flexiblecylinder.cpp (modified) (3 diffs)
• multishell.cpp (modified) (2 diffs)

sansmodels/src/sans/models/c_models/CMultiShellModel.cpp

@@ -97 +97 @@
         self->model->core_radius.dispersion->accept_as_source(visitor, self->model->core_radius.dispersion, disp_dict);
         PyDict_SetItemString(self->dispersion, "core_radius", disp_dict);
+        disp_dict = PyDict_New();
+        self->model->s_thickness.dispersion->accept_as_source(visitor, self->model->s_thickness.dispersion, disp_dict);
+        PyDict_SetItemString(self->dispersion, "s_thickness", disp_dict);
+        disp_dict = PyDict_New();
+        self->model->w_thickness.dispersion->accept_as_source(visitor, self->model->w_thickness.dispersion, disp_dict);
+        PyDict_SetItemString(self->dispersion, "w_thickness", disp_dict);
 
 
@@ -162 +168 @@
     disp_dict = PyDict_GetItemString(self->dispersion, "core_radius");
     self->model->core_radius.dispersion->accept_as_destination(visitor, self->model->core_radius.dispersion, disp_dict);
+    disp_dict = PyDict_GetItemString(self->dispersion, "s_thickness");
+    self->model->s_thickness.dispersion->accept_as_destination(visitor, self->model->s_thickness.dispersion, disp_dict);
+    disp_dict = PyDict_GetItemString(self->dispersion, "w_thickness");
+    self->model->w_thickness.dispersion->accept_as_destination(visitor, self->model->w_thickness.dispersion, disp_dict);
 
         
@@ -226 +236 @@
     disp_dict = PyDict_GetItemString(self->dispersion, "core_radius");
     self->model->core_radius.dispersion->accept_as_destination(visitor, self->model->core_radius.dispersion, disp_dict);
+    disp_dict = PyDict_GetItemString(self->dispersion, "s_thickness");
+    self->model->s_thickness.dispersion->accept_as_destination(visitor, self->model->s_thickness.dispersion, disp_dict);
+    disp_dict = PyDict_GetItemString(self->dispersion, "w_thickness");
+    self->model->w_thickness.dispersion->accept_as_destination(visitor, self->model->w_thickness.dispersion, disp_dict);
 
         
@@ -283 +297 @@
     if (!strcmp(par_name, "core_radius")) {
         self->model->core_radius.dispersion = dispersion;
+    } else    if (!strcmp(par_name, "s_thickness")) {
+        self->model->s_thickness.dispersion = dispersion;
+    } else    if (!strcmp(par_name, "w_thickness")) {
+        self->model->w_thickness.dispersion = dispersion;
     } else {
             PyErr_SetString(CMultiShellModelError,

sansmodels/src/sans/models/c_models/flexiblecylinder.cpp

@@ -52 +52 @@
  */
 double FlexibleCylinderModel :: operator()(double q) {
-        double dp[5];
+        double dp[6];
 
         // Fill parameter array for IGOR library
@@ -64 +64 @@
 
         // Get the dispersion points for the length
-        vector<WeightPoint> weights_length;
-        length.get_weights(weights_length);
+        vector<WeightPoint> weights_len;
+        length.get_weights(weights_len);
 
         // Get the dispersion points for the kuhn_length
-        vector<WeightPoint> weights_kuhn_length;
-        kuhn_length.get_weights(weights_kuhn_length);
+        vector<WeightPoint> weights_kuhn;
+        kuhn_length.get_weights(weights_kuhn);
 
         // Get the dispersion points for the radius
-        vector<WeightPoint> weights_radius;
-        radius.get_weights(weights_radius);
+        vector<WeightPoint> weights_rad;
+        radius.get_weights(weights_rad);
 
         // Perform the computation, with all weight points
@@ -80 +80 @@
 
         // Loop over semi axis A weight points
-        for(int i=0; i< (int)weights_length.size(); i++) {
-                dp[1] = weights_length[i].value;
+        for(int i=0; i< (int)weights_len.size(); i++) {
+                dp[1] = weights_len[i].value;
 
                 // Loop over semi axis B weight points
-                for(int j=0; j< (int)weights_kuhn_length.size(); j++) {
-                        dp[2] = weights_kuhn_length[j].value;
+                for(int j=0; j< (int)weights_kuhn.size(); j++) {
+                        dp[2] = weights_kuhn[j].value;
 
                         // Loop over semi axis C weight points
-                        for(int k=0; k< (int)weights_radius.size(); k++) {
-                                dp[3] = weights_radius[k].value;
+                        for(int k=0; k< (int)weights_rad.size(); k++) {
+                                dp[3] = weights_rad[k].value;
 
-                                sum += weights_length[i].weight
-                                        * weights_kuhn_length[j].weight * weights_radius[k].weight * FlexExclVolCyl(dp, q);
-                                norm += weights_length[i].weight
-                                        * weights_kuhn_length[j]*weights_radius[k].weight;
+                                sum += weights_len[i].weight
+                                        * weights_kuhn[j].weight*weights_rad[k].weight * FlexExclVolCyl(dp, q);
+                                norm += weights_len[i].weight
+                                        * weights_kuhn[j].weight*weights_rad[k].weight;
+                        }
                 }
         }

sansmodels/src/sans/models/c_models/multishell.cpp

@@ -66 +66 @@
         core_radius.get_weights(weights_core_radius);
 
+        // Get the dispersion points for the s_thickness
+        vector<WeightPoint> weights_s_thickness;
+        s_thickness.get_weights(weights_s_thickness);
+
+        // Get the dispersion points for the w_thickness
+        vector<WeightPoint> weights_w_thickness;
+        w_thickness.get_weights(weights_w_thickness);
+
         // Perform the computation, with all weight points
         double sum = 0.0;
@@ -73 +81 @@
         for(int i=0; i< (int)weights_core_radius.size(); i++) {
                 dp[1] = weights_core_radius[i].value;
+                for(int j=0; j< (int)weights_s_thickness.size(); j++){
+                        dp[2] = weights_s_thickness[j].value;
+                        for(int k=0; k< (int)weights_w_thickness.size(); k++){
+                                dp[3] = weights_w_thickness[k].value;
 
-                sum += weights_core_radius[i].weight
-                        * MultiShell(dp, q);
-                norm += weights_core_radius[i].weight;
+                                sum += weights_core_radius[i].weight*weights_s_thickness[j].weight
+                                        *weights_w_thickness[k].weight* MultiShell(dp, q);
+                                norm += weights_core_radius[i].weight*weights_s_thickness[j].weight
+                                        *weights_w_thickness[k].weight;
+                        }
+                }
         }
         return sum/norm + background();