15.12.2012 Views

scipy tutorial - Baustatik-Info-Server

scipy tutorial - Baustatik-Info-Server

scipy tutorial - Baustatik-Info-Server

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

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

Notes<br />

This is a LU factorization routine written for Scipy.<br />

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

Compute pivoted LU decomposition of a matrix.<br />

The decomposition is:<br />

A = P L U<br />

where P is a permutation matrix, L lower triangular with unit diagonal elements, and U upper triangular.<br />

See Also:<br />

Parameters<br />

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

Matrix to decompose<br />

overwrite_a : boolean<br />

Whether to overwrite data in A (may increase performance)<br />

Returns<br />

lu : array, shape (N, N)<br />

Matrix containing U in its upper triangle, and L in its lower triangle. The unit diagonal<br />

elements of L are not stored.<br />

piv : array, shape (N,)<br />

Pivot indices representing the permutation matrix P: row i of matrix was interchanged<br />

with row piv[i].<br />

lu_solve<br />

solve an equation system using the LU factorization of a matrix<br />

Notes<br />

This is a wrapper to the *GETRF routines from LAPACK.<br />

lu_solve((lu, piv), b, trans=0, overwrite_b=0)<br />

Solve an equation system, a x = b, given the LU factorization of a<br />

Parameters<br />

(lu, piv) :<br />

b : array<br />

Factorization of the coefficient matrix a, as given by lu_factor<br />

Right-hand side<br />

trans : {0, 1, 2}<br />

Returns<br />

x : array<br />

Type of system to solve:<br />

trans system<br />

0 a x = b<br />

1 a^T x = b<br />

2 a^H x = b<br />

234 Chapter 3. Reference

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

Saved successfully!

Ooh no, something went wrong!