> I need to get incrementally the shear elastic potential energy for HM. > For the normal direction I can work out the integral, so no problem > but in the shear direction it has to be done incrementally. Bruno, I > see you used the OpenMPAccumulator to get plastic dissipation > incrementally (I followed more or less this thread > http://www.mail-archive.com/[email protected]/msg03490.html). Is > this really necessary? Would you suggest to follow this way?
If you store it in the interaction (i.e. InteractionPhysics/InteractionGeometry), there is no need to use the accumulator -- it serves for accumulating the values concurrently from different threads, but interactions are isolated. You can just have a normal Real/Vector3r as you need. HTH, Vaclav _______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp

