Ši column-by-column using the Olsen method, which is described below. The columns of the matrix Ši are used to define a projected coordinate system, and in this section vectors and matrices in this space are decorated with a grave accent (e.g. v). The Algorithm 1 Iterative diagonalization algorithm evaluate yi = ȞL ši through forward finite difference (see equation 12) 13: