Sorry for the late reply. My yahoo email was messed up by their change of format. I will describe the details of the problem with the replica exchange as Axel asked directly to the list:
The number of partitions is four (and there’s four temperatures) but when I
run it (mpirun -n 4 ~/lmp_openmpi -p 4 -cuda off -in in.run), it just says
it’s running on one partition and it does nothing and exits. With the gpu
it does the same thing but it throws a CUDA error like “Cuda driver error 4
in call at file ‘geryon/nvd_device.h’ in line 116.” sometimes before and
sometimes after it prints the information about 1 partition each time I try
it.
Everything’s below and I would greatly appreciate any help you could
provide. I will take your advice and switch to lj/sdk.
Thanks,
Luke Czapla
The restart2data output looked like this (I noticed that even when the
restart2data is from the same source tree as my binary, the version numbers
are different, but this restart file is from running with an earlier build):
luke@…4543…:~/random/q0.5_replica$ ./restart2data run4.rest DATA.FILE
in5.run
Reading restart file …
WARNING Restart file version does not match restart2data version
restart2data version = 23 July 2013
Restart file version = 31 May 2013
Ntimestep = 25000000
Nprocs = 4
Natoms = 40000
Nbonds = 3732
Nangles = 3599
Unit style = real
Atom style = full
Pair style = cg/cmm/coul/long
Bond style = harmonic
Angle style = cg/cmm
Xlo xhi = -76.2836 76.2836
Ylo yhi = -76.2836 76.2836
Zlo zhi = -76.2836 76.2836
Periodicity = 1 1 1
Boundary = 0 0, 0 0, 0 0
Writing data file …
Writing input file …
ERROR: Cannot write pair_style cg/cmm/coul/long to input file
The input file looks like this:
package gpu force/neigh 0 1 -1
units real
dimension 3
atom_style full
read_restart run2.rest
#read_data DATA.FILE
#include PARM.FILE
kspace_style pppm 1.0e-5
kspace_modify mesh 32 32 32 order 3
reset_timestep 0
neighbor 2.0 bin
#original neigh_modify delay 5
#neigh_modify exclude molecule sphere
neigh_modify delay 0 every 2 check yes
timestep 20.0
#run_style verlet
run_style respa 2 2 bond 1 angle 1 pair 2 kspace 2
#velocity all create 303.0 87285 dist gaussian
#fix 1 all npt temp 303.0 303.0 1000.0 iso 1.0 1.0 1000.0
#fix 2 all momentum 1 linear 1 1 1 angular
dump 1 all dcd 1000 run3.dcd
thermo 1000
thermo_style multi
variable t world 310.0 330.0 350.0 370.0
fix myfix all nvt $t $t 100.0
reset_timestep 0
temper 20000000 500 $t tempfix 37 69
write_restart run3.rest