cvec TestEseNonGrayFixture::conv(const cmat &a, const cvec &b) {
// check a.rows()==b.length()
int len = a.cols()+b.length()-1;
cvec ret = zeros_c(len);
for(int i=0; i<a.cols(); i++) {
cvec tmp = zeros_c(len);
tmp.set_subvector(i, elem_mult(a.get_col(i), b));
ret += tmp;
}
return ret;
}
