oops

Solver/TESTCASES/Stommel.lua

@@ -8,7 +8,6 @@ groups = dgGroupCollection(model, dimension, order)
 -- groups:split...
 groups:buildGroupsOfInterfaces()
 solution = dgDofContainer(groups, claw:getNbFields())
---[[
 a = fullMatrix(3,1);
 a:set(0,0, 1.);
 a:set(1,0, 2.);
@@ -16,9 +15,9 @@ a:set(2,0, 3.);
 f = functionConstant(a):getName()
 solution:L2Projection(f)
 solution:exportMsh('output/init')
-]]--
 
-for i=1,00000 do
+
+for i=1,100000 do
   norm = dg:RK44(150*(3/(2.*order+1)))
   if ( i%100 ==0 ) then
     print ('iter ', i, norm)