Re: [Yade-dev] sign normal forces

2011-03-01 Thread Vincent Richefeu
Hi, I think that: (1) if you consider the negative/positive force decomposition, so the negative viscous forces are to be taken into account like any other negative force. So just test the sign and it's ok (in my opinion, but it depends on why you want to do this decomposition) (2) if you cons

Re: [Yade-dev] time step with viscous damping

2011-02-25 Thread Vincent Richefeu
Le 25 févr. 2011 à 11:38, Chiara Modenese a écrit : > Yes, I agree. I think I was confused by Vincent's words such as: "let say > that the critical time step determined for non-viscous contact will be > smaller that the effectively minimum time step". Vincent, maybe you could > clarify this po

Re: [Yade-dev] time step with viscous damping

2011-02-25 Thread Vincent Richefeu
Le 23 févr. 2011 à 17:43, Bruno Chareyre a écrit : > p.s. Vincent, you defined viscosity in order to keep critical time-step below > a given value IIRC. What was the reasoning behind? No. I (and other people like Sergei) defined the viscosity so that it is lower than the critical viscosity. If

Re: [Yade-dev] time step with viscous damping

2011-02-24 Thread Vincent Richefeu
Just a pdf file I found on the web: http://www.m-hikari.com/ces/ces2009/ces1-4-2009/teufelsbauerCES1-4-2009.pdf Maybe it helps !? Ciao Vincent Le 23 févr. 2011 à 19:01, Chiara Modenese a écrit : > Hi Bruno and Vincent, > > thanks for your answers. > > I have PFC manual on my desk (though I ne

Re: [Yade-dev] time step with viscous damping

2011-02-22 Thread Vincent Richefeu
Hi, Choosing an appropriate time step requires a certain expertise. Maybe a good reference could be the book of Allen and Tildesley (?) Keep in mind that the smallest time step depends also on the integration scheme. For instance, bigger time step can be used with velocity Verlet,Runge-Kutta or

Re: [Yade-dev] subdomains

