Skip to content
Snippets Groups Projects
Commit 937fdff2 authored by Jean-François Remacle's avatar Jean-François Remacle
Browse files

No commit message

No commit message
parent a09ff1c5
No related branches found
No related tags found
No related merge requests found
...@@ -20,7 +20,7 @@ end ...@@ -20,7 +20,7 @@ end
Example of a lua program driving the DG code Example of a lua program driving the DG code
--]] --]]
order = 5 order = 4
xyz = functionCoordinates.get() xyz = functionCoordinates.get()
FS = functionLua(4, 'free_stream', {xyz}) FS = functionLua(4, 'free_stream', {xyz})
...@@ -54,19 +54,20 @@ DG:exportSolution('output/cyl_0') ...@@ -54,19 +54,20 @@ DG:exportSolution('output/cyl_0')
print'*** solve ***' print'*** solve ***'
CFL = 5.1; CFL = 4.2;
dt = CFL * DG:computeInvSpectralRadius(); dt = CFL * DG:computeInvSpectralRadius();
print('DT = ',dt) print('DT = ',dt)
T = 0; T = 0;
x = os.clock()
for i=1,100000 do for i=1,100000 do
dt = CFL * DG:computeInvSpectralRadius(); dt = CFL * DG:computeInvSpectralRadius();
norm = DG:RK44(dt) norm = DG:RK44(dt)
T = T + dt T = T + dt
if (i % 10 == 0) then if (i % 10 == 0) then
print('*** ITER ***',i,norm,dt,T) print('*** ITER ***',i,norm,dt,T,os.clock() - x)
end end
if (i % 100 == 0) then if (i % 100 == 0) then
DG:exportSolution(string.format("output/cyl-%06d", i)) -- DG:exportSolution(string.format("output/cyl-%06d", i))
-- DG:saveSolution(string.format("output/cyl-%06d.bin", i)) -- DG:saveSolution(string.format("output/cyl-%06d.bin", i))
end end
end end
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment