Open Access. The submitted version of this work has been reproduced here according to the publisher's polices.
NOTICE: this is the author's version of a work that was submitted for publication in Applied Numerical Mathematics. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. A definitive version of this article was published in Applied Numerical Mathematics, 109: 123-137, 2016, and can be found at http://dx.doi.org/10.1016/j.apnum.2016.06.010 (Please note: access via this link may only be available with a subscription)
Copyright (2016) IMACS. Published by Elsevier B.V.