better way to do things is therefore if we build the matrix $\hat A$
as the vector Laplace matrix $\hat A_{ij} = (\nabla {\mathbf v}_i,
\eta \nabla {\mathbf v}_j)$ and then apply the same boundary condition
-as we applied to <i>A</i>. If this is Dirichlet boundary conditions all
+as we applied to <i>A</i>. If this is a Dirichlet boundary conditions all
around the domain, the $\hat A$ will decouple to three diagonal blocks
as above, and if the boundary conditions are of the form $\mathbf u
\cdot \mathbf n = 0$ then this will introduce a coupling of degrees of
step-17 and step-18, we use interfaces to the <a
href="http://trilinos.sandia.gov">Trilinos</a> library (see the
deal.II README file for installation instructions) in this program. Trilinos
-is a very large collection of
+is a very large collection of
everything that has to do with linear and nonlinear algebra, as well as all
sorts of tools around that (and looks like it will grow in many other
directions in the future as well).