'1 set 2d x, 1set 2D y, 2 param_sets, NonLinearFitter'

@pytest.mark.skipif(u'not    HAS_SCIPY')
def test_nonlinear_lsqt_Nset_2d(self):

    with pytest.raises(ValueError):
        g2 = models.Gaussian2D([10, 10], [3, 3], [4, 4], x_stddev=[0.3, 0.3], y_stddev=[0.2, 0.2], theta=[0, 0], n_models=2)
        z = g2(self.x.flatten(), self.y.flatten())
        gfit = fitting.LevMarLSQFitter()
        model = gfit(g2, self.x, self.y, z)
