On 03/31/2016 02:49 PM, Adolfo Rodriguez wrote:
How did you create the region? Using Domain.create_region() or in some other way?
An example script would be handy.
r.
Robert,
I am using
frac0_region = Region.from_vertices(frac0_vertices,domain,'Fracture0', 'facet')
where frac0_vertices is an array containing the nodes I want to include in the region.
Then you need to tell the Domain instance about that region - try appending it to the domain.regions list. Basically, you need to manually do all the things Domain.create_region() does after the visit_stack() call.
So, in the end, it might be easier to just call Domain.create_region() :) Maybe with the 'vertices by function' selector, where the function would just return the vertices.
But I see that some shortcut for creating a fully initialized region might be handy - maybe Region.finalize() could be changed to do just that.
r.