boundlab.ibp.matmul#
Interval matrix products.
The bilinear form splits over components:
\[(c_1 + e_1)(c_2 + e_2) = c_1 c_2 + c_1 e_2 + e_1 c_2 + e_1 e_2 .\]
The three terms with a constant factor are linear maps (exact); only the error-error product needs an interval bound, \([-a, a] \cdot [-b, b] \subseteq [-(|a| \cdot |b|),\, |a| \cdot |b|]\).
Classes
Peel the |
|
Interval bound for the error-error product: \([-a, a] @ [-b, b] \subseteq \pm(a @ b)\) for \(a, b \ge 0\), with the reasons blended by mass. |