- [X] Faster linear solves with `Lineax`
- [X] Support for both `vals` and `coeffs` as fields inputs to differential and rhs operators
- [X] Fixed the `args` inputs to construct the local matrix for the `nodal_div_grad`: (if array, else, etc.)
- [X] Added colorbar to animate fields
- [X] Implemented the Darcy flow problem in demos