yade-dev team mailing list archive

[Bruno Chareyre <bruno.chareyre@xxxxxxxxxxx>]

> > Why does it return a wrong value? spheresVolume is internally correct
> > (it is used to define porosity), wrapping changes the value?!
> >      
> I don't believe that this is even possible.
>
>    
It is apparently...
> If I read the code right, then spheresVolume is calculated once, at
> start. And then it is being modified in that line 210:
>
> 	spheresVolume = spheresVolume * std::pow ( multiplier,3 );
>
>    
Correct.
> Why are you doing this, btw?
>
>    
If all sizes are multiplied by "multiplier", then the volume is 
multiplied by multiplier^3.

> TriaxialStressController.cpp line 165. How many times that line is called?
>
>    
Each time the size of spheres is modified (i.e. probably 0 times in your 
case, since you didn't even realize it could be changed ;-)).


> now, question is if I read the code right... And why porosity is
> correct if spheresVolume is not??
>
>    
That is the interesting question. spheresVolume is correct, really. 
Would python convert it to int for some reason?

Bruno