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

git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@10047 f3b2605a-c512-4ea7-a41b-209d697bcdaa

Browse Source
This commit is contained in:
sjplimp
2013-06-06 20:18:08 +00:00
parent 51814ba9bd
commit aa1ebc8444
3 changed files with 12 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
src/fix_langevin.cpp
View File

@ -160,7 +160,7 @@ FixLangevin::FixLangevin(LAMMPS *lmp, int narg, char **arg) :
if (gjfflag) {
nvalues = 3;
grow_arrays(atom->nlocal);
grow_arrays(atom->nmax);
atom->add_callback(0);
// initialize franprev to zero
@ -308,10 +308,10 @@ void FixLangevin::post_force(int vflag)
if (tbiasflag == BIAS)
if (rmass)
if (zeroflag) post_force_templated<1,1,1,1,1,1>();
else post_force_templated<1,1,1,1,1,0>();
else post_force_templated<1,1,1,1,1,0>();
else
if (zeroflag) post_force_templated<1,1,1,1,0,1>();
else post_force_templated<1,1,1,1,0,0>();
else post_force_templated<1,1,1,1,0,0>();
else
if (rmass)
if (zeroflag) post_force_templated<1,1,1,0,1,1>();
@ -884,9 +884,7 @@ double FixLangevin::memory_usage()
void FixLangevin::grow_arrays(int nmax)
{
if (!gjfflag) return;
memory->grow(franprev,nmax,3,"fix_langevin:franprev");
array = franprev;
}
/* ----------------------------------------------------------------------
@ -896,7 +894,7 @@ void FixLangevin::grow_arrays(int nmax)
void FixLangevin::copy_arrays(int i, int j, int delflag)
{
for (int m = 0; m < nvalues; m++)
array[j][m] = array[i][m];
franprev[j][m] = franprev[i][m];
}
/* ----------------------------------------------------------------------
@ -905,7 +903,7 @@ void FixLangevin::copy_arrays(int i, int j, int delflag)
int FixLangevin::pack_exchange(int i, double *buf)
{
for (int m = 0; m < nvalues; m++) buf[m] = array[i][m];
for (int m = 0; m < nvalues; m++) buf[m] = franprev[i][m];
return nvalues;
}
@ -915,7 +913,6 @@ int FixLangevin::pack_exchange(int i, double *buf)
int FixLangevin::unpack_exchange(int nlocal, double *buf)
{
for (int m = 0; m < nvalues; m++) array[nlocal][m] = buf[m];
for (int m = 0; m < nvalues; m++) franprev[nlocal][m] = buf[m];
return nvalues;
}