/** Get boundary matrix 'bm' of 'order', compute reduces boundary
* matrix 'rbm', and the transformation matrix 'tm' such that rbm = bm *
* tm. */
/** Get boundary matrix 'bm' of 'order', compute reduces boundary
* matrix 'rbm', and the transformation matrix 'tm' such that rbm = bm *
* tm. */