Write some instructions.

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);