From 9a20545e798a68ebd73f92ca0f19a2cfbebf5983 Mon Sep 17 00:00:00 2001 From: Erich Eckner Date: Tue, 4 Aug 2015 15:55:25 +0200 Subject: alle Runge-Kuttas in .inc ausgelagert, minimales dT auf dX gesetzt (sonst Bug z.B. in AY) --- rk108.inc | 64 +++++++++++++++++++++++++++++++-------------------------------- 1 file changed, 32 insertions(+), 32 deletions(-) (limited to 'rk108.inc') diff --git a/rk108.inc b/rk108.inc index 7635367..0f98175 100644 --- a/rk108.inc +++ b/rk108.inc @@ -3,9 +3,9 @@ felders[aktuelleFelder], 0.1 * dT ); - felders[2].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + felders[2].berechneAbleitungen(dX,iDX,pDNMax); - if pruefeMaxDT(2,mDT,dT,dTMin) then + if pruefeMaxDT(2,dTMax,dT,dTMin) then continue; felders[3].liKo( @@ -15,9 +15,9 @@ -0.91517656137529144052001501927534215431895138766436972056466 * dT, 1.45453440217827322805250021715664459117622483736537873607016 * dT ); - felders[3].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + felders[3].berechneAbleitungen(dX,iDX,pDNMax); - if pruefeMaxDT(3,mDT,dT,dTMin) then + if pruefeMaxDT(3,dTMax,dT,dTMin) then continue; felders[4].liKo( @@ -27,9 +27,9 @@ 0.202259190301118170324681949205488413821477543637878380814562 * dT, 0.606777570903354510974045847616465241464432630913635142443687 * dT ); - felders[4].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + felders[4].berechneAbleitungen(dX,iDX,pDNMax); - if pruefeMaxDT(4,mDT,dT,dTMin) then + if pruefeMaxDT(4,dTMax,dT,dTMin) then continue; felders[5].liKo( @@ -41,9 +41,9 @@ 0.197966831227192369068141770510388793370637287463360401555746 * dT, -0.0729547847313632629185146671595558023015011608914382961421311 * dT ); - felders[5].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + felders[5].berechneAbleitungen(dX,iDX,pDNMax); - if pruefeMaxDT(5,mDT,dT,dTMin) then + if pruefeMaxDT(5,dTMax,dT,dTMin) then continue; felders[6].liKo( @@ -55,9 +55,9 @@ 0.410459702520260645318174895920453426088035325902848695210406 * dT, 0.482713753678866489204726942976896106809132737721421333413261 * dT ); - felders[6].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + felders[6].berechneAbleitungen(dX,iDX,pDNMax); - if pruefeMaxDT(6,mDT,dT,dTMin) then + if pruefeMaxDT(6,dTMax,dT,dTMin) then continue; felders[7].liKo( @@ -71,9 +71,9 @@ 0.48966295730945019284450701113589820117801547843379009721079 * dT, -0.0731856375070850736789057580558988816340355615025188195854775 * dT ); - felders[7].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + felders[7].berechneAbleitungen(dX,iDX,pDNMax); - if pruefeMaxDT(7,mDT,dT,dTMin) then + if pruefeMaxDT(7,dTMax,dT,dTMin) then continue; felders[8].liKo( @@ -87,9 +87,9 @@ 0.0325402621549091330158899334391231259332716675992700000776101 * dT, -0.0595780211817361001560122202563305121444953672762930724538856 * dT ); - felders[8].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + felders[8].berechneAbleitungen(dX,iDX,pDNMax); - if pruefeMaxDT(8,mDT,dT,dTMin) then + if pruefeMaxDT(8,dTMax,dT,dTMin) then continue; felders[9].liKo( @@ -103,9 +103,9 @@ 0.32176370560177839010089879904987890408140436860307712925111 * dT, 0.51048572560806303157775901228512341674467213703175235406759 * dT ); - felders[9].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + felders[9].berechneAbleitungen(dX,iDX,pDNMax); - if pruefeMaxDT(9,mDT,dT,dTMin) then + if pruefeMaxDT(9,dTMax,dT,dTMin) then continue; felders[10].liKo( @@ -121,9 +121,9 @@ 0.499269229556880061353316843969978567860276816592673201240332 * dT, 0.50950460892968610423609869004538625398664323235298960218506 * dT ); - felders[10].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + felders[10].berechneAbleitungen(dX,iDX,pDNMax); - if pruefeMaxDT(10,mDT,dT,dTMin) then + if pruefeMaxDT(10,dTMax,dT,dTMin) then continue; felders[11].liKo( @@ -141,9 +141,9 @@ 0.0107558956873607455550609147441477450257136782823280838547024 * dT, -0.327769124164018874147061087350233395378262992392394071906457 * dT ); - felders[11].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + felders[11].berechneAbleitungen(dX,iDX,pDNMax); - if pruefeMaxDT(11,mDT,dT,dTMin) then + if pruefeMaxDT(11,dTMax,dT,dTMin) then continue; felders[12].liKo( @@ -163,9 +163,9 @@ -0.0782942710351670777553986729725692447252077047239160551335016 * dT, -0.114503299361098912184303164290554670970133218405658122674674 * dT ); - felders[12].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + felders[12].berechneAbleitungen(dX,iDX,pDNMax); - if pruefeMaxDT(12,mDT,dT,dTMin) then + if pruefeMaxDT(12,dTMax,dT,dTMin) then continue; felders[13].liKo( @@ -191,9 +191,9 @@ -3.25479342483643918654589367587788726747711504674780680269911 * dT, -2.16343541686422982353954211300054820889678036420109999154887 * dT ); - felders[13].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + felders[13].berechneAbleitungen(dX,iDX,pDNMax); - if pruefeMaxDT(13,mDT,dT,dTMin) then + if pruefeMaxDT(13,dTMax,dT,dTMin) then continue; felders[14].liKo( @@ -221,9 +221,9 @@ -1.71561608285936264922031819751349098912615880827551992973034 * dT, -0.0244036405750127452135415444412216875465593598370910566069132 * dT ); - felders[14].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + felders[14].berechneAbleitungen(dX,iDX,pDNMax); - if pruefeMaxDT(14,mDT,dT,dTMin) then + if pruefeMaxDT(14,dTMax,dT,dTMin) then continue; felders[15].liKo( @@ -241,9 +241,9 @@ 0.0910895662155176069593203555807484200111889091770101799647985 * dT, 0.777333643644968233538931228575302137803351053629547286334469 * dT ); - felders[15].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + felders[15].berechneAbleitungen(dX,iDX,pDNMax); - if pruefeMaxDT(15,mDT,dT,dTMin) then + if pruefeMaxDT(15,dTMax,dT,dTMin) then continue; felders[16].liKo( @@ -255,9 +255,9 @@ -0.157178665799771163367058998273128921867183754126709419409654 * dT, 0.157178665799771163367058998273128921867183754126709419409654 * dT ); - felders[16].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + felders[16].berechneAbleitungen(dX,iDX,pDNMax); - if pruefeMaxDT(16,mDT,dT,dTMin) then + if pruefeMaxDT(16,dTMax,dT,dTMin) then continue; felders[17].liKo( @@ -293,9 +293,9 @@ -0.34275815984718983994222055341385087174233873470395891993726 * dT, -0.675 * dT ); - felders[17].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + felders[17].berechneAbleitungen(dX,iDX,pDNMax); - if pruefeMaxDT(17,mDT,dT,dTMin) then + if pruefeMaxDT(17,dTMax,dT,dTMin) then continue; felders[1-aktuelleFelder].liKo( -- cgit v1.2.3-70-g09d2