anchor.diff:  call load_dp_array1(vec5,5,vec2,2,5,'vec2',ierr)
   call test_value(ierr,1,'load_dp_array1: overflow')
   write(stdout,*) 'subroutine load_dp_array1() passed unit test'
  -
  +write(stdout,*)
   
   !---- Test: loading 2-D double precision array
   mat5x4 = 0.0d0
  
positive.code: call test_value(vec2(1),-3.90d0,'solve_dgels: vec2(1)')
  call test_value(vec2(2), 1.98d0,'solve_dgels: vec2(2)')
  write(stdout,*) 'subroutine solve_dgels() (least-squares solver) passed unit test'
  #endif
  
  !---- Test: non-negative least squares solver
  
negative.code: !     write(0,*) vec3(i)
  ! enddo
  write(stdout,*) 'subroutine solve_dgssv() (sparse linear equation solver) passed unit test'
  write(stdout,*)
  #endif
  
  !---- Test: least squares solver
  