Expose solve method with new system matrix in pyamgcl