yade-dev team mailing list archive

[Václav Šmilauer <eudoxos@xxxxxxxx>]

> Hi Vaclav, thanks for suggestion, it seems, I have fixed it. 
Hi, what I suggested was wrong, as I saw from results of
examples/concrete/uniax.py. I fixed that in the last commit. The results
seem to be the same. There are still some things to wrap in python, run
scripts/regression-tests.py, but I take care of that.

The uniax.py simulation runs about 10x slower with eigen, probably due
to lots of construction/descruction of temporary objects.

Please test some of your simulations as well, then we will ask people to
test theirs.

Once this is done, we might either migrate yade code to eigen, or (which
would be smoother and would let us compare performance), first
progressively adapt miniWm3 interface to plain eigen (without our
wrapper types) by adding extra methods, then changing sources of yade.
We would create a few extra functions for things that have no direct
correspondence in eigen, or ask upstream for including those.

Long way ahead, but with solid basis of nice work of yours, it will be
feasible, I think.

Cheers, Vaclav


> 
> When I ran some tests with old variant and a new one (with null-vector
> checks), I got a little bit different results with regular-hexa. Does
> it mean, that in some cases geometry was wrongly calculated
> previously?
> 
> 
> 2010/3/15 <noreply@xxxxxxxxxxxxx>
>         ------------------------------------------------------------
>         revno: 2079
>         committer: Anton Gladky <gladky.anton@xxxxxxxxx>
>         branch nick: trunk
>         timestamp: Mon 2010-03-15 12:24:35 +0100
>         message:
>          Previous commit correction. Thanks to Vaclav.
>         modified:
>          pkg/dem/DataClass/InteractionGeometry/Dem3DofGeom_SphereSphere.cpp
>         
>         
>         --
>         lp:yade
>         https://code.launchpad.net/~yade-dev/yade/trunk
>         
>         Your team Yade developers is subscribed to branch lp:yade.
>         To unsubscribe from this branch go to
>         https://code.launchpad.net/~yade-dev/yade/trunk/+edit-subscription.
>         
>         === modified file
>         'pkg/dem/DataClass/InteractionGeometry/Dem3DofGeom_SphereSphere.cpp'
>         ---
>         pkg/dem/DataClass/InteractionGeometry/Dem3DofGeom_SphereSphere.cpp  2010-03-12 21:20:32 +0000
>         +++
>         pkg/dem/DataClass/InteractionGeometry/Dem3DofGeom_SphereSphere.cpp  2010-03-15 11:24:35 +0000
>         @@ -45,9 +45,10 @@
>                        Vector3r axis=planeNormal.Cross(planePt);
>         axis.Normalize();
>                        Real angle=planePt.Length()/radius;
>                        normal2pt.FromAxisAngle(axis,angle);
>         +               Quaternionr ret; return
>         ret.Align(Vector3r::UNIT_X,normal2pt*planeNormal);
>         +       } else {
>         +               return Quaternionr::IDENTITY;
>                }
>         -       Quaternionr ret;
>         ret.Align(Vector3r::UNIT_X,normal2pt*planeNormal);
>         -       return ret;
>          }
>         
>         
>         
>         
>         _______________________________________________
>         Mailing list: https://launchpad.net/~yade-dev
>         Post to     : yade-dev@xxxxxxxxxxxxxxxxxxx
>         Unsubscribe : https://launchpad.net/~yade-dev
>         More help   : https://help.launchpad.net/ListHelp
>         
> 
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-dev
> Post to     : yade-dev@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-dev
> More help   : https://help.launchpad.net/ListHelp