[gmx-developers] Extracting i-j forces from nb_generic.c

Mark Abraham mark.j.abraham at gmail.com
Mon Feb 3 02:28:42 CET 2014


Yeah, that's what I meant. Multiply thing by -1 or take some other thing to
the power -1, it's all the same right? :-P

Mark


On Sun, Feb 2, 2014 at 8:51 PM, Erik Lindahl <erik.lindahl at scilifelab.se>wrote:

> Hi,
>
> On 02 Feb 2014, at 20:11, Mark Abraham <mark.j.abraham at gmail.com> wrote:
> >
> > Do take care - components of fscal are not actually the scalar force,
> but the negative of it (because mostly it saves a flop).
>
> For the record, it can be even worse than this. In many (most) places we
> need to multply the force with the interparticle vector to get the
> vectorial force. However, it would be expensive (and stupid) computationall
> to normalize this vector when we already have 1/r available. Thus, fscal is
> often the scalar force *divided by r*, so we can then multiply it by the
> non-normalized interparticle vector.
>
> Cheers,
>
> Erik
>
> --
> Gromacs Developers mailing list
>
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List before
> posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-developersor send a mail to
> gmx-developers-request at gromacs.org.
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-developers/attachments/20140203/fa5a77e6/attachment.html>


More information about the gromacs.org_gmx-developers mailing list