15.12.2012 Views

scipy tutorial - Baustatik-Info-Server

scipy tutorial - Baustatik-Info-Server

scipy tutorial - Baustatik-Info-Server

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

solveh_banded(ab, b, overwrite_ab=0, overwrite_b=0, lower=0)<br />

Solve equation a x = b. a is Hermitian positive-definite banded matrix.<br />

The matrix a is stored in ab either in lower diagonal or upper diagonal ordered form:<br />

SciPy Reference Guide, Release 0.8.dev<br />

ab[u + i - j, j] == a[i,j] (if upper form; i = j)<br />

Example of ab (shape of a is (6,6), u=2):<br />

upper form:<br />

* * a02 a13 a24 a35<br />

* a01 a12 a23 a34 a45<br />

a00 a11 a22 a33 a44 a55<br />

lower form:<br />

a00 a11 a22 a33 a44 a55<br />

a10 a21 a32 a43 a54 *<br />

a20 a31 a42 a53 * *<br />

Cells marked with * are not used.<br />

Parameters<br />

ab : array, shape (M, u + 1)<br />

Banded matrix<br />

b : array, shape (M,) or (M, K)<br />

Right-hand side<br />

overwrite_ab : boolean<br />

Discard data in ab (may enhance performance)<br />

overwrite_b : boolean<br />

Discard data in b (may enhance performance)<br />

lower : boolean<br />

Is the matrix in the lower form. (Default is upper form)<br />

Returns<br />

c : array, shape (M, u+1)<br />

Cholesky factorization of a, in the same banded format as ab<br />

x : array, shape (M,) or (M, K)<br />

The solution to the system a x = b<br />

det(a, overwrite_a=0)<br />

Compute the determinant of a matrix<br />

Parameters<br />

a : array, shape (M, M)<br />

Returns<br />

det : float or complex<br />

Determinant of a<br />

3.7. Linear algebra (<strong>scipy</strong>.linalg) 219

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!