Dear Dr. Case,
Thank you very much for the kind explanation. My research project requires
multiplying a three- or four-body terms to a LJ pairwise term; therefore, it
is critical for me to understand how the indices of atoms and their x y z
components are assigned near a pair of atoms.
So far I have not identified the patterns in this code. Could you please
help me? Thank you.
Sincerely,
Antonios
On Fri, Jun 11, 2010 at 7:48 PM, case <case.biomaps.rutgers.edu> wrote:
> On Fri, Jun 11, 2010, Antonios Samiotakis wrote:
>
> > I don't understand what takes place in the ew_directp.h file.
> >
> >
> > ! prologue: gather the data and put it in temporary arrays.
> >
> > n=ipairs(m)
> > itran=ishft(n,-27)
> > n = iand(n,mask27)
> > j = bckptr(n)
> >
> > delx = imagcrds(1,n) + xktran(1,itran)
> > dely = imagcrds(2,n) + xktran(2,itran)
> > delz = imagcrds(3,n) + xktran(3,itran)
> > delr2 = delx*delx + dely*dely+delz*delz
>
> Once you have delx, dely, delz, it should be straightforward to compute
> any function of r, and its derivatives. The "epilogue" file shows the
> actual
> compution of evdw (the contribution of the LJ 12/6 terms to the energy) and
> the update to the force array.
>
> It may be more difficult than this, depending on exactly what you want to
> do,
> but some of the internal manipulations are probably not relevant (?)
>
> ....dac
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
--
Antonios Samiotakis
Department of Physics
University of Houston
4800 Calhoun Rd
SR1 Rm 629G, Houston
TX, 77204
e-mail: asamiotakis2_at_uh.edu
http://www.phys.uh.edu/~mscheung/Antonios.CV/Resume10.2.website.htm
Tel: 713-743-3725
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sat Jun 12 2010 - 08:30:04 PDT
