alle Runge-Kuttas in .inc ausgelagert,

minimales dT auf dX gesetzt (sonst Bug z.B. in AY)

1 files changed, 47 insertions, 0 deletions

diff --git a/rk3_8.inc b/rk3_8.inc
new file mode 100644
index 0000000..3419dab
--- /dev/null
+++ b/rk3_8.inc
@@ -0,0 +1,47 @@
+  felders[2].liKo(
+    felders[aktuelleFelder],
+    felders[aktuelleFelder],
+    1/3 * dT
+  );
+  felders[2].berechneAbleitungen(dX,iDX,pDNMax);
+
+  if pruefeMaxDT(2,dTMax,dT,dTMin) then
+    continue;
+
+  felders[3].liKo(
+    felders[aktuelleFelder],
+    felders[aktuelleFelder],
+    felders[2],
+    -1/3 * dT,
+    1 * dT
+  );
+  felders[3].berechneAbleitungen(dX,iDX,pDNMax);
+
+  if pruefeMaxDT(3,dTMax,dT,dTMin) then
+    continue;
+
+  felders[4].liKo(
+    felders[aktuelleFelder],
+    felders[aktuelleFelder],
+    felders[2],
+    felders[3],
+    1 * dT,
+    -1 * dT,
+    1 * dT
+  );
+  felders[4].berechneAbleitungen(dX,iDX,pDNMax);
+
+  if pruefeMaxDT(4,dTMax,dT,dTMin) then
+    continue;
+
+  felders[1-aktuelleFelder].liKo(
+    felders[aktuelleFelder],
+    felders[aktuelleFelder],
+    felders[2],
+    felders[3],
+    felders[4],
+    1/8 * dT,
+    3/8 * dT,
+    3/8 * dT,
+    1/8 * dT
+  );