Changeset 56b0435 in sasmodels

Timestamp:

Dec 20, 2016 3:07:15 PM (7 years ago)

Author:

Paul Kienzle <pkienzle@…>

Branches:

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

Children:

a0168e8

Parents:

abb6e4f

Message:

round to nearest Nlayers rather than truncate

File:

• sasmodels/models/lamellar_stack_caille.c (modified) (3 diffs)

sasmodels/models/lamellar_stack_caille.c

@@ -3 +3 @@
 */
 
-double Iq(double qval,
-      double del,
-      double Nlayers, 
-      double dd,
-          double Cp, 
-      double sld,
-      double solvent_sld);
-
-double Iq(double qval,
-      double del,
-      double Nlayers, 
-      double dd,
-          double Cp, 
-      double sld,
-      double solvent_sld)
+static double
+Iq(double qval,
+   double del,
+   double fp_Nlayers,
+   double dd,
+   double Cp,
+   double sld,
+   double solvent_sld)
 {
-  double contr,NN;   //local variables of coefficient wave
+  int Nlayers = (int)(fp_Nlayers+0.5);    //cast to an integer for the loop
+  double contr;   //local variables of coefficient wave
   double inten,Pq,Sq,alpha,temp,t2;
   //double dQ, dQDefault, t1, t3;
-  int ii,NNint;
   // from wikipedia 0.577215664901532860606512090082402431042159335
   const double Euler = 0.577215664901533;   // Euler's constant, increased sig figs for new models Feb 2015
@@ -28 +21 @@
   //dQ = dQDefault; // REMOVED UNUSED dQ calculations for new models Feb 2015
 
-  NN = trunc(Nlayers);    //be sure that NN is an integer
-  
   contr = sld - solvent_sld;
 
   Pq = 2.0*contr*contr/qval/qval*(1.0-cos(qval*del));
 
-  NNint = (int)NN;    //cast to an integer for the loop
-  ii=0;
   Sq = 0.0;
   // the vital "=" in ii<=  added March 2015
-  for(ii=1;ii<=(NNint-1);ii+=1) {
+  for (int ii=1; ii<=Nlayers-1; ii++) {
 
     //fii = (double)ii;   //do I really need to do this? - unused variable, removed 18Feb2015
@@ -48 +37 @@
     //t3 = dQ*dQ*dd*dd*ii*ii;
 
-    temp = 1.0-ii/NN;
+    temp = 1.0 - (double)ii / (double)Nlayers;
     //temp *= cos(dd*qval*ii/(1.0+t1));
     temp *= cos(dd*qval*ii);