A computationally efficient and mechanically compatible multi-phase-field model applied to coherently stressed three-phase solids