Maybe a more convenient way of getting global indices of subdomain objects

Maybe we’ve talked about it in the past. I think getting global indices of subdomain objects should be more convenient than doing:

network.map_pores(geom_A.Ps, origin=geom_A)

[jgostick]

Yes, but the only problem is that you must have labelled the pores on obj. We could enhance the pores method to accept an array as well as string. If string -> look up pores, if array -> just use them directly. Would be about 3 lines to add.

The only down side is that accepting such a wide variety of arguments starts to make the methods unpredictable and complex. For instance, we used to accept objects handles AND names in many functions, but decided to go with handles everywhere for consistency.

I’d be willing to make an exception this time though.

[ma-sadeghi]

I see, it’s just that obj.pores([1,3], target=network) to me seems more readable and self-descriptive than network.map_pores([1, 3], origin=obj), though it might as well just be my preference.

But, I see your point about consistency.