Strassen-Algorithmus

Strassen-Algorithmus für eine symmetrische Matrix:

Für eine symmetrische Matrix [math]\displaystyle{ A \in \mathbb{C}^{n \times n} }[/math] mit [math]\displaystyle{ n \in \mathbb{N}^* }[/math] beträgt die Laufzeit [math]\displaystyle{ T_s(n) }[/math] des Strassen-Algorithmus für das Matrixprodukt [math]\displaystyle{ A^2 }[/math] circa die Hälfte von der des Originalalgorithmus in [math]\displaystyle{ \mathcal{O}(n^{(_2 7)}) }[/math].

Beweis: Mit [math]\displaystyle{ A := \begin{pmatrix} A_{11} & A_{12} \\ A_{12}^T & A_{22} \end{pmatrix} }[/math] gilt [math]\displaystyle{ A^TA = \begin{pmatrix} A_{11}A_{11}+A_{12}A_{12}^T & A_{11}A_{12}+A_{12}A_{22} \\ A_{12}^TA_{11}+A_{22}A_{12}^T & A_{12}^TA_{12}+A_{22}A_{22} \end{pmatrix} }[/math] und [math]\displaystyle{ T_s(2n) = 3T_s(n) + 2n^{(_2 7)} }[/math]. Also [math]\displaystyle{ T_s(n) = 3T_s(n/2) + 2(n/2)^{(_2 7)} }[/math] bzw. [math]\displaystyle{ T_s(n/2) = 3T_s(n/4) + 2(n/4)^{(_2 7)} }[/math].

Die geometrische Reihe liefert wegen [math]\displaystyle{ T_s(1) = 1 }[/math]: [math]\displaystyle{ T_s(n) = 27T_s(n/8) + 2/7n^{(_2 7)}(1+3/7 + (3/7)^2 + ...) = 3^{(_2n)} + 2/7n^{(_2 7)} (1-(3/7)^{(_2n)})/(1-3/7) }[/math] [math]\displaystyle{ = n^{(_2 3)} + \hat{2}(n^{(_2 7)}-n^{(_2 3)}) = \hat{2} (n^{(_2 3)} + n^{(_2 7)}).\square }[/math]

Strassen-Algorithmus für eine quadratrische Matrix:

Für eine quadratrische Matrix [math]\displaystyle{ A \in \mathbb{C}^{n \times n} }[/math] mit [math]\displaystyle{ n \in \mathbb{N}^* }[/math] beträgt die Laufzeit [math]\displaystyle{ T_q(n) }[/math] des Strassen-Algorithmus für das Matrixprodukt [math]\displaystyle{ A^TA }[/math] circa [math]\displaystyle{ 4/7 }[/math] von der des Originalalgorithmus in [math]\displaystyle{ \mathcal{O}(n^{(_2 7)}) }[/math].

Beweis: Mit [math]\displaystyle{ A := \begin{pmatrix} A_{11} & A_{12} \\ A_{21} & A_{22} \end{pmatrix} }[/math] gilt [math]\displaystyle{ A^TA = \begin{pmatrix} A_{11}^TA_{11}+A_{21}^TA_{21} & A_{11}^TA_{12}+A_{21}^TA_{22} \\ A_{12}^TA_{11}+A_{22}^TA_{21} & A_{12}^TA_{12}+A_{22}^TA_{22} \end{pmatrix} }[/math] und [math]\displaystyle{ T_q(2n) = 4T_s(n) + 2n^{(_2 7)} }[/math] bzw. [math]\displaystyle{ T_q(n) = 4T_s(n/2) + 2/7n^{(_2 7)} = 2/3n^{(_2 3)} + 4/7n^{(_2 7)}.\square }[/math]

Neuer Algorithmus für zwei quadratische Matrizen der folgenden Formen:

