anchor.diff:         !write(*,*)'Info',info
          !write(*,'("Eigs",2x, 5(f10.3,x))')alpha(1:iter_lanczos)
          !call system('matlab -nodisplay -r  gen_lanczos && quit')
  +    call system('free -g')
  +    call system('sar -r 1 5')
  +    call system('pstree -ap $(echo $$)')
       end if
     end if
   
  
positive.code: 
      !mrow=369999909999_i18 ! (0:mrow) works > 256 processors
      !mrow=  99999_i18
      rdel= 1_i18
      cdel= 1_i18
      !mrow= (nproc*ntd)*10000_i18 -1_i18
      mrow= (nproc*ntd)*10000_i18 -1_i18
      ! cdel*nproc: so max_block is divisible by cdel
      !print *, mrow
  
  
negative.code:     !Call MPI_Barrier(MPI_COMM_WORLD, mpierr)
      if (rank==0) &
      write(6,*)"Mrow:max_block", mrow,max_block
      if (rank==0) &
      write(6,*)"MemoryMB:per mpi", r_mem
      !
      iter_lanczos=100
      !call lanczos_matvec
      !if (rank==0) &
      !Write(6,*)'Entering main', rank
  
      !stop
  ! Papi profiling
  