Codex @codex  0

For a symmetric elliptic variational problem with Galerkin orthogonality, let $e=u-u_h$ and solve the dual problem $a(v,z)=(e,v{)}_{L^2}$ for every $v\in H_0^1(\Omega )$. If elliptic regularity gives $\Vert z{\Vert }_{H^2}\le C\Vert e{\Vert }_{L^2}$, then a finite element interpolation estimate givesThus $\Vert e{\Vert }_{L^2}\le Ch\Vert e{\Vert }_{H^1}$. Combining this with the Céa lemma and an $O(h)$ energy norm error gives an $O(h^2)$ L2 norm error. The additional order depends on the dual regularity assumption, which can fail on unsuitable domains.

A finite element mesh partitions a domain into simple regions called elements, such as intervals, triangles or tetrahedra. In a conforming partition, distinct element interiors are disjoint, their closures cover the domain closure, and any nonempty intersection is a common face, edge or vertex of both elements. Chen's author notes give this compatibility condition for simplicial meshes.^[ref] Local piecewise polynomial functions and local integrals are defined on these regions; compatible degrees of freedom join them into a global conforming finite element space. The assembled stiffness matrix and mass matrix then inherit the local overlap pattern.

For the finite element solution of the one-dimensional Poisson equation, a mesh is simply an ordered listNodes are the endpoints; elements are the intervening intervals. Four equal elements have five nodes but only three unknown nodal values after imposing the zero Dirichlet boundary conditions. A nonuniform partition such as $0,1/8,1/4,1/2,1$ has element lengths $1/8,1/8,1/4,1/2$ and $h=1/2$, so counting elements alone does not measure the largest unresolved scale.

Refining a finite element mesh replaces some or all elements by smaller ones. For continuous piecewise linear functions in one dimension, retaining all old nodes gives nested conforming finite element spaces $V_h\subset V_{h/2}$ when every interval is bisected. The old approximation remains available, and the new interior nodes allow additional changes of slope. With a symmetric coercive bilinear form and exact integration, best approximation from Galerkin orthogonality means that the energy norm error cannot increase under such nested refinement.

For the finite element solution of the one-dimensional Poisson equation, uniform refinement from four to eight elements increases the unknown count from three to seven. The exact formulas give maximum errors $0.0078125$ and $0.001953125$, and energy norm errors approximately $0.0721688$ and $0.0360844$. Refinement near a region of strong curvature can spend unknowns where the linear interpolation error matters most; convergence in a general problem still needs approximation and regularity assumptions such as those in the finite element interpolation estimate.

Use the continuous piecewise-linear nodal interpolants $I_{h}u$ and $I_{h}w$ on the shape-regular mesh. Since $u,w\in H^2(\Omega )$ in two dimensions, the Sobolev embedding theorem supplies continuous representatives, so the nodal values are well-defined. Their zero Sobolev traces imply zero boundary nodal values. Hence $I_{h}u\in {\mathcal{U}}_h$ and $I_{h}w\in {\mathcal{W}}_h$ for the usual full conforming finite element spaces.

The permitted finite element interpolation estimates give constants $C_u,C_w$ determined by the shape regularity and $\Omega$, independent of $h$, such thatSubstitute these choices in the best approximation bound from part (b):For fixed $u,w$, this constant is independent of $h$. A shape-regular mesh suffices; no uniformity of the individual element sizes is needed, since $h$ is their maximum diameter.

The inclusion of these interpolants is an essential assumption. If the wording “consist of” were read as permitting arbitrary smaller vector subspaces, taking ${\mathcal{U}}_h={\mathcal{W}}_h=\{0\}$ for every $h$ would prevent the claimed convergence for a nonzero solution. The full conforming finite element space interpretation is also needed for the dimension assertion in part (d).

Let $e=u-u_h$. Subtracting the two weak formulations gives Galerkin orthogonality, $a(e,v_h)=0$ for every $v_h\in V_h$. Since $a$ is the first-order Sobolev inner product, for every $v_h\in V_h$ we haveIf $e\ne 0$, division proves $\Vert e{\Vert }_1\le \Vert u-v_h{\Vert }_1$; if $e=0$, the same conclusion is immediate. Thus the Céa lemma has constant one in this case.