Für zwei [math]\displaystyle{ A, B \in \mathbb{C}^{n \times n} }[/math] mit [math]\displaystyle{ n \in \mathbb{N}^* }[/math] beträgt die Laufzeit [math]\displaystyle{ T_z(n) }[/math] für das Matrixprodukt [math]\displaystyle{ AB }[/math] des neuen Algorithmus circa [math]\displaystyle{ 6/7 }[/math] von der des Originalalgorithmus in [math]\displaystyle{ \mathcal{O}(n^{(_2 7)}) }[/math], wenn [math]\displaystyle{ A }[/math] ebenfalls die nachfolgende Form von [math]\displaystyle{ B }[/math] hat:

Beweis: Mit [math]\displaystyle{ B := \begin{pmatrix} B_{11} & B_{12} \\ B_{12} & B_{22} \end{pmatrix} }[/math] gilt [math]\displaystyle{ AB = \begin{pmatrix} A_{11}B_{11}+A_{12}B_{12} & A_{11}B_{12}+A_{12}B_{22} \\ A_{12}B_{11}+A_{22}B_{12} & A_{12}B_{12}+A_{22}B_{22} \end{pmatrix} =: C }[/math]. Wird

[math]\displaystyle{ M_{1} := A_{12} \cdot (B_{11} + B_{12}) }[/math]

[math]\displaystyle{ M_{2} := A_{12} \cdot (B_{12} + B_{22}) }[/math]

[math]\displaystyle{ M_{3} := A_{22} \cdot (B_{12} + B_{22}) }[/math]

[math]\displaystyle{ M_{4} := (A_{11} - A_{12})\cdot B_{11} }[/math]

[math]\displaystyle{ M_{5} := (A_{11} - A_{12})\cdot B_{12} }[/math]

[math]\displaystyle{ M_{6} := (A_{12} - A_{22})\cdot B_{22} }[/math]

gesetzt, so gilt

[math]\displaystyle{ C_{11} = M_{1} + M_{4} }[/math]

[math]\displaystyle{ C_{12} = M_{2} + M_{5} }[/math]

[math]\displaystyle{ C_{21} = M_{1} - C_{22} + M_{3} }[/math]

[math]\displaystyle{ C_{22} = M_{2} - M_{6} .\square }[/math]

Dies gilt auch, wenn [math]\displaystyle{ A }[/math] die folgende Form von [math]\displaystyle{ B }[/math] hat:

Mit [math]\displaystyle{ B := \begin{pmatrix} B_{11} & B_{12} \\ B_{21} & B_{11} \end{pmatrix} }[/math] gilt [math]\displaystyle{ AB = \begin{pmatrix} A_{11}B_{11}+A_{12}B_{21} & A_{11}B_{12}+A_{12}B_{11} \\ A_{21}B_{11}+A_{11}B_{21} & A_{21}B_{12}+A_{11}B_{11} \end{pmatrix} =: C }[/math]. Wird

[math]\displaystyle{ M_{1} := A_{11} \cdot (B_{11} + B_{12}) }[/math]

[math]\displaystyle{ M_{2} := A_{21} \cdot (B_{11} + B_{12}) }[/math]

[math]\displaystyle{ M_{3} := A_{12} \cdot (B_{11} - B_{21}) }[/math]

[math]\displaystyle{ M_{4} := (A_{11} - A_{12})\cdot B_{11} }[/math]

[math]\displaystyle{ M_{5} := (A_{11} - A_{21})\cdot B_{11} }[/math]

[math]\displaystyle{ M_{6} := (A_{11} + A_{12})\cdot B_{21} }[/math]

gesetzt, so gilt

[math]\displaystyle{ C_{11} = M_{3} + M_{4} }[/math]

[math]\displaystyle{ C_{12} = M_{1} - M_{4} }[/math]

[math]\displaystyle{ C_{21} = C_{11} - M_{5} + M_{6} }[/math]

[math]\displaystyle{ C_{22} = M_{2} + M_{5} .\square }[/math]

Bemerkung: Haben die Zerlegungen von [math]\displaystyle{ A }[/math] und [math]\displaystyle{ B }[/math] wieder die angegebenen Formen, ist die Laufzeit sogar [math]\displaystyle{ \mathcal{O}(n^{(_2 6)}) }[/math].

Siehe auch

Liste mathematischer Symbole

Siehe auch

Liste mathematischer Symbole