Current Status:

• 1. Works with Periodic and Isolated (Dirichlet) boundary conditions.
• 2. Use work array to work with MG-level Structure.