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@2040 f3b2605a-c512-4ea7-a41b-209d697bcdaa

Browse Source
This commit is contained in:
sjplimp
2008-08-08 14:28:23 +00:00
parent dfa5c73928
commit c48558f297
1 changed files with 34 additions and 33 deletions
Download Patch File Download Diff File Expand all files Collapse all files

67
tools/restart2data.cpp
View File

@ -2488,10 +2488,10 @@ void Data::write(FILE *fp, FILE *fp2)
fprintf(fp,"\n");
fprintf(fp,"%g %g xlo xhi\n",xlo,xhi);
fprintf(fp,"%g %g ylo yhi\n",ylo,yhi);
fprintf(fp,"%g %g zlo zhi\n",zlo,zhi);
if (triclinic) fprintf(fp,"%g %g %g xy xz yz\n",xy,xz,yz);
fprintf(fp,"%-1.16e %-1.16e xlo xhi\n",xlo,xhi);
fprintf(fp,"%-1.16e %-1.16e ylo yhi\n",ylo,yhi);
fprintf(fp,"%-1.16e %-1.16e zlo zhi\n",zlo,zhi);
if (triclinic) fprintf(fp,"%-1.16e %-1.16e %-1.16e xy xz yz\n",xy,xz,yz);
// write ff styles to input file
@ -2981,7 +2981,8 @@ void Data::write(FILE *fp, FILE *fp2)
fprintf(fp,"\n");
} else {
fprintf(fp,"%d %d %g %g %g",tag[i],type[i],x[i],y[i],z[i]);
fprintf(fp,"%d %d %-1.16e %-1.16e %-1.16e",
tag[i],type[i],x[i],y[i],z[i]);
for (int k = 1; k <= style_hybrid; k++) {
if (k == style_angle) write_atom_angle_extra(fp,i);
if (k == style_atomic) write_atom_atomic_extra(fp,i);
@ -3017,7 +3018,7 @@ void Data::write(FILE *fp, FILE *fp2)
fprintf(fp,"\n");
} else {
fprintf(fp,"%d %g %g %g",tag[i],vx[i],vy[i],vz[i]);
fprintf(fp,"%d %-1.16e %-1.16e %-1.16e",tag[i],vx[i],vy[i],vz[i]);
for (int k = 1; k <= style_hybrid; k++) {
if (k == style_angle) write_vel_angle_extra(fp,i);
if (k == style_atomic) write_vel_atomic_extra(fp,i);
@ -3072,62 +3073,62 @@ void Data::write(FILE *fp, FILE *fp2)
void Data::write_atom_angle(FILE *fp, int i, int ix, int iy, int iz)
{
fprintf(fp,"%d %d %d %15.12g %15.12g %15.12g %d %d %d",
fprintf(fp,"%d %d %d %-1.16e %-1.16e %-1.16e %d %d %d",
tag[i],molecule[i],type[i],x[i],y[i],z[i],ix,iy,iz);
}
void Data::write_atom_atomic(FILE *fp, int i, int ix, int iy, int iz)
{
fprintf(fp,"%d %d %15.12g %15.12g %15.12g %d %d %d",
fprintf(fp,"%d %d %-1.16e %-1.16e %-1.16e %d %d %d",
tag[i],type[i],x[i],y[i],z[i],ix,iy,iz);
}
void Data::write_atom_bond(FILE *fp, int i, int ix, int iy, int iz)
{
fprintf(fp,"%d %d %d %15.12g %15.12g %15.12g %d %d %d",
fprintf(fp,"%d %d %d %-1.16e %-1.16e %-1.16e %d %d %d",
tag[i],molecule[i],type[i],x[i],y[i],z[i],ix,iy,iz);
}
void Data::write_atom_charge(FILE *fp, int i, int ix, int iy, int iz)
{
fprintf(fp,"%d %d %15.12g %15.12g %15.12g %15.12g %d %d %d",
fprintf(fp,"%d %d %-1.16e %-1.16e %-1.16e %-1.16e %d %d %d",
tag[i],type[i],q[i],x[i],y[i],z[i],ix,iy,iz);
}
void Data::write_atom_dipole(FILE *fp, int i, int ix, int iy, int iz)
{
fprintf(fp,"%d %d %15.12g %15.12g %15.12g %15.12g %15.12g %15.12g %15.12g %d %d %d",
fprintf(fp,"%d %d %-1.16e %-1.16e %-1.16e %-1.16e %-1.16e %-1.16e %-1.16e %d %d %d",
tag[i],type[i],q[i],x[i],y[i],z[i],mux[i],muy[i],muz[i],ix,iy,iz);
}
void Data::write_atom_dpd(FILE *fp, int i, int ix, int iy, int iz)
{
fprintf(fp,"%d %d %15.12g %15.12g %15.12g %d %d %d",
fprintf(fp,"%d %d %-1.16e %-1.16e %-1.16e %d %d %d",
tag[i],type[i],x[i],y[i],z[i],ix,iy,iz);
}
void Data::write_atom_ellipsoid(FILE *fp, int i, int ix, int iy, int iz)
{
fprintf(fp,"%d %d %15.12g %15.12g %15.12g %15.12g %15.12g %15.12g %15.12g %d %d %d",
fprintf(fp,"%d %d %-1.16e %-1.16e %-1.16e %-1.16e %-1.16e %-1.16e %-1.16e %d %d %d",
tag[i],type[i],x[i],y[i],z[i],
quatw[i],quati[i],quatj[i],quatk[i],ix,iy,iz);
}
void Data::write_atom_full(FILE *fp, int i, int ix, int iy, int iz)
{
fprintf(fp,"%d %d %d %15.12g %15.12g %15.12g %15.12g %d %d %d",
fprintf(fp,"%d %d %d %-1.16e %-1.16e %-1.16e %-1.16e %d %d %d",
tag[i],molecule[i],type[i],q[i],x[i],y[i],z[i],ix,iy,iz);
}
void Data::write_atom_granular(FILE *fp, int i, int ix, int iy, int iz)
{
fprintf(fp,"%d %d %15.12g %15.12g %15.12g %15.12g %15.12g %d %d %d",
fprintf(fp,"%d %d %-1.16e %-1.16e %-1.16e %-1.16e %-1.16e %d %d %d",
tag[i],type[i],2.0*radius[i],density[i],x[i],y[i],z[i],ix,iy,iz);
}
void Data::write_atom_molecular(FILE *fp, int i, int ix, int iy, int iz)
{
fprintf(fp,"%d %d %d %15.12g %15.12g %15.12g %d %d %d",
fprintf(fp,"%d %d %d %-1.16e %-1.16e %-1.16e %d %d %d",
tag[i],molecule[i],type[i],x[i],y[i],z[i],ix,iy,iz);
}
@ -3150,29 +3151,29 @@ void Data::write_atom_bond_extra(FILE *fp, int i)
void Data::write_atom_charge_extra(FILE *fp, int i)
{
fprintf(fp," %15.12g",q[i]);
fprintf(fp," %-1.16e",q[i]);
}
void Data::write_atom_dipole_extra(FILE *fp, int i)
{
fprintf(fp," %15.12g %15.12g %15.12g %15.12g",q[i],mux[i],muy[i],muz[i]);
fprintf(fp," %-1.16e %-1.16e %-1.16e %-1.16e",q[i],mux[i],muy[i],muz[i]);
}
void Data::write_atom_dpd_extra(FILE *fp, int i) {}
void Data::write_atom_ellipsoid_extra(FILE *fp, int i)
{
fprintf(fp," %15.12g %15.12g %15.12g %15.12g",quatw[i],quati[i],quatj[i],quatk[i]);
fprintf(fp," %-1.16e %-1.16e %-1.16e %-1.16e",quatw[i],quati[i],quatj[i],quatk[i]);
}
void Data::write_atom_full_extra(FILE *fp, int i)
{
fprintf(fp," %d %15.12g",molecule[i],q[i]);
fprintf(fp," %d %-1.16e",molecule[i],q[i]);
}
void Data::write_atom_granular_extra(FILE *fp, int i)
{
fprintf(fp," %15.12g %15.12g",2.0*radius[i],density[i]);
fprintf(fp," %-1.16e %-1.16e",2.0*radius[i],density[i]);
}
void Data::write_atom_molecular_extra(FILE *fp, int i)
@ -3187,54 +3188,54 @@ void Data::write_atom_molecular_extra(FILE *fp, int i)
void Data::write_vel_angle(FILE *fp, int i)
{
fprintf(fp,"%d %15.12g %15.12g %15.12g",tag[i],vx[i],vy[i],vz[i]);
fprintf(fp,"%d %-1.16e %-1.16e %-1.16e",tag[i],vx[i],vy[i],vz[i]);
}
void Data::write_vel_atomic(FILE *fp, int i)
{
fprintf(fp,"%d %15.12g %15.12g %15.12g",tag[i],vx[i],vy[i],vz[i]);
fprintf(fp,"%d %-1.16e %-1.16e %-1.16e",tag[i],vx[i],vy[i],vz[i]);
}
void Data::write_vel_bond(FILE *fp, int i)
{
fprintf(fp,"%d %15.12g %15.12g %15.12g",tag[i],vx[i],vy[i],vz[i]);
fprintf(fp,"%d %-1.16e %-1.16e %-1.16e",tag[i],vx[i],vy[i],vz[i]);
}
void Data::write_vel_charge(FILE *fp, int i)
{
fprintf(fp,"%d %15.12g %15.12g %15.12g",tag[i],vx[i],vy[i],vz[i]);
fprintf(fp,"%d %-1.16e %-1.16e %-1.16e",tag[i],vx[i],vy[i],vz[i]);
}
void Data::write_vel_dipole(FILE *fp, int i)
{
fprintf(fp,"%d %15.12g %15.12g %15.12g",tag[i],vx[i],vy[i],vz[i]);
fprintf(fp,"%d %-1.16e %-1.16e %-1.16e",tag[i],vx[i],vy[i],vz[i]);
}
void Data::write_vel_dpd(FILE *fp, int i)
{
fprintf(fp,"%d %15.12g %15.12g %15.12g",tag[i],vx[i],vy[i],vz[i]);
fprintf(fp,"%d %-1.16e %-1.16e %-1.16e",tag[i],vx[i],vy[i],vz[i]);
}
void Data::write_vel_ellipsoid(FILE *fp, int i)
{
fprintf(fp,"%d %15.12g %15.12g %15.12g %15.12g %15.12g %15.12g",
fprintf(fp,"%d %-1.16e %-1.16e %-1.16e %-1.16e %-1.16e %-1.16e",
tag[i],vx[i],vy[i],vz[i],angmomx[i],angmomy[i],angmomz[i]);
}
void Data::write_vel_full(FILE *fp, int i)
{
fprintf(fp,"%d %15.12g %15.12g %15.12g",tag[i],vx[i],vy[i],vz[i]);
fprintf(fp,"%d %-1.16e %-1.16e %-1.16e",tag[i],vx[i],vy[i],vz[i]);
}
void Data::write_vel_granular(FILE *fp, int i)
{
fprintf(fp,"%d %15.12g %15.12g %15.12g %15.12g %15.12g %15.12g",
fprintf(fp,"%d %-1.16e %-1.16e %-1.16e %-1.16e %-1.16e %-1.16e",
tag[i],vx[i],vy[i],vz[i],omegax[i],omegay[i],omegaz[i]);
}
void Data::write_vel_molecular(FILE *fp, int i)
{
fprintf(fp,"%d %15.12g %15.12g %15.12g",tag[i],vx[i],vy[i],vz[i]);
fprintf(fp,"%d %-1.16e %-1.16e %-1.16e",tag[i],vx[i],vy[i],vz[i]);
}
// ---------------------------------------------------------------------
@ -3251,14 +3252,14 @@ void Data::write_vel_dpd_extra(FILE *fp, int i) {}
void Data::write_vel_ellipsoid_extra(FILE *fp, int i)
{
fprintf(fp," %15.12g %15.12g %15.12g",angmomx[i],angmomy[i],angmomz[i]);
fprintf(fp," %-1.16e %-1.16e %-1.16e",angmomx[i],angmomy[i],angmomz[i]);
}
void Data::write_vel_full_extra(FILE *fp, int i) {}
void Data::write_vel_granular_extra(FILE *fp, int i)
{
fprintf(fp," %15.12g %15.12g %15.12g",omegax[i],omegay[i],omegaz[i]);
fprintf(fp," %-1.16e %-1.16e %-1.16e",omegax[i],omegay[i],omegaz[i]);
}
void Data::write_vel_molecular_extra(FILE *fp, int i) {}