'This case covers: N param sets , N sets 1D x --> N N sets 1D y data'

def test_non_linear_NXYset(self):

    g1 = models.Gaussian1D([10, 10], [3, 3], [0.2, 0.2], n_models=2)
    xx = np.array([self.x1, self.x1])
    y1 = g1(xx)
    assert_allclose(y1[:, 0], y1[:, 1], atol=(10 ** (-12)))