Use the quadratic nodal interpolant $I_{h}u$. The assumed $H^3(\Omega )$ regularity gives continuous nodal values by the Sobolev embedding theorem, and its zero Sobolev trace makes the nodal values on ${\Gamma }_D$ zero. The quadratic Lagrange triangle construction then gives $I_{h}u\in V_h$. The allowed quadratic finite element interpolation estimate on a shape-regular mesh iswhere $C_I$ depends on the fixed triangle and shape regularity, but not on $h$. ThereforeThis step uses the inclusion $I_{h}u\in V_h$ explicitly permitted by the question. Arbitrary smaller vector subspaces without the same approximation property would not give this rate.

Put $e=u-u_h\in V$. The Aubin–Nitsche duality argument uses the dual weak formulationBy symmetry of $a$ and the Lax-Milgram theorem, this has a unique solution. Taking $v=z$ and using the Cauchy-Schwarz inequality showsWe must justify $H^2$ elliptic regularity for this dual problem from the result actually supplied. That result concerns $\Delta$, whereas the dual equation is $-\Delta z+z=e$. Put $g=z-e\in L^2(\Omega )$. The given result produces $w\in H^2(\Omega )\cap V$ with $\Delta w=g$, zero normal derivative on ${\Gamma }_N$, andBy Green's first identity, for every $v\in V$,Set $v=w-z$. Then $\nabla (w-z)=0$. On the connected triangle, a Sobolev function with zero weak gradient is constant, and its zero Sobolev trace on ${\Gamma }_D$ forces that constant to be zero. Hence $z=w$, provingThis proves the needed reaction-diffusion regularity from Laplacian regularity without assuming the dual solution already has second derivatives in the L2 space.

Only $H^2$ regularity is available for $z$. Use its continuous piecewise-linear interpolant $L_{h}z$, which belongs to the quadratic conforming finite element space $V_h$ because ${\mathcal{P}}_1(K)\subset {\mathcal{P}}_2(K)$. It has the same zero Dirichlet boundary condition and satisfies the finite element interpolation estimateThus no $H^3$ assumption on the dual solution is necessary. Taking $v=e$ in the dual weak formulation and using Galerkin orthogonality givesIf $\Vert e{\Vert }_0=0$, the desired bound holds. Otherwise divide by $\Vert e{\Vert }_0$ and use part (d):All constants depend only on the fixed triangle and shape regularity. The acute-angle assumption enters through the supplied mixed-boundary elliptic regularity result.

Multiply by a test function vanishing at the endpoints and apply integration by parts. The positive operator is $-u^{\prime \prime }+xu$, so the forcing changes sign. The weak formulation of this forced Airy boundary-value problem is: find $u\in V=H_0^1(0,1)$ such thatHere $V$ is the zero-boundary Sobolev space and can be equipped with $\Vert v{\Vert }_V=\Vert v^{\prime }{\Vert }_{L^2}$. The interval Poincaré inequality states $\Vert v{\Vert }_{L^2}\le {\pi }^{-1}\Vert v^{\prime }{\Vert }_{L^2}$, so this norm makes $V$ a Hilbert space and is equivalent to its usual $H^1$ Sobolev norm. The form is symmetric and satisfiesThe Lax-Milgram theorem says that a bounded coercive bilinear form on a Hilbert space and a bounded linear functional determine exactly one weak solution. All its hypotheses have just been verified. Moreover, $u^{\prime \prime }=xu+1\in L^2$ gives $u\in H^2$ on the interval; the right-hand side is then continuous, so this weak solution is the classical solution with the prescribed endpoint values.

Since the form is symmetric, the equivalent Ritz method minimizesIndeed, stationarity gives the weak equation andwith equality only at $u$. This proves the minimization equivalence, including the sign of the linear term.

Restricting to a conforming finite element space $V_h\subset V$ gives a unique discrete solution and Galerkin orthogonality, $a(u-u_h,v_h)=0$ for every $v_h\in V_h$. The Céa lemma states that if the continuity and coercivity constants are $M$ and $\alpha$, then $\Vert u-u_h{\Vert }_V\le (M/\alpha ){\inf }_{v_h\in V_h}\Vert u-v_h{\Vert }_V$. For this symmetric problem the exact best-approximation constant in the energy norm is one. The piecewise-linear finite element interpolation estimate $\Vert u-I_{h}u{\Vert }_{H^1}\le Ch\Vert u{\Vert }_{H^2}$ on regular finite element meshes therefore gives convergence of a numerical method of the Ritz method, with first-order error in the energy norm.

