From 5112c3bcf9348c8ce2808f7a64991d61b8f48229 Mon Sep 17 00:00:00 2001 From: Jonathan Lambrechts <jonathan.lambrechts@uclouvain.be> Date: Fri, 8 Jan 2010 15:20:12 +0000 Subject: [PATCH] correct comment --- Common/LuaBindings.cpp | 3 ++- Solver/TESTCASES/WavePulse.lua | 8 +++----- Solver/dgGroupOfElements.cpp | 3 ++- Solver/dgSystemOfEquations.cpp | 2 +- 4 files changed, 8 insertions(+), 8 deletions(-) diff --git a/Common/LuaBindings.cpp b/Common/LuaBindings.cpp index c829911287..c75451330a 100644 --- a/Common/LuaBindings.cpp +++ b/Common/LuaBindings.cpp @@ -123,7 +123,8 @@ void binding::interactiveSession(){ add_history(expansion); report_errors(L, luaL_dostring(L, expansion)); } - free(expansion); + if(expansion) + free(expansion); } clear_history(); #else diff --git a/Solver/TESTCASES/WavePulse.lua b/Solver/TESTCASES/WavePulse.lua index 982cfeb34d..a9dee53840 100644 --- a/Solver/TESTCASES/WavePulse.lua +++ b/Solver/TESTCASES/WavePulse.lua @@ -1,5 +1,3 @@ - - --[[ Function for initial conditions --]] @@ -55,9 +53,9 @@ print('DT = ',dt) for i=1,N do -- norm = DG:multirateRK43(dt) norm = DG:RK44(dt) - print('*** ITER ***',i,norm) - if (i % 10 == 0) then - DG:exportSolution(string.format("output/solution-%04d", i)) + if (i % 100 == 0) then + print('*** ITER ***',i,norm) + DG:exportSolution(string.format("output/solution-%04d", i)) end end diff --git a/Solver/dgGroupOfElements.cpp b/Solver/dgGroupOfElements.cpp index f93e120e72..42d9044416 100644 --- a/Solver/dgGroupOfElements.cpp +++ b/Solver/dgGroupOfElements.cpp @@ -88,6 +88,7 @@ dgGroupOfElements::dgGroupOfElements(const std::vector<MElement*> &e, int polyOr for (int l=0;l<nbPsi;l++) { imass(k,l) += f[k]*f[l]*weight*detjac; } + // (iQP*3+iXYZ , iFace*NPsi+iPsi) (*_dPsiDx)(j*3 ,i*nbPsi+k) = g[k][0]*ijac[0][0]+g[k][1]*ijac[0][1]+g[k][2]*ijac[0][2]; (*_dPsiDx)(j*3+1,i*nbPsi+k) = g[k][0]*ijac[1][0]+g[k][1]*ijac[1][1]+g[k][2]*ijac[1][2]; (*_dPsiDx)(j*3+2,i*nbPsi+k) = g[k][0]*ijac[2][0]+g[k][1]*ijac[2][1]+g[k][2]*ijac[2][2]; @@ -266,7 +267,7 @@ void dgGroupOfFaces::init(int pOrder) { getElementLeft(i)->getJacobian ((intLeft)(j,0), (intLeft)(j,1), (intLeft)(j,2), jac); inv3x3(jac,ijac); //compute dPsiLeftDxOnQP - //(iPsi*3+iXYZ,iQP+iFace*NQP); + // (iQP*3+iXYZ , iFace*NPsi+iPsi) int nPsi = _fsLeft->coefficients.size1(); for (int iPsi=0; iPsi< nPsi; iPsi++) { (*_dPsiLeftDxOnQP)(j*3 ,i*nPsi+iPsi) = g[iPsi][0]*ijac[0][0]+g[iPsi][1]*ijac[0][1]+g[iPsi][2]*ijac[0][2]; diff --git a/Solver/dgSystemOfEquations.cpp b/Solver/dgSystemOfEquations.cpp index a4d074086c..8616b5bf53 100644 --- a/Solver/dgSystemOfEquations.cpp +++ b/Solver/dgSystemOfEquations.cpp @@ -43,7 +43,7 @@ static dgSystemOfEquations *myConstructorPtr(GModel* gm){ return new dgSystemOfEquations(gm); } -void dgSystemOfEquations::registerBindings(binding *b){ +void dgSystemOfEquations::registerBindings(binding *b) { classBinding *cb = b->addClass<dgSystemOfEquations>("dgSystemOfEquations"); cb->setDescription("a class to rule them all :-) -- bad description, this class will be removed anyway"); cb->setConstructor<dgSystemOfEquations,GModel*>(); -- GitLab