anchor.diff:        complex cx(*), cy(*), c, cdotc
         integer n, incx, incy
         
  -      cdotc(n, cx, incx, cy, incy, c)
  -      wzdotc = zdotc
  +      c = cdotc(n, cx, incx, cy, incy)
  +      wcdotc = 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:       double complex function wzdotc(n, zx, incx, zy, incy)
        double complex zx(*), zy(*), z
        double complex zdotc
        integer n, incx, incy
        
        z = zdotc(n, zx, incx, zy, incy)
        wzdotc = z
        
        end
        
  