template <class VecX, class VecY> inline
void do_add(VecX& x, VecY& y)
{
  for (int i = 0; i < N; ++i) {
    x[i] = i;
    y[i] = i + 1;
  }
  blais_vv::add<N>(x.begin(), y.begin());
}
  int ix[N], iy[N];
  external_vec<int> x1(ix, N);
  external_vec<int> y1(iy, N);
  do_add(x1, y1);

  dense1D<int> x2(N);
  dense1D<int> y2(N);
  do_add(x2, y2, N);