2011-02-09 Thread Vincent Richefeu
Le 9 févr. 2011 à 17:45, Bruno Chareyre a écrit : > >> Oh simply cubic. Say that you divide the main cell into 8 cubic parts >> and then you get the average of stress in each of them (I simply want >> to compare them). > Here is the algorithm : > > foreach sphere: >cubeStress+=bodyStress(s

Re: [Yade-dev] subdomains

2011-02-09 Thread Vincent Richefeu
This is exactly what JJ Moreau did (so I think one can have some confidence with that) for spheres into a spherical probe (subvolume). Since the intersected volume can be analytically computed, he summed the moment tensor of each sphere weighted by the ratio "insidePartVolume" / "sphereVolume",

Re: [Yade-dev] weak/strong contacts

2011-02-07 Thread Vincent Richefeu
Le 7 févr. 2011 à 18:03, Bruno Chareyre a écrit : > , >> Say that you have cohesion between particles, then if you project the >> normal contact force along the normal to the plane you can see that >> the sign of the force could be either negative or positive (depending >> also on the convention)

[Yade-dev] Fwd: Another stress in bodies

2011-01-20 Thread Vincent Richefeu
> Note that the integral is computed on deformed contour. > In , is the contact point, closer to center > than radius. Hence, dividing by a volume of the Voronoi sort would not be > consistent with the integration domain. I don't see what other volume could > be used than the one of the sphere

Re: [Yade-dev] Another stress in bodies

2011-01-20 Thread Vincent Richefeu
>> You're right! Some people consider granular matter as a continuous media. I >> don't. >> > The integral is for one grain alone. Considering one grain as a > continuum doesn't sound so absurd, does it? > We are far above the atomic scale. As I said, I see no problem for the computation. I jus

Re: [Yade-dev] Another stress in bodies

2011-01-20 Thread Vincent Richefeu
Le 20 janv. 2011 à 13:01, Bruno Chareyre a écrit : > Hi Vincent, > > I didn't see the need to find a reference for contour integral of Cauchy > stress, since it is something classical in continuum mechanics courses. The > demonstration is in function's documentation. You're right! Some people

Re: [Yade-dev] Another stress in bodies

2011-01-19 Thread Vincent Richefeu
Hi, Have a look to this paper: http://www.springerlink.com/content/l133vr26m7142101/ I can send you the pdf if necessary. The stress computation at particle level is explained in section 5. Basically, it correspond from the dimensional point of view to the moment tensor divided by a volume. The

Re: [Yade-dev] mgpost

2010-12-15 Thread Vincent Richefeu
Le 15 déc. 2010 à 14:46, Václav Šmilauer a écrit : > Hi (Vincent and possibly others), > > would you mind if I move mgpost to a separate branch (still under the > project) at launchpad? I am trying to reduce size of sources, and mgpost is > really a separate project. Yes, you can move it to s

Re: [Yade-dev] twist plastic moment

2010-12-10 Thread Vincent Richefeu
Le 8 déc. 2010 à 18:58, Václav Šmilauer a écrit : > In the end, I wrote MatchMaker, which does just that: you can define value > for any given couple of materials, and in addition a fallback function/value > (such as average, minimum, maximum etc) if that particular combination is not > handle

Re: [Yade-dev] twist plastic moment

2010-12-08 Thread Vincent Richefeu
Hi everybody, Just a few comments (as usual, do not feel obliged to agree with me!) I understand (I think) what Chiara means. This is something that bothers me a lot (mentally) and, unfortunately, it leads me to continue using my poor bits of code instead of yade (please, do not be offended by

Re: [Yade-dev] energy conservation?

2010-11-14 Thread Vincent Richefeu
Hi, Some comments (useful or not): - If it's possible, try using velocity-Verlet integrator so that the the value of velocity is computed at time step (not mid-step) - Time integration (trapezoidal rule) involves approximations. Vaclav is right in believing that this error accumulation should de

Re: [Yade-dev] periodic weighted voronoi tesselation library (c++, open-source)

2010-10-29 Thread Vincent Richefeu
Hi, I personally use voro++. I think it is more appropriate (simple) if you only need the cells with circumcenter-vertices (instead of particle-center-vertices). In voro++ there is not iterator to loop over edges, faces or something else. Instead, you need to know the internal structure. The

Re: [Yade-dev] [Yade-users] Seminar of Vaclav this Friday at lab 3S-R, Grenoble, France

2010-03-10 Thread Vincent Richefeu
> Why is it so late to announce ;) ? Sorry for that! I did not think to make the announcement also on yade-dev and yade-users since I didn't do that for the other seminars. Jerome Duriez (thank you!) advised me to make this announcement here by remarking that some people could come (in particul

[Yade-dev] Seminar of Vaclav this Friday at lab 3S-R, Grenoble, France

2010-03-10 Thread Vincent Richefeu
Hi everybody, A announcement for those who are in Grenoble. Vaclav is here since this week (even if I have not seen him yet!). He will give a seminar on Friday 12th at 14:00 about Yade. Here is the announcement (in French): http://www.3s-r.hmg.inpg.fr/3sr/spip.php?article468 Vincent

Re: [Yade-dev] EigenDecomposition wrapper problem

2010-02-08 Thread Vincent Richefeu
ust let it be? > __ > > Anton Gladkyy > > > 2010/2/8 Vincent Richefeu > if numerical recipes, numpy and eigen give the same result, we can have some > confidence with this result. > It means that Wm3 result is "less good". > Howev

Re: [Yade-dev] EigenDecomposition wrapper problem

2010-02-08 Thread Vincent Richefeu
new values" it will not work? > __ > > Anton Gladkyy > > > 2010/2/8 Vincent Richefeu > Hi, > > At first glance, it seems the eigenvalues (and thus the eigenvectors) are not > ordered with Wm3. > The difference can also be due to the fact that th

Re: [Yade-dev] EigenDecomposition wrapper problem

2010-02-08 Thread Vincent Richefeu
Hi, At first glance, it seems the eigenvalues (and thus the eigenvectors) are not ordered with Wm3. The difference can also be due to the fact that the solution is not converged. If needed, I can send a peace of code (based on "jacobi" in numerical recipes) for doing that task VR Le 8 févr.

Re: [Yade-dev] [yade dev] About blocked rotation

2010-01-28 Thread Vincent Richefeu
Hi Vaclac, Have you think about vector for this purpose? As you certainly know it's a specialized template that hold bits (instead of bools). In fact, I have not seen the code, but in the c++ point of view, something like : std::vector blockedDOFs; blockedDOFs.reserve(3); // Actually, I do

Re: [Yade-dev] KLA

2010-01-12 Thread Vincent Richefeu
Hi Bruno, I have just commit a new version of KinematicLocalisationAnalyser::DefToFile. The code were NOT checked and I don't know if it compile (I can not compile on my apple laptop). Please send me an image from paraview when it will work... For everybody: If it crash (I sorry in advance) comm

Re: [Yade-dev] [Branch ~yade-dev/yade/trunk] Rev 1932: - merge stiffness/inertia control in PeriTriaxEngine and remove periEngine

2010-01-04 Thread Vincent Richefeu
Le 4 janv. 2010 à 13:36, Václav Šmilauer a écrit : >>> - with log, there is no upper limit to strain (except the grain >>> size/period size condition), you can compress 100, 200, 300%. With your >>> definition, e=100% has no meaning. >> e=-100% ;-) > (This BTW is not really "valid" argument...

Re: [Yade-dev] periodic cell

2009-12-22 Thread Vincent Richefeu
> I was surprised by the elegance of using unsheared coordinates without > changing the collider. It limits maximum shear, though, as the 1/cos(φ) > term grows without bounds. Salut Vaclav (désolé pour l'absence d'accent sur le a !), Un but de l'espace virtuel 's' dont te parlait Gael est justem

Re: [Yade-dev] PeriTriaxController

2009-12-10 Thread Vincent Richefeu
Bruno, I am not as confident as you. Imagine spheres that do not overlap on a regular grid. If you impose a shear rate at constant volume (without gravity) by making the displacement of each sphere proportional to their height (the total height of the cell remain constant), I think that the sph

[Yade-dev] Fwd: [Yade-users] From Blender2YADE

2009-12-10 Thread Vincent Richefeu
Début du message réexpédié : > De : Vincent Richefeu > Date : 10 décembre 2009 09:43:19 HNEC > À : Sergei D. > Objet : Rép : [Yade-users] From Blender2YADE > > Hi (Anton), > > Could you please briefly specify what is imported into what? > I suppose that onl

Re: [Yade-dev] some smaller changes

2009-12-09 Thread Vincent Richefeu
Le 9 déc. 2009 à 16:49, Václav Šmilauer a écrit : > Hi, I am in the process of writing programmer's manual to yade and some > idiosyncracies are showing up. > > 1. I will rename StandAloneEngine -> GlobalEngine (since it describeds > what it really does, globally affecting the scene). But what

[Yade-dev] [Yade-users] Questions on equations (contact laws, velocity, libraries)

2009-12-09 Thread Vincent Richefeu
> De : Vincent Richefeu > Date : 9 décembre 2009 09:35:42 HNEC > À : boon chiaweng > Objet : Rép : [Yade-users] Questions on equations (contact laws, velocity, > libraries) > > > Le 9 déc. 2009 à 03:06, boon chiaweng a écrit : > >> Q3. The position/velo

Re: [Yade-dev] PeriTriaxController

2009-12-08 Thread Vincent Richefeu
Le 8 déc. 2009 à 15:18, yade a écrit : > Bruno Chareyre a écrit : >> Václav Šmilauer a écrit : I'm not sure I understand completely. First reaction : Why not change "goal" with time, in comp=M*goal? It would give a strain rate. >>> Goal is user-defined and shouldn't

Re: [Yade-dev] new wiki & DNS

2009-10-11 Thread Vincent Richefeu
I prefer yade-dem.org - yade-open-dem -> to long - yade-project and yade-code -> we don't see explicitly what the project/code is for Vincent Le 11 oct. 2009 à 15:30, Bruno Chareyre a écrit : Very nice! There is not more than 1 vote per name for now! :) I voted yade-open-dem, but I like ya

[Yade-dev] Fwd: [Yade-users] New Yade wiki

2009-09-17 Thread Vincent Richefeu
Here is a mail that Rémi sent to me. I added an attempt of translation for people who do not understand French. Début du message réexpédié : De : Rémi Cailletaud Date : 16 septembre 2009 18:13:22 HAEC À : Vincent Richefeu Objet : Rép : Fwd: [Yade-users] New Yade wiki Salut, Vincent

[Yade-dev] Questions on the integration method

2009-09-16 Thread Vincent Richefeu
Hi everybody, This morning I had a look at the file LeapFrogPositionIntegrator.cpp and I a little scarred of what I saw (I exaggerate of course). Can someone enlighten me on that point: If I translate the c++ code into human readable operations, I've got: v(t+dt) = v(t) + a(t)*dt p(t+dt) = p(

[Yade-dev] Questions about Leap-Frog position integrator.

2009-09-16 Thread Vincent Richefeu
(I re-send this mail because it didn't appear in yade-dev list... sorry if you receive it for a second time) Hi everybody, This morning I had a look at the file LeapFrogPositionIntegrator.cpp and I a little scarred of what I saw (I exaggerate of course). Can someone enlighten me on that po

[Yade-dev] Move the wiki: script wiki_import

2009-09-04 Thread Vincent Richefeu
Hi, I address this message to yade-wikia administrator (Vaclav?). As you know, mediaWiki is installed at Grenoble and the ports will be soon open for access from outside. It's time to think about importing the database. Can you link up with Remi Cailletaud (remi.caillet...@hmg.inpg.fr) for

Re: [Yade-dev] ld: duplicated symbol (Mac OS X)

2009-09-02 Thread Vincent Richefeu
$ scons scons: Reading SConscript files ... @@@ Using profile default (scons.profile-default) @@@ Yade version is `vr', installed files will be suffixed with `-vr'. All intermediary files will be in `/Users/richefeu/Documents/yade- local/build-vr'. Checking whether c++ compiler "g++" works...no

Re: [Yade-dev] ld: duplicated symbol (Mac OS X)

2009-09-01 Thread Vincent Richefeu
I have no idea on that, and googling that didn't find anything really useful (besides extern inline you mentioned, which is not the case here). I think you have problem with miniWm3 just because it is the first lib that you link (right?) and it would be the same anywhere else. No, for instance

[Yade-dev] ld: duplicated symbol (Mac OS X)

2009-09-01 Thread Vincent Richefeu
Hi, I currently try (and try again) compiling yade on Mas OS X. Lot of problems have been resolved but I still have troubles to link miniWm3... I've seen with google that the problem could come from extern inline function, but in the present case, inline functions are not extern... I can n

Re: [Yade-dev] ld: new problems on Mac OS X

2009-08-29 Thread Vincent Richefeu
Le 29 août 09 à 14:21, Václav Šmilauer a écrit : Does runtimePREFIX could replace the origin option? No. Besides that, I think runtimePREFIX will not be needed soon as we use $ORIGIN now. Let's not confuse linking stage (at compilation) and running (when libraries are being found also). If $O

Re: [Yade-dev] ld: new problems on Mac OS X

2009-08-29 Thread Vincent Richefeu
maybe he could try compiling monolithic version? Thank you Janek for your idea... Would be the same IMO. Unfortunately Vaclav is right Or even the static one? $ scons-h ## :-( ok thank you, I will try scons LIBPATH=... Does runtimePREFIX could replace the origin option?

[Yade-dev] funny...

2009-08-28 Thread Vincent Richefeu
Have you ever heard of Yade X? This is Yet Another Dvd eXtractor: http://www.netfox2.net/modules/wfdownloads/singlefile.php?cid=45&lid=455 VR ___ Mailing list: https://launchpad.net/~yade-dev Post to : yade-dev@lists.launchpad.net Unsubscribe : h

[Yade-dev] ld: new problems on Mac OS X

2009-08-28 Thread Vincent Richefeu
Hi, I still have problem on linkage stage on my apple laptop. I have tried lot of thinks without results... Please give me new ideas!!! here is the log: $ scons pretty=0 scons: Reading SConscript files ... @@@ Using profile default (scons.profile-default) @@@ Yade version is `vr', installed fil

Re: [Yade-dev] scons and linkage rules

2009-08-27 Thread Vincent Richefeu
Oh, if apples use different extension for shared object, have a look Omega.cpp:193, where we check what files are plugins. It will compile just fine, but will not load any plugins at startup. You can therefore put some other checks there, the point is to not fix in advance where we have plu

Re: [Yade-dev] Fwd: Re: scons and linkage rules

2009-08-27 Thread Vincent Richefeu
"_FSGetCatalogInfo", referenced from: Wm3::System::GetPath(char const*, char const*)in miniWm3.os Wm3::System::GetPath(char const*, char const*)in miniWm3.os Wm3::System::GetPath(char const*, char const*)in miniWm3.os "_GetProcessBundleLocation", referenced from: Wm3::System::GetPath(c

Re: [Yade-dev] scons and linkage rules

2009-08-27 Thread Vincent Richefeu
g++ -o /Users/richefeu/Documents/yade-local/build-vr/lib/ libminiWm3.dylib -rdynamic -z origin -dynamiclib -rdynamic /Users/ richefeu/Documents/yade-local/build-vr/lib/miniWm3.os -L/Users/ richefeu/Documents/YADE/lib/yade-vr/extra -L/Users/richefeu/ Documents/YADE/lib/yade-vr/gui -L/Users/riche

[Yade-dev] scons and linkage rules

2009-08-27 Thread Vincent Richefeu
I assume you mean CLOCK_MONOTONIC -- it is clock source that ensures that it is, well, monotonic (i.e. resistent to leap seconds, changing hw clock etc). OK I don't know what works on leopard, try man clock_gettime. STFW gives http://www.wand.net.nz/~smr26/wordpress/2009/01/19/monotonic-ti

[Yade-dev] CLOCK_MONOLITIC portability

2009-08-27 Thread Vincent Richefeu
Hi, I'm currently trying to compile yade on Mac OS X (Leopard). for the moment, all optional packages are disabled and also qt3. My problem is about clock_gettime(CLOCK_MONOLITIC) that seem to be not support by this OS. It appears only one time in the file core/Timing.hpp coded by Vaclav. Do

[Yade-dev] Strange Behavior: an empty simulation auto-starts!

2009-08-26 Thread Vincent Richefeu
Hi everybody, I checked out the version 1736 with bzr and compiled it. I report here a strange behavior: when running yade without arguments (./yade-version), an "empty simulation" begins to run and the control buttons are disabled! Do you observe the same behavior? If yes, is this normal?

Re: [Yade-dev] WIKIA advertisements

2009-07-30 Thread Vincent Richefeu
> It seems there are a lot of advertisements on WIKIA last days. Is it > not? > > Yesterday I even could not click the right link because of that: first > I needed to click ads :( > They even started to insert their ads directly to our articles! I didn't see any advertisements. Ever heard of "

Re: [Yade-dev] yade crash launching xml

2009-07-15 Thread Vincent Richefeu
boost::shared_ptr::operator->() const [with T = Engine]: L'assertion « px != 0 » a échoué. Hi Emanuele, I imagine the assert(px != 0) is into your class FlowEngine. If it is the case, try to check the reason why you've got px == 0. If it is not the case, my comment will not help. Ciao, V

Re: [Yade-dev] Material and body State

2009-07-15 Thread Vincent Richefeu
Le 15 juil. 09 à 17:01, Václav Šmilauer a écrit : As I said before, the solution described here before seems to be nice (although the question of random-values dispatching remains, cf. Bruno's comments). As written, cf. http://yade.wikia.com/wiki/ SharedMaterialProperties, it is enhanced i

Re: [Yade-dev] Material and body State

2009-07-15 Thread Vincent Richefeu
Le 14 juil. 09 à 09:52, Václav Šmilauer a écrit : Hi, based on discussion with Vincent and also (older) proposal of Janek, and for other reasons of simplicity, I would like to create new Material class that would hold material properties of bodies. There would be only a few instances of

Re: [Yade-dev] Material and body State

2009-07-14 Thread Vincent Richefeu
Hi, Your proposal sounds good for me. Actually I have not studied it in details. I will re-answer your post in a few days (There's a lot of informations!). Just a remark that comes to me about naming: frictionAngle does not only relate to 'material' but also to 'surface state'. Ciao, V

Re: [Yade-dev] strided collider speedups

2009-07-11 Thread Vincent Richefeu
Actually, Vincent was suprised yesterday when he realised this strategy was not yet used in Yade. Yes, this is a very common and basic way to save a lot time in DEM simulations! The underlying idea is that no computation is performed when it is not necessary (and there is no better optimiza

Re: [Yade-dev] [deprecated list] Special serialization treatment

2009-07-08 Thread Vincent Richefeu
Vincent, to me, the GroupRelationData concept is not very clear. Can I have some questions? 1. Why not have just one table inside, and every engine that needs the data could create as much instances as necessary, one for every parameter? It would seem more logical than having that inside MetaBo

Re: [Yade-dev] [svn] r1843 - trunk/core

2009-07-07 Thread Vincent Richefeu
Does it mean, that if I have 3 particles from 3 different materials I can make the next "database": P1-P2 - Kn(1-2); Kt(1-2)...; P2-P3 - Kn(2-3); Kt(2-3)...; P1-P3 - Kn(1-3); Kt(1-3)...; P1-P1 - Kn(1-1); Kt(1-1)...; P2-P2 - Kn(2-2); Kt(2-2)...; P3-P3 - Kn(3-3); Kt(3-3)...; Particle interaction

Re: [Yade-dev] [deprecated list] Special serialization treatment

2009-07-06 Thread Vincent Richefeu
When a string is serialized, all space characters become   Is this normal? How can I do to avoid this? thanks, VR ___ Mailing list: https://launchpad.net/~yade-dev Post to : yade-dev@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-dev

Re: [Yade-dev] [deprecated list] Special serialization treatment

2009-07-06 Thread Vincent Richefeu
Le 6 juil. 09 à 12:11, Janek Kozicki a écrit : Are you trying to serialize a C array? It's not supported. And won't ever be. You don't want to add support for it. Easiest if you switch to std::vector<> and it will work straight away. The class I want to serialize is GroupRelationData. The pro

[Yade-dev] Special serialization treatment

2009-07-03 Thread Vincent Richefeu
Janek, Could you please remind me where I could include the code to perform special treatment for the serialization? I hope my question is claire enough. Practically, I would like to do something like this: if (the type of my object is 'GroupRelationData') { I write what I want in xml for

Re: [Yade-dev] [svn] r1822 - in trunk: pkg/dem/Engine/StandAloneEngine scripts scripts/test

2009-06-29 Thread Vincent Richefeu
1. Rename ef2_Dem3Dof_Elastic_ElasticLaw to Law2_Dem3Dof_Elastic_Elastic (in scripts as well) Hi, Some news about naming convention? VR ___ Mailing list: https://launchpad.net/~yade-dev Post to : yade-dev@lists.launchpad.net Unsubscribe : https:/

Re: [Yade-dev] [Yade-users] YADE documentation

2009-06-24 Thread Vincent Richefeu
to have 1 website with all YADE stuff. PPS domain name yade.info is available __ [ENG] Best Regards [GER] Mit freundlichen Grüßen [RUS] С наилучшими пожеланиями [UKR] З найкращими побажаннями Anton Gladkyy 2009/6/24 Vincent Richefeu Hi, If no objection, I propose

Re: [Yade-dev] [Yade-users] YADE documentation

2009-06-24 Thread Vincent Richefeu
Hi, If no objection, I propose to create a second wiki (on wikidot.com) exclusively dedicated to yade documentation. I personally use wikidot and I find it much more user-friend than wikia. I could "design" and propose (in a first stage) a structure for the site. Then, any contributor will

[Yade-dev] Fwd: [Yade-users] YADE documentation

2009-06-24 Thread Vincent Richefeu
Début du message réexpédié : De : Vincent Richefeu Date : 24 juin 2009 11:23:06 HAEC À : yade-us...@lists.launchpad.net Objet : Rép : [Yade-users] YADE documentation Hi, If no objection, I propose to create a second wiki (on wikidot.com) exclusively dedicated to yade documentation. I

Re: [Yade-dev] build crash

2009-06-23 Thread Vincent Richefeu
Hello I have this error in yade.cpp (related to log4cxx), after an svn update, anyone got this before? Thanks Bruno C core/yade.cpp core/yade.cpp:45: error: wrong number of arguments specified for 'constructor' attribute l.45 is : __attribute__((constructor(1000))) void initLog4cxx(

[Yade-dev] Fwd: Another script with clumps

2009-06-18 Thread Vincent Richefeu
-- Forwarded message -- From: Vincent Richefeu Date: 2009/6/18 Subject: Another script with clumps To: eudo...@arcig.cz Another script... VR # Parameters Kn = 1e6 Ks = 0.8e6 FrictionAngle = 0.5

[Yade-dev] Fwd: python script with clumps

2009-06-18 Thread Vincent Richefeu
-- Forwarded message -- From: Vincent Richefeu Date: 2009/6/18 Subject: python script with clumps To: eudo...@arcig.cz The attached file is a yade-python script that makes use of clumps. (my python style should be funny for you!) # Parameters

[Yade-dev] Total moment at the mass center of clumps

2009-06-18 Thread Vincent Richefeu
Hello everybody, I try to use clumps for sometimes and it seems to me that there is a small "bug" (maybe I'm wrong...). When clumps are used, they are part of the bodyContainer (like common bodies). But when interaction forces are computed, it seems that the resulting moment is computed at

Re: [Yade-dev] [svn] r1788 - in trunk: core gui/py lib/serialization pkg/dem pkg/dem/Engine/StandAloneEngine

2009-06-03 Thread Vincent Richefeu
2. Replace vector by Vector2 and make it properly serializable (should be faster) Hi, Quick question: Is it possible to serialize a vector? VR ___ Mailing list: https://launchpad.net/~yade-dev Post to : yade-dev@lists.launchpad.net Unsubscribe

Re: [Yade-dev] [svn] r1784 - trunk/extra/SpherePadder

2009-05-29 Thread Vincent Richefeu
Add a function to SpherePadder_wrapper to pass the result back to python. As Vaclav suggested, it is a list of (x,y,z,R) tuples. (Complet integration into yade is not done yet) You can keep your module independent of yade, it can be loaded through the built-in python interpreter just fine and

Re: [Yade-dev] Exposing SpherePadder to python

2009-05-28 Thread Vincent Richefeu
Easy. You just need to link with boost_python library and the shared object should be called something.so (not libsomething.so): g++ src1.cpp src2.cpp -shared -o something.so -lboost_python - lCGAL ... I fill that I begin to love python... Sorry for my poor english. I wanted to say 'I fee

Re: [Yade-dev] Exposing SpherePadder to python

2009-05-27 Thread Vincent Richefeu
Easy. You just need to link with boost_python library and the shared object should be called something.so (not libsomething.so): g ++ src1.cpp src2.cpp -shared -o something.so -lboost_python -lCGAL ... I fill that I begin to love python... If you want, I could make scons detect CGAL (in SCons

Re: [Yade-dev] Exposing SpherePadder to python

2009-05-27 Thread Vincent Richefeu
3. SpherePadder should be added to scons (as an optional part), since it will also install the python module tothe appropriate place. OK, thank you for the clear answer. I found some tutorial on boost::python on the web, but I'm not clear with the compilation procedure (and I would like to av

[Yade-dev] Exposing SpherePadder to python

2009-05-27 Thread Vincent Richefeu
Hi, Thank you Vaclav for your previous answer. I've got an other question about boost::python and dynamic libraries widely used in yade. I would like to expose the class SpherePadder (and maybe TetraMesh) to python. I think that it could be possible (?) to add something like: class_()

[Yade-dev] Doxygen docs

2009-05-27 Thread Vincent Richefeu
Hi, Does someone can tell me if the doxygen docs are regularly generated? And where I can find it? Thanks, VR ___ Mailing list: https://launchpad.net/~yade-dev Post to : yade-dev@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-d

Re: [Yade-dev] Relations between (m,tc,en,es) and (kn,cn,ks,cs)

2009-05-06 Thread Vincent Richefeu
we have tc and en: (*) tc = pi 2 m_eff / sqrt( 4 kn m_eff - cn^2 ) (note your delta in brackets) en = exp( -cn tc / 2 / meff ) where we have kn and cn: kn = m_eff / tc^2 * ( pi^2 + ln^2(en) ) cn = -2 * m_eff / tc * ln(en) Then, for pair spheres with mass m1 we have m1_eff = m1/2, en1, tc1 =>

Re: [Yade-dev] Relations between (m,tc,en,es) and (kn,cn,ks,cs)

2009-05-06 Thread Vincent Richefeu
Yes, en and es are normal and tangential coefficients of restitution (ratio of velocities after and before the 'impact': vf/vi). I wanted to say that the parameters in SimpleViscoelasticBodyParameters could be the stiffnesses and restitution coefficients. Then, the effective mass together with the

Re: [Yade-dev] Relations between (m,tc,en,es) and (kn,cn,ks,cs)

2009-05-06 Thread Vincent Richefeu
Thanks a lot. Now, I understand a bit more... The implanted law is for sphere-sphere contacts. I suggest something (tell me your opinion please): what about computing cn and cs for each contact as a function of the corresponding effective mass? In this way, the parameters would be kn,en,ks, and es

Re: [Yade-dev] Relations between (m,tc,en,es) and (kn,cn,ks,cs)

2009-05-05 Thread Vincent Richefeu
Vincent Richefeu ?: I suppose this means that you would like to learn about me (?) So briefly, I'm French and I work in the laboratory 3S-R since a few months. My PhD is about unsaturated granular materials simulated by DEM. You can find more with google if you want... Hi,

[Yade-dev] Relations between (m,tc,en,es) and (kn,cn,ks,cs)

2009-05-05 Thread Vincent Richefeu
he current law, but I need more explanations, so please... Thanks, Vincent Richefeu ___ Mailing list: https://launchpad.net/~yade-dev Post to : yade-dev@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.lau

Re: [Yade-dev] Python scripts and clump creation

2009-05-05 Thread Vincent Richefeu
said that python syntax is obscure but I know that the problem come from me. I have to think different... Thank for adding the clump wrapper. It works fine. Regards, Vincent Richefeu ___ Mailing list: https://launchpad.net/~yade-dev Post to :

[Yade-dev] Python scripts and clump creation

2009-05-04 Thread Vincent Richefeu
to access the method 'createOneClump' (for example) from a python script. Does anyone (Vaclav ?) can help me? Thanks, Vincent Richefeu ___ Mailing list: https://launchpad.net/~yade-dev Post to : yade-dev@lists.launchpad.net

Re: [Yade-dev] [Yade-users] Be carefull with "normalize"

2009-02-05 Thread Vincent Richefeu
Hi all, I've got just a comment: a simple solution to avoid this (without modifying [mini]wm3) is to change the user units. I think yade computation can be done with any system of units as they are consistent. For example, one can use centimeters instead of meters and deca- Newtons instead o

[Yade-dev] Fwd: inherited PhysicalParameters for InteractionPhysicsEngineUnit

2008-11-27 Thread Vincent Richefeu
UP ? I've got a problem with 'SimpleViscoelasticRelationships'. This 2d-functor needs 2 'PhysicalParameters' that are casted into 'SimpleViscoelasticBodyParameters' in the method 'go'. I want to use BcpConnection that inherit from SimpleViscoelasticBodyParameters for one of the parameters.

Re: [Yade-dev] inherited PhysicalParameters for InteractionPhysicsEngineUnit

2008-11-21 Thread Vincent Richefeu
I've got a problem with 'SimpleViscoelasticRelationships'. This 2d-functor needs 2 'PhysicalParameters' that are casted into 'SimpleViscoelasticBodyParameters' in the method 'go'. I want to use BcpConnection that inherit from SimpleViscoelasticBodyParameters for one of the parameters. The probl

[Yade-dev] inherited PhysicalParameters for InteractionPhysicsEngineUnit

2008-11-20 Thread Vincent Richefeu
Hello everybody, I've got a problem with 'SimpleViscoelasticRelationships'. This 2d-functor needs 2 'PhysicalParameters' that are casted into 'SimpleViscoelasticBodyParameters' in the method 'go'. I want to use BcpConnection that inherit from SimpleViscoelasticBodyParameters for one of the parame

Re: [Yade-dev] (x,y,z) convention

2008-11-19 Thread Vincent Richefeu
Hi everybody, I have just a brief comment concerning this discussion. I have write a visualizer for DEM simulations (mgpost, see google). This tool could be incorporated 'into' Yade in the future. It may be interesting to known that I used the 'vertical Y' convention. Vincent Le 18 nov. 08 à

Re: [Yade-dev] Damping at the contact level

2008-11-07 Thread Vincent Richefeu
Without going into sterile discussions, I would like to know if something like that exists in Yade. Otherwise, I will need to add the class 'NormalViscousInteraction' (and/or 'NormalViscousShearInteraction') that store the damping parameter. I will need also to add a relative normal velocity p

[Yade-dev] Damping at the contact level

2008-11-07 Thread Vincent Richefeu
Hi everybody, I am very surprised that no viscosity is taken into account at contact level in Yade. Maybe I'm wrong, but I have not seen anything like it! I do not usually use a global damping ('NonViscousDamping' in Yade) and I prefer to use local (at the contact level) damping that can be

Re: [Yade-dev] Incomprehensible crash

2008-10-30 Thread Vincent Richefeu
PS. You really _want_ to use those ugly names of classes? OMG... at least skip the _make part and write just _BssSweptSpheres..., no? I agree with you. I usually use short names but these long names aim to get a quick identification of the task performed by the plugin. Admittedly, your nam

Re: [Yade-dev] Incomprehensible crash

2008-10-30 Thread Vincent Richefeu
Re-enabling this line makes me run without crash :-) Thank you Vaclav. It was also necessary to remove the engine 'PhysicalActionContainerReseter' in order to have no crash... It's a pleasure to work with you because you are very reactive! PS. You really _want_ to use those ugly names of c

Re: [Yade-dev] Incomprehensible crash

2008-10-29 Thread Vincent Richefeu
The following is the beginning of the console log: It's better if you send the whole message. It seems hard to believe but some of us succeed in understanding something from it ! Jerome OK, sorry. here is the complet log: In [1]: SIGSEGV/SIGABRT handler called; gdb batch file is `/home/

[Yade-dev] Incomprehensible crash

2008-10-29 Thread Vincent Richefeu
Hello, I still got problems with 'MembraneTest'. It crash at the first time step. I really can't understand why! If someone has an idea, I'm taker... The files MembraneTest.* have been added in the svn repository into the folder pkg/dem/PreProcessor. The following is the beginning of the console

[Yade-dev] commit not reported!

2008-10-29 Thread Vincent Richefeu
Hi everybody, I added and committed new plugins. The commit is not reported in the mailling list. Instead, I received the following message from yade-dev-bounces: >>> Your mail to 'Yade-dev' with the subject [svn] r1564 - in trunk/pkg: common common/DataClass/GeometricalModel common/Dat

Re: [Yade-dev] Mandatory engines

2008-10-29 Thread Vincent Richefeu
Hi Vaclav, Thanks a lot for explanations. I will write the needed physicalAction (s) after. For the moment, I would simply like to test (with display) two new functors: - the one make interactingGeometry (BssSweptSphereLineSegment) from GeometricalModel (BshTube), - the other make AABB from

[Yade-dev] Mandatory engines

2008-10-27 Thread Vincent Richefeu
Hello everybody, I'm Vincent, a new yade-dev member. I am currently in a learning phase in yade. * HERE IS THE STORY: Say I want to implement a new geometry for bodies. I code pluggins that I needed (GeometriclModel, InteractingGeometry,Functors, GLDraw*,Preprocessor...) after discussing