Lorentzian 3d gravity with wormholes via matrix models