Add scalar regression test
@@ -492,6 +492,7 @@ def test_ecpl_integrate():
# regression test to check the numerical integration for small energy bins
ecpl = ExpCutoffPowerLawSpectralModel()
value = ecpl.integral(1 * u.TeV, 1.1 * u.TeV)
+ assert value.isscalar
assert_quantity_allclose(value, 8.380714e-14 * u.Unit("s-1 cm-2"))
