yade-users team mailing list archive

[luc scholtes <lscholtes63@xxxxxxxxx>]

Hi all,

A remark since I tried to compute the coordination number for a given
simulation.

When I use those c++ lines:

Real inter=0;

FOREACH(const shared_ptr<Interaction>& I, *ncb->interactions){

      if( !I->isReal() ) continue;
      inter += 1;
}

Real coordN = 0;
coordN = 2* ( inter/nbElt );

it gives me a lower average number of interactions K per body than the
python script defined in utils.plotDirections().

I suspected an exclusion of the sphere/wall interactions but, even in that
case, I should have a greater K.

Thanks for your advices


  Luc