This test is written like a class. A function at the module level is probably sufficient, something like:

```
def test_gaussian2D():
    m = models.Gaussian2DModel(2., 3., 4., x_stddev=1., y_stddev=5., theta=30.)
    assert self.m(3., 4.) == 2.
```