void
dswap(size_t n,
      double *x, ptrdiff_t incX,
      double *y, ptrdiff_t incY);