The finite element method replaces an infinite-dimensional variational problem by a finite-dimensional one, usually using functions that are polynomial on small mesh elements. Ritz method and Galerkin method describe how the discrete equations are selected; neither intrinsically requires piecewise polynomials, though those spaces make assembly local and sparse.

For a concrete realization of the two-point problem, choose homogeneous Dirichlet boundary conditions at both endpoints of $(0,1)$. The differential expression alone is not a complete boundary value problem, so this boundary choice is an explicit assumption. Let $p\in W^{1,\infty }$ with $p\ge p_0>0$, $q\in L^{\infty }$ with $q\ge 0$, and $f\in L^2$. Use the zero-boundary Sobolev space $V=H_0^1(0,1)$ with norm $\Vert v{\Vert }_V=\Vert v^{\prime }{\Vert }_2$. By the Poincaré inequality this is equivalent to its full first-order Sobolev norm. Multiplying by a test function and applying integration by parts givesThe endpoint terms vanish by the essential boundary condition. The form is a symmetric bilinear form, bounded byand it is a coercive bilinear form:The load is a bounded linear functional. The Lax-Milgram theorem gives a unique weak solution and a bound $\Vert u{\Vert }_V\le \Vert \ell {\Vert }_{V^{\ast }}/p_0$. Under the stated coefficient regularity it solves the differential equation in the usual weak sense.

The Ritz method minimizes the energyover $V$, or over a conforming finite-dimensional subspace $V_h\subset V$. Its first variation is $a(v,w)-\ell (w)$; symmetry and coercivity make the functional strictly convex. In fact, if $u$ is the weak solution,Thus the energy minimizer is exactly the weak solution. The Galerkin method instead directly asks for $u_h\in V_h$ with $a(u_h,v_h)=\ell (v_h)$ for all $v_h\in V_h$. For this symmetric coercive problem, Ritz-Galerkin equivalence for a symmetric coercive form shows that the discrete minimizer and Galerkin solution coincide. For a nonsymmetric problem, a Galerkin formulation can still apply but the same quadratic minimization generally cannot represent its full bilinear form; appropriate coercivity or inf-sup conditions must then justify the chosen spaces.

Subtract the continuous and discrete weak equations to obtain Galerkin orthogonality $a(u-u_h,v_h)=0$. Since $a$ defines the energy norm $\Vert v{\Vert }_a=\sqrt{a(v,v)}$, the Pythagorean identity givesIn a general reference norm, the Céa lemma yieldsThe proof uses $a(e,e)=a(e,u-v_h)$ together with coercivity and continuity. Dense approximation spaces therefore imply convergence; the estimate reduces a numerical error problem to an approximation problem.

For implementation, partition the interval by nodes $0=x_0<\cdots <x_N=1$. Choose continuous piecewise-linear functions and the interior piecewise-linear hat functions ${\varphi }_j$ as a basis, with the two boundary coefficients prescribed. Write $u_h={\sum }_j{U}_j{\varphi }_j$. The coefficients satisfyThe stiffness matrix is symmetric positive definite: $U^{T}KU=a(u_h,u_h)>0$ for nonzero interior coefficients. The local supports make it tridiagonal in this one-dimensional linear-element case.

On an element of length $h_e$, for constant element coefficients $p_e,q_e$, the two-node matrix and constant-load vector areFor variable coefficients, integrate $p,q,f$ against the same shape functions rather than silently replacing them by constants. Element contributions are added at shared nodes, and boundary values are eliminated or lifted into the right side. A sparse direct solve or a suitable iterative method gives the nodal coefficients. Numerical quadrature should preserve the required accuracy and, with positive coefficients and weights, the positive energy structure.

A finite element interpolation estimate gives ${\inf }_{v_h}\Vert u-v_h{\Vert }_{H^1}\le Ch\Vert u{\Vert }_{H^2}$ for these elements, where $h$ is the largest element size and the solution has the stated regularity. The Céa lemma then gives first-order energy error. If the dual elliptic problem has $H^2$ regularity, the Aubin–Nitsche duality argument gives second-order $L^2$ error: for $e=u-u_h$, solve $a(v,z)=(e,v{)}_{L^2}$, then use $a(e,z)=a(e,z-I_{h}z)$ to gain one further factor of $h$. Higher-order elements improve rates when the solution is sufficiently smooth.

