anchor.diff:      else
         iy = 0
       endif
  -    ! ..
       do i = 1, N
         stemp = stemp + ( sx( ix ) * sy( iy ) )
         ix = ix + strideX
  
positive.code: !
  ! @param {integer} N - number of values to use from each array
  ! @param {Array<real>} sx - first array
  ! @param {integer} strideX - stride length on the first array
  ! @param {Array<real>} sy - second array
  ! @param {integer} strideY - stride length of the second array
  ! @returns {real} the dot product of sx and sy
  !<
  real function sdot( N, sx, strideX, sy, strideY )
    implicit none
  
negative.code:   endif
    sdot = stemp
    return
  end function sdot