'Polynomial evaluation of multiple data sets with different domain'

def test_nset_domain(self):

    xx = np.array([self.x1, self.x1]).T
    xx[(0, 0)] = 100
    xx[(1, 0)] = 100
    xx[(2, 0)] = 99
    p1 = models.Polynomial1D(5, n_models=2)
    yy = p1(xx, model_set_axis=(-1))
    x1 = xx[:, 0]
    x2 = xx[:, 1]
    p1 = models.Polynomial1D(5)
    assert_allclose(p1(x1), yy[0, :], atol=(10 ** (-12)))
    p1 = models.Polynomial1D(5)
    assert_allclose(p1(x2), yy[1, :], atol=(10 ** (-12)))
