yade-dev team mailing list archive

[Bruno Chareyre <bruno.chareyre@xxxxxxxxxxx>]

Václav Šmilauer a écrit :
> > > 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.
> > >       
>
> As said, refR{1,2}<0 is effective for facets, where it makes perfect
> sense. It will have no effect if you use ScGeom, so please keep those
> things there.
>
>   
Ok.

Just to be sure, did you double check the equation? It says :
if (R1<0) R1 = 2*R2;
if (R2<0) R2 = R1;

The factor two is for the double-radius sphere?
Is it correct to multiply only in R1=... and not in R2=...?

Bruno