mult_grid function

Kronecker matrix product functionality