[ase-users] NEB code
Peterson, Andrew
Andrew_Peterson at brown.edu
Fri Aug 23 21:42:53 CEST 2019
You probably saw the line
<https://gitlab.com/ase/ase/blob/master/ase/neb.py#L315> f = forces[i - 1],
and then f is modified. That first line returns a “view” of the original
array, so when f is modified so is the original forces array. It's just a
numpy way of dealing with slices.
On Tue, Aug 20, 2019 at 11:37 AM Alaithan, Zainab via ase-users <
ase-users at listserv.fysik.dtu.dk> wrote:
> I am examining NEB code because I am facing some convergence problems when
> I increase the number of images from 3. I noticed that NEB get.forces
> method returns the real forces instead of the modified forces.
>
>
> The modified forces are calculated as 'f'
> but the method returns 'forces' rather than 'f'
>
>
> am I missing something?
>
>
> Thank you,
>
>
> Regards,
> _______________________________________________
> ase-users mailing list
> ase-users at listserv.fysik.dtu.dk
> https://listserv.fysik.dtu.dk/mailman/listinfo/ase-users
--
Andrew Peterson
Associate Professor
Brown University School of Engineering
Barus & Holley 247
184 Hope Street
Providence, RI 02912
(401) 863-2153
http://brown.edu/go/catalyst
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://listserv.fysik.dtu.dk/pipermail/ase-users/attachments/20190823/9ccd0bf2/attachment.html>
More information about the ase-users
mailing list