anchor.diff:        complex cx(*), cy(*), c, cdotu
         integer n, incx, incy
         
  -      cdotu(n, cx, incx, cy, incy, c)
  -      wzdotu = zdotu
  +      c = cdotu(n, cx, incx, cy, incy)
  +      wcdotu = c
         
         return
         end
  
positive.code:       double complex zx(*), zy(*), z, zdotu
        integer n, incx, incy
        
        zdotu(n, zx, incx, zy, incy, z)
        wzdotu = zdotu
        
        return
        end
  
negative.code:       complex cx(*), cy(*), c, cdotc
        integer n, incx, incy
        
        c = cdotc(n, cx, incx, cy, incy)
        wcdotc = c
        
        return
        end
  