'N param_sets, N 2D data sets, Poly2d'

def test_p2_nset_npset(self):

    p2 = models.Polynomial2D(5, n_models=2)
    xx = np.array([self.x, self.x])
    yy = np.array([self.y, self.y])
    z = p2(xx, yy)
    assert (z.shape == (2, 10, 10))
