From 4edf3df52942514999ae22f36e87293a028b58dc Mon Sep 17 00:00:00 2001
From: Emilie Marchandise <emilie.marchandise@uclouvain.be>
Date: Thu, 7 Jan 2010 10:16:34 +0000
Subject: [PATCH]
---
Solver/TESTCASES/ForwardFacingStep.lua | 12 +++++-------
1 file changed, 5 insertions(+), 7 deletions(-)
diff --git a/Solver/TESTCASES/ForwardFacingStep.lua b/Solver/TESTCASES/ForwardFacingStep.lua
index d1bbe4aeab..c1bb142ad0 100644
--- a/Solver/TESTCASES/ForwardFacingStep.lua
+++ b/Solver/TESTCASES/ForwardFacingStep.lua
@@ -25,7 +25,7 @@ end
--[[
Example of a lua program driving the DG code
--]]
-
+sv
order = 1
print'*** Loading the mesh and the model ***'
myModel = GModel ()
@@ -55,14 +55,12 @@ print'*** export ***'
DG:exportSolution('output/solution_0')
print'*** solve ***'
+CFL = 2.0;
-LC = 0.1*.1
-dt = .3*LC/(SOUND+V);
-print('DT=',dt)
-
-for i=1,10 do
+for i=1,1000 do
+ dt = CFL * DG:computeInvSpectralRadius();
norm = DG:RK44_limiter(0.1*dt)
- print('*** ITER ***',i,norm)
+ print('*** ITER ***',i,dt, norm)
if (i % 1 == 0) then
DG:exportSolution(string.format("output/solution-%06d", i))
end
--
GitLab