Solving linear system with sum of kronecker products

For 2D Poisson, I know how to deal with Dirichlet BCs, but could you provide simple explanation for how to deal with Neumann BCs with this tensor trick? Thank you.

Partial Neumann BCs (only one side in each direction) are easy but I do not know how to cope with full Neumann and Tensor Trick (explaining why I use GMG in
the fluid simulation example)…

Here is an example with partial Neumann (DDNN)

I did update the repo to fix issuez with PrettyTable and Makie evolutions.

2 Likes