  We present an explicit formula for the expected value of a product of several
independent symplectically invariant matrices in which the trace and real part
function may be applied, possibly to different subexpressions. This takes the
form of a topological expansion; however, each term has two topologies: one for
the trace, and another for the real part.
  The traces and real parts can always be written in terms of index
contraction, but in some cases, it is possible to write the expression as a
product in which the two functions are applied to bracketed intervals in a
legal bracket diagram. We present the conditions under which this may be done,
and an algorithm to construct such an expression given the contracted indices
when possible.
  The summands in the topological expansion are written in terms of matrix
cumulants. We compute the matrix cumulants of quaternionic Ginibre, Gaussian
symplectic, quaternionic Wishart, and Haar-distributed symplectic matrices,
which allow direct computation of an expression constructed from several
independent ensembles of any of these matrices.
