1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
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
|