Calculation of tunneling matrix elements in rigid systems: mixed-valence dithiaspirocyclobutane molecules