def test_6():
# Test that _c_fitkappa is consistent with the two-step python implementation
np.random.seed(42)
vm = VonMisesHMM(n_states=13)
kappas = np.random.randn(13, 7)
posteriors = np.random.randn(100, 13)
obs = np.random.randn(100, 7)
means = np.random.randn(13, 7)
vm.kappas_ = kappas
vm._c_fitkappas(posteriors, obs, means)
c_kappas = np.copy(vm._kappas_)
vm._py_fitkappas(posteriors, obs, means)
py_kappas = np.copy(vm._kappas_)
np.testing.assert_array_almost_equal(py_kappas, c_kappas)
