jianmu
/
lammps
forked from lijiext/lammps
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

try to fix code processing norms for minimizers that looked incorrect

This commit is contained in:
Axel Kohlmeyer 2019-10-14 16:40:24 +02:00
parent 29f2569ee2
commit 661238f627
No known key found for this signature in database
GPG Key ID: D9B44E93BF0C375A
4 changed files with 42 additions and 20 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
src/min.cpp
View File

@ -852,20 +852,22 @@ double Min::fnorm_max()
for (int m = 0; m < nextra_atom; m++) {
fatom = fextra_atom[m];
n = extra_nlen[m];
for (i = 0; i < n; i+=3)
fdotf = fvec[i]*fvec[i]+fvec[i+1]*fvec[i+1]+fvec[i+2]*fvec[i+2];
for (i = 0; i < n; i+=3) {
fdotf = fatom[i]*fatom[i]+fatom[i+1]*fatom[i+1]+fatom[i+2]*fatom[i+2];
local_norm_max = MAX(fdotf,local_norm_max);
}
}
}
double norm_max = 0.0;
MPI_Allreduce(&local_norm_max,&norm_max,1,MPI_DOUBLE,MPI_MAX,world);
if (nextra_global)
for (i = 0; i < n; i+=3)
fdotf = fvec[i]*fvec[i]+fvec[i+1]*fvec[i+1]+fvec[i+2]*fvec[i+2];
if (nextra_global) {
for (i = 0; i < nextra_global; i+=3) {
fdotf = fextra[i]*fextra[i];
norm_max = MAX(fdotf,norm_max);
}
}
return norm_max;
}
@ -906,14 +908,14 @@ double Min::total_torque()
double Min::inf_torque()
{
double fmsq,fmaxsqone,fmaxsqall;
double fmaxsqone,fmaxsqall;
int nlocal = atom->nlocal;
double hbar = force->hplanck/MY_2PI;
double tx,ty,tz;
double **sp = atom->sp;
double **fm = atom->fm;
fmsq = fmaxsqone = fmaxsqall = 0.0;
fmaxsqone = fmaxsqall = 0.0;
for (int i = 0; i < nlocal; i++) {
tx = fm[i][1]*sp[i][2] - fm[i][2]*sp[i][1];
ty = fm[i][2]*sp[i][0] - fm[i][0]*sp[i][2];

14
src/min_fire.cpp
View File

@ -251,10 +251,16 @@ int MinFire::iterate(int maxiter)
// sync across replicas if running multi-replica minimization
if (update->ftol > 0.0) {
if (normstyle == MAX) fdotf = fnorm_max(); // max force norm
else if (normstyle == INF) fdotf = fnorm_inf(); // inf force norm
else if (normstyle == TWO) fdotf = fnorm_sqr(); // Euclidean force 2-norm
else error->all(FLERR,"Illegal min_modify command");
if (normstyle == MAX) {
fdotfloc = fnorm_max(); // max force norm
MPI_Allreduce(&fdotfloc,&fdotf,1,MPI_DOUBLE,MPI_MAX,universe->uworld);
} else if (normstyle == INF) {
fdotfloc = fnorm_inf(); // inf force norm
MPI_Allreduce(&fdotfloc,&fdotf,1,MPI_DOUBLE,MPI_MAX,universe->uworld);
} else if (normstyle == TWO) {
fdotf = fnorm_sqr(); // Euclidean force 2-norm
MPI_Allreduce(&fdotfloc,&fdotf,1,MPI_DOUBLE,MPI_SUM,universe->uworld);
} else error->all(FLERR,"Illegal min_modify command");
if (update->multireplica == 0) {
if (fdotf < update->ftol*update->ftol) return FTOL;
} else {

14
src/min_quickmin.cpp
View File

@ -216,10 +216,16 @@ int MinQuickMin::iterate(int maxiter)
// sync across replicas if running multi-replica minimization
if (update->ftol > 0.0) {
if (normstyle == MAX) fdotfloc = fnorm_max(); // max force norm
else if (normstyle == INF) fdotfloc = fnorm_inf(); // inf force norm
else if (normstyle == TWO) fdotfloc = fnorm_sqr(); // Euclidean force 2-norm
else error->all(FLERR,"Illegal min_modify command");
if (normstyle == MAX) {
fdotfloc = fnorm_max(); // max force norm
MPI_Allreduce(&fdotfloc,&fdotf,1,MPI_DOUBLE,MPI_MAX,universe->uworld);
} else if (normstyle == INF) {
fdotfloc = fnorm_inf(); // inf force norm
MPI_Allreduce(&fdotfloc,&fdotf,1,MPI_DOUBLE,MPI_MAX,universe->uworld);
} else if (normstyle == TWO) {
fdotfloc = fnorm_sqr(); // Euclidean force 2-norm
MPI_Allreduce(&fdotfloc,&fdotf,1,MPI_DOUBLE,MPI_SUM,universe->uworld);
} else error->all(FLERR,"Illegal min_modify command");
if (update->multireplica == 0) {
if (fdotf < update->ftol*update->ftol) return FTOL;
} else {

16
src/min_sd.cpp
View File

@ -17,6 +17,7 @@
#include "update.h"
#include "output.h"
#include "timer.h"
#include "universe.h"
using namespace LAMMPS_NS;
@ -77,11 +78,18 @@ int MinSD::iterate(int maxiter)
return ETOL;
// force tolerance criterion
// sync across replicas if running multi-replica minimization
if (normstyle == MAX) fdotf = fnorm_max(); // max force norm
else if (normstyle == INF) fdotf = fnorm_inf(); // infinite force norm
else if (normstyle == TWO) fdotf = fnorm_sqr(); // Euclidean force 2-norm
else error->all(FLERR,"Illegal min_modify command");
if (normstyle == MAX) {
fdotfloc = fnorm_max(); // max force norm
MPI_Allreduce(&fdotfloc,&fdotf,1,MPI_DOUBLE,MPI_MAX,universe->uworld);
} else if (normstyle == INF) {
fdotfloc = fnorm_inf(); // infinite force norm
MPI_Allreduce(&fdotfloc,&fdotf,1,MPI_DOUBLE,MPI_MAX,universe->uworld);
} else if (normstyle == TWO) {
fdotfloc = fnorm_sqr(); // Euclidean force 2-norm
MPI_Allreduce(&fdotfloc,&fdotf,1,MPI_DOUBLE,MPI_MAX,universe->uworld);
} else error->all(FLERR,"Illegal min_modify command");
if (fdotf < update->ftol*update->ftol) return FTOL;
// set new search direction h to f = -Grad(x)