From 6b53eb31ef95a8b9612627344263ddc684da3124 Mon Sep 17 00:00:00 2001
From: jpekkila <johannes.pekkila@aalto.fi>
Date: Tue, 6 Aug 2019 19:29:40 +0300
Subject: [PATCH] Errors with forcing now down from 3 to 1 after switching from
 fast & inaccurate trig functions to more accurate ones

---
 src/core/kernels/integration.cuh | 11 ++++++-----
 1 file changed, 6 insertions(+), 5 deletions(-)

diff --git a/src/core/kernels/integration.cuh b/src/core/kernels/integration.cuh
index 3503372..fec3b95 100644
--- a/src/core/kernels/integration.cuh
+++ b/src/core/kernels/integration.cuh
@@ -59,17 +59,18 @@ create_rotz(const AcReal radians)
 }
 
 #if AC_DOUBLE_PRECISION == 0
+/*
+// Fast but inaccurate
 #define sin __sinf
 #define cos __cosf
 #define exp __expf
+*/
+#define sin sinf
+#define cos cosf
+#define exp expf
 #define rsqrt rsqrtf // hardware reciprocal sqrt
 #endif               // AC_DOUBLE_PRECISION == 0
 
-/*
-typedef struct {
-    int i, j, k;
-} int3;*/
-
 /*
  * =============================================================================
  * Level 0 (Input Assembly Stage)