function [op_l, op_r] = assemble_boundary_op(op,m)
l = length(op);

op_l = sparse(m,1);
op_l(1:l) = op;

op_r = sparse(m,1);
op_r(m-l+l:m) = flipr(op);
end
