reformat; increase ncells; add @time.

Author: alsam

examples/scalar_law/PROGRAM1/consdiff.jl

@@ -23,7 +23,8 @@
 
 function consdiff(fic::Int,lac::Int,fif::Int,laf::Int,fix::Int,lax::Int,
                   ifirst::Int,ilast::Int,
-                  x::OffsetArray{Float64,1},flux::OffsetArray{Float64,1},
+                  x::OffsetArray{Float64,1},
+                  flux::OffsetArray{Float64,1},
                   conservd::OffsetArray{Float64,1})
 #   ******************************************************************
 #   update conservd to new time

examples/scalar_law/PROGRAM1/linaddmain.jl

@@ -44,15 +44,15 @@ include("consdiff.jl")
 include("upwind.jl")
 
 function main()
+    #const ncells = 100
     const ncells = 10000
 
     u    = OffsetArray(Float64, -2:ncells+1)
     x    = OffsetArray(Float64,  0:ncells)
     flux = OffsetArray(Float64,  0:ncells)
     dfdu = OffsetArray(Float64, -2:ncells+1)
 
-    const nsteps   =  100
-    #const nsteps   =  10000
+    const nsteps   =  10000
     const tmax     =  0.8
     const cfl      =  0.9
 
@@ -87,16 +87,17 @@ function main()
     end
 
     # write final results (plot later)
-    #@unsafe for ic=0:ncells-1
-    for ic=0:ncells-1
+    @unsafe for ic=0:ncells-1
+    #for ic=0:ncells-1
         xc = (x[ic]+x[ic+1])*0.5
         uc = u[ic]
-        @printf("%e %e\n",xc,uc)
+        #@printf("%e %e\n",xc,uc)
     end
 
 
 end # main
 
-main()
+@time main()
+@time main()
 
 

examples/scalar_law/PROGRAM1/riemprob.jl

@@ -22,7 +22,8 @@
 #***********************************************************************
 
 function initsl(ncells::Int, fc::Int, lc::Int, fx::Int, lx::Int, ifirst::Int, ilast::Int,
-                conservd::OffsetArray{Float64,1}, x::OffsetArray{Float64,1})
+                conservd::OffsetArray{Float64,1},
+                x::OffsetArray{Float64,1})
 #   ******************************************************************
 #   interior values only; others defined by boundary conditions
 #   ******************************************************************
@@ -64,9 +65,9 @@ function bcmesh(fim::Int,lam::Int,ncells::Int, x::OffsetArray{Float64,1})
 end
 
 function bccells(fic::Int,lac::Int,ncells::Int, conservd::OffsetArray{Float64,1})
-#     ******************************************************************
-#     right side 
-#     ******************************************************************
+#   ******************************************************************
+#   right side 
+#   ******************************************************************
     if lac >= ncells
         # println("treating right side")
         @unsafe for ic=ncells:lac