Solvers support matrix-free operations
Addresses: #2 (closed), but does not fully implement it, yet as the operators still need to become matrix-free
Addresses: #2 (closed), but does not fully implement it, yet as the operators still need to become matrix-free