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

MatmulBiased

Peel the Bias centers off both operands.

MatmulNoise

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.