An Operator-integration-factor splitting method for time-dependent problems: Application to incompressible fluid flow