yade-dev team mailing list archive

[Anton Gladky <gladky.anton@xxxxxxxxx>]

spheresVolume=1 is in constructor. Maybe it is not updating later?

http://bazaar.launchpad.net/%7Eyade-dev/yade/trunk/annotate/head%3A/pkg/dem/Engine/PartialEngine/TriaxialCompressionEngine.hpp#L109


2010/7/7 Bruno Chareyre <bruno.chareyre@xxxxxxxxxxx>

> The data was duplicated in the inheriting CompressionEngine. Commit in a
> few minutes.
>
> B.
>
>
> On 07/07/10 19:54, Bruno Chareyre wrote:
>
>>
>>  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
>>
>>
>>
>>
>> _______________________________________________
>> Mailing list: https://launchpad.net/~yade-dev<https://launchpad.net/%7Eyade-dev>
>> Post to     : yade-dev@xxxxxxxxxxxxxxxxxxx
>> Unsubscribe : https://launchpad.net/~yade-dev<https://launchpad.net/%7Eyade-dev>
>> More help   : https://help.launchpad.net/ListHelp
>>
>>
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-dev<https://launchpad.net/%7Eyade-dev>
> Post to     : yade-dev@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-dev<https://launchpad.net/%7Eyade-dev>
> More help   : https://help.launchpad.net/ListHelp
>