From c48558f2972738d9a91eb8087978ca256e7259eb Mon Sep 17 00:00:00 2001 From: sjplimp Date: Fri, 8 Aug 2008 14:28:23 +0000 Subject: [PATCH] git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@2040 f3b2605a-c512-4ea7-a41b-209d697bcdaa --- tools/restart2data.cpp | 67 +++++++++++++++++++++--------------------- 1 file changed, 34 insertions(+), 33 deletions(-) diff --git a/tools/restart2data.cpp b/tools/restart2data.cpp index 9770392673..cdb61f4264 100644 --- a/tools/restart2data.cpp +++ b/tools/restart2data.cpp @@ -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) {}