# Setting inhomogeneous Dirichlet boundary conditionsΒΆ

A mesh stores boundary elements, which know the *bc* name given in the geometry.
The Dirichlet boundaries are given as a regular expression of these names to the finite element space:

```
V = FESpace(mesh,order=3,dirichlet="top|back")
u = GridFunction(V)
```

The BitArray of free (i.e. unconstrained) dof numbers can be obtained via

```
freedofs = V.FreeDofs()
print (freedofs)
```

Inhomogeneous Dirichlet values are set via

```
u.Set(x*y, definedon=mesh.Boundaries("top|back"))
```

This function performs an element-wise L2 projection combined with arithmetic averaging of coupling dofs.

As usual we define biform a and liform f. Here, the full Neumann matrix and non-modified right hand sides are stored.

Boundary constraints are treated by the preconditioner. For example, a Jacobi preconditioner created via

```
c = Preconditioner(a, "local")
```

inverts only the unconstrained diagonal elements, and fills the remaining values with zeros.

The boundary value solver keeps the Dirichlet-values unchanged, and solves only for the free values

```
BVP(bf=a,lf=f,gf=u, pre=c).Do()
```

A do-it-yoursolve version of homogenization is:

```
res = f.vec.CreateVector()
res.data = f.vec - a.mat * u.vec
u.vec.data += a.mat.Inverse(v.FreeDofs()) * res
```