Changes in / [39c2de9:0168844] in sasmodels

File:

• sasmodels/kernel_iq.c (modified) (2 diffs)

sasmodels/kernel_iq.c

@@ -84 +84 @@
   out_spin = clip(out_spin, 0.0, 1.0);
   // Previous version of this function took the square root of the weights,
-  // under the assumption that
+  // under the assumption that 
   //
   //     w*I(q, rho1, rho2, ...) = I(q, sqrt(w)*rho1, sqrt(w)*rho2, ...)
@@ -188 +188 @@
     QACRotation *rotation,
     double qx, double qy,
-    double *qab_out, double *qc_out)
+    double *qa_out, double *qc_out)
 {
+    const double dqc = rotation->R31*qx + rotation->R32*qy;
     // Indirect calculation of qab, from qab^2 = |q|^2 - qc^2
-    const double dqc = rotation->R31*qx + rotation->R32*qy;
-    const double dqab_sq = -dqc*dqc + qx*qx + qy*qy;
-    //*qab_out = sqrt(fabs(dqab_sq));
-    *qab_out = dqab_sq > 0.0 ? sqrt(dqab_sq) : 0.0;
+    const double dqa = sqrt(-dqc*dqc + qx*qx + qy*qy);
+
+    *qa_out = dqa;
     *qc_out = dqc;
 }