def compute_rmsd(a, b):
return np.sqrt(((a-b)**<NUM_LIT:2>).mean())<EOL>
Compute the root-mean-square deviation between two arrays

       Arguments:
         a, b  --  Two numpy arrays with the same shape