(fenicsx) subdomain integral fails depending on form of meshtags

dokken · January 10, 2023, 1:28pm

Try using dtype=np.int32 or np.intc, as this is what form supports.
See for instance

FEniCS/dolfinx/blob/ddcad2e3167858157d413f74781a2f458d34afb0/python/test/unit/mesh/test_meshtags.py#L48


      
          
          def test_ufl_id():
              """Test that UFL can process MeshTags (tests ufl_id attribute)"""
              comm = MPI.COMM_WORLD
              mesh = create_unit_cube(comm, 6, 6, 6)
              tdim = mesh.topology.dim
              marked_facets = locate_entities(mesh, tdim - 1, lambda x: np.isclose(x[1], 1))
              f_v = mesh.topology.connectivity(tdim - 1, 0).array.reshape(-1, 3)
          
              entities = create_adjacencylist(f_v[marked_facets])
              values = np.full(marked_facets.shape[0], 2, dtype=np.int32)
              ft = meshtags_from_entities(mesh, tdim - 1, entities, values)
              ds = Measure("ds", domain=mesh, subdomain_data=ft, subdomain_id=(2, 3))
              a = 1 * ds
              assert isinstance(a.subdomain_data(), dict)

and

github.com

FEniCS/dolfinx/blob/2be366d20632377e0688df585fcd7e201d5dc86c/python/test/unit/fem/test_assemble_domains.py#L57


      
          mesh = create_unit_square(MPI.COMM_WORLD, 10, 10, ghost_mode=mode)
          V = FunctionSpace(mesh, ("Lagrange", 1))
          u, v = ufl.TrialFunction(V), ufl.TestFunction(V)
          
          # Prepare a marking structures
          # indices cover all cells
          # values are [1, 2, 3, 3, ...]
          cell_map = mesh.topology.index_map(mesh.topology.dim)
          num_cells = cell_map.size_local + cell_map.num_ghosts
          indices = np.arange(0, num_cells)
          values = np.full(indices.shape, 3, dtype=np.intc)
          values[0] = 1
          values[1] = 2
          marker = meshtags_factory(mesh, mesh.topology.dim, indices, values)
          dx = ufl.Measure('dx', subdomain_data=marker, domain=mesh)
          w = Function(V)
          w.x.array[:] = 0.5
          
          # Assemble matrix
          a = form(w * ufl.inner(u, v) * (dx(1) + dx(2) + dx(3)))
          A = assemble_matrix(a)

Topic		Replies	Views
V0.6.0 container subdomain multiple variable General	1	172	September 26, 2023
Integrating on subdomain	30	2355	November 11, 2022
Subdomain based on cell index	9	794	October 19, 2022
Subdomains, problems following Tutorial General	4	210	March 30, 2024
Transfer meshtags to submesh in dolfinx mesh	9	1091	December 22, 2023

(fenicsx) subdomain integral fails depending on form of meshtags

Related topics