Improve RMSD function to take molecular symmetry into consideration?

[ravila4]

Just wondering what the algorithm is behind the RMSD calculation.

[rasbt]

It should be the "usual" RMSD equation, i.e.

(https://en.wikipedia.org/wiki/Root-mean-square_deviation_of_atomic_positions)

Here's a link to the actual code:
https://github.com/rasbt/biopandas/blob/master/biopandas/pdb/pandas_pdb.py#L152

and the relevant lines are

        total = ((df1['x_coord'] - df2['x_coord'])**2 +
                 (df1['y_coord'] - df2['y_coord'])**2 +
                 (df1['z_coord'] - df2['z_coord'])**2)
        rmsd = round((total.sum() / df1.shape[0])**0.5, 4)
        return rmsd