Let $\Phi(t)$ be an $n\times n$ complex matrix whose columns are (independent) solutions of the system of ordinary differential equations (ODE): $\frac{d}{dt}y= A(t) y$, where $A(t)$ is a $t$-dependent $n\times n$ complex matrix.

Then, the Liouville formula provides a very simple relation between the determinant of $\Phi(t)$ at different $t$: $\det(\Phi(t)) = e^{\int_0^t Tr A(s) ds} \det(\Phi(0))$.

Now, consider the alternative system of ODE:

$\frac{d}{dt}y= A(t) y + B(t) \bar{y}$,

where $\bar{x}$ represents the complex conjugate of $x$. (Actually, I am interested only in the case where $A(t)$ is diagonal and $B(t)$ is skew-symmetric, if that helps.)

Does it exist any generalization of the Liouville formula applicable to this case? or any other expression for $\det(\Phi(t))$, which is simpler (to compute numerically) than finding the explicit solutions for all the starting conditions in $\Phi(0)$ and then computing the determinant explicitly?

The kind of manipulations used by the classical proof of Liouville formula do not seem very helpful to me, in this slightly modified case, unfortunately.

The solution of this problem would be very beneficial to the method described in arXiv:1205.3996, and to its very relevant physical applications.