Add bounds checking for sncosmo models

[jeremykubica]

Instead of crashing when given wavelengths out of bounds, return zero flux in those wavelengths.

[github-actions]

Before [0386b25]	After [4f37069]	Ratio	Benchmark (Parameter)
2.03±0.7ms	1.32±1ms	~0.65	benchmarks.TimeSuite.time_fnu_to_flam
16.5±0.1μs	17.1±0.2μs	1.04	benchmarks.TimeSuite.time_sample_x1_from_hostmass
5.11±0.1ms	5.24±0.2ms	1.03	benchmarks.TimeSuite.time_evaluate_salt3_passbands
1.40±0.01s	1.42±0.01s	1.01	benchmarks.TimeSuite.time_make_x1_from_hostmass
10.5±0.5ms	10.5±0.9ms	1.00	benchmarks.TimeSuite.time_evaluate_salt3_model
8.83±0.07ms	8.82±0.1ms	1.00	benchmarks.TimeSuite.time_load_passbands
28.8±0.2μs	28.9±0.4μs	1.00	benchmarks.TimeSuite.time_make_new_salt3_model
131±1μs	128±0.4μs	0.98	benchmarks.TimeSuite.time_sample_x0_from_distmod
1.52±0.1ms	1.47±0.1ms	0.96	benchmarks.TimeSuite.time_apply_passbands
5.96±0.9ms	5.45±1ms	0.91	benchmarks.TimeSuite.time_chained_evaluate

Click here to view all benchmarks.

[mi-dai]

This is good for now. We should allow users to fill the out-of-bound values the way they want in the future.