include SimIO in solver class, residuum calculation