login  home  contents  what's new  discussion  bug reports     help  links  subscribe  changes  refresh  edit

LUDecomposition (LUD)

LU decomposition for ordinary matrices.

fricas
(1) -> )expose LUDecomposition
LUDecomposition is now explicitly exposed in frame initial A:=matrix [[subscript('a,[10*i+j]) for i in 1..3] for j in 1..3]

\label{eq1}\left[ 
\begin{array}{ccc}
{a_{11}}&{a_{21}}&{a_{31}}
\
{a_{12}}&{a_{22}}&{a_{32}}
\
{a_{13}}&{a_{23}}&{a_{33}}
(1)
Type: Matrix(Polynomial(Integer))
fricas
diagProduct(x) == reduce(*,[x(i,i) for i in 1..nrows(x)])
Type: Void
fricas
B:=LUDecomp A;
Type: Record(LU: Matrix(Fraction(Polynomial(Integer))),Perm: Vector(Integer),Pivots: List(Fraction(Polynomial(Integer))))
fricas
B.LU

\label{eq2}\left[ 
\begin{array}{ccc}
{a_{11}}&{a_{21}}&{a_{31}}
\
{\frac{a_{12}}{a_{11}}}&{\frac{{{a_{11}}\ {a_{22}}}-{{a_{12}}\ {a_{21}}}}{a_{11}}}&{\frac{{{a_{11}}\ {a_{32}}}-{{a_{12}}\ {a_{3
1}}}}{a_{11}}}
\
{\frac{a_{13}}{a_{11}}}&{\frac{{{a_{11}}\ {a_{23}}}-{{a_{13}}\ {a_{21}}}}{{{a_{11}}\ {a_{22}}}-{{a_{12}}\ {a_{21}}}}}&{\frac{{{\left({{a_{11}}\ {a_{22}}}-{{a_{12}}\ {a_{21}}}\right)}\ {a_{33}}}+{{\left(-{{a_{11}}\ {a_{23}}}+{{a_{13}}\ {a_{21}}}\right)}\ {a_{32}}}+{{\left({{a_{12}}\ {a_{23}}}-{{a_{13}}\ {a_{22}}}\right)}\ {a_{31}}}}{{{a_{11}}\ {a_{22}}}-{{a_{12}}\ {a_{21}}}}}
(2)
Type: Matrix(Fraction(Polynomial(Integer)))
fricas
B.Perm

\label{eq3}\left[ 1, \: 2, \: 3 \right](3)
Type: Vector(Integer)
fricas
B.Pivots

\label{eq4}\begin{array}{@{}l}
\displaystyle
\left[{\frac{{{\left({{a_{11}}\ {a_{22}}}-{{a_{12}}\ {a_{21}}}\right)}\ {a_{33}}}+{{\left(-{{a_{11}}\ {a_{23}}}+{{a_{13}}\ {a_{21}}}\right)}\ {a_{32}}}+{{\left({{a_{12}}\ {a_{23}}}-{{a_{13}}\ {a_{22}}}\right)}\ {a_{31}}}}{{{a_{11}}\ {a_{22}}}-{{a_{12}}\ {a_{21}}}}}, \right.
\
\
\displaystyle
\left.\:{\frac{{{a_{11}}\ {a_{22}}}-{{a_{12}}\ {a_{21}}}}{a_{1
1}}}, \:{a_{11}}\right] 
(4)
Type: List(Fraction(Polynomial(Integer)))
fricas
diagProduct(B.LU)=determinant A
fricas
Compiling function diagProduct with type Matrix(Fraction(Polynomial(
      Integer))) -> Fraction(Polynomial(Integer))

\label{eq5}\begin{array}{@{}l}
\displaystyle
{
\begin{array}{@{}l}
\displaystyle
{{\left({{a_{11}}\ {a_{22}}}-{{a_{12}}\ {a_{21}}}\right)}\ {a_{33}}}+{{\left(-{{a_{11}}\ {a_{23}}}+{{a_{13}}\ {a_{21}}}\right)}\ {a_{32}}}+ 
\
\
\displaystyle
{{\left({{a_{12}}\ {a_{23}}}-{{a_{13}}\ {a_{22}}}\right)}\ {a_{31}}}
(5)
Type: Equation(Fraction(Polynomial(Integer)))
fricas
%::Boolean

\label{eq6} \mbox{\rm true} (6)
Type: Boolean




  Subject:   Be Bold !!
  ( 15 subscribers )  
Please rate this page: