You are working with Boost so why not use its built-in matrices (and simply pass their “references” between your routines)?