  The operator function (A,B)\to\tr f(A,B)(K^*)K, defined on pairs of bounded
self-adjoint operators in the domain of a function f of two real variables, is
convex for every Hilbert Schmidt operator K, if and only if f is operator
convex. As a special case we obtain a new proof of Lieb's concavity theorem for
the function (A,B)\to\tr A^pK^*B^{q}K, where p and q are non-negative numbers
with sum p+q\le 1. In addition, we prove concavity of the operator function
(A,B)\to \tr(A(A+\mu_1)^{-1}K^* B(B+\mu_2)^{-1}K) on its natural domain
D_2(\mu_1,\mu_2), cf. Definition 4.1
