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

restart bug fix in new fix ttm/grid

Browse Source
This commit is contained in:
Steve Plimpton
2021-08-26 10:08:34 -06:00
parent a0dfae9876
commit 0fb03072b1
1 changed files with 5 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
src/EXTRA-FIX/fix_ttm_grid.cpp
View File

@ -593,10 +593,12 @@ void FixTTMGrid::restart(char *buf)
for (ix = nxlo_in; ix <= nxhi_in; ix++) {
iglobal = nygrid*nxgrid*iz + nxgrid*iy + ix;
T_electron[iz][iy][ix] = rlist[n+iglobal];
if (ix == 0 && iy == 0 & iz == 0)
printf("READRESTART 000 n+iglobal %d %d value %g\n",n,iglobal,
T_electron[ix][iy][iz]);
}
// communicate new T_electron values to ghost grid points
gc->forward_comm(GridComm::FIX,this,1,sizeof(double),0,
gc_buf1,gc_buf2,MPI_DOUBLE);
}
/* ----------------------------------------------------------------------