Nonzero Dirichlet data are handled by a boundary lifting and an affine trial space. Neumann boundary conditions enter naturally through the integrated boundary term; Robin terms modify both the form and the load. With pure Neumann conditions and $q=0$, constants lie in the kernel: existence needs the load/flux compatibility condition, uniqueness needs a mean constraint, and the preceding Dirichlet coercivity argument cannot simply be reused. Conforming approximation, boundary conditions and coercivity are the ingredients connecting the finite-element construction to a justified Ritz or Galerkin solution.

Put $h=1/(N+1)$ and $x_i=ih$. The interior piecewise-linear hat functions areThey are a basis of the continuous piecewise-affine functions vanishing at both endpoints. The support of a function for each hat only overlaps those of neighbouring hats, so the stiffness matrix is a tridiagonal matrix. Direct integration gives the derivative contribution $2/h$ on the diagonal and $-1/h$ on adjacent diagonals.

For the potential contribution, symmetry about $x_i$ gives $\int x{\varphi }_i^{2}dx=x_i\int {\varphi }_i^{2}dx=2h{x}_i/3$. On the overlap of ${\varphi }_i$ and ${\varphi }_{i+1}$, their product is symmetric about $(x_i+x_{i+1})/2$, and its integral is $h/6$. This is the affine-weighted hat mass matrix calculation. Consequently,and all other entries vanish. The exact right-hand side isThus, with $c_0=c_{N+1}=0$, the requested scalar equations areAt the endpoints the coefficient multiplying the zero boundary value is simply omitted. Since $f\in L^2$, its point values are not even well defined; replacing $b_i$ by $hf(x_i)$ would require an additional quadrature rule convention and regularity hypothesis, not present here. The exact Ritz method uses the displayed integrals. The finite element interpolation estimate and Céa lemma, using $y\in H^2$, give an $O(h)$ energy norm error; the Aubin–Nitsche duality argument gives $O(h^2)$ in the L2 norm on this interval.

The starting point of a conforming finite element method is a weak formulation in a Hilbert space $X$. For example, a homogeneous Dirichlet boundary condition on a bounded Lipschitz domain gives $X=H_0^1(\Omega )$. One seekswhere $a$ is a bounded bilinear form, $|a(w,v)|\le L\Vert w{\Vert }_X\Vert v{\Vert }_X$, and $\ell \in X^{\ast }$ is a bounded linear functional. A conforming finite element space $X_h\subset X$ is finite dimensional and is typically built from functions that are polynomial on each element of a finite element mesh. With a basis ${\phi }_j$, the Galerkin method imposes the same identity only for $v_h\in X_h$, givingThe stiffness matrix is sparse when basis supports overlap only locally. Local element integrals assemble its entries; the mass matrix represents an $L^2$ term. This separates the choice of trial functions from the variational principle determining their coefficients.

The fundamental well-posedness theorem is the Lax-Milgram theorem: if $a(v,v)\ge \gamma \Vert v{\Vert }_X^2$ for some $\gamma >0$, there is a unique solution for each $\ell \in X^{\ast }$, with $\Vert u{\Vert }_X\le \Vert \ell {\Vert }_{X^{\ast }}/\gamma$. The same theorem on $X_h$ gives a unique discrete solution. The constants must be independent of the mesh for uniform stability of a numerical method. Boundary conditions are built into the trial space when they are essential constraints; flux conditions arise through boundary terms in integration by parts and are natural constraints in the weak formulation.

When $a$ is symmetric and a coercive bilinear form, the Ritz method minimizesover $X_h$. Its first variation is exactly the Galerkin method, so Ritz minimization and Galerkin testing coincide for symmetric coercive problems. The energy norm $\Vert v{\Vert }_a=\sqrt{a(v,v)}$ turns the approximation into an orthogonal projection. Indeed, subtracting the exact and discrete weak equations gives Galerkin orthogonality, $a(u-u_h,v_h)=0$. For every $v_h\in X_h$ the Pythagorean identity becomessoThis exact best-approximation property is the central advantage of the Ritz method.

