From 9c3c4ad068c0322d5e7aff464ee732fd42af9c2a Mon Sep 17 00:00:00 2001
From: niklas <niklas.birk.nb@gmail.com>
Date: Wed, 22 Nov 2023 19:51:30 +0100
Subject: [PATCH] Initial

---
 uebung_01/exercise_1.py                       |  22 ++++++
 uebung_01/exercise_2a.py                      |  24 +++++++
 uebung_01/exercise_2b.py                      |  22 ++++++
 uebung_01/exercise_3.py                       |  68 ++++++++++++++++++
 uebung_02/main_matrix.c                       |  64 +++++++++++++++++
 uebung_02/main_openmp.c                       |  19 +++++
 uebung_02/matrix_mul                          | Bin 0 -> 16328 bytes
 uebung_02/openmp_hello                        | Bin 0 -> 16448 bytes
 uebung_04/exercise_01/exercise_01.c           |  17 +++++
 uebung_04/exercise_01/exercise_01.out         | Bin 0 -> 16208 bytes
 uebung_04/exercise_01/exercise_01.script      |  13 ++++
 .../uebung_04__exercise_01.e662620            |   0
 .../uebung_04__exercise_01.o662620            |   8 +++
 uebung_04/exercise_03/exercise_03.c           |  26 +++++++
 uebung_04/exercise_03/exercise_03.out         | Bin 0 -> 16320 bytes
 uebung_04/exercise_03/exercise_03.script      |  13 ++++
 uebung_04/exercise_05/exercise_05.c           |  33 +++++++++
 uebung_04/exercise_05/exercise_05.out         | Bin 0 -> 16392 bytes
 18 files changed, 329 insertions(+)
 create mode 100644 uebung_01/exercise_1.py
 create mode 100644 uebung_01/exercise_2a.py
 create mode 100644 uebung_01/exercise_2b.py
 create mode 100644 uebung_01/exercise_3.py
 create mode 100644 uebung_02/main_matrix.c
 create mode 100644 uebung_02/main_openmp.c
 create mode 100755 uebung_02/matrix_mul
 create mode 100755 uebung_02/openmp_hello
 create mode 100644 uebung_04/exercise_01/exercise_01.c
 create mode 100755 uebung_04/exercise_01/exercise_01.out
 create mode 100644 uebung_04/exercise_01/exercise_01.script
 create mode 100644 uebung_04/exercise_01/uebung_04__exercise_01.e662620
 create mode 100644 uebung_04/exercise_01/uebung_04__exercise_01.o662620
 create mode 100644 uebung_04/exercise_03/exercise_03.c
 create mode 100755 uebung_04/exercise_03/exercise_03.out
 create mode 100644 uebung_04/exercise_03/exercise_03.script
 create mode 100644 uebung_04/exercise_05/exercise_05.c
 create mode 100755 uebung_04/exercise_05/exercise_05.out

