Matrix class is inefficient #10
Description
The Matrix class implemented in #8 uses a NumPy array to represent its data and it uses numpy.linalg.solve to solve linear systems. We would be better using petsc4py.
{{ message }}