From 70ae3adaff2cb0005e66dcd52c669ee003361f6b Mon Sep 17 00:00:00 2001
From: Emilie Marchandise <emilie.marchandise@uclouvain.be>
Date: Tue, 15 Dec 2009 09:48:01 +0000
Subject: [PATCH]
---
Solver/TESTCASES/AdvectionDiffusion.lua | 4 ++--
Solver/TESTCASES/WavePulse.lua | 8 ++++----
2 files changed, 6 insertions(+), 6 deletions(-)
diff --git a/Solver/TESTCASES/AdvectionDiffusion.lua b/Solver/TESTCASES/AdvectionDiffusion.lua
index 847b2e9e2e..90a85807b8 100644
--- a/Solver/TESTCASES/AdvectionDiffusion.lua
+++ b/Solver/TESTCASES/AdvectionDiffusion.lua
@@ -1,6 +1,6 @@
model = GModel ()
model:load ('square.geo')
-model:mesh(2)
+model:load ('square.msh')
dg = dgSystemOfEquations (model)
dg:setOrder(5)
@@ -40,7 +40,7 @@ dg:L2Projection(createFunction.lua(1,'initial_condition','XYZ'))
dg:exportSolution('output/Advection_00000')
-- main loop
-for i=1,10 do
+for i=1,10000 do
norm = dg:RK44(0.001)
if (i % 1 == 0) then
print('iter',i,norm)
diff --git a/Solver/TESTCASES/WavePulse.lua b/Solver/TESTCASES/WavePulse.lua
index a4d0074051..815eb769e9 100644
--- a/Solver/TESTCASES/WavePulse.lua
+++ b/Solver/TESTCASES/WavePulse.lua
@@ -8,8 +8,8 @@ function initial_condition( _x , _f )
X = xyz:get(i,0) - .5
Y = xyz:get(i,1) - .5
Z = xyz:get(i,2)
- VALUE = math.exp(-40*(X*X+Y*Y+Z*Z));
- f:set(i,0,7)
+ VALUE = math.exp(-40*(X*X+Y*Y+Z*Z));
+ f:set(i,0,VALUE)
f:set(i,1,0.0)
f:set(i,2,0.0)
end
@@ -21,8 +21,8 @@ end
print'*** Loading the mesh and the model ***'
myModel = GModel ()
-myModel:load('square.geo')
-myModel:mesh(2)
+myModel:load('box.geo')
+myModel:load('box.msh')
print'*** Create a dg solver ***'
DG = dgSystemOfEquations (myModel)
--
GitLab