Distributed and vectorized method of characteristics for fast transient simulations in water distribution systems