You can always define the sparse matrix dimension in this way:<br>scipy.sparse.coo_matrix(([5.], [[7.],[7.]]), dims=(16, 16))<br><br>L.<br><br><div class="gmail_quote">On Dec 14, 2007 8:04 PM, Nils Wagner &lt;<a href="mailto:nwagner@iam.uni-stuttgart.de">
nwagner@iam.uni-stuttgart.de</a>&gt; wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi all,<br><br>I use the coo constructor to build sparse matrices.
<br>A_i = coo_matrix( (vals_i, (rows_i, cols_i)) )<br><br>I read the data for each matrix from a different file.<br><br>One of the matrices is rank-one with only non-zero entry.<br><br>(7, 7) &nbsp; &nbsp; &nbsp; &nbsp;5.0<br><br>Another matrix has the entries
<br><br> &nbsp; (0, 0) &nbsp; &nbsp; &nbsp; &nbsp;80000000.0<br> &nbsp; (0, 1) &nbsp; &nbsp; &nbsp; &nbsp;-120000000.0<br> &nbsp; (1, 0) &nbsp; &nbsp; &nbsp; &nbsp;-120000000.0<br> &nbsp; (1, 1) &nbsp; &nbsp; &nbsp; &nbsp;480000000.0<br> &nbsp; (0, 2) &nbsp; &nbsp; &nbsp; &nbsp;40000000.0<br> &nbsp; (2, 0) &nbsp; &nbsp; &nbsp; &nbsp;40000000.0<br> &nbsp; (2, 2) &nbsp; &nbsp; &nbsp; &nbsp;160000000.0
<br> &nbsp; (1, 3) &nbsp; &nbsp; &nbsp; &nbsp;-240000000.0<br> &nbsp; (3, 1) &nbsp; &nbsp; &nbsp; &nbsp;-240000000.0<br> &nbsp; (2, 3) &nbsp; &nbsp; &nbsp; &nbsp;-120000000.0<br> &nbsp; (3, 2) &nbsp; &nbsp; &nbsp; &nbsp;-120000000.0<br> &nbsp; (3, 3) &nbsp; &nbsp; &nbsp; &nbsp;480000000.0<br> &nbsp; (1, 4) &nbsp; &nbsp; &nbsp; &nbsp;120000000.0<br> &nbsp; (4, 1) &nbsp; &nbsp; &nbsp; &nbsp;120000000.0
<br> &nbsp; (2, 4) &nbsp; &nbsp; &nbsp; &nbsp;40000000.0<br> &nbsp; (4, 2) &nbsp; &nbsp; &nbsp; &nbsp;40000000.0<br> &nbsp; (4, 4) &nbsp; &nbsp; &nbsp; &nbsp;160000000.0<br> &nbsp; (3, 5) &nbsp; &nbsp; &nbsp; &nbsp;-240000000.0<br> &nbsp; (5, 3) &nbsp; &nbsp; &nbsp; &nbsp;-240000000.0<br> &nbsp; (4, 5) &nbsp; &nbsp; &nbsp; &nbsp;-120000000.0<br> &nbsp; (5, 4) &nbsp; &nbsp; &nbsp; &nbsp;-120000000.0
<br> &nbsp; (5, 5) &nbsp; &nbsp; &nbsp; &nbsp;480000000.0<br> &nbsp; (3, 6) &nbsp; &nbsp; &nbsp; &nbsp;120000000.0<br> &nbsp; (6, 3) &nbsp; &nbsp; &nbsp; &nbsp;120000000.0<br> &nbsp; (4, 6) &nbsp; &nbsp; &nbsp; &nbsp;40000000.0<br> &nbsp; : &nbsp; &nbsp; :<br> &nbsp; (9, 11) &nbsp; &nbsp; &nbsp; -240000000.0<br> &nbsp; (11, 9) &nbsp; &nbsp; &nbsp; -240000000.0<br> &nbsp; (10, 11) &nbsp; &nbsp; &nbsp;-
120000000.0<br> &nbsp; (11, 10) &nbsp; &nbsp; &nbsp;-120000000.0<br> &nbsp; (11, 11) &nbsp; &nbsp; &nbsp;480000000.0<br> &nbsp; (9, 12) &nbsp; &nbsp; &nbsp; 120000000.0<br> &nbsp; (12, 9) &nbsp; &nbsp; &nbsp; 120000000.0<br> &nbsp; (10, 12) &nbsp; &nbsp; &nbsp;40000000.0<br> &nbsp; (12, 10) &nbsp; &nbsp; &nbsp;40000000.0<br> &nbsp; (12, 12) &nbsp; &nbsp; &nbsp;
160000000.0<br> &nbsp; (11, 13) &nbsp; &nbsp; &nbsp;-240000000.0<br> &nbsp; (13, 11) &nbsp; &nbsp; &nbsp;-240000000.0<br> &nbsp; (12, 13) &nbsp; &nbsp; &nbsp;-120000000.0<br> &nbsp; (13, 12) &nbsp; &nbsp; &nbsp;-120000000.0<br> &nbsp; (13, 13) &nbsp; &nbsp; &nbsp;480000000.0<br> &nbsp; (11, 14) &nbsp; &nbsp; &nbsp;120000000.0<br> &nbsp; (14, 11) &nbsp; &nbsp; &nbsp;
120000000.0<br> &nbsp; (12, 14) &nbsp; &nbsp; &nbsp;40000000.0<br> &nbsp; (14, 12) &nbsp; &nbsp; &nbsp;40000000.0<br> &nbsp; (14, 14) &nbsp; &nbsp; &nbsp;160000000.0<br> &nbsp; (13, 15) &nbsp; &nbsp; &nbsp;120000000.0<br> &nbsp; (15, 13) &nbsp; &nbsp; &nbsp;120000000.0<br> &nbsp; (14, 15) &nbsp; &nbsp; &nbsp;40000000.0<br> &nbsp; (15, 14) &nbsp; &nbsp; &nbsp;
40000000.0<br> &nbsp; (15, 15) &nbsp; &nbsp; &nbsp;80000000.0<br><br><br>How can I circumvent the<br>ValueError: shape mismatch: objects cannot be broadcast to<br>a single shape<br><br>if I try to add the matrices ?<br><br>Nils<br><br><br>_______________________________________________
<br>SciPy-user mailing list<br><a href="mailto:SciPy-user@scipy.org">SciPy-user@scipy.org</a><br><a href="http://projects.scipy.org/mailman/listinfo/scipy-user" target="_blank">http://projects.scipy.org/mailman/listinfo/scipy-user
</a><br></blockquote></div><br>