[SciPy-Dev] Sylvester equation solver for complex matrices
Wed Feb 15 15:44:56 CST 2012
On Wed, 15 Feb 2012 12:42:04 -0600
Warren Weckesser <firstname.lastname@example.org> wrote:
> On Sat, Feb 11, 2012 at 11:24 AM, Warren Weckesser <
> email@example.com> wrote:
>> On Sat, Feb 11, 2012 at 4:45 AM, Nils Wagner
>> > wrote:
>>> Hi all,
>>> I run a short test on the recent solver for Sylvester
>>> In case of complex input matrices the method returns
>>> ** On entry to ZTRSYL parameter number 2 had an
>> Thanks, Nils. The new solve_sylvester function didn't
>> matrices correctly. I have submitted a pull request (
>> https://github.com/scipy/scipy/pull/155) that should fix
> Nils, the fix has been committed to the master branch.
> Let us know if you
> find any other problems.
>>> import numpy as np
>>> from scipy import linalg as la
>>> a = np.random.rand(20,20)+1j*np.random.rand(20,20)
>>> b = np.random.rand(20,20)+1j*np.random.rand(20,20)
>>> q = np.random.rand(20,20)+1j*np.random.rand(20,20)
>>> x = la.solve_sylvester(a,b,q)
>>> res = la.norm(np.dot(a,x)+np.dot(x,b)-q)
>>> SciPy-Dev mailing list
Great. It works fine for me.
BTW, is the Sylvester equation solver part of the next
scipy release ?
How about Lyapunov and Riccati equation solvers ?
More information about the SciPy-Dev