Point(2) = {0, -0.5, 0, .03};
Point(3) = {0, 0.5, 0, .03};
Point(4) = {-0.5, 0.5, 0, .03};

Point(5) = {0.5, -0.5, 0, .03};
Point(6) = {0.5, 0.5, 0, .03};

Line(1) = {4, 3};
Line(2) = {3, 2};
Line(3) = {2, 1};
Line(4) = {1, 4};

Line(7) = {3, 6};
Line(8) = {6, 5};
Line(9) = {5, 2};

Line Loop(5) = {2, 3, 4, 1};
Line Loop(6) = {8, 9, -2, 7};
Plane Surface(7) = {5};
Plane Surface(8) = {6};
Transfinite Line {4,8} = 20 ;
Transfinite Line{2}=20;
Transfinite Surface {7};
Recombine Surface {7};
Transfinite Line {3,1,7,9} = 10;
Transfinite Surface{8};

Physical Line("Border") = {1, 3, 4, 7, 8, 9};
Physical Surface("Inside1") = {7};
Physical Surface("Inside2") = {8};

-This line, and those below, will be ignored--

M    Common/LuaBindings.cpp
M    Geo/MVertex.cpp
M    Geo/GVertex.h
M    Geo/MVertex.h
M    Solver/dgAlgorithm.cpp
M    Solver/TESTCASES/ForwardFacingStep.lua
M    Solver/dgSlopeLimiter.cpp
M    Solver/dgSystemOfEquations.cpp
M    Solver/dgConservationLawPerfectGas.cpp

No commit message

GAMMA = 1.4;
U = 3.0
V = 0.0 
RHO  = 1.4;

PRES = RHO*U*U/(GAMMA*MACH*MACH)
--PRES = 1;
--PRES = ./(MACH*RHO*RHO*GAMMA*GAMMA) 

SOUND = math.sqrt(U*U+V*V)/MACH

--[[ 
     Function for initial conditions
--]]
function free_stream( XYZ, FCT )
  for i=0,XYZ:size1()-1 do
    FCT:set(i,0,RHO) 
    FCT:set(i,1,RHO*U) 
    FCT:set(i,2,RHO*V) 
    FCT:set(i,3, 0.5*RHO*(U*U+V*V)+PRES/(RHO*GAMMA-1)) 
  end
end

--[[ 
     Example of a lua program driving the DG code
--]]
order = 1
print'*** Loading the mesh and the model ***'
myModel   = GModel  ()
myModel:load ('step.geo')
myModel:load ('step.msh')

print'*** Create a dg solver ***'
DG = dgSystemOfEquations (myModel)
DG:setOrder(order)
FS = functionLua(4, 'free_stream', {'XYZ'}):getName()

law=dgPerfectGasLaw2d()
DG:setConservationLaw(law)

law:addBoundaryCondition('Walls',law:newWallBoundary())
law:addBoundaryCondition('LeftRight',law:newOutsideValueBoundary(FS))

DG:setup()

print'*** setting the initial solution ***'

DG:L2Projection(FS)
DG:limitSolution()

print'*** export ***'

DG:exportSolution('output/solution_0')

print'*** solve ***'
CFL = 2.0;

for i=1,1000 do
    dt = CFL * DG:computeInvSpectralRadius();
    norm = DG:RK44_limiter(0.1*dt)
    print('*** ITER ***',i,dt, norm)
    if (i % 1 == 0) then 
       DG:exportSolution(string.format("output/solution-%06d", i)) 
    end
end

print'*** done ***'

No commit message

No commit message

Added Navier Stokes & 2 test cases in DG
Added binary restart capabilities in DG
Added mesh renumbering option (-renumber)@

No commit message

No commit message

No commit message

though. doc interface is present but no doc is generated

No commit message

No commit message

spectral radius computed
quads work

No commit message

No commit message