New question #670610 on Yade: https://answers.launchpad.net/yade/+question/670610
Hi, I am trying to tessellate particle geometry generated with Yade. For this purpose I am using the TwoPhaseFlowEngine. When I run the TwoPhaseFlowEngine on geometry generated directly from makeCloud everything runs smoothly. However, when I try to run the TwoPhaseFlowEngine on geometry which I have generated to simulate an aggregate. I receive a segmentation fault. In order to track down this issue, I have written the code below. This code reads a packing of 113 spheres. If I run the TwoPhaseFlowEngine at this point, there is no error, however looking at the exported pore network, only 8 pores are generated by the tesselation, and they are all located at the same location (see here: http://pasteall.org/1012264 ) In order to try and dig into the issue. I wrote the following code. First I load the sphere geometry. I try to run the tesselation. I then save the individual sphere locations and radii and add them back in one at a time to the simulation. At every step I try to run the tesselation. In this case, when I add the 5th particle I receive a Segmentation fault (core dumped) and the contents of PoreBodyLocation.txt is a list of 0,0,0 's . from yade import export,ymport sp = ymport.textClumps('/tmp/clump0.txt') print('Number of particles = ' , len(O.bodies)) TwoPhase=TwoPhaseFlowEngine() TwoPhase.initialization() TwoPhase.savePoreNetwork('networkTest') #raw_input('Pause') spheres = [] for b in O.bodies: if b.isClump == False: spheres.append([b.state.pos, b.shape.radius]) O.resetThisScene() for pos,r in spheres: O.bodies.append(sphere(pos,r)) print('Number of particles = ' , len(O.bodies)) #raw_input('Pause') TwoPhase.initialization() TwoPhase.savePoreNetwork('networkTest') try: from yade import qt qt.View() except: pass and clump0.txt at http://pasteall.org/1012270 Does anyone know the reason for this segmentation fault, and or the reason for the pore locations = 0,0,0 ? many thanks in advance Jesse -- You received this question notification because your team yade-users is an answer contact for Yade. _______________________________________________ Mailing list: https://launchpad.net/~yade-users Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
