2010/2/11 Luc Sibille <[email protected]> > 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 ); >> > > How do you compute "nbElt", do you take into account the elements not > involved in the contact force network (i.e. elements involved in less than 2 > contacts)? Maybe element not involved in the contact force network are not > taken into account in the python script. >
In my case (C++ code), nbElt is the total number of elements (spheres only) contained in the simulation. I did not make anything to distinguish spheres with less than 2 contacts... I checked the python script in _utils.cpp:147 and I cannot see why there is a difference (I didn't see anything to avoid spheres with less than 2 contacts in the script), but, I must confess that the code in there is not very meaningful to me... I have to add that when I look into the xml file, my C++ code tends to be accurate (just divide the total number of interactions by the number of particles and multiply by 2). I have to work on the understanding of the python script, but if you have any comments, do not hesitate. Luc
_______________________________________________ Mailing list: https://launchpad.net/~yade-users Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
