bayespecon._ops.kron_solve_matrix¶

bayespecon._ops.kron_solve_matrix(Lo, Ld, B, n)[source]¶

Solve \((L_o \otimes L_d)\,H = B\) for a matrix RHS via batched two-step solve.

Applies the same Kronecker algorithm as kron_solve_vec() to all k columns of B simultaneously using a single \(L_d\) factorisation and a single \(L_o^\top\) factorisation (both of size \(n \times n\)).

Parameters:¶

Lo : scipy.sparse.csr_matrix, shape (n, n)¶

Ld : scipy.sparse.csr_matrix, shape (n, n)¶

B¶

n : int¶

Returns:¶

H

Return type:¶

ndarray, shape (N, k)