hpa.fun.toy(x)  Internal function pdm.maker() shows how the B matrix is
  obtained from the various elements of input argument x.  Note
  how, in this simple example, the B matrices are diagonal, but
  generalizing to non-diagonal matrices should be straightforward (if
  you can guarantee that they remain positive definite).