  We compute the soft anomalous dimension (SAD) matrices for production of
massive quarks $Q$ and $\bar{Q}$ in association with a gluon jet, from massless
quarks $q$ and antiquarks $\bar{q}$: $q\bar{q}\rightarrow Q\bar{Q}g$, and in
the gluon scattering $gg\rightarrow Q\bar{Q}g$. To analyze the behaviour of the
eigenvalues of SAD matrices we perform numerical studies of their eigensystems
at two special kinematical configurations.
