<br><br><div class="gmail_quote">On Sat, Feb 11, 2012 at 4:45 AM, Nils Wagner <span dir="ltr"><<a href="mailto:nwagner@iam.uni-stuttgart.de">nwagner@iam.uni-stuttgart.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

Hi all,<br>
<br>
I run a short test on the recent solver for Sylvester<br>
equations.<br>
In case of complex input matrices the method returns<br>
  ** On entry to ZTRSYL parameter number  2 had an illegal<br>
value<br>
<br>
Nils<br>
<br></blockquote><div><br><br>Thanks, Nils.  The new solve_sylvester function didn't handle complex matrices correctly. I have submitted a pull request (<a href="https://github.com/scipy/scipy/pull/155">https://github.com/scipy/scipy/pull/155</a>) that should fix this.<br>

<br>Warren<br><br><br></div><blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
import numpy as np<br>
from scipy import linalg as la<br>
np.random.seed(10)<br>
a = np.random.rand(20,20)+1j*np.random.rand(20,20)<br>
b = np.random.rand(20,20)+1j*np.random.rand(20,20)<br>
q = np.random.rand(20,20)+1j*np.random.rand(20,20)<br>
<br>
x = la.solve_sylvester(a,b,q)<br>
res = la.norm(np.dot(a,x)+np.dot(x,b)-q)<br>
_______________________________________________<br>
SciPy-Dev mailing list<br>
<a href="mailto:SciPy-Dev@scipy.org">SciPy-Dev@scipy.org</a><br>
<a href="http://mail.scipy.org/mailman/listinfo/scipy-dev" target="_blank">http://mail.scipy.org/mailman/listinfo/scipy-dev</a><br>
</blockquote></div><br>