anchor.diff:  call test_value(vec5(5),  0.0d0,'load_dp_array_constraints: vec5(5)')
   call test_value(ncols,2,'load_dp_array_constraints: ncols output')
   write(stdout,*) 'subroutine load_dp_array_constraints() passed unit test'
  +write(stdout,*)
   
   
   !---- Test: linear equation solver
  
positive.code: call load_dp_array2(mat5x4,5,4,mat2x3,2,3,3,3,'mat2x3',ierr)
  call test_value(ierr,1,'load_dp_array2: overflow')
  write(stdout,*) 'subroutine load_dp_array2() passed unit test'
  
  
  !---- Test: loading 2-D integer array
  
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,*)
  
  
  
  