hi

in order to get a difference approximation to a time derivative I need to have my compute called at one step before the thermo output step (as well as the output step). I thought that adding this in the ::compute function:

int next_compute = output->next_thermo + output->thermo_every -1;

this->addstep(next_compute);

and this in the constructor

timeflag = 1; // stores list of timesteps to compute

would do the trick but then I looked at Verlet.cpp and it looks like computes are only called at thermo output times.

So, is it possible to call a compute at steps other than the thermo output steps?

thanks,

Reese