Is there a standard solution for Gauss elimination in Python?

by

I finally found, that it can be done using LU decomposition. Here the U matrix represents the reduced form of the linear system.

from numpy import array
from scipy.linalg import lu

a = array([[2.,4.,4.,4.],[1.,2.,3.,3.],[1.,2.,2.,2.],[1.,4.,3.,4.]])

pl, u = lu(a, permute_l=True)

Then u reads

array([[ 2.,  4.,  4.,  4.],
       [ 0.,  2.,  1.,  2.],
       [ 0.,  0.,  1.,  1.],
       [ 0.,  0.,  0.,  0.]])

Depending on the solvability of the system this matrix has an upper triangular or trapezoidal structure. In the above case a line of zeros arises, as the matrix has only rank 3.

More Related Contents:

  1. Difference between numpy.array shape (R, 1) and (R,)
  2. What are the differences between numpy arrays and matrices? Which one should I use?
  3. Very large matrices using Python and NumPy
  4. numpy matrix vector multiplication [duplicate]
  5. Numpy matrix to array
  6. how does multiplication differ for NumPy Matrix vs Array classes?
  7. Bin elements per row – Vectorized 2D Bincount for NumPy
  8. Convolve2d just by using Numpy
  9. Numpy ‘smart’ symmetric matrix
  10. Convert a 1D array to a 2D array in numpy
  11. numpy get index where value is true
  12. Deprecation status of the NumPy matrix class
  13. Efficiently Calculating a Euclidean Distance Matrix Using Numpy
  14. Block tridiagonal matrix python
  15. Iterating over Numpy matrix rows to apply a function each?
  16. how to perform max/mean pooling on a 2d array using numpy
  17. How to get element-wise matrix multiplication (Hadamard product) in numpy?
  18. Why does numpy.linalg.solve() offer more precise matrix inversions than numpy.linalg.inv()?
  19. NumPy array/matrix of mixed types
  20. Replace sub part of matrix by another small matrix in numpy
  21. How to convert a column or row matrix to a diagonal matrix in Python?
  22. How to find linearly independent rows from a matrix
  23. Convert row vector to column vector in NumPy
  24. Numpy select rows based on condition
  25. How to find the pairwise differences between rows of two very large matrices using numpy?
  26. Multiple matrix multiplication
  27. A numpy array unexpectedly changes when changing another one despite being separate
  28. How to print the full NumPy array, without truncation?
  29. Load numpy array in google-cloud-ml job
  30. How can I serialize a numpy array while preserving matrix dimensions?

Leave a Comment