parameter(nsize=1000)
      dimension a(nsize+1,nsize),  b(nsize,1), ipiv(nsize)
c     call set_d_stream(1)
c     call set_i_stream(1)
      write(6,*) 'test of sgesv'
      do ibid=1,10
      do i=1,nsize
      b(i,1)=i
      do j=1,nsize
      a(i,j)=0.
      enddo
      enddo
      do i=1,nsize
      a(i,i)=2.
      enddo
      lda=nsize+1
      ldb=nsize
      nrhs=1
      n=nsize
      t0=second()
      call SGESV( N, NRHS, A, LDA, IPIV, B, LDB, INFO )
      t1=second()
      tdel=t1-t0
      write(6,*) 'test number ',ibid,' nsize =',nsize
      write(6,*) 'mflops =',.666667*nsize*nsize*nsize/1.e6/tdel
      write(6,*)
      end do
      end