diff --git a/uebung_01/exercise_1.py b/uebung_01/exercise_1.py
new file mode 100644
index 0000000..5ee8612
--- /dev/null
+++ b/uebung_01/exercise_1.py
@@ -0,0 +1,22 @@
+import numpy as np
+
+for n in [3,5,10,100,1000]:
+    A = np.zeros((n,n))
+    x = np.zeros((n,1))
+    y_c = np.zeros((n,1))
+
+    for i in range(0,n):
+        x[i] = i+1
+        y_c[i] = (i+1) * n;
+
+        for j in range(0,n):
+            A[i,j] =(i+1) / (j+1)
+
+    y = A@x
+    
+    if n == 3 or n == 5:
+        print(x)
+        print(A)
+        print(y)
+
+    print("Fuer n =",n,":",np.linalg.norm(y - y_c))
diff --git a/uebung_01/exercise_2a.py b/uebung_01/exercise_2a.py
new file mode 100644
index 0000000..a9d3653
--- /dev/null
+++ b/uebung_01/exercise_2a.py
@@ -0,0 +1,24 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+# x = np.linspace(0, 2 * np.pi, 100)
+x_1 = np.arange(0, 2 * np.pi, 1)
+x_0_5 = np.arange(0, 2 * np.pi, 0.5)
+x_0_1 = np.arange(0, 2 * np.pi, 0.1)
+
+y_s_1 = np.sin(x_1)
+y_s_0_5 = np.sin(x_0_5)
+y_s_0_1 = np.sin(x_0_1)
+y_c = np.cos(x_1)
+
+plt.plot(x_1,y_s_1,"purple", label="sin(x) with stepsize 1")
+plt.plot(x_0_5,y_s_0_5,"pink", label="sin(x) with stepsize 0.5")
+plt.plot(x_0_1,y_s_0_1,"red", label="sin(x) with stepsize 0.1")
+plt.plot(x_1,y_c,"green", label="cos(x)")
+
+plt.title("A2.a) Sine and Cosine")
+plt.xlabel("x")
+plt.ylabel("y")
+plt.legend()
+
+plt.show()   
diff --git a/uebung_01/exercise_2b.py b/uebung_01/exercise_2b.py
new file mode 100644
index 0000000..ed7c67c
--- /dev/null
+++ b/uebung_01/exercise_2b.py
@@ -0,0 +1,22 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+plt.rcParams['text.usetex'] = True
+
+# x = np.linspace(0, 2 * np.pi, 100)
+x = np.arange(-2, 2, 0.1)
+
+f = np.exp(-x**2)
+g = np.sin(x**2)
+h = np.sin(1 / (x**3 + 9))
+
+plt.plot(x, f, "purple", label=r"$f(x) = \exp(-x^2)$")
+plt.plot(x, g, "red", label=r"$g(x) = \sin(x^2)$")
+plt.plot(x, h, "green", label=r"$h(x) = \frac{1}{x^3 + 9}$")
+
+plt.title("A2.b)")
+plt.xlabel("x")
+plt.ylabel("y")
+plt.legend()
+
+plt.show()   
diff --git a/uebung_01/exercise_3.py b/uebung_01/exercise_3.py
new file mode 100644
index 0000000..da1d328
--- /dev/null
+++ b/uebung_01/exercise_3.py
@@ -0,0 +1,68 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+def runge(x):
+    return 1 / (1 + 25 * x**2)
+
+def divided_differences(x, y):
+    n = len(y)
+    a = np.zeros((n,n))
+    
+    a[:,0] = y
+    
+    for i in range(1,n):
+        for j in range(1,i+1):
+            a[i,j] = (a[i,j-1] - a[i-1,j-1]) / (x[i] - x[i-j]) # (links daneben - links darüber) / (x[zeile] - x[zeile - spalte])
+    
+    return a
+    
+def newton_interpolation(a, data, x):
+    n = len(a)
+    
+    p = np.zeros(len(x))
+    
+    for i in range(1,n+1):
+        p = a[n-i] + (x - data[n-i]) * p # Horner Schema; '-' and '*' are overloaded for numpy arrays: so at least one argument 'data' or 'x' has to be numpy array 
+    
+    return p
+
+####################################################################################################################################################################
+n = 12
+
+x = np.linspace(-1, 1, 200)
+x_e = np.linspace(-1, 1, n) # equidistant grid points
+
+# Chebyshev grid points
+x_c = np.zeros(n)
+for i in range(0,n):
+    x_c[i] = np.cos((2 * i + 1) * np.pi / (2 * n))
+    
+f = runge(x)
+y_e = runge(x_e) # values for grid points for interpolation with equidistant grid points
+y_c = runge(x_c) # values for grid points for interpolation with Chebyshev grid points
+
+# Interpolation with equidistant grid points and evaluation of interpolated values at x
+a_e = np.diag(divided_differences(x_e, y_e))
+p_e = newton_interpolation(a_e, x_e, x)
+
+# Interpolation with Chebyshev grid points and evaluation of interpolated values at x
+a_c = np.diag(divided_differences(x_c, y_c))
+p_c = newton_interpolation(a_c, x_c, x)
+       
+       
+# Plotting of Runge and the two interpolated polynomials
+plt.rcParams['text.usetex'] = True
+
+plt.plot(x, f, label=r"$f(x) = \frac{1}{1 + 25 x^2}$")
+plt.plot(x, p_e, label=r"$p_n(x)$ equidistant")
+plt.plot(x, p_c, label=r"$p_n(x)$ Chebyshev")
+plt.plot([], [], ' ', label="$n = {}$".format(n))
+
+plt.title("A3)")
+plt.xlabel("x")
+plt.ylabel("y")
+plt.legend()
+
+plt.show()
+
+# braendel@math.tu-freiberg.de
diff --git a/uebung_02/main_matrix.c b/uebung_02/main_matrix.c
new file mode 100644
index 0000000..9380c05
--- /dev/null
+++ b/uebung_02/main_matrix.c
@@ -0,0 +1,64 @@
+#include <stdio.h>
+#include <stdlib.h>
+
+int compute(int N) {
+
+  int i_0, i_1;
+  // claim space
+  double **A = (double **)malloc(sizeof(double *) * N);
+  for (i_0 = 0; i_0 < N; i_0++) {
+    A[i_0] = (double *)malloc(sizeof(double) * N);
+  }
+  double *x = (double *)malloc(sizeof(double) * N);
+  double *y = (double *)malloc(sizeof(double) * N);
+  double *y_sol = (double *)malloc(sizeof(double) * N);
+
+  // fill matrix and vector
+  for (i_0 = 0; i_0 < N; i_0++) {
+    x[i_0] = i_0 + 1;
+    y_sol[i_0] = (i_0 + 1) * N;
+    for (i_1 = 0; i_1 < N; i_1++) {
+      A[i_0][i_1] = (i_0 + 1.0) / (i_1 + 1.0);
+    }
+  }
+
+  // for (i_0 = 0; i_0 < N; i_0++){
+  // 	printf("[ ");
+  // 	for(i_1 = 0; i_1 < N; i_1++){
+  // 		printf("%.3f ",A[i_0][i_1]);
+  // 	}
+  // 	printf("]");
+  // 	printf("[ %.1f ] = [ %.1f ][ %.1f ]\n",x[i_0],y[i_0],y_sol[i_0]);
+  // }
+
+  // compute matrix-vector product
+  for (i_0 = 0; i_0 < N; i_0++) {
+    y[i_0] = 0;
+    for (i_1 = 0; i_1 < N; i_1++) {
+      y[i_0] += A[i_0][i_1] * x[i_1];
+    }
+  }
+
+  // compare result to solution
+  double sum = 0.0;
+  for (i_0 = 0; i_0 < N; i_0++) {
+    sum += (y[i_0] - y_sol[i_0]) * (y[i_0] - y_sol[i_0]);
+  }
+
+  return sum;
+}
+
+int main(int argc, char *argv[]) {
+
+  int N;
+  if (argc != 2) {
+    N = 1000;
+  } else {
+    N = atoi(argv[1]); // Convert the argument to an integer
+  }
+
+  double err = compute(N);
+  printf("Error for N = %d: %.8f\n", N, err);
+
+  return 0;
+}
diff --git a/uebung_02/main_openmp.c b/uebung_02/main_openmp.c
new file mode 100644
index 0000000..92be7c1
--- /dev/null
+++ b/uebung_02/main_openmp.c
@@ -0,0 +1,19 @@
+#include <omp.h> // load Open-MP library
+#include <stdio.h>
+#include <stdlib.h>
+
+int main() {
+  int N = 5;
+  int num;
+
+#pragma omp parallel for private(num) // start Open-MP section
+  // here the program (1 thread) is split into multiple threads
+  // each thread has its own copy of i (loop) and num (private)
+  // all threads share N (by default, because defined outside)
+  for (int i = 0; i < N; i++) {
+    num = omp_get_thread_num();
+    printf("Thread %d does iteration %d\n", num, i);
+  }
+
+  return 0;
+}
diff --git a/uebung_02/matrix_mul b/uebung_02/matrix_mul
new file mode 100755
index 0000000000000000000000000000000000000000..f0a36d91830c9c674aceebc08977cca57d7e5e01
GIT binary patch
literal 16328
zcmeHOeQZ<L6~A@@1WIl5wL&_;v@nHsc_a{`z?xx)IJ{KhV<c=;1zw%lCT<-&wVz8;
zRg5ZYH4jqJNn5*5(f!dPP1Be(O{|^TL;(t0w?Qjauu2=-3JPm!srvzJDsSi9ch2MG
zVYgbhe<r;vJ)e7i=lssS_q_M)`~2?xdZeSx>+uLK72?YRwF3bY36Y_yi7FsMqDkBg
zp9{qTF&XqqiJA70DZpx$)82`)itq}M$eW9m67usVELd_6Q6jHg8Nm)gVJT=F$eV-)
z(w(xI`oWUzR(nAtLkY@HOKw63C@k$f?8?DfX%+{YC(2&ROFv2Gp^Ci2#5+tpmY1kK
zmYh$_2_8)}4yUXG8y1#MUN7;Svel#sme>ZC&U1R~v$VH~cr`&8MmXhNCM;NTdtU(_
zom%qvWWlF-GqqRlJZz=>v!wKjWvOI$W5cpktTvTQ=Lc#Bni^{x8<bo|SuPvKeg)yc
zYiiwwEds)fAVz5v-aV$y@z6G#?C4)O`R%73nKIh>;_uGyTKZZ?<@x;g4xkVA8#(BM
z1>^B0f;!gm+@X(|k4K?qVkl3=|ImlBTv-f;nI(+(0~RcSKVAlJ1iVz7`DJh%S1JGR
zl!=eyyA6&Til+4pM72aR9Ze-4jtfoe&nDA)LhI>$K=egZsZ0-GaPoSzoF2{UT3<Ao
z7CAkuXHo)TJcW&7T}OM%T5Y+qLTRw-%axT@eMN)N+Pl_kv3NGVGnv!l*{=0#Q<-$U
zE83mHR(JMg($s=xx|Vc=?t9@}dEh)&!2OBRELg)^o3B%SUUDkl!RrXiSEF7bI!XWK
zli!8Y?G>6yx4uvKKHf9KG#sMu3&ul)SA|SF<E)=vW8!!%Va~anwc+;bZ`_9SSdeqc
zhEL%jAg|bPZX3-bgrhHB_gs$KaQS^rhLbkjyzZIdN;d;;2HXs|8E`Y;X28vWn}Po%
z1MmA6{8JrzyGq?V>C7@A)W=42uW?cxI#M+z!x$^S1#r<=vKNGZp+Ns$v=lG=)i8|x
zQpbs7@xm!f#|dNc!YJ#&9+vD@A^D(39D1`uUuoCtFSZv>?NAGUQ-|KXw5cmnKT>}}
z9d53MW{f$r!1IHI;$Qd;?3Dm~6(H*Hq~$19UD4;j&^}Vkp-nZ;`WNm&zcG^FKUMm#
z^rBdLO)XqfN8eeaj$WxyJ;&5j*Yw#C;3x&CGR`La3t=40$MNrJMzhG@xkVjn{s~sq
z!a04qy0`gRP{yxbGmP;VB;(kmAA#=K4%`X+3%kIGfx=%6ZwYs8=?oV>?&!JPK6<{Q
zqi}ir4Dch5UGzV3JLFj{990Wrt%dQ9;n-1u=I6nz7An-Cv8s+$f6u>-*s~HF!|D&9
z+Fm$<)%_BExwY_*NMZa%%tFy0m|yJ^Vz|wd@V`4!^xpx;jLP?=5+?1ky0_y)aTxsp
ze-)fwFZ#E>cAc<wEa9*DZPCBuIFP50szW0lb=ArIpVgsbo``2$UG;kY$FjKzflff<
zm7_)fJ(Xh+?l`k96vGhY1LPv{6>j$g;DO2`MgO+Smw^QRyv}*E=>OvAz1zdv!rQ|;
zpfm5Dc<UnOv-5><ct#mAyt-2O%lP+U_}uBjpVZ-%l{X9V5?;H9@nihaj}1e%cJG+F
zcX!qE(*Ts+9$yER!q{6c0DxiT{bJoA0J_D^fSUm~18xS~47eF^GvH>x&48N$Hv|6%
z8Sub6rARiL$p#Z}YzQ_7m&8^FmncmMeEi{ZX^rRBiZ9_iS$y|$<z2(*2ibJVFrEgf
zf;<Fr0OVGX(f17FILKcV4WkMlt)Sq$Pb{8?Hw(`|z;o;LDOLNyhB&_Cp9|bCphN2v
z`ATdtR2$&f18s#s3SXejH}`J;O}nc0h&9!#?^v-|z8*yWHaMOJKYo{n?uhGf><7FK
z0V+6#eQ-Pp?a6ocZN9)0-WK27;ffYta9?HESM&I!u&-`taz&$Ony)783xcu57ijTS
zg>Qm>0{0a-o`ErTffNy6V92{`nrTJAEp7(f47eF^GvH>x&48N$Hv?`4+zhxG`0p~n
z>-%^;A3avdHG~*>I)|r+Ab7PUOWiC}&Lf@IA>T?mugR+>o!1ZIGZ_}Pe|XKvpnixJ
z{qfqoLzRf~T2vuUlbxR)e{PxuJUYpO&rew9n8Nxbnm9=7{CLfwN*?$;hlSS`-eRg^
zfEEKXewaGO{T?AdSWGDkx4(dN_WPG%jt}xGA$V6P$2&mr;0>+RIq&fJDE0fuewB0{
z2ZjHi310v03|L3!J4A95NvGeZ>AW6v-P*OQgEd>a^JzUFtglxZl)Bpbyj1IVFJG?I
zH7q6glOtlDD#d$-(+kgo2Q2$2Q9*LR%e0g~o+mG?l0E4d7xEAnw7wNoSysR&;(VVF
zEZU*s5!E6PG7&zHxXRiiJ}3CPH?2ZE0u_wI_Z{Qkl>T$9^|y>a2bEI&|5X{h+<1@3
z__IaFB93gza`Dkxs@+!s_lR4B^SXIMhMQ%ryO%0_yd(YRTK6C0m5``X>$VpG9%Lt*
z)FckyHX$1b)n>p0VxE<MwzdL3(>p`(_Y~uI13trhV?Vc39B2Id%lPjGJjg+9av1PZ
z<2oSY%oBV+Fy&{6BZGp`*Pi%InK*-h2jF^}8{!D2?2-NfYk!UW&kIa{B>#ZJ{|)j#
zL!FiXIf4IKiJLL)pH~5&;k_}>li&-rR6gecJ`;Lcp!-PvCk6o*mGN%`JZNu_>YD(+
ztpxr6;4{7O>*Tn9MQ2Y|uU9g#ls?*>)b!|10YFaAClX4JxW4RO)BCg@TzH=Y&sav=
znaXrWQ(8>VWOG_HKOlNCef_Dp9*-$?r9E)3KvIilv(Z5<p4PL2B9V>u#kE+zuWt}2
zHcbOheS&Kykw|5tG42;`CV(9UdNw%#z0|a}&EfSCEwZ6igG~Uf_iPBSZ(j>AY~Ii!
zDp}Rm%|ct(v9TrGp>1qy>x^`1UE!9F2r!)c6hhXvgiqOi5RXOmsE{LrDAY;qp(x}e
zo{j;8c5{1z&7R1od)Q&Z{snV;gw;B3gn$sajMf`X$8ftt`$h;BgR$px@fdBhaJYBp
za@3&Qabe$=VR_<i3r&Na6r3`9*Wnt#trd7a^Ez&^;9+66$td?^l-y`Baa)Mhq8tjI
zRUwqzU>`IA5;lWK*~>N9KoZaP3niV=<I2u-Ug^(f`r}!B&<1tqldu^m86zm%(q4-f
zt#t2==6Zz^8%#qJrqr{h<)L^sm&~LmC>q$Z@l+HERO?UaLXq86psMW5fQpq|yhkW{
zd;nCuiorW0&%6@vrR%mg2GLkEW0)7P$>2IP741txFp~}^ODK@(KDZ%D+_awo!uce8
z^}uhYS)BX&cs@yH>C1JLx5I&RIOdnzzlOe8c;04eLg0QRJKwpF4c{l&3x9tzhVuhE
zpXVnmgX~Ds$^RPQAHf`u`_J=amfSu~S8wQl6!<ttW<JlOElKTSA2H87exA=J0mCU3
z^Lbv*lIP*zA`8Z_mj^+?IWF^gUdxhyE4lyd$8tCDan8&-&zo6>s4-{%2^Pap!JIOm
z=kY8@$dTK3j{hm*w~!w{7qC=qe9W1X{|sQ56OM2FUeWIsKQB=j*;lTk{2?56evo)9
zL+E5H{By|iIk4IJ_mP1m`*8@?Ssrrm`^mtvj)^3h$MV+>{sA(u9Oe__;B(lQ9DJV7
zvE=7b_IKVtM~QDW6C_WPHKelO^T+z{p@R8ieofG{u)NDutxsqBI4`#IdxNGc%W;co
zR-Emh0~i+4TKoI>Ih*&LI{jH^`5tinHa<U(9I6En3j3Mbb)E5#!GPB;KY#MPpZ5{-
zeqp@-urQC$|4#wK@iU*F!w)V)Zy~mt!un)B%gI2t^F@eU>rBb*Z~*49J{`!2ar-<E
z4iMkzk8Rn@3~Ha;Efe$av=Fnxe7odp$L<=UfK7D4%m*X>7s37E@879_ao_Fyfd*51
Z&=H*5WnY|rm*QW*A2uk}4uV4!{|5UI4GI7N

literal 0
HcmV?d00001

diff --git a/uebung_02/openmp_hello b/uebung_02/openmp_hello
new file mode 100755
index 0000000000000000000000000000000000000000..00099eb1a04e743d029b09271e9149f4722d6f6a
GIT binary patch
literal 16448
zcmeHOYiu0V6~4Q+!xF;UAqhB1OBQfwDD}opOqu{8yN>OwkwenNDJmsQ#{2ND*oV72
z3l31~%7C&BR$ZxTLshk+0zpM3sMMnLB}@}uLaj*YkAxQ0AdsBWrok;Ofs*ZaX3iN;
zhFxf*{+ug~zk9y(xOdLn8Q<Bt_ba^veSV)$Noi2mDUzmJEGEK=vGrU5MAUk9KK`yy
z7pr-Y8%@sAM=XHWDo6ciVu1MuAo;aYr3QY&qLrZZkdXZ9m8o+qjiBQ5AipM3vFyf&
zxu1ffx7!PmGfK$(Yhw=?pa?qsC^!4sO|{y$(S_XJB%}C|NPY*{?*RJ=DzRsjd{R#6
zSkL)%V+V9p1l@jV_H$#8<)#Fw4Jtd|G!kdi-gfqD3z=b*8?Uu!B`EEE4t|uE|F+0w
z-ofqFyAGqAM?t=vYIQa<w61G)HqoBV<jbS&qwCkTuj>kz3gI<o7pPw$95kl3+^~x-
z%M@{(M&%!~WQj+*#jGd$N48vl%T1U3>Yh^@9(erv$)<-sb<L&UCmZTF`H&42iYJ$x
zWU9}{K{i%CzKWW~5gtkYHb+>md>e`jYA7jQs0Kd>u9p1|>+r*M_)mfl;c%C)fv6S#
zGvGsRW5|6tYT4fj{?Z!uO6wooswI-e<Zz~BB#VPvH)jj^<e)y3O)5<r&K2@n$<T|2
zrdg&n4YzFDx?LO5i+VPj%&J0eL>o>TnvpIh^@Nr$=WGgfrb~z(AJw%~Ca-5RcTi&^
z#Z2BvY4P-}5HP;t+;UFO<iVgD@moo&rSwcz!C)9I(U9vTAx^b*%Ir#4xI<BwI)PgT
z`nxx4Yr<>7SK0No;V!#AS6kaf3G!p?G~ij$0BpiBPthaO2ma?k8h`YRc<HjsF&yWq
zBg{qGNN!NCK&IkDgy(>ApOu$o&Q*x?I&;QSlAF{CmKSF-=hGAFJ(hp+jc;N2`c)$h
zJQR68%5!OIp~X+~^H2CO<^vH+FTBV{H(5O08ywmsrX9TVewuOc*8ODppLFoj25H`P
z@N^$ad4hS$r91;jchJGp{Viq6!JGe=Sn<4rr&v;6a_~s3SG){(8SpaTWx&gTmjN#W
zXPtq+H(&f-Z2T{Q*#4$ttCfn~J7xH*FU7_W2M(Drt5@9tdbYaqPC)YtMfPdZ%pQNO
zTCGl+GEKB*j~}sRnqbWyf7+I5LNt5)$F@upq1odPij3$vGkT9Q<3nEJlFxU^O)mL^
zF!gTE{~f2?`Pcr+ky~Pw*JI;voZLRx+d0+wLTqB=CUm>H<Yid?BNc95aSu9$1War^
z1tvDpG)u7IU1JGS^CZ=gnscku%`3*p?hqsFUoh>j`Zd87r(>0qv8PXLiaq^qL(KPV
z?8s?jF#_Dn0Rq+Ol-W=5r~Hj=B(*AExhpolafqt1%A3Z5*#3<v2s6h{SF1A#B<tCx
z-H?5|;cK@~{>N{Jk}A<%(ZO9i>2mI0G4ngPRi2r73=pjxt{i>fUmsMSsT_%J`2Ahi
zRo;p=KmJGGi@sMky!GXk0~771FgdRr>#4jxP&x4wUCx6z)XewbJbuWx;ngy7cR29H
z-O-z)yAMsIeLo}hTL-C8G}oGl?QdMU3RJAp(-!EQdakFfk!H2A38Ssiw54|oO71Rn
zb>>EdjqQhJB6g+aSDtS+h(=lu%><z4$bq?$NCCX!Wx&gTmjN#WUIx4jcp30A;AOze
zfR_O;18187`tCDme#Z)}OoS4JWGR%vXENQ$6!KuqPaMAG4cE};fBHW1@4r>6S>Vkl
ztJQA;v%nt#e+i^-g%7`7tsVw0o2^#q1I%L7src^Lp?sq)zU2$f2~09i-}d*w?{YHe
zImi5z)P^d3rrL*g_5f9|r7zfeee-#@2gcN<WmkWCZJYU-kK*2pV+wZi4T{W(H*g#T
zznus!)c^Z&97B6$;`)Lu6aMaC>-`PUVCcTaXmHhdQ%|tN|M~gBRbZmQ*6v_ScQ6n=
zPqo7Qc^ntu+-@S@STN)tm~ZMx^oo}OF9Ti%ybO35@G{_Kz{`M_0WSky2L6`}$a*(f
z-$viysF>f)fFjeHJ1Pq-uovLr`<+=1@Z!j&EX%sO<t)ouwq-2K`Z{_SqayluPFD*g
zKf#OHWR2Q%BeAmfRH>t^m$h)s+|lzaGxObvs!LcFU^jJu*R{!-yiqow_c$uD2JRwD
zR52!m-^Yvgq~B9)Cu`*B-Hpn{7TD+ae$_fZbk}snCr9&q0~`;ZDNR}Wh1W(?KFNtV
z#j>0ShyTBUz24guv4dZBBaGV__b`qyy4*fil=YKaHgCQ<v})H-Id7Chot<H<G;QxJ
zn_}nPYu1E2x~^dQqeId#XjE@k-Bxs6Xzl3zoVXAhECSZDr|ag&D$^ET_sQle#C*iX
zuVr5DBeCyC#iy33mWajhFjch5+NUm5avxY)rS3w7+LilK_<g4R5_`R>@Q<KUtN%~d
z;psOLivI%l-~3GiT)jNJ44n_ZZbU9Y#Zl_DOD47Y^BQ=cx=6Xl*PqRpi|qCIriecu
znD(vqeJcDSB(T;x><;iDF>>HTCXc65#6hAu0KP>nwXch4cY?1KKTo`Wp_1pE*c)~1
zzYcz(f00@h5=#qz0G`@)_vi6CabBzw=TG2Wq0I(PnD$GR+*i{2`^3|ku@gT|)hpmz
zFg{x&Vrk){;4gq?n(u4#w=$^a<H1DsGM>fjLhxjN9VgHHy$Z_ZCNEZV*ad!}|HJu@
znEZ$HHdu$xfDg@uwQJ+xYhB0hfWOqXnydVj?EU!f=6;Utj2Dg0Z~^P&^`VSr=)(#`
z$tb5%;kcSx->(@tElw-^Y5jeopbcjWLwZ(A7=>a<)61hOjy3h!q>)Tu>3mHKTCtze
z^kPxpt0nVBaj!}h^;}X*lykYgaB)Z)ER8cv3#nALpeLkXw6b3-j3o0|#~)WDg*Dne
zkjk5__w9&o?bUj3=%LLAJzuyXy0w2ZTEN~9tv9Aw+E~vHrEMA5)*T(tw)OSx>>bnw
zqum3&@Nn-Th{#p|EcE}#Z4b$WZs_)&g}H4E5!q8<<^T<xxh27lkMJ6wsdk%=#3u6z
zP!!eL(qPAS?1z{OZu_Rn`M5Zrxf4T6lnPo}&nIY?M*lWMO(2uyQj*#kDwVhy-Z?S1
zvBK^m?SRlU?DmlK+Kr;a7(Uka4><|yh>)~j<l{GiMC^SSXKwSbJ3!kp%<-XAxU@Hi
z7J)|5!nD+|AtYHGQQ>^SNQQ^=<?u)m!=-5Kbx=d)47MR<5==$A``c;wo91b~lvd%y
z-aMMH&?s7(+mgi+e#1XQ(4Z?OvpP9&Z6s@`Fb!iBMm0QK;7TbOS79SL3X$#ySQpGe
z6Hcc24oD{uO(ZLZbt_mNQb$vIE`wl}JFcAy;~vc6)|$KP&%Q?UT0A20ztk%3{e&{F
zWwmK5b;8SV&|Fgd>+PE(TPiZov?Qh0q2ly+@1LXZ8O|cld*Nta;`En!kzhy+8QuOj
zg8wJxqSAkvuM0~1Jgxmu|8@A&TwVNSJ}(%c<f2d?sVe6$^Xd#Zn#zg)G<Q%?-Unb}
z77iz(`Vwk1*A{=7UkhHymePN*6TBOGn!AfE^LN29_Lm$>|Cv@3s8E;1U*0DKr`S;1
z7k@c_+D}LR-E1fCCxY@mL;mFB_J0r@<wWA!|M&R+pS-uh#4Ma%5!H85bNYwaPjHWI
zX;tLks^B9o|1nkwik*ZJS?~#$|0F8}w~HgA_zBYDP-p+ASs{3kl#U{PLOq3w(_iL=
zg7RK0_U`;W!~SypRfsJaFJV!@?f)AnsD1Ho3t5(e*IJ_e*X{oYXq^6gLYArE8@8ua
zk?|w+o2UeF$a|=~@5_Eyx4p=MZ^O6Q;V<t&lkFf;#LklDWZ~)kfX1%0FYgtye_0Yk
z_a7DUlk4~xI68mZpMMWHu$ruuN?Tz66@S5baCZ8u2%C0TDD6l9@e_FgoQae6Wj@`)
z{%(6}%UKq3`)qEj#Q#d0u`1&4L>W6y(+CGl^MIlKv-BH^^hcgQ=Yo^IJN+N+vZP<Q
Xf=jz%ORt}`{NtasES9+(T%!69ogxcz

literal 0
HcmV?d00001

diff --git a/uebung_04/exercise_01/exercise_01.c b/uebung_04/exercise_01/exercise_01.c
new file mode 100644
index 0000000..ca5fb6f
--- /dev/null
+++ b/uebung_04/exercise_01/exercise_01.c
@@ -0,0 +1,17 @@
+#include "mpi.h"
+#include <stdio.h>
+
+int main(int argc, char **argv) {
+  int rank, size, err, n;
+
+  MPI_Init(&argc, &argv);
+
+  MPI_Comm_size(MPI_COMM_WORLD, &size);
+  MPI_Comm_rank(MPI_COMM_WORLD, &rank);
+
+  printf("rank=%d size=%d\n", rank, size);
+
+  MPI_Finalize();
+
+  return 0;
+}
diff --git a/uebung_04/exercise_01/exercise_01.out b/uebung_04/exercise_01/exercise_01.out
new file mode 100755
index 0000000000000000000000000000000000000000..1421d5eb8be0fbfa0b98fb70aacc5149a628a7de
GIT binary patch
literal 16208
zcmeHOeQXrR6`woXaFh?96DR}{I3)p*QZHvP6sQ~y+h^954G!39E0HYcyS2UYedTVi
z*r-h`(I6~i=}6nOQW`0=e}pJSqEaiR{i8LdrJ+(QNZW|2A9Aar+7fAu3z6F<T;H30
zZ@u2$(Ws63*N(OKoA;ad-psz8-JA8jc{bY7=J)vomn!kNkW!t+M5v+o3DyuIB6f-8
z@VP;(7mI+dh9g4FZoKZd8W;nlR|SZ6Yp_xVAGBz}$Rk8VyGrHiVylHw1T7iuYOqSO
z8^1#1WaM_czYrQj3CT|z-y@OHX@{{HskanHnigV_`a8!bL*qgk?TTbqBs)eyc8nZP
zbc4q(il-aH(1wN4Z8t!6Zfv#81S9r=<?szN@!#~<MRsdK(v5IqqeTlw?(Ye(!?=9V
z!k76F^;hXSj8Hro@p@y~nlyTMHf~KO>XSygFjPOZYiIq=MkSw7w#jxJR|pQAQ~M9}
z;N_S=j>0JXBbLnmux%x^qyCWt-@Mhe`}|WsI$krFf3^4g(vHp-Q3uD3HmHLI{qZG7
z8S6{oKpiU{CxnPt+%n*<>Ay&wS^AlGp<1^50NDG=w4lFG8NL^|a{4&Va{3>yz_(TK
z(_Mi-Ux9zN0{>zK{xt9p!{IK2K-7t~A{3z(KcS}er!r|RZ^m+_rU|XRyHiW(IlbS=
zn|iLhvn82H>)o;5q;9p9ZRzZ4*IF{El$JM+>%x|DvGkyn+S7(9GO4Vg;Sop8W0_ns
zAsgF_bS#N>nie05X?^n|%*i;7EEO}-pa8>+4{Gs&L9H)lBt<r7q|H77df*o)YhxIN
zI97Lx{T=PiE!s9^d!t>QuWeIyU<CXyv*!P*pk7T}m3R}{`M&zjnHOOCRg2#dx6&|{
z<81v4$%VDX67+TpWGwvr<LAyKeZQQd=Nt1c6Ca33Gr>H|*Y{dH&V9aQ$V@tTypOm{
zIe7W}r-m5^kM|RoIS2m{CINY#pbxyyxSS(>{Mx~Lh|2{B&(8#Gd)>iX_pfC!;o#{j
zOg6mX;5io9aMi&>puOTnz>9zv0WShx1iT1%5%|Cn_*?Dzf2b#J2Gp^dD_ez7pPMlK
zrK{?R9|tZ;AEg~Rpl3@P^8mFQ1nLi9%k1<YOQlj#%DC{EoxWzvxG<ZYzG%z10Ggft
zwk_iVZFc&sE#snUcKQs<&|8n}tyReWKz?7r-sj~zF8OVje8VOG-X;IqDL1^`K6>pb
zb@Wf_i9gSDbw?W}8m_40O+SNCl~(-{G;jAQwHv+&A_Qo>={LZr<27$0?4C1MLFit<
zI)rXXX|i_12<lxT1pQZ~{*G@TJa$(dol!6Tb+3AHu1fV?R<GSPSAv5l$w8nr*(cAF
z?J>S1O@%PHENtmfPc-cTO`(q7Fqf%gO>H1dUAbE-O(h_Bmuosf&-WzQ+Wn*b^f73}
zN)On<#SV>r;ifuTm`tc+)f+cUchR|q33a?@QXP*%(_8B3Pp6&-w;<kDM=!(K-iH52
zAbZ+kXDs#)#GV9ptmpPvbVe$_2Co93{5rT6Q+sfD&%uWj)v=zrgt}#=?ANTBIsy$>
zV5;Hp<UF)z9vXYai+~paF9Kc!ya;#^@FL(vz>9zv0WShx1U{??_~2I*{8_4LV<Ln<
zRsqVNJ$&n`K8e4<;;VUJrc}ar{vp63;3>d!fKz~7fY)x7N)v!s@Hb8@zT=04Zz$+n
zw`_5s2>%~I9)HU(g3TsWXk8pccoYuo=WFkjN?9aU1lv}u`Bd%FV}TK|ckP~!ZGZG(
zq|yEe9Fw3Y4UvR`ya~s3;Ew`w!M^afza=oXmynFE2%hw}EMMgRGK>fbuXqvgBH%^9
zi+~paF9Kc!ya;#^@FL(v;Qt!|&V%E;H_jpZ2@Qi~%)P^6=bq6^5d6|6%PC6Ie1K%m
zw_8Utzw)dlne*=OJqZi9-?>}LpnQRn*6{rgOOPs@vv!@@IY;h&OHJmcVs#D4s|Xp$
zOcW`9jdS>lq=4^pSU3l6wIzxO5zG(JKzQ8ONsn{kxc~LUv)(^T*7-qPB?Uh}%Jbb!
z{@|BVDRbQ6bydnMsQvFG^Ldc_4-M=*Z<og~UGE5C7vWLDEa3>D%ltozfrzEDzolhQ
zXmd|*A#D~y4Gl`860UD3NU`CWZQGP^<6}hM+a<=YTHGqRweb32PU&e1V2Bj%2P~(L
z*Ut}Gr(wzqC3O)Z{-}%JPCVaFtlth5ALjcuTMS=c%(C{0j|#pYSpV}-!G8JvWd4-Y
zUuEaBGXDxx%8kF$d3{&<St;g1>|q|u`MD&U@3-^iS^F)h_~35i`-%CRl3y*{^C<-<
z9un}Bh_I%GcL2)ec|Q&M_dy&4#d#B%1*D&aK2TbX`u;k>&m-1&Ncv%E9ye9+(+qqa
zjwod5R+GsS74!#z50*KvA>iddkmkFTwJ!l*?tISzU*~^71S1O#&J)k~3AcW~f}b~l
z55jzN`}tJ`{dX$x?~$L)G;g_IpTB(m`YPZt4(@&*1isw)?v{K|Ow;w1|0{z|9!K5}
z-!|^&$>Dg;Y)~?gdmigGG&9yOfXJJLzCI-`<`d;LGo{5bRXz`ziHz2t%=E^RTEfia
z@>;AgBp~TMo77D`p@hq8VB))>#d5jWu%@TY+_31&#ZtPKD5O%uVB(N8&@>mQ>O*=i
zZsc_>+yKMUw6;V0I-^?jKr5#7Yp}UOi>gGbt%roRzvE!@z7FkRTifAix7NL{xg!dS
zt)D%xud}^{(&{6eKaZ*M?rj1%*!k4{d+!=J2x<EXBdrzkdg30N2q11TGiK%cLx0Y_
z01Jj<Z!eIW4Kxk*DA;G3$Y->HSULf_6>K%R3qjL*^Lg9%!fgb$_eR%#1U;Q7yK6zy
zailOt*9HaK5$rgyHZ0(U5<<xjr_5L{pqaC9fNQXSLC<A{lFpdA(w{CU*<2>8=geUT
z)mt!NUyhL=YF~4EJ<da^J`l?f2qiI`hAu2Lb5_ge^jzM^q!$PpwB__<3=OE3O`1ZH
z!&IQE^k+cCN?wl(#ngvDl-F9xW#maK`T$MEfdqJC$?{>%Ny~!k&{ZsDfHTV+E`v~T
z2MEk@zPbMg#PuFLGT^(*D(?MtyxyZ`smpc5N8rFU7TZ_aM~1pscwJ>lLhOW!)84(G
z4FBKnEc{$&4!@_J_PkDD46!1i+y2wQzYA+L9zU;R88?$1E$9Ai{4r?9{UB`5>s-cN
z7+fegMy&Gr^Lo<&4wp7;&+BGJUMGW!EaZ+w@if$M4a)Ys?quZkJCC3B7@q-qT+_15
z>srPL_2eEu(PA7b7*n?Abv7ffv$=ow`G1w{o5^7WJA%Skv_!dg4P)lEKMNejg#Fv!
zKlDAs?=$2^j+N_(Ux&kKA0j)(tWz24)NJD`F8d-WFtR?oVVUuQ%l>6*U_8o(glxz7
zJ(oSNyBJTpCJfuL+aJ2@c^$^c?~APOzJKt&6T9H+FG8eA=zgDb+y5NYo%U<Uj<J#4
zF}Ur2>$2}6d&Y0EAtBEnrf)!H1swExVXxD9|EF7@WyV`zTkEjr_kvUPKtf?XOPZIN
z$NL-SF89ytZ{E+#fu-BR_Q%)pHL%6_+4l6FT-*vIl&l5zC-=*^2#lTfB0{QR3%MV5
zz;-Mz17qa4|4_tg36g!lo){J%vcUeVqy9<V=GlIW%~%y2D;6iZx{&P=ayUZ+dk9!u
o&+_>B`HSzDJa%bFA7Pq!tuE)>FYDG-us89G4a!=VflCzs0u)KU_W%F@

literal 0
HcmV?d00001

diff --git a/uebung_04/exercise_01/exercise_01.script b/uebung_04/exercise_01/exercise_01.script
new file mode 100644
index 0000000..649503a
--- /dev/null
+++ b/uebung_04/exercise_01/exercise_01.script
@@ -0,0 +1,13 @@
+#!/bin/bash
+
+## parameters
+#PBS -N uebung_04__exercise_01
+#PBS -q entry_teachingq
+#PBS -l select=1:ncpus=8:mpiprocs=8
+#PBS -l walltime=00:05:00
+
+## environment 
+cd ~/pwr/uebung_04/exercise_01
+
+## execute
+mpiexec -n 8 exercise_01.out
diff --git a/uebung_04/exercise_01/uebung_04__exercise_01.e662620 b/uebung_04/exercise_01/uebung_04__exercise_01.e662620
new file mode 100644
index 0000000..e69de29
diff --git a/uebung_04/exercise_01/uebung_04__exercise_01.o662620 b/uebung_04/exercise_01/uebung_04__exercise_01.o662620
new file mode 100644
index 0000000..89bfac5
--- /dev/null
+++ b/uebung_04/exercise_01/uebung_04__exercise_01.o662620
@@ -0,0 +1,8 @@
+rank=2 size=8
+rank=4 size=8
+rank=5 size=8
+rank=6 size=8
+rank=7 size=8
+rank=0 size=8
+rank=1 size=8
+rank=3 size=8
diff --git a/uebung_04/exercise_03/exercise_03.c b/uebung_04/exercise_03/exercise_03.c
new file mode 100644
index 0000000..2adcaf3
--- /dev/null
+++ b/uebung_04/exercise_03/exercise_03.c
@@ -0,0 +1,26 @@
+#include "mpi.h"
+#include <stdio.h>
+
+int main(int argc, char **argv)
+{
+	int rank, size, err, n;
+
+	MPI_Init(&argc, &argv);
+
+	MPI_Comm_size(MPI_COMM_WORLD,&size);
+	MPI_Comm_rank(MPI_COMM_WORLD,&rank);
+
+	printf("rank=%d size=%d\n", rank, size);
+
+	err = MPI_Barrier(MPI_COMM_WORLD);
+
+	n = (rank + 1) * 4711;
+
+	err = MPI_Bcast(&n, 1, MPI_INT, 1, MPI_COMM_WORLD);
+
+	printf("Received =%d\n", n);
+
+	MPI_Finalize();
+
+	return 0;
+}
diff --git a/uebung_04/exercise_03/exercise_03.out b/uebung_04/exercise_03/exercise_03.out
new file mode 100755
index 0000000000000000000000000000000000000000..dedb0d5375f17aadda8f0431c1f08ff4086d01b0
GIT binary patch
literal 16320
zcmeHOeQX>@6`woXiJLDwZA{aot(Vp*NrN|bVq2!j<ot2Yx^g~9Y$&31d%j!S+nn#*
zyS>D&+Hxgrg>fu7pz<NL$VH%nC?ZI$lp-jOLPG;8a)eY>q7n`w!j@8SOjArtbG$eE
z-g<9)7bKv6AninZZ{F{{-^}dW+4XGp%@@Pnoxwmra48X=5~zjN+enBF&0GPSXclYX
zw_4mPRsdWIM~KW``9{z-uq-FP5|GH-h?N5JW*ZhPd4wpDSFB7@2P{RU?TNfHtdj1P
zjWiCHY<K!okPJnUzm6<L2PiDvJdDk7tu%|_mV7xyaXE!fSZG`bBa|gxmUt`$@mO*^
zF(!C4(|LGh9oVq2^zsIY=an5cO|Zm1upGK*#r~824iK+Kk<AFNY_MU$lKcBK@HkZD
za!tXP`5^UI>^e+RJXzA^6g%Q(UsJ=5c&s*VCNd+nBh5{<O$}-~sqU2h;<ywza8B*n
z*Nc~98Zio^2xe`a+e6z1vZH_1Bem_%{HAnc{B-qn_G8_5R}2O-=!4@%4*Foh_V|*s
zz7{IzW5?r|5N$!YAB9*4M+F?`jo%KPSbfhs=9^plK`mJB1N8a=mJ3C2>`y7+KT^aW
z&!dq4$BW<xinQ}U5xlbqK2rpLstEpa5&Xx1-v)=bJO;peu}LUxWTA`m9mYfXxX|^1
zp=3f&TalEd>q779*{jEllrdnYEhE*lw>_Rr7(J1`xMAB0TJ|33(%Yhuv?Z|Co*Wv|
z)8<2l)Y>Aclxd`vEvZQ2u++K|rX`X?!={c$6j?`;sd!8pJIzESj%_#&X2OEeMn@ug
zKihO2Zi6TdekfulzyZ!AdRUJR9@hIKW?T%XfY&d;4?4m7qoEF!C{~-qp6;%;c73P1
ztHG%**LJFn5Uc>q??B$a63|K=v_xD1FPu+UUO4e}nD(V22X(|Ym?lj960r&<DC(I_
z=5=_-Yyce#KX>_gr_iJyq30^&v-J0_JS1s?an|41ZR2>)@jZ^ntP7WaUu-uq=faW0
zWx<8xJ<Daug|A@{q!O|fY#Z-wE<FDb$NQekH>sVK08j{|Y08Dul$H4RT{u2}xSVw1
z*cO*lE}Ub8wrLj*7r-g6y6{yFD#VNnU+u#EsviMA0)7Pi2>223Bk<uz;BOVTzN;O(
zRIZJeo!cRV_QbRm%*|-Wep!A-cAjhO19&mFB?_XVTA=?RS{COo<Z`)5spCRpasIrc
z;{s!G{<Ndxf@yL7CytH_jm7zA9UT`ji}NQO9T!ZC^N%??E<hINA7mXy(<n#NA>?>K
ze@(^FE$gp%^q+h57d`rqJo@)M`ZFH=DYst#R@d11`?Rq?YRCSxaG)n#KV5%Ln`rqp
zj5=5KHSm19U#+Nq98?JKM9a4T(I(2iisBPXRuzQ)GgybvUzeM$sLrC_8Is^{O8>?_
z6x*(9V+-2pzwFjdFO_J4v)cKq)&^+c7HXh8H`_1cz<j(u*_K9pY0Yf!)sD5iTMJcf
z?4q?s8*f<ww7DN#&E@7|5Ztq6Yr!w@1>idUBY*xV7_rg|JeVwXV~@P4jb&zH+IZ=f
zt+KiBQvI|x(L19}gu(Q>HukHzFF~`QUed<S!r8unlMB#?ZFJ5?zfI^7K*xJ8jfdx?
z^I=#SgYzk9TFl)k8D(j5T@C73_-$}KHXSgh2Vgxv7G8jk7P=?O_JPw&I2*Nzsx~;$
zmt`Hg-NLK|*}-Oa2j||0^REF1*`GGvyA;#5Z;&TumCd~e)^lg%Jd^i{b7$b-{NhKz
zkANQmKLUOP{0R6F@FU<yz>k0*f&ZTf;MW)YVAQfDrr;+fkgLQ&BWjpOjF=)d5!hVv
z$vQYCd`<u7Uvs$_$m0vS+;NakgZwVYn=j{b-5|>rbGe^^#BZHg0uLP&fsx9<<~1wJ
zC*f@jar~zLD)8<=hmMu<(@qUk@tf}?v_FaXy2{RV8}F`IeY89)c5k}t&RsRPA&mTc
z;aCDcNkk9|;uahWfM-E+!M<>utuVf;2!_{HJ{fFZyCQf6MudQ0{0R6F@FU<yz>k0*
z0Y3tM1pEm25%44M;Y5J*&p6)<KTu(rvxSqhg*k6n*4V;HfP)WLvTVJ<mN${kd2X9Y
z=U1Ceq;tL-z7Jtx`<1J?B<d$AsSV%pu-s1-&Rv@&JLhmIG;q${+Gx{66-gFE6Isgl
z;@q=oa=`aCESzI@qpgY=0vMmAf$+E&$d7Z)xc^%TXTSGy_W8kZB?O=O<oQlfd+=pb
z>U<sGwN&a^vj2m0J`ZaCe<qx~ZcmGK^l}>_d4S~oB!@|6NlucSBI#-F->RJFxTn4S
zE@f+PUnXH?l=^zLL9MH;&q%fY(VaWhx`u58U%N?+XQ{ZH^LpX+!tBoD6i|g6a5@+A
z$Lkpsm13CkGs#^jgim?!U4--f#r|DT31HrChmG*{?v?mQ1m7R*|1earU%p=%|FZP2
za`H+Ue-0{z#$W8bUX<;??EfgYuq+F;b4HqPaPsBZ`*o-U#75`-V*HZCZxr77h0hgA
zfn`@IpM-!Y0-h$JJXqk?kV5`nk$9!EPvj1wVcxBh--Og|NB`h@!Otyr*#@}Grw$Q9
z4~Qo8XD6534Y=acx%zq84lu}AN`?oD_&*MKWr6cO0l55+CT^JI^MDr`&rbngAG}Fa
zhVlt#3FrHYt-mYM&IQ0LVSaksd8>&3-;3baL;MPjvl?)0e=E&@9#0M6h4R^(rT+)d
zYj=_M@0EC^*hZbp|7(LzB8Ug!ubB6_bSRp#>eVD9lt=nZ-HHqd0Mb^bzh8}t<pg=%
z8q%YfBA*7&SW+K|C;K9CJ!U0SX+4q|5s=nC95*Z@rq&hqz-)F?ha~yPsBR>z)TroB
zMTQJLmKhov1&T}4!PClfHAak7)Jz+C-7XlGu6G`6-5b`!`#La{Ux!T^dRQY+>o_R%
zJ>C1;TD$fAot=lmJ$g@TTXz^7J3hCsb#GUDUb?=WHy`2zdCZ0PZWX|R=cngm_P@;P
zzxFnQ{g9*|GZT6yZN#pzs{jHTu_AVMKXmTS{LdR#{+<I}hfNDk2x94^J{U>FVEcmO
zCO07HdS5zCZKK8BzTn+H;PlgwzbOH(RBZXwU0WA)1IG-b^=wRV8iJh$_C^J~mO`lM
z(IG3+2hvK}a*%7VZ^1|n3pJ6n40RxpQHN8>VIyUYx}d&{30r{77(uOVU9~t>rTbtc
zJt)-JXac&hrIoTRM~qb3OeXRa9c(Eh9zg=thU1n{<uFyKssl+-v642TLbZ$$P~{X*
zQ%QN!sxe5jbub34v1YenPg|S8b?7QGWI|&$9WIDaF_v%}@csYaN?dQkqX_<1+Qqw{
zkJp=ImcCp^c{?1q#$$f5eQW58h1YeqCPWid+<fo;HT<sNF8p_aF<ck8`MfS+sjwqS
zFaNWEzXNMT9zU<6S=JGcmXjYCe+2BfhGssmvsq>_xKMD6SmpEQ^{xpRE~S{y>von3
z@xetFk18Gn6xYDa=XEbj#ij`!Kl`zK6zsUBW_^krS%!%39Y4Wh0xB3&=FbwJ<pMc!
z|K9We2Jzd-kKY?uPP_ORGcW&Hz%VA<zVmmA{$BC>6g4Bq%5{`a!{O#D#A7+>c7{5c
z9r>Jxf07(nvOhP&I?IzD{xlg_W|>Hmc`RS@@I$noX8C8&gkc^x`wI`B*LN)WeU<&a
z_s=WD=j$&N@+67R2Y2!EUjri!k@+=>?ai{mR-Ioj|My^V^Y2k?SC-E^OuNGKhvAD*
zSqBH+gIM@IocEi0{aI&u88{ViFrVK;MryIO<@suB%R2Ml2LsMs=JR@=_YZSm@czTX
zJid-U1PssL;nRD1b_aS3k(Iy8Kg?&j0?2N@2$64{Ex8|VfqAU20WxCTzY?-7mBcUi
zTrkBpJHP9xe{y$l=5KcpyMklI;+8yjY>WLHq6V9(|62gYr7e%2pTGD%%VW2h_U{Ny
bya7+++%Nm$?N*5Y2mHeZWs`^CQN_OiZA>Ev

literal 0
HcmV?d00001

diff --git a/uebung_04/exercise_03/exercise_03.script b/uebung_04/exercise_03/exercise_03.script
new file mode 100644
index 0000000..34e10b3
--- /dev/null
+++ b/uebung_04/exercise_03/exercise_03.script
@@ -0,0 +1,13 @@
+#!/bin/bash
+
+## parameters
+#PBS -N uebung_04__exercise_03
+#PBS -q entry_teachingq
+#PBS -l select=1:ncpus=8:mpiprocs=8
+#PBS -l walltime=00:05:00
+
+## environment 
+cd ~/pwr/uebung_04/exercise_03
+
+## execute
+mpiexec -n 8 exercise_03.out
diff --git a/uebung_04/exercise_05/exercise_05.c b/uebung_04/exercise_05/exercise_05.c
new file mode 100644
index 0000000..c236687
--- /dev/null
+++ b/uebung_04/exercise_05/exercise_05.c
@@ -0,0 +1,33 @@
+#include "mpi.h"
+#include <stdio.h>
+
+int main(int argc, char **argv) {
+  int rank, size, err, n, sum;
+
+  MPI_Init(&argc, &argv);
+
+  MPI_Comm_size(MPI_COMM_WORLD, &size);
+  MPI_Comm_rank(MPI_COMM_WORLD, &rank);
+
+  printf("rank=%d size=%d\n", rank, size);
+
+  n = rank + 1;
+
+  sum = 4710;
+
+  MPI_Barrier(MPI_COMM_WORLD);
+
+  MPI_Reduce(&n, &sum, 1, MPI_INT, MPI_SUM, 0, MPI_COMM_WORLD);
+  printf("Received in reduction=%d\n", sum);
+
+  MPI_Barrier(MPI_COMM_WORLD);
+
+  sum = 3710;
+
+  MPI_Allreduce(&n, &sum, 1, MPI_INT, MPI_SUM, MPI_COMM_WORLD);
+  printf("Received in allreduction =%d\n", sum);
+
+  MPI_Finalize();
+
+  return 0;
+}
diff --git a/uebung_04/exercise_05/exercise_05.out b/uebung_04/exercise_05/exercise_05.out
new file mode 100755
index 0000000000000000000000000000000000000000..24df63b1c1b98a1d02c4ba7aede2dba6bbf10a5c
GIT binary patch
literal 16392
zcmeHOeQX@n5r1dLiGd`(G;s+IO_q;Tp{X}^;<$mz&DqX5YibAM*rlIr&Ub5j#eJB)
zJ+KL+t_&2*SPF>JRBA=|1BldCN>z|jq^geiMk;coibgFhP?gG7LkVs`p@>{(-kVvk
zmpfXa{#EU>_U6rRKHlu?*IVzr@5Z}2Dk2fVtx|kNAa`JmLSiH^c2NgNOzaY?;CGYQ
zD3$}iOkzqNQv{@{ovv6)s|l|JigxR<Q-WTuV8N6lM2dFhns$+rFcne7qg@qti4W4}
zC=RA9_r}xE7@8)3J=%l{(3tw|FgL^XQY?mBmeOOCmzmWHXFG(^Zj|gs$&RTYJEojZ
z%n5G0Xg`9q0c6;i2JME)E=cz%Gr<(cz_$M_H}z@8>mj=yO}Y_5+N5B?l*hXf>@Y9?
zZR5k-OXHP04o4`TOzCinZE1U;xoKNERiC!Ah0*%aUCs5)O?p13Z<irpTpDb6PVL>-
zkB4IhF-oJT7*{;|L)ki#qyEr_N1Fe9*YEf3Klt*&>DM3r^<#rupF<st8*NYr8~Woz
z&ODwnsH5^RA+)v%ff(B2WjkxVlDTnm<K+)r^7Z%a*R~7-RJP?n@|jXn7@!o6<7y@J
zCjl#^|6my$W7S{_wl4yp1wr6DU@PVSU>Urx4E|0T{K+!-pUU8`0DdKG!Da%mM$`#y
z4@u}kJIZ();dD?>G9Fus3)39R<g#YoO*pP;3bV7X+e}%GHDu>q%jxTGPv^2$Ut%C_
zDOpKLcTcC;p37v+ynTlykn6Qlg`_2NnGxH>EjME33mGYDOE`{gIkGQ}oJ95($#rIJ
zS9V&{X(!+yiSBRDIq8(t=&-YiG!EvCU}s&}ndE4~9Auek!qt_eeaj^5EGWQ^CT}s5
z!?&1&2|F!D9IzV{pa)UnW!KbzOcJ}zVsBSxTf4bk-_hiC7kk_FotVoATsM)xUnTI%
zJX~p2igl62$HJ2jJPeoBGO>}cb+)|<ufl7HFVxxagsBiafXBx7b-tg*=p-GX`!?fz
z@>a*BnP8mx)4LUp*EHWl5P8Fgqm0{}4`0C`(6c_AA1x?8@55n-yr#ue7w=>A$d`ad
z=PBZRFGT^L*NE33K<zO3$5<#eQgXzH(<v|Uqdpw(d)$usa2$i%aUagPL)i%*PLF<R
zGULN9W+y;j^x>EI@L3<e(uarLFau!*!VH8N2s7~in1O%RZ2Xroai-dss(NXg5XL<-
zZbfm{n0U7O8M(W~o!<m_p}6^5Kx#G#)E`F4!u*@XV)3x#@qxcEf6C+WLBBA6(&KSL
zvoQY)kH-i8!u+Eij}Q8V`3F26CpruB-}88!R4mNj>hU-MSeQ>R57BLw(ZQS>;{pCD
z9iv_3e;(i;3-CV)@IMIf(*gdT0Do72FZlV!w>u|K-E2($#hCc(*`B_5<4ogA#&pZe
zuouPJJZPR9)N3}~4J;gh>6S4-jOi)|$+hR*S~xV#*oQ-NQSptMP2;Hd3{lYUl=?f{
zkX&=Ym^^Eo{KszN<oQY?^1N~Cg1ZhJd<6eaKnK;uHwI;1*dFsW-coxVxGZe#Hzr#C
zfd6<fCf{;b8dEKQ17_|=7mCHX6eRoks<*F&+1+=*)*Bz~=WhogcKX2%PNs&*yIwOU
z3$rO>YT4#3(p~&~<BT!gKWj|KLG-FI`NEt7Zh<{xOg;~L`w*UdfZnOlIfWh~^fo}J
z`p-<o=cMviVFnAzC&0Cst3z3QZek`f89)0%{M;bSz`<cVng0DF9Hwv%hB_-x^y%6y
zAnI-X6<T*qSJgG5e(hS!(B!G7@m$70)?y$N&qk->=j51kFe~)Nd|W!MdL2^^y59#D
zMn1F}SQ@Hpdeupgg8!!*5Uu?syxIEwPk;^ke%$Z-dJxP#1wV54jH&+fDP!w8nMb#3
z?tT!w1eYeJTkei`rQ+^{EzCfefiMGM2Eq)483;2FW+2Q!n1L_@<uia^N$@*Q%jT4Z
zABlioEPAb^W#49{G&`%w?>(-a%YICfppQW))T9)KiEOC+GJd<oSNhp^ibeci!N^(o
zM*+}rpvQnl-z^q<fPMgU1}HZCHi|8BN3V#CMk5<muBbi??_!8!4ZsPo{X8n{Ss_11
zTn$|u=Kzdx6sTAo?O473`kG5_uO1h>>#qC4j;pUk8126e+iB2~h6qAK+=UHe#Mrpu
zSog!$2(e}GT`;~ndT&Mhs^yib2!zCHAq_JSW+2Q!n1L_@VFtnsgc%4k5N06EK-mm%
zT^QGa;j0C<71UuKYtgX57a`T?r3Jot$#!C`qL&fR^=BK1=Qor(;<*lu>!?}&-i2Zg
zRgP1!8P{4Jp$^xY&5@jIsCH9)TsyU%#;GOBgv`V!)m3rr*J)C~_b+T*yLPE!g+T!_
z9-|J&T}|V2trw5Kk?dLT!=l<hP?r#V&6E4f1Hu<e$&)()ua}ZPO*{Sp@w^}8{(nZi
zdhUS720Gp`qCG?p5FH^pPV_L*qePDr4Y>V}EY~6KZEwF$+tNQ!$hrlsu~Bc*8|oVi
zl5IS+eY@V!bPd6u>=N_3OuSnRYT<Fi8qE7Cu^K7hX<bSmk86d9io^cmM8Gt{X9D;R
z!udJ_R<=&~jbL4Fm*V-jW0iG8d{*#vqNGCH2^}2QyKYHvpVY7Q>ONWj3FwrHzubO3
zCH<@u4Og&-MOw<wGg7?PtCwf(SD_OT>%HrY@iP*?R0Pkn_W&<dNBs~aL=$jViLnPo
zacM{?{UZ{Odh0~w|7O4h?_b(d7m~XQ^()p0zJIX_*44^-=@=n&qi9BbR`OB2|CGYJ
z0oQzb+_eE;gBB!QMuO2Y`VRsgEs=*u0GIzzV#h?WUZ7MwF9N=%;xZA9Ej9QZ;e5TZ
z^vyDU-T^!c=WWo>hh_AaL;gy|b43|k2ORxxq3eR<*#WpbE^UMWk8LmGP+KD#-EeOL
z9P<{8^Jc(H#dA>NQPD|4`M)-BIF5J)JpF?AvHeNMZPas6DV`XxO*b(l0LZ(A!9hJK
z77OD|H)AHTG(HcSshl~K&J85eX3EVuc{5QM6;Q7}lD1qcr8kt;z%q5)gj)E-m}zBQ
zXG{z_iHv2Y3Yp9pnD{soG~Fet)~Mwq?Yw0+?1ZpPv!l1QJ8s7J?ZFa$6Bc!taf3i(
zPp>fdcHPj{+GXC*(Xl_?XZE$Wb;Ut(&o}n9c6YYRf_kpXR|WAgu6)NL_uvu%418(b
zJy!geDBl0%#RNA%Eq=<*nuWZT`h>LwkluuwP(}O@iod{rNt8?1AebgBUGNf;%ID1C
zL^cJB7(6w(>cBJy@_F*@Tec8fCg2U!v~=+S9JthCCVWd2yj2Ob?!mNhCm{U5f&|Yo
ztT0du6!1_Bq36dkZejqa>nJ+RJy@?`IU_>P=3Gl3$`<qyCpTg_?wAi6DA=%~$W9T|
z+SXZ*XR=fuPUMG$o*K)-5Q@5vlH6uFdH5gu62=6XW2F;lK)sQ)D|8vA4qbgH2P}5-
zR#ND$H43aeLv$x6cUrfG>0};Gfj8!q4|P^63+}^EiHr@-lsOzBp<{00DxsUfrzbdn
zg<Bdtb5#qj3*`ALiKQ<0kzNHGP8Hd{-1;}v#m4hQ#R<_29lw2W9UOjd@Hc*rFoyFB
zzdg@ym};y@G-!Vl;P1nnljG-kInySx<9T^7{shQzPR{l`Z)bWOlM4-F#4g+O{MiN!
zr(|r;?*UAC-VZ9WkvnF^t<b|cHQVz%nknax<7Yjlhrk}^=*-JN(K`kVEhc|I9zVfi
z8akL$wvVwLQGRdW@q_#SL$Yrphx{JGl;2y>9&;A7e-tpx3H$e+U-TT~_cKtDEqMHY
z45;5;BRi%?m6`X;k4L6Y1ng%>z?Aja0rO0c2kd7_!1OR15@kE4SWW4V{{YR)nSPBA
zOu(Mq{w840^F*fn{>%Eo{5?nZ-Y^<z6743INznecAi((9zDHA<O!@j#s5;lUkMxhQ
z`R(u15EX*y&lKzZ^7+HSThM{>QE%-Rzt{7+)u2A}Oy33D8lOGC&y3Uq2#xg=x5zX8
z5eV?yrPng=y@A&ub7JUn^8E4f|0&qw{(JWH-aonxKxp?V;{EctOqYYP-(JK>wLwuH
zhaIpT^DDs^F&<xwDM^&<t6@s%Z<?Qh{u&ye)IBNNZ}kw>!C0~RDW5w^qJCrKaF|ZG
s&jE~6UXGvdzZU_<u`5gUdyG!JdjrmST-L?Qt(3j~C8bdpFbJ^X-&2~2ssI20

literal 0
HcmV?d00001