add a method to presolve iterative system

make this the default because multiplication is much faster than
elimination/decomposition