anchor.diff:    
     contains
     
  +  subroutine randsvd(nu, nv, m, n, x, s, u, vt, info)
     subroutine randsvd(method, jobu, jobvt, m, n, a, k, q, s, u, vt, info)
       ! in/out
       character*1         method, jobu, jobvt
  
positive.code:     
      return
    end
  
  
  
  
    ! stage a from the paper, with the orthonormalization step
    subroutine randsvd_stage_a(m, n, k, q, a, y, u, omega, work, lwork, info)
      ! in/out
  
negative.code:     ! general
      do j = 2, q
        call randsvd_subspaceiter(m, n, k, q, a, u, y, tau, work, lwork, info)
        if (info /= 0) exit
      end do
      
      
  