aboutsummaryrefslogtreecommitdiff
path: root/charf.c
diff options
context:
space:
mode:
Diffstat (limited to 'charf.c')
-rw-r--r--charf.c16
1 files changed, 13 insertions, 3 deletions
diff --git a/charf.c b/charf.c
index 56e7b74..d1dbd33 100644
--- a/charf.c
+++ b/charf.c
@@ -40,13 +40,17 @@ Floating point without error bounds are the default.*/
typedef unsigned long index_t;
-/*maximum length of the circle*/
+/*maximum length of the circle
+Adjust as desired.*/
#define N 36
-/*maximal order of derivative*/
+/*maximal order of derivative
+Adjust as desired.*/
#define K 64
-/*maximal number of exponents*/
+/*maximal number of exponents
+Adjust as desired.
+Which exponents are considered is set using the array exponents in the beginning of the main routine.*/
#define P 5
/*Given function f on domain [0,D-1], compute derivative of order k and store result in df.*/
@@ -523,6 +527,12 @@ void* handle_threads(void* arguments){
}
int main() {
+ /*
+ Array of exponents to be considered.
+ Default: 1,2,4,8,∞.
+ Adjust as desired.
+ Note, that the macro P defined in the beginning of the file has to be adjusted in case fewer or more exponents are to be considered.
+ */
EXPTYPE exponents[P];
exponents[0] = int_to_exptype(1);
exponents[1] = int_to_exptype(2);