yade-users team mailing list archive

Thread
Date

Re: [Question #686131]: flow.nCells() != "cell number" recorded by flow.saveVtk()

To: yade-users@xxxxxxxxxxxxxxxxxxx
From: Robert Caulk <question686131@xxxxxxxxxxxxxxxxxxxxx>
Date: Sat, 23 Nov 2019 13:17:48 -0000
Reply-to: question686131@xxxxxxxxxxxxxxxxxxxxx
Sender: bounces@xxxxxxxxxxxxx

Question #686131 on Yade changed:
https://answers.launchpad.net/yade/+question/686131

    Status: Open => Answered

Robert Caulk proposed the following answer:
Thanks for the clarification.

> how many finite cells are produced by the triangulation, 383 or 242?

242 withoutBoundaries (default). Try using
flow.saveVTK(withBoundaries=True), the numbers should match.

>the particle ids comprising the cell id=0 are 38, 41, 37, 39. Am I
wrong?

Yes, you are wrong :-) Those are not the particle IDs. Those are
redirected IDs for paraview since it needs zero based vertex IDs [1]

[1]https://gitlab.com/yade-
dev/trunk/blob/master/lib/triangulation/FlowBoundingSphere.ipp#L1690

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.