13 


14 
\section{The Einstein Summation Convention} 
\section{The Einstein Summation Convention} 
15 


16 

The Einstein Summation Convention (ESC) is a notational convention that is prefered by the \ESCRIPT developers. It is a condensed and practical way to deal with multidimensional and convoluted PDEs. By suppressing the need to write out the many terms of each problem it is possible to increase efficiency and reduce the number of errors created through poor working. According to the convention, when an index variable appears twice in a single term, it implies that we are summing over all of its possible values. 
17 

So we have; 
18 

\begin{equation} 
19 

a_1\frac{\partial^2 f}{\partial x_1^2} + a_2\frac{\partial^2 f}{\partial x_2^2} = a_i\frac{\partial^2 f}{\partial x_i^2} 
20 

\end{equation} 
21 


22 

For a scalar function $f(x_1,x_2,..x_i)$ and a vector $\mathbf{u}(u_1,u_2,..u_i)$ with $u_i(x_1,x_2,..x_i)$, we have the following notation: 
23 

\begin{equation} 
24 

\mathbf{u}=\sum_{i}u_ie^i = u_ie^i 
25 

\end{equation} 
26 

\begin{equation} 
27 

\mathbf{grad}(f) = \mathbf{\nabla}(f) = \sum_{i}\frac{\partial f}{\partial x_i}e^i = (\partial_i f)e^i = f_{,i}e^i 
28 

\end{equation} 
29 

\begin{equation} 
30 

div(\mathbf{u}) = \mathbf{\nabla}.\mathbf{u} = \sum_{i}\frac{\partial u_i}{\partial x_i} = \partial_i u_i = u_{i,i} 
31 

\end{equation} 
32 

\begin{equation} 
33 

div(\mathbf{grad}(f)) = \nabla^2 f = \Delta f = \sum_{i}\frac{\partial^2 f}{\partial x_i^2} = f_{,ii} 
34 

\end{equation} 