Journal of Scientific Computing, 32(3) (2007), 411-424 http://dx.doi.org/10.1007/s10915-007-9141-0