cwpearson/astaroth
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Errors with forcing now down from 3 to 1 after switching from fast & inaccurate trig functions to more accurate ones

Browse Source
This commit is contained in:
jpekkila
2019-08-06 19:29:40 +03:00
parent d7e26e8f21
commit 6b53eb31ef
1 changed files with 6 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

11
src/core/kernels/integration.cuh
View File

@@ -59,17 +59,18 @@ create_rotz(const AcReal radians)
} }
#if AC_DOUBLE_PRECISION == 0 #if AC_DOUBLE_PRECISION == 0
/*
// Fast but inaccurate
#define sin __sinf #define sin __sinf
#define cos __cosf #define cos __cosf
#define exp __expf #define exp __expf
*/
#define sin sinf
#define cos cosf
#define exp expf
#define rsqrt rsqrtf // hardware reciprocal sqrt #define rsqrt rsqrtf // hardware reciprocal sqrt
#endif // AC_DOUBLE_PRECISION == 0 #endif // AC_DOUBLE_PRECISION == 0
/*
typedef struct {
int i, j, k;
} int3;*/
/* /*
* ============================================================================= * =============================================================================
* Level 0 (Input Assembly Stage) * Level 0 (Input Assembly Stage)