Changeset f12357f in sasmodels for sasmodels/models/lib/Si.c

Timestamp:

Feb 3, 2016 6:49:57 AM (8 years ago)

Author:

krzywon

Branches:

master, core_shell_microgels, costrafo411, magnetic_model, release_v0.94, release_v0.95, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests

Children:

a17fe40

Parents:

8696a87

Message:

Modified the pearl_necklace model for better double precision accuracy
and added unit tests.

File:

• sasmodels/models/lib/Si.c (modified) (2 diffs)

sasmodels/models/lib/Si.c

@@ -13 +13 @@
 
                 // Explicitly writing factorial values triples the speed of the calculation
-                out_cos = 1/x - 2/pow(x,3) + 24/pow(x,5) - 720/pow(x,7) + 40320/pow(x,9);
-                out_sin = 1/x - 6/pow(x,4) + 120/pow(x,6) - 5040/pow(x,8) + 362880/pow(x,10);
+                out_cos = 1./x - 2./pow(x,3) + 24./pow(x,5) - 720./pow(x,7);
+                out_sin = 1./pow(x,2) - 6./pow(x,4) + 120./pow(x,6) - 5040./pow(x,8);
 
                 out -= cos(x) * out_cos;
@@ -22 +22 @@
 
         // Explicitly writing factorial values triples the speed of the calculation
-        out = x - pow(x, 3)/18 + pow(x,5)/600 - pow(x,7)/35280 + pow(x,9)/3265920;
+        out = x - pow(x, 3)/18. + pow(x,5)/600. - pow(x,7)/35280. + pow(x,9)/3265920. - pow(x,11)/439084800.;
 
-        //printf ("Si=%g %g\n", x, out);
         return out;
 }