anchor.diff:      write(*,*) "-------------------------------------------------------------"
       write(*,*) "Statsitcs:"
       write(*,*) "-------------------------------------------------------------"
  -    write(*,*) "Time to compute = ", endTime-startTime
  +    write(*,*) "Time to compute = ", elapsedTime
       write(*,*) ""
       write(*,*) "Avg Iters for iterating betn. soln. =", avgIters
       write(*,*) "Max Iters for iterating betn. soln. =", maxIters
       write(*,*) "Avg Iters for linear solver =", avgNit
       write(*,*) "Max Iters for linear solver =", maxNit  
       write(*,*) "Writing statistics to file"
  -    call reportStats(avgIters,maxIters,avgNit,maxNit,endTime-startTime)
  +    call reportStats(avgIters,maxIters,avgNit,maxNit,elapsedTime)
       
       write(*,*) "Execution completed"
     
  
positive.code:         
  
      write(*,*) "Starting Solver"
      call system_clock(startTime, clock_rate, clock_max)
      call solveOrder2(tEnd,nOuts,tDel,n,row,col,M,C,yLen,xDel,yDel,&
         c0,m0,y,nu,gam,alpha,beta,k,delta,ndiag,e1,e2,nit,esoln,dSelect,&
         fSelect,gSelect,avgIters,maxIters,avgNit,maxNit)
      call system_clock(endTime, clock_rate, clock_max)
      
  !    call solveOrder1(tEnd,nOuts,tDel,n,row,col,M,C,xLen,yLen,xDel,yDel&
  !       ,c0,m0,nu,gam,alpha,beta,k,delta,ndiag,e1,e2,nit,dSelect,fSelect)
    
  
negative.code: subroutine reportStats(avgIters,maxIters,avgNit,maxNit,time)
    implicit none
    real,intent(in)::avgIters,maxIters,avgNit,maxNit
    real,intent(in)::time
    integer :: stat
  
    open(UNIT = 125, IOSTAT = stat, FILE = "statReport.dat", STATUS = "old")
  