@cached<EOL><INDENT>def inertia_tensor(self):<DEDENT>
result = np.zeros((<NUM_LIT:3>,<NUM_LIT:3>), float)<EOL>for i in range(self.size):<EOL><INDENT>r = self.coordinates[i] - self.com<EOL>result.ravel()[::<NUM_LIT:4>] += self.masses[i]*(r**<NUM_LIT:2>).sum()<EOL>result -= self.masses[i]*np.outer(r,r)<EOL><DEDENT>return result<EOL>
the intertia tensor of the molecule