def METHOD_NAME (self):
    fcisolvers = [fci.solver (molsym, symm=True, singlet=False) for i in range (2)]
    fcisolvers[0].nroots = fcisolvers[1].nroots = 2
    fcisolvers[0].wfnsym = 'A1'
    fcisolvers[1].wfnsym = 'B1'
    mc = mcscf.addons.state_average_mix (mcscf.CASSCF (msym, 4, 4), fcisolvers, [0.25,]*4)
    mo = mc.sort_mo([4,5,6,10], base=1)
    mc.kernel(mo)
    self.assertAlmostEqual (mc.e_tot, mc_ref.e_tot, 8)
    for e1, e0 in zip (numpy.sort (mc.e_states), mc_ref.e_states):
        self.assertAlmostEqual (e1, e0, 5)