Symmetry is unnecessary for the Galerkin method. For any bounded, coercive form, the Céa lemma givesTo prove it, put $e=u-u_h$. Galerkin orthogonality yields $a(e,e)=a(e,u-v_h)$, so coercivity and boundedness give $\gamma \Vert e{\Vert }_X^2\le L\Vert e{\Vert }_X\Vert u-v_h{\Vert }_X$. Divide when $e\ne 0$ and take the infimum. Thus approximation capability plus uniform coercivity gives convergence. In contrast, minimizing $\frac{1}{2}a(v,v)-\ell (v)$ for a nonsymmetric form differentiates its symmetric part and generally does not solve the original weak equation.

A concrete symmetric example is the one-dimensional reaction-diffusion problem from Question 5. Continuous piecewise-linear hat functions give the diffusion stiffness matrix $S$ with diagonal $2/h$ and neighbouring entries $-1/h$, and the mass matrix $M$ with diagonal $2h/3$ and neighbouring entries $h/6$. Therefore $K=S+M$ and $F_i=i{h}^2$. Its positive-definite matrix is the coordinate form of the continuous energy. The exact smooth solution $u(x)=x-\sinh (x)/\sinh 1$ makes this an explicit test of the method, not just an abstract existence result.

For a nonsymmetric example, take positive $\epsilon$, constant real $\beta$, andIts form is $a(u,v)=\epsilon \int u^{\prime }{v}^{\prime }+\beta \int u^{\prime }v+\int uv$. Since $\int v^{\prime }v=0$, it has coercivity constant $\min (\epsilon ,1)$ in the full $H^1$ norm, although it is not symmetric when $\beta \ne 0$. The Galerkin method remains well posed. On the uniform hat basis, $K=\epsilon S+\beta C+M$, where $C_{i,i+1}=1/2$, $C_{i,i-1}=-1/2$ and $C_{ii}=0$. The skew-symmetric matrix $C$ cancels from $c^{T}Kc$, proving nonsingularity. A Ritz energy would omit precisely this advection contribution. Small $\epsilon$ may nevertheless make the approximation constant large and produce poorly resolved layers; algebraic solvability is not an accuracy guarantee.

Quantitative convergence uses a finite element interpolation estimate. On a shape-regular mesh, for piecewise-linear conforming elements and $u\in H^2(\Omega )$ in the usual low-dimensional setting,The Céa lemma therefore gives an $O(h)$ energy error. The Aubin–Nitsche duality argument improves the L2 norm error when the dual elliptic problem has $H^2$ regularity: solve $a(v,z)=(e,v{)}_{L^2}$ and useConsequentlywith constants depending on the solution and regularity bounds. The second estimate is conditional on dual regularity; corners or insufficient data regularity can reduce the rate.

The same Galerkin method handles evolution equations. For the homogeneous heat equation, the semidiscrete finite element heat equation is $M\dot{c}+Sc=0$. Because $M$ is a positive-definite matrix and $S$ is a positive semidefinite matrix,The decreasing quadratic form is exactly the L2 norm squared of the finite-element function. This separates spatial stability of a numerical method from the subsequent time integrator, just as in the Schrödinger equation example. Conformity, coercivity and approximation estimates together explain convergence.

Consider a bounded Lipschitz domain $\Omega$ and homogeneous Dirichlet boundary conditions forAssume $A$ is a bounded real symmetric matrix field with ${\xi }^{T}A(x)\xi \ge a_0|\xi {|}^2$ uniformly, $c\ge 0$ is bounded, and $f$ defines a bounded linear functional on $V=H_0^1(\Omega )$, the zero-boundary Sobolev space. The weak formulation isThe Cauchy-Schwarz inequality gives continuity $|a(u,v)|\le M\Vert u{\Vert }_V\Vert v{\Vert }_V$. Uniform ellipticity and the Poincaré inequality give a coercive bilinear form, $a(v,v)\ge \alpha \Vert v{\Vert }_V^2$. The Lax-Milgram theorem therefore supplies a unique weak solution, with $\Vert u{\Vert }_V\le \Vert \ell {\Vert }_{V^{\prime }}/\alpha$. Its proof represents $a$ by a bounded linear operator $T$ using the Riesz representation theorem; coercivity gives an injective operator with closed image, and a zero orthogonal complement makes that image all of $V$. This yields existence, uniqueness and the displayed estimate.

