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