On Wed 18 Dec 2013 12:39:05 GMT, Jaime wrote:
> In scipy/signal/signaltools.py:
>
> > @@ -1609,6 +1609,10 @@ def lfilter_zi(b, a):
> >      a = np.atleast_1d(a)
> >      if a.ndim != 1:
> >          raise ValueError("Denominator a must be rank 1.")
> > +
> > +    typeout = (a[0] * b[0]).dtype
>
> It is probably better to call np.result_type
> <http://docs.scipy.org/doc/numpy/reference/generated/numpy.result_type.html>
> to do type promotion.
>
> —
> Reply to this email directly or view it on GitHub
> <https://github.com/scipy/scipy/pull/3159/files#r8434451>.
>

Good suggestion - I haven't come across that function before