inline
void
spop_repmat::apply(SpMat<typename T1::elem_type>& out, const SpOp<T1, spop_repmat>& in)
{
arma_extra_debug_sigprint();
typedef typename T1::elem_type eT;
const unwrap_spmat<T1> U(in.m);
const SpMat<eT>& X = U.M;
const uword X_n_rows = X.n_rows;
const uword X_n_cols = X.n_cols;
const uword copies_per_row = in.aux_uword_a;
const uword copies_per_col = in.aux_uword_b;
// out.set_size(X_n_rows * copies_per_row, X_n_cols * copies_per_col);
//
// const uword out_n_rows = out.n_rows;
// const uword out_n_cols = out.n_cols;
//
// if( (out_n_rows > 0) && (out_n_cols > 0) )
// {
// for(uword col = 0; col < out_n_cols; col += X_n_cols)
// for(uword row = 0; row < out_n_rows; row += X_n_rows)
// {
// out.submat(row, col, row+X_n_rows-1, col+X_n_cols-1) = X;
// }
// }
SpMat<eT> tmp(X_n_rows * copies_per_row, X_n_cols);
if(tmp.n_elem > 0)
{
for(uword row = 0; row < tmp.n_rows; row += X_n_rows)
{
tmp.submat(row, 0, row+X_n_rows-1, X_n_cols-1) = X;
}
}
// tmp contains copies of the input matrix, so no need to check for aliasing
out.set_size(X_n_rows * copies_per_row, X_n_cols * copies_per_col);
const uword out_n_rows = out.n_rows;
const uword out_n_cols = out.n_cols;
if( (out_n_rows > 0) && (out_n_cols > 0) )
{
for(uword col = 0; col < out_n_cols; col += X_n_cols)
{
out.submat(0, col, out_n_rows-1, col+X_n_cols-1) = tmp;
}
}
}
