Computes the missing gradients for QM/MM linkatoms that is bonded terms on QM atoms, LA forces projection on QM and MM atoms. To obtain the correct forces in output, qmg should already contain the QM forces, so that LA forces could be projected on QM and MM force vectors
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
type(ommp_qm_helper), | intent(in) | :: | qm | |||
type(ommp_system), | intent(inout) | :: | mm | |||
real(kind=rp), | intent(inout) | :: | qmg(3,qm%qm_top%mm_atoms) | |||
real(kind=rp), | intent(inout) | :: | mmg(3,mm%top%mm_atoms) | |||
real(kind=rp), | intent(in) | :: | original_qmg(3,qm%qm_top%mm_atoms) |