'Scalar parameters should broadcast with an array input to result in an
array output of the same shape as the input.'

def test_scalar_parameters_3d_array_input(self):

    t = TModel_1_2(1, 10, 1000)
    (y, z) = t((np.arange(12).reshape(2, 3, 2) * 100))
    assert isinstance(y, np.ndarray)
    assert isinstance(z, np.ndarray)
    assert (np.shape(y) == np.shape(z) == (2, 3, 2))
    assert np.all((y == [[[11, 111], [211, 311], [411, 511]], [[611, 711], [811, 911], [1011, 1111]]]))
    assert np.all((z == (y + 1000)))
