leeping / forcebalance

Systematic force field optimization.
Other
145 stars 75 forks source link

Computing FD gradient of angle with OpenMM 7.7.0 #258

Open SimonBoothroyd opened 2 years ago

SimonBoothroyd commented 2 years ago

Description

OpenMM 7.7.0 introduced a number of checks to validate user inputs, including that all harmonic angle force angles should be between 0 and pi. This causes an error however when computing the finite difference gradient of an angle parameter with an angle of 180.0 as 180.0 + h > pi.