Ilhan Polat

Is there any progress on this feature? Anything along these lines would be beneficial for my use case. For what is worth, my use case is a little more convoluted...

Since it is pinged let me put an amateurish solution I've did for the example above. Since things are logarithmically plotted it's a bit annoying to track the plot but...

@X-Raym Would you care to open a Pull Request for a more general audience? Then interested people can review and leave comments. Also your code can be tested against the...

I have basically two comments one is do we really want to add this. Because typically blowups are related to the stability of matrix and not the norm bounds and...

Typically for evolution in time you might want to use a state space approach that is what simulink does; that is pick a timestep and do `dot(x) = Ax` and...

Yes because your array is not Hurwitz meaning you have eigenvalues with positive real parts. So this is not something a code can fix but as you said this is...

> Would you be able to help me implement this for solve_triangular? The idea was to finish the cursed #12824 and add nD support to all linalg, but, for some...

The point of expm overhaul was to gain speed from the internal algorithm and nD array capability was a bonus, because scipy.sparse.linalg.expm was used for the dense input too. The...

If it was a wrapper issue we would probably see a segfault rather than precision problem And precision violation is ever so little increased . So I am cautiously looking...

We need to provide Python wrappers for this which would take some time but in the meantime you can use `scipy.linalg.cython_lapack.dtgsyl` within Cython environment if you know your way around...