I made a typing error in the previous post, the first code segment is this one instead:<div><br></div><div><div style="border-collapse: collapse; font-family: arial, sans-serif; font-size: 13px; ">Asym = matrix(dot(A.T,A))</div>

<div style="border-collapse: collapse; font-family: arial, sans-serif; font-size: 13px; ">bsym = matrix(dot(A.T,b))</div><div style="border-collapse: collapse; font-family: arial, sans-serif; font-size: 13px; ">sol = cg(Asym,bsym,tol = 1e-10,maxiter=30)</div>

<div style="border-collapse: collapse; font-family: arial, sans-serif; font-size: 13px; "><br></div><div style="border-collapse: collapse; font-family: arial, sans-serif; font-size: 13px; "><span class="Apple-style-span" style="font-family: arial; font-size: small; border-collapse: separate; "><br>

</span></div><div style="border-collapse: collapse; font-family: arial, sans-serif; font-size: 13px; "><span class="Apple-style-span" style="font-family: arial; font-size: small; border-collapse: separate; ">On Thu, Dec 9, 2010 at 9:12 PM, Raul Acuña <span dir="ltr">&lt;<a href="mailto:raultron@gmail.com">raultron@gmail.com</a>&gt;</span> wrote:</span></div>

<div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Hi,<div><br></div><div>Am using the iterative methods of scipy.sparse.linalg for solving a linear system of equations Ax = b. My matrix A is non symmetric. I&#39;ve been using the <a href="http://scipy.sparse.linalg.cg" target="_blank">scipy.sparse.linalg.cg</a>() function multiplying both matrix &quot;A&quot; and &quot;b&quot; with the transpose of A so the matrix will become symmetric:</div>


<div><br></div><div>Asym = matrix(dot(A.T,A))</div><div>bsym = matrix(dot(A.T,b))</div><div>sol = cg(A,b,tol = 1e-10,maxiter=30)</div><div><br></div><div>Also i&#39;ve been reading about the biconjugate gradient method, and if i am not mistaken the literature says that this method works on non-symmetric matrix, but when i try to use  scipy.sparse.linalg.bicg() it wont work:</div>


<div><br></div><div><div>sol = bicg(A,b,tol = 1e-10,maxiter=30)</div><div><br></div><div>   File &quot;C:\Python26\lib\site-packages\scipy\sparse\linalg\isolve\iterative.py&quot;, line 74, in bicg</div><div>          A,M,x,b,postprocess = make_system(A,M,x0,b,xtype)</div>


<div>   File &quot;C:\Python26\lib\site-packages\scipy\sparse\linalg\isolve\utils.py&quot;, line 65, in make_system</div><div>          raise ValueError(&#39;expected square matrix (shape=%s)&#39; % shape)</div><div>          NameError: global name &#39;shape&#39; is not defined</div>


</div><div><br></div><div><div>Any help will be greatly appreciated, am comparing this methods for my data with an important emphasis on speed for my master thesis, so any discrepancy with the theory would be a great problem for me.
</div></div><div><br></div><div>Thanks in advance,</div><div><br></div><div><br></div><div>Raúl Acuña.</div>
</blockquote></div><br><br clear="all"><br>-- <br><div>Ing. Raúl Acuña</div><div><br></div>Profesor @Universidad Simón Bolívar <div><div>Grupo de Mecatrónica</div><div>Departamento de Electrónica y Circuitos<br><div>Tel. +58-212-4121983 / Cel +58-412-5840317</div>

<div><br></div></div></div><br>
</div>