Spectral IDM: speed up by vectorizing and using matrices

We can use odespec also for matrix ODEs thus we should be able to remove the loop over the rigid coordinates. And maybe also over the elastic coordinates.