Changeset 5efe850 in sasmodels for sasmodels/kernel_template.c

Timestamp:

May 20, 2016 12:53:41 PM (8 years ago)

Author:

Paul Kienzle <pkienzle@…>

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:

c98c772, 4bfbca2

Parents:

33af590

Message:

Use tinycc if available; support float32 models in tinycc

File:

• sasmodels/kernel_template.c (modified) (3 diffs)

sasmodels/kernel_template.c

@@ -12 +12 @@
 // Note: if using a C++ compiler, then define kernel as extern "C"
 #ifndef USE_OPENCL
+// Use SAS_DOUBLE to force the use of double even for float kernels
+#  define SAS_DOUBLE dou ## ble
 #  ifdef __cplusplus
       #include <cstdio>
@@ -38 +40 @@
      #if defined(__TINYC__)
          #include <math.h>
-         inline double trunc(double x) { return x>=0?floor(x):-floor(-x); }
-         inline double fmin(double x, double y) { return x>y ? y : x; }
-         inline double fmax(double x, double y) { return x<y ? y : x; }
          // TODO: test isnan
          inline double _isnan(double x) { return x != x; } // hope this doesn't optimize away!
          #undef isnan
          #define isnan(x) _isnan(x)
+         // Defeat the double->float conversion since we don't have tgmath
+         inline SAS_DOUBLE trunc(SAS_DOUBLE x) { return x>=0?floor(x):-floor(-x); }
+         inline SAS_DOUBLE fmin(SAS_DOUBLE x, SAS_DOUBLE y) { return x>y ? y : x; }
+         inline SAS_DOUBLE fmax(SAS_DOUBLE x, SAS_DOUBLE y) { return x<y ? y : x; }
          #define NEED_EXPM1
          #define NEED_TGAMMA
@@ -70 +73 @@
 
 #if defined(NEED_EXPM1)
-   static double expm1(double x) {
+   static SAS_DOUBLE expm1(SAS_DOUBLE x_in) {
+      double x = (double)x_in;  // go back to float for single precision kernels
       // Adapted from the cephes math library.
       // Copyright 1984 - 1992 by Stephen L. Moshier