anchor.diff:  #endif
   
   
  -!---- Test: loading integer array
  -!---- Test: loading double precision array
  +!---- Test: loading 1-D double precision array
  +vec5 = 0.0d0
  +vec2(1) =  8.0d0
  +vec2(2) = -3.0d0
  +call load_dp_array1(vec5,5,vec2,2,3,'vec2',ierr)
  +call test_value(vec5(1), 0.0d0,'load_dp_array1: vec5(1)')
  +call test_value(vec5(2), 0.0d0,'load_dp_array1: vec5(2)')
  +call test_value(vec5(3), 8.0d0,'load_dp_array1: vec5(3)')
  +call test_value(vec5(4),-3.0d0,'load_dp_array1: vec5(4)')
  +call test_value(vec5(5), 0.0d0,'load_dp_array1: vec5(5)')
  +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'
  +
  +
  +!---- Test: loading 2-D double precision array
   mat5x4 = 0.0d0
   do i = 1,2
       do j = 1,3
  
positive.code: 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'
  
  !---- Test: linear equation solver
  ! dgesv
  
negative.code: 
  implicit none
  
  integer :: i, j, ierr, nrows, ncols, imat5x4(5,4), imat2x3(2,3)
  double precision :: vec2(2), vec5(5), vec4_constraint(4)
  double precision :: mat2x3(2,3), mat5x2(5,2), mat5x4(5,4)
  logical :: isConst(4)
  