trans_coef.T.dot(solve(H.dot(V).dot(H.T)), trans_coef)
trans_coef.T @ solve(H @ V @ H.T, trans_coef)
trans_coef' * ((H * V * H') \ trans_coef)