StructuredPenaltyOperator.quad_form()#

StructuredPenaltyOperator.quad_form(x)[source]#

Efficiently computes a quadratic form x.T @ P @ x.

P is the overall precision matrix of the tensor product. The quadratic form is computed without actually materializing the full precision matrix by exploiting its structured composition of marginal precisions.