p1 = astmodels.Polynomial1D(1)
p2 = astmodels.Polynomial1D(1)
p3 = astmodels.Polynomial1D(1)
p4 = astmodels.Polynomial1D(1)
m1 = p1 & p2
m2 = p3 & p4
m1.inverse = m2
return m1
test_models = [
astmodels.Identity(2),
astmodels.Polynomial1D(2, c0=1, c1=2, c2=3),
astmodels.Polynomial2D(1, c0_0=1, c0_1=2, c1_0=3),
astmodels.Shift(2.),
astmodels.Hermite1D(2, c0=2, c1=3, c2=0.5),
astmodels.Legendre1D(2, c0=2, c1=3, c2=0.5),
astmodels.Chebyshev1D(2, c0=2, c1=3, c2=0.5),
astmodels.Chebyshev2D(1, 1, c0_0=1, c0_1=2, c1_0=3),
astmodels.Legendre2D(1, 1, c0_0=1, c0_1=2, c1_0=3),
astmodels.Hermite2D(1, 1, c0_0=1, c0_1=2, c1_0=3),
astmodels.Scale(3.4),
astmodels.RotateNative2Celestial(5.63, -72.5, 180),
astmodels.Multiply(3),
astmodels.Multiply(10 * u.m),
astmodels.RotateCelestial2Native(5.63, -72.5, 180),
astmodels.EulerAngleRotation(23, 14, 2.3, axes_order='xzx'),
astmodels.Mapping((0, 1), n_inputs=3),
astmodels.Shift(2. * u.deg),
astmodels.Scale(3.4 * u.deg),
astmodels.RotateNative2Celestial(5.63 * u.deg, -72.5 * u.deg, 180 * u.deg),
c0=2,
c1=3,
c2=0.5,
domain=(0.0, 1.0),
window=(1.5, 2.5)),
astropy_models.Chebyshev2D(1, 1, c0_0=1, c0_1=2, c1_0=3),
astropy_models.Chebyshev2D(1,
1,
c0_0=1,
c0_1=2,
c1_0=3,
x_domain=(1.0, 2.0),
y_domain=(3.0, 4.0),
x_window=(5.0, 6.0),
y_window=(7.0, 8.0)),
astropy_models.Hermite1D(2, c0=2, c1=3, c2=0.5),
astropy_models.Hermite2D(1, 1, c0_0=1, c0_1=2, c1_0=3),
astropy_models.Legendre1D(2, c0=2, c1=3, c2=0.5),
astropy_models.Legendre2D(1, 1, c0_0=1, c0_1=2, c1_0=3),
astropy_models.Polynomial1D(2, c0=1, c1=2, c2=3),
astropy_models.Polynomial2D(1, c0_0=1, c0_1=2, c1_0=3),
# astropy.modeling.powerlaws
astropy_models.BrokenPowerLaw1D(amplitude=10,
x_break=0.5,
alpha_1=2.0,
alpha_2=3.5),
astropy_models.ExponentialCutoffPowerLaw1D(10, 0.5, 2.0, 7.),
astropy_models.LogParabola1D(
amplitude=10,
x_0=0.5,