Compute the the center of mass (COM) of the current configuration (stored in <tt>SR</tt>).
The result (although not accessible from the user) will be used in subsequent calls of functions
<tt>translate_com</tt>, <tt>rotate_com</tt>.

