Hello, I’m working with an internal boundary integration as the right side of a variational system: V = FunctionSpace(mesh, 'CG', 1) n = FacetNormal(mesh) S = 2.0*mu*inner(e(u),e(u)) + lmbda*div(u)**2 theta,xi = [TrialFunction(V), TestFunction(V)] left_side = assemble((inner(grad(theta),grad…

Hi Nate! Here it is. Thank you! from dolfin import * import numpy as np lx,ly = [2.0,1.0]; Nx, Ny = [62, 31] def phi(x): return (x[0] - 1.0)**2 + (x[1]-0.5)**2 - 0.2**2 class Omega(SubDomain): def inside(self, x, on_boundary): return phi(x) > 0 mesh = RectangleMesh(Point(0.0,…

You need to choose on which side of the facet to evaluate xi, e.g. xi("+").

Thank you! This works, but I’m still a bit confused on what this restriction means, since I’m interested in the solution th defined on the whole domain. If I choose (’+’), will I somehow be concentrating the solution on “this side”? Also, could you maybe give me some insight on the element dS? Is …

Discontinuous type Argument must be restricted

Errors

nate July 9, 2021, 6:29pm 2

If you make this a MWE, we can help debug it.

Topic		Replies	Views
What does the "must be restricted" error mean? dolfinx	11	59	April 7, 2026
Boundary conditions for internal facets dolfinx boundary-conditions	1	305	January 21, 2024
Discontinuous type Coefficient must be restricted type error General	3	237	May 30, 2024
Integrate a flux over an internal curve General	2	78	May 25, 2025
Applying loading on interior facets in L-form variational formulation dolfinx	6	94	August 16, 2024

Discontinuous type Argument must be restricted

Related topics