Question #698778 on Yade changed: https://answers.launchpad.net/yade/+question/698778
Robert Caulk posted a new comment: Hey Zoheir, I don't think you are missing anything. It seems that the pseudocode I wrote here is incomplete. Looking a bit closer at the source, I am unsure we can get the facet fluxes in python with the current wrapping. Do you *really* need this in python? The only reason I suggested this pseudocode in python was so that you could double-check the function that I wrote in C++ [1]. Is [1] not working for you? Apparently you said "Now, the influx and outflux are the same." which tells me there is no need to be double-checking this in python. Please clarify your needs. If you for some reason need the facet fluxes to be wrapped in python, then I can do that - but it needs justification first. For some reason, I started chasing geese again on this problem, before realizing that this was in fact solved, and not a problem, when you said that you got influx and outflux to be the same using flow.getCavityFlux()....oh well, that is my fault for not rereading the entire thread. Anyways, during that goose chase, I visualized your cavity with flow.saveVTK(). The cavity appears to be scattered and incorrectly assigned, based on the VTK. Although, perhaps you have an updated script where you made sure the cavity was properly defined. If you did not change the cavity assignment, I would not expect any cavity fluxes to be properly estimated. -rc [1]https://gitlab.com/yade- dev/trunk/-/blob/master/lib/triangulation/FlowBoundingSphere.ipp#L1310 Cheers, Robert -- 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
