2 calls to
THETA_AND_THETAM_LBC_ONLY

SOLVE_EM
SOLVE_EM