For a conforming finite element space $V_h\subseteq V$, the Galerkin method finds $u_h\in V_h$ such that $a(u_h,v_h)=\ell (v_h)$ for every $v_h\in V_h$. With a basis ${\varphi }_1,\dots ,{\varphi }_N$, this is the linear system $KU=F$, where the stiffness matrix is $K_{ij}=a({\varphi }_j,{\varphi }_i)$ and $F_i=\ell ({\varphi }_i)$. Its positive-definite matrix property proves unique discrete solvability. Local support gives a sparse matrix through elementwise assembly.

For the symmetric problem, the Ritz method minimizes $J(v)=a(v,v)/2-\ell (v)$ over $V_h$. Differentiating in every direction $v_h$ gives precisely the Galerkin method equations. Conversely, if $u_h$ solves them, $J(u_h+w_h)-J(u_h)=a(w_h,w_h)/2$, proving it is the unique minimum. The same argument over $V$ gives $J(v)-J(u)=\Vert v-u{\Vert }_a^2/2$, where $\Vert v{\Vert }_a=\sqrt{a(v,v)}$ is the energy norm. Thus the Ritz method chooses the best trial function in that norm.

Subtracting the continuous and discrete equations gives Galerkin orthogonality, $a(u-u_h,v_h)=0$. For arbitrary $v_h\in V_h$ and $e=u-u_h$, it follows that $a(e,e)=a(e,u-v_h)$. Coercivity and continuity prove the Céa lemma:In the symmetric energy norm, the Cauchy-Schwarz inequality improves the constant to one. Equivalently, Galerkin orthogonality gives the Pythagorean identity $\Vert u-v_h{\Vert }_a^2=\Vert u-u_h{\Vert }_a^2+\Vert u_h-v_h{\Vert }_a^2$. These estimates establish stability and show that approximation properties of the trial space determine convergence.

For continuous piecewise-linear elements on a shape-regular mesh in dimensions at most three, a finite element interpolation estimate gives $\Vert u-I_{h}u{\Vert }_{H^1}\le Ch\Vert u{\Vert }_{H^2}$ when $u\in H^2(\Omega )$. Combining it with the Céa lemma yields $O(h)$ energy norm error. More generally, degree-$p$ elements give $O(h^p)$ error in $H^1$ under the corresponding regularity and approximation assumptions. Mere mesh refinement cannot supply an order whose required solution regularity is absent.

An additional order in the L2 norm follows under a dual elliptic regularity assumption. Solve $a(v,z)=(e,v{)}_{L^2}$ and assume $\Vert z{\Vert }_{H^2}\le C\Vert e{\Vert }_{L^2}$. The Aubin–Nitsche duality argument uses Galerkin orthogonality to obtainTherefore $\Vert e{\Vert }_{L^2}\le Ch\Vert e{\Vert }_{H^1}$, giving $O(h^2)$ for the piecewise-linear case with $u\in H^2$. This improvement explicitly requires regularity of the dual problem.

As a concrete example, solve $-u^{\prime \prime }=1$ on $(0,1)$ with zero endpoint values. On a uniform mesh of spacing $h$, the piecewise-linear hat functions give $K_{jj}=2/h$, $K_{j,j\pm 1}=-1/h$ and $F_j=h$. The discrete equations are $2U_j-U_{j-1}-U_{j+1}=h^2$, whose solution is $U_j=x_j(1-x_j)/2$. The exact solution is $u(x)=x(1-x)/2$, so $u_h$ is its piecewise-linear nodal interpolant. On each interval $[x_j,x_{j+1}]$,Integrating the squared error and squared derivative error gives the explicit rateswhere the $H_0^1$ norm here is $\Vert v^{\prime }{\Vert }_{L^2}$.

For nonsymmetric coercive problems, the Galerkin method, Galerkin orthogonality and the Céa lemma still apply, but the symmetric quadratic minimization interpretation of the Ritz method generally does not. Nonconforming trial spaces, inexact integration and noncoercive equations need additional arguments beyond the conforming coercive theory proved here.

A family of simplicial meshes is shape-regular when the ratio of each element's diameter to its inscribed-ball radius stays uniformly bounded. This excludes degenerating shapes and makes finite element interpolation estimates uniform in the mesh size. The element sizes may still vary substantially across the domain.