anchor.diff:  call load_int_array2(imat5x4,5,4,imat2x3,2,3,3,3,'imat2x3',ierr)
   call test_value(ierr,1,'load_int_array2: overflow')
   write(stdout,*) 'subroutine load_int_array2() passed unit test'
  +write(stdout,*)
   
   
   !---- Test: load equation constraints
  
positive.code: !     write(0,*) vec3(i)
  ! enddo
  write(stdout,*) 'subroutine solve_dgssv() (sparse linear equation solver) passed unit test'
  #endif
  
  !---- Test: least squares solver
  
negative.code: call test_value(vec2(1), 4.10d0,'solve_nnls: vec2(1)')
  call test_value(vec2(2), 1.98d0,'solve_nnls: vec2(2)')
  write(stdout,*) 'subroutine solve_nnls() (non-negative least-squares solver) passed unit test'
  write(stdout,*)
  
  
  
  