Closed kjs73 closed 9 years ago
Just to be sure: Is it intended that some of the harmonics here compute the actual gradient and others -g/r?
Just to be sure: Is it intended that some of the harmonics here compute the actual gradient and others -g/r?
I don't understand what you mean. In the pair interaction formulation, the interaction classes define a function energy_gradient which should return g/r. A BasePotential class of course needs to compute the actual gradient.
Great, thank you for the explanation.
OK, I think that is as good (or bad) as it gets for now.
thanks. I'll merge it then.
Mostly adds some tests. Also replaces some explicit for loops by pele::dot function.