ICODE-ADC/lammps
4
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Update minimization documentation to reflect line quadratic default

Browse Source 
git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@12533 f3b2605a-c512-4ea7-a41b-209d697bcdaa
This commit is contained in:
athomps
2014-09-24 00:48:02 +00:00
parent fb705efa84
commit e1d881b11c
1 changed files with 15 additions and 11 deletions
Download Patch File Download Diff File Expand all files Collapse all files

26
doc/min_modify.html
View File

@ -49,8 +49,20 @@ highly overlapped atoms from being moved long distances (e.g. through
another atom) due to large forces. another atom) due to large forces.
</P> </P>
<P>The choice of line search algorithm for the <I>cg</I> and <I>sd</I> minimization <P>The choice of line search algorithm for the <I>cg</I> and <I>sd</I> minimization
styles can be selected via the <I>line</I> keyword. The default styles can be selected via the <I>line</I> keyword.
backtracking search is robust and should always find a local energy The default <I>quadratic</I> line search algorithm starts out using
the robust backtracking method described below. However, once
the system gets close to a local
minimum and the linesearch steps get small, so that the energy
is approximately quadratic in the step length, it uses the
estimated location of zero gradient as the linesearch step,
provided the energy change is downhill.
This becomes more efficient than backtracking
for highly-converged relaxations. The <I>forcezero</I>
line search algorithm is similar to <I>quadratic</I>.
It may be more efficient than <I>quadratic</I> on some systems.
</P>
<P>The backtracking search is robust and should always find a local energy
minimum. However, it will "converge" when it can no longer reduce the minimum. However, it will "converge" when it can no longer reduce the
energy of the system. Individual atom forces may still be larger than energy of the system. Individual atom forces may still be larger than
desired at this point, because the energy change is measured as the desired at this point, because the energy change is measured as the
@ -58,14 +70,6 @@ difference of two large values (energy before and energy after) and
that difference may be smaller than machine epsilon even if atoms that difference may be smaller than machine epsilon even if atoms
could move in the gradient direction to reduce forces further. could move in the gradient direction to reduce forces further.
</P> </P>
<P>By contrast, the <I>quadratic</I> line search algorithm tries to
reduce the forces to zero, while guaranteeing that the energy
changes is not positive (uphill). For some systems, it may also
be more efficient than the backtracking algorithm by
requiring fewer energy/force evaluations. The <I>forcezero</I>
line search algorithm is similar to <I>quadratic</I>.
It may be more efficient than <I>quadratic</I> on some systems.
</P>
<P><B>Restrictions:</B> none <P><B>Restrictions:</B> none
</P> </P>
<P><B>Related commands:</B> <P><B>Related commands:</B>
@ -74,6 +78,6 @@ It may be more efficient than <I>quadratic</I> on some systems.
</P> </P>
<P><B>Default:</B> <P><B>Default:</B>
</P> </P>
<P>The option defaults are dmax = 0.1 and line = backtrack. <P>The option defaults are dmax = 0.1 and line = quadratic.
</P> </P>
</HTML> </HTML>