Re: [Yade-users] Are our systems conservatives?

2010-05-07 Thread chiara modenese
On 7 May 2010 16:40, Bruno Chareyre wrote: > > >> I tried the sample example with the two balls and the time step estimated >> by PWave. I attach the figure. >> > Mmmm... It is not "that" better. I expected better matching with corrected > kinetic (as in the simple mass-spring test). "Shifting"

Re: [Yade-users] Are our systems conservatives?

2010-05-07 Thread chiara modenese
On 7 May 2010 16:49, Bruno Chareyre wrote: > > >> So, you need in fact to __substract__ accel*dt/2, and you will get at time >> "t" the energy of time "t-dt". Sorry if I wrote something different in >> previous mails. >> > Not sure it was a good idea, since nobody will understand why kineticEnerg

Re: [Yade-users] Are our systems conservatives?

2010-05-07 Thread Bruno Chareyre
So, you need in fact to __substract__ accel*dt/2, and you will get at time "t" the energy of time "t-dt". Sorry if I wrote something different in previous mails. Not sure it was a good idea, since nobody will understand why kineticEnergy computed a time t refers to time t-dt. The best is to

Re: [Yade-users] Are our systems conservatives?

2010-05-07 Thread chiara modenese
On 7 May 2010 16:40, Bruno Chareyre wrote: > > >> I tried the sample example with the two balls and the time step estimated >> by PWave. I attach the figure. >> > Mmmm... It is not "that" better. I expected better matching with corrected > kinetic (as in the simple mass-spring test). > > > Now I

Re: [Yade-users] Are our systems conservatives?

2010-05-07 Thread Bruno Chareyre
I tried the sample example with the two balls and the time step estimated by PWave. I attach the figure. Mmmm... It is not "that" better. I expected better matching with corrected kinetic (as in the simple mass-spring test). Now I suppose that if the kinetic energy is computed correctly (is

Re: [Yade-users] Are our systems conservatives?

2010-05-07 Thread chiara modenese
On 7 May 2010 15:51, chiara modenese wrote: > > > On 7 May 2010 13:31, Bruno Chareyre wrote: > >> >>> Ok, can I change something in Shop for the dt/2 shift or should I save >>> velocities and accelerations and then interpolate? >>> >>> >> You can change in shop. >> > I changed in shop the kineti

Re: [Yade-users] Are our systems conservatives?

2010-05-07 Thread chiara modenese
On 7 May 2010 13:31, Bruno Chareyre wrote: > >> Ok, can I change something in Shop for the dt/2 shift or should I save >> velocities and accelerations and then interpolate? >> >> > You can change in shop. > I changed in shop the kinetic as ret+=.5*(b->state->mass*(b->state->vel+b->state->accel*d

Re: [Yade-users] Experiments on a rotating cylinder

2010-05-07 Thread Sergei D.
> > I also added > https://www.yade-dem.org/sphinx/prog.html#local-documentation with > instructions how to generate documentation locally. > When I was generating the local documentation, I get an error: INFO yade.Omega /home/sega/work/yade/trunk/core/Omega.cpp:329 saveSimulation: Saving fi

Re: [Yade-users] Are our systems conservatives?

2010-05-07 Thread Bruno Chareyre
Ok, can I change something in Shop for the dt/2 shift or should I save velocities and accelerations and then interpolate? You can change in shop. Another note : your remark on different dt for different problems is correct. PWave and GSTS both consider a dense packing. If you have problems

Re: [Yade-users] Are our systems conservatives?

2010-05-07 Thread chiara modenese
On 7 May 2010 12:08, Bruno Chareyre wrote: > > Perhaps I guess that dealing with _quasi-static conditions_, energy >> considerations become superfluous (that is maybe why nobody here has >> approached the problem?). On the other hand, if the problem is _dynamic_ I >> think they matter. >> >> Is

Re: [Yade-users] Are our systems conservatives?

2010-05-07 Thread Bruno Chareyre
Perhaps I guess that dealing with _quasi-static conditions_, energy considerations become superfluous (that is maybe why nobody here has approached the problem?). On the other hand, if the problem is _dynamic_ I think they matter. Is not supefluous under quasistatic conditions, and I approac

Re: [Yade-users] Are our systems conservatives?

2010-05-07 Thread Bruno Chareyre
I am only applying an initial velocity to one particle, hence the total energy of the system is only given by (kinetic + potential elastic). Look at the two attached figures. The first one (dt_crit) is the result obtained using the critical time step computed by Pwave in Yade. The second one

Re: [Yade-users] Experiments on a rotating cylinder

2010-05-07 Thread Václav Šmilauer
> Vaclav, I added doc in bzr2207. > Please, review it and tell me if I done some wrong (I mean format, etc...). Hi Sega, great, thanks! It's almost fine, I fixed syntax, no problem. I already observed different behavior for syntax in pure python and c++ docstrings, so I put the question on sphinx

Re: [Yade-users] Experiments on a rotating cylinder

2010-05-07 Thread Sergei D.
> Sega, it should be documented (in py/_utils.cpp), at the end. Please > name arguments (py::arg("name")) and also put formuas that are used into > the docstring (either in $...$ or ".. math::"), please. > Hello! Vaclav, I added doc in bzr2207. Please, review it and tell me if I done some wrong