script for generating data from integrator program

1 files changed, 68 insertions, 0 deletions

diff --git a/run.sh b/run.sh
new file mode 100755
index 0000000..ca72d9c
--- /dev/null
+++ b/run.sh
@@ -0,0 +1,68 @@
+#!/bin/bash
+
+# Copyright 2019 Wojciech Kosior
+
+# This is free and unencumbered software released into the public domain.
+
+# Anyone is free to copy, modify, publish, use, compile, sell, or
+# distribute this software, either in source code form or as a compiled
+# binary, for any purpose, commercial or non-commercial, and by any
+# means.
+
+# In jurisdictions that recognize copyright laws, the author or authors
+# of this software dedicate any and all copyright interest in the
+# software to the public domain. We make this dedication for the benefit
+# of the public at large and to the detriment of our heirs and
+# successors. We intend this dedication to be an overt act of
+# relinquishment in perpetuity of all present and future rights to this
+# software under copyright law.
+
+# THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
+# EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
+# MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.
+# IN NO EVENT SHALL THE AUTHORS BE LIABLE FOR ANY CLAIM, DAMAGES OR
+# OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,
+# ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR
+# OTHER DEALINGS IN THE SOFTWARE.
+
+# For more information, please refer to <http://unlicense.org/>
+
+if [ -z "$1" ]; then echo "please pass number of images"; exit; fi
+
+COMMAND="cafrun -np $1 ./integrator"
+
+QUADRATURE_TYPES="newton-cotes gauss"
+
+FUNCTIONS="exp sin poly"
+
+# for simplicity always use the same
+IBEG="-0.5"
+IEND="1"
+
+SUBINTERVALS="1 4 16 64 128"
+
+valid_poly_orders() {
+    if [ "$1" = newton-cotes ]; then seq 0 2; fi
+    if [ "$1" = gauss ];        then seq 1 3; fi
+}
+
+echo "\
+# results of numerical integration of three different functions
+# (exp, sin and some 10 degree polynomial) on [$IBEG; $IEND] interval
+#
+#FUNC ___QUAD_TYPE POLY_ORDER SUBINTERVALS _______COMPUTED_VALUE __________VALUE_ERROR"
+
+
+for FUNC in $FUNCTIONS; do
+    echo -n "# function $FUNC;  exact integral:"
+    $COMMAND analytical $FUNC $IBEG $IEND 
+    
+    for TYPE in $QUADRATURE_TYPES; do
+	for N in $SUBINTERVALS; do
+	    for P in `valid_poly_orders $TYPE`; do
+		printf "%5s %12s %10s %12s " $FUNC $TYPE $P $N
+		$COMMAND $TYPE $FUNC $IBEG $IEND $P $N
+	    done
+	done
+    done
+done