yade-dev team mailing list archive

[Bruno Chareyre <bruno.chareyre@xxxxxxxxxxxxxxx>]

Hi Jérôme,
Thanks for update. Not sure I understand the pdf. Why is the new code
faster than the old code with -j1?
Taking the difference between -j1 and -j4 with the new code suggests
that the speedup of the capillary law is about x1.15, correct?
Before commiting, could you show the diff?
Thanks
Bruno

On 11/03/15 18:28, Jerome Duriez wrote:
> Hi guys,
>
> I finally reconsidered parallelization of
> Law2_ScGeom_CapillaryPhys_Capillarity, after our passed discussion (*).
>
> I noticed only after re-starting this work that implementations of
> e.g. remove() and resize() in in-house class OpenMPVector would be
> necessary to do it completely right (using OpenMPVector objects in
> BodiesMenisciiList class). Then, I only coded yet parallel loops for
> the "addForce" block, at the end of the Law2....
>
> Obviously, the results (in terms of capillary stresses as given by
> getCapillaryStress() function) are unaffected by the changes : see the
> four attached pictures, showing sXXcap or sYYcap for two triaxial
> loadings with 20000 particles, repeted before and after changing code.
> One loading is with 10 kPa capillary pressure : meniscii fusion occur
> in the sample, affecting the simulation (with fusionDetection =
> 1,binaryFusion=0 attributes of the Law2). For the other one, the
> suction value of 300 kPa prevents any meniscii fusion.
>
> As for the speed up, *maybe* such changes in the code allow a
> significant speed up (of capillary Law2) to be obtained with parallel
> computations... See the data sheet "NBl2cp" (in ods or pdf format).
> Anyway, I hope it is clear this little change in code does not harm.
> Since it might still be little useful, and exists already on my
> computer, I plann to commit it, after time for discussion.
>
> Jerome
>
> PS : these results might still be of interest for you if you are
> curious of results scattering through parallel computations, or
> capillary law2 and NewtonIntegrator (NI) timings... I obtain here e.g.
> a speed up of NI about 2.3 from j1 to j4
>
>
> (*)
> https://www.mail-archive.com/yade-dev@xxxxxxxxxxxxxxxxxxx/msg10842.html)
>
>
> _______________________________________________
> 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


-- 
_______________
Bruno Chareyre
Associate Professor
ENSE³ - Grenoble INP
Lab. 3SR
BP 53
38041 Grenoble cedex 9
Tél : +33 4 56 52 86 21
Fax : +33 4 76 82 70 43
________________