plusQS Usage: C = plusQS(A, B [, bfac, opts]) adds one QSpace to another as in C = A+bfac*B ([bfac=1]) block-diagonal identity matrix. Options '-v' verbose mode (e.g. on remaining outer multiplicity) (C) AW : Jun 2006 ; Oct 2014