yade-dev team mailing list archive

[Christian Jakob <jakob@xxxxxxxxxxxxxxxxxxx>]

Hi,

I want to do some updates/improvements of PFV method.

First thing I wanna do is a simulation with an increasing pore  
pressure in all cells to investigate liquefaction of a loose sand at  
different stress conditions. I tested imposePressure(), but for now it  
is not possible to use it within a loop over PFV cells. Actually it  
needs position of the cells (or at least a point within the cell), but  
I want it to do something like:

for i in range(flow.ncells()):
  flow.imposePressure(i,p)

I had a short look at the source code, it is quite confusing. I  
searched code for imposePressure. It finds 6 lines in  
pkg/dem/FlowEngine.hpp, but where is corresponding code in .cpp file?  
I noticed there is a macro called TPL. How can edit the code of this  
macro?

Regards,

Christian