146
146
double boxyhi = *((double *) lammps_extract_global(info->lmp,"boxyhi"));
147
147
double boxzlo = *((double *) lammps_extract_global(info->lmp,"boxzlo"));
148
148
double boxzhi = *((double *) lammps_extract_global(info->lmp,"boxzhi"));
149
double boxxy = *((double *) lammps_extract_global(info->lmp,"xy"));
150
double boxxz = *((double *) lammps_extract_global(info->lmp,"xz"));
151
double boxyz = *((double *) lammps_extract_global(info->lmp,"yz"));
150
153
double xprd = (boxxhi-boxxlo)/ANGSTROM_per_BOHR;
151
154
double yprd = (boxyhi-boxylo)/ANGSTROM_per_BOHR;
152
155
double zprd = (boxzhi-boxzlo)/ANGSTROM_per_BOHR;
156
double xy = boxxy/ANGSTROM_per_BOHR;
157
double xz = boxxz/ANGSTROM_per_BOHR;
158
double yz = boxyz/ANGSTROM_per_BOHR;
154
160
if (info->me == 0) {
155
161
sprintf(boxlines[0],"%g %g %g\n",xprd,0.0,0.0);
156
sprintf(boxlines[1],"%g %g %g\n",0.0,yprd,0.0);
157
sprintf(boxlines[2],"%g %g %g\n",0.0,0.0,zprd);
162
sprintf(boxlines[1],"%g %g %g\n",xy,yprd,0.0);
163
sprintf(boxlines[2],"%g %g %g\n",xz,yz,zprd);
160
166
// xlines = x for atoms on each proc converted to text lines