Went back to old interface for the matvec routine that does the real work, will make it easier for JIT syntax later