else if(d3dg){Da=d3dg->refR1>0?d3dg->refR1:2*d3dg->refR2; Db=d3dg->refR2>0?d3dg->refR2:d3dg->refR1;

Apparently, the equations are not symmetric and the radius of the virtual sphere will be 2*r "sometimes"... I suggest using 1*r in all cases, but this is the sort of thing that will modify results (at least slightly) for everybody.
Oooops, no, it was just a typo, D=2*r anyway... Its fixed now.

*B*

_______________________________________________
Mailing list: https://launchpad.net/~yade-dev
Post to     : [email protected]
Unsubscribe : https://launchpad.net/~yade-dev
More help   : https://help.launchpad.net/ListHelp

Reply via email to