\section{Implementierung}
-
-
\subsection{Datenstruktur}
Für die Implementierung in \Matlab~und C++ wollen wollen wir eine einheitliche Datenstruktur einführen.
Die für die Triangulierung $\mathcal{T}_{\ell} = \{T_1\ldots T_M\}$ benötigen Knoten $\mathcal{K}_{\ell} = \{C_1\ldots C_N\}$ stellen wir in einer $ N \times 3$ Matrix dar. Dabei enthält die $j$-te Zeile die Koordinaten des Knoten $C_j$ im $\R^3$.
(Siehe Abb.:\ref{exmpl3:nei:part})
\begin{figure}[ht]
\caption{Nachbarschaftsrelationen Element 4 aus Abb.\ref{exmpl3}}
-\label{exmpl3:nei:part}
+% \label{exmpl3:nei:part}
\centering
\subfloat[Lage]{\includegraphics{fig/Net_Neigh}}
\subfloat[Nachbarn]{\input{fig/exmpl3_nei_part}}
Index & n1 & n2 & n3 & n4 & n5 & n6 & n7 & n8\\
4 & 9 & 7 & 12 & 0 & 0 & 0 & 14 & 0
\end{tabular}
-\label{exmpl3:nei}
+\label{exmpl3:nei:part}
function export_mesh(coo, ele, nei, f2s, file)
-plotShape(coo,ele,'db');view(2);
+plotShape(coo,ele,'db');
print('-r600','-depsc',['../doc/fig/' file '_ref.eps'])
+system(['epstopdf ../doc/fig/' file '_ref.eps'])
%% Koordinaten
fid = fopen(['../doc/fig/' file '_coo.tex'],'w');
steps = 30;
%Art der Berechnungen
-type = [1 3];
+type = [1];
%LShape adaptiv anisotrop
-compute('exmpl_2DQuad', steps, 0.7, type, 0.5, 0.5, 0, 'testAA_')
+compute('exmpl_2DQuad', steps, 0.7, type, 0.5, 0, 0, 'testAA_')
% compute('exmpl_2DQuad', steps, 0.7, type, 0.5, 0.5, 1, 'testAAvcon_')
% %LShape adaptiv isotrop