I have a question about lines 115-120 of angle_harmonic.cpp. I would have thought the gradients would have been a series of cross products of the bond vectors: r12 x (r32 x r12) appropriately normalized multiplied by the force constants etc. . However this is something very different. I tried checking the manual for a reference but there was none. What is going on here?