]> git.leopard-lacewing.eu Git - bacc.git/commitdiff
[doc] Beispiele Eingebunden
authortreecity <treecity@26120e32-c555-405d-b3e1-1f783fb42516>
Tue, 13 Mar 2012 12:43:00 +0000 (12:43 +0000)
committertreecity <treecity@26120e32-c555-405d-b3e1-1f783fb42516>
Tue, 13 Mar 2012 12:43:00 +0000 (12:43 +0000)
git-svn-id: https://drops.fb12.tu-berlin.de/svn/bacc/trunk@102 26120e32-c555-405d-b3e1-1f783fb42516

25 files changed:
doc/doc1.pdf
doc/doc1.tex
doc/doc2.tex [new file with mode: 0644]
doc/fig/Net_Neigh.eps [new file with mode: 0644]
doc/fig/Net_Neigh.pdf [new file with mode: 0644]
doc/fig/Net_Neigh.svg [new file with mode: 0644]
doc/fig/exmpl1_coo.tex [new file with mode: 0644]
doc/fig/exmpl1_ele.tex [new file with mode: 0644]
doc/fig/exmpl1_f2s.tex [new file with mode: 0644]
doc/fig/exmpl1_nei.tex [new file with mode: 0644]
doc/fig/exmpl1_ref.eps [new file with mode: 0644]
doc/fig/exmpl1_ref.pdf [new file with mode: 0644]
doc/fig/exmpl2_coo.tex [new file with mode: 0644]
doc/fig/exmpl2_ele.tex [new file with mode: 0644]
doc/fig/exmpl2_f2s.tex [new file with mode: 0644]
doc/fig/exmpl2_nei.tex [new file with mode: 0644]
doc/fig/exmpl2_ref.eps [new file with mode: 0644]
doc/fig/exmpl2_ref.pdf [new file with mode: 0644]
doc/fig/exmpl3_coo.tex [new file with mode: 0644]
doc/fig/exmpl3_ele.tex [new file with mode: 0644]
doc/fig/exmpl3_f2s.tex [new file with mode: 0644]
doc/fig/exmpl3_nei.tex [new file with mode: 0644]
doc/fig/exmpl3_nei_part.tex [new file with mode: 0644]
doc/fig/exmpl3_ref.eps [new file with mode: 0644]
doc/fig/exmpl3_ref.pdf [new file with mode: 0644]

index 3aa4b9652e3a120cc21810b4f9a2308f82ef4db7..e3e58fc992eca5b30ec10541e95725766beaae6c 100644 (file)
Binary files a/doc/doc1.pdf and b/doc/doc1.pdf differ
index 7b9937bee226264c04e3b025f59ee4d2a9a3fe56..b8eb38734512d5eb014396a02442a408e2ad7196 100644 (file)
@@ -70,7 +70,7 @@ Berechnet werden soll:
 
 \section{Gitter}
 Alle Gitter bestehen nur aus achsenorientierten Rechtecken.
-Entscheidend für die Tests sind Gitter in einer Ebene und Ganze Körper...
+Entscheidend für die Tests sind Gitter in einer Ebene und Körper...
 
 \chapter{Analytische Berechnung der Integrale}
 
diff --git a/doc/doc2.tex b/doc/doc2.tex
new file mode 100644 (file)
index 0000000..0545be6
--- /dev/null
@@ -0,0 +1,167 @@
+\documentclass[a4paper,10pt,fleqn]{report}
+\usepackage{fullpage}  %Seiten etwas Größer
+\usepackage{amsmath,amssymb}   %Mathematische Symbole
+%\usepackage{moreverb}
+\usepackage{graphicx,psfrag,subfig}    %Grafiken einbinden/Texte ersetzen/Bilder nebeneinander
+%\usepackage{ifthen}
+%\usepackage{showkeys}
+%\usepackage[final,numbered]{mcode}
+\usepackage{colortbl}  %Einfache Färbungen in Tabellen
+\usepackage{color}     %Farben benutzen und Definieren
+\usepackage{subfig}    %mehrere Figuren in einer
+\usepackage{hyperref}  %Links im Inhaltsverzeichnis
+\hypersetup{linkbordercolor={1 1 1},citebordercolor={1 1 1},urlbordercolor={1 1 1}}
+
+
+\usepackage[ngerman]{babel}    %Sprachpacket für Überschriften
+\usepackage[utf8]{inputenc}    %Eingabekodierung
+\usepackage{fixltx2e}  %Deutschsprach Bugs
+
+\definecolor{gray}{gray}{.95}
+
+\def\todo#1{\textcolor{red}{#1}}
+\def\why#1{\textcolor{blue}{#1}}
+\def\Matlab{{\sc Matlab}}
+\def\q{\Q}
+
+\def\Ea{$E$}
+\def\Eb{$\tilde E$}
+\def\v{\boldsymbol{v}}
+\def\D{\mathcal{D}}
+\def\Q{\mathcal{Q}}
+\def\G{\mathcal{G}}
+\def\L{\mathcal{L}}
+\def\T{\mathcal{T}}
+\def\oder{\vee}
+\def\und{\wedge}
+
+\newcommand{\dif}[2]{\frac{\partial#1}{\partial#2}}
+\newcommand{\enorm}[1]{\left| \! \left| \! \left|#1\right| \! \right| \! \right|}
+\newcommand{\norm}[1]{\lVert#1\rVert}
+\newcommand{\abs}[1]{\lvert#1\rvert}
+\newcommand{\showMesh}[2][]{\begin{figure}[ht]
+\caption{#1}
+\label{#2}
+\centering
+\subfloat[Lage]{\includegraphics[width=0.5\textwidth]{fig/#2_ref}}
+\subfloat[Koordinaten]{\input{fig/#2_coo}}\\
+\subfloat[Elemente]{\input{fig/#2_ele}}
+\subfloat[Nachbarn]{\input{fig/#2_nei}}
+\end{figure}}
+
+\def\N{\mathbb{N}}
+\def\R{\mathbb{R}}
+
+\author{P. Schaefer}
+
+\begin{document}
+\tableofcontents
+\clearpage
+\chapter{Implementierung}
+
+
+\section{Datenstruktur}
+Alle Koordinaten werden in einer $ N \times 3,N \in \N$ Matrix abgespeichert.
+\begin{displaymath}
+ COO[i,1:3] = C_i := (x_1,x_2,x_3)^{-1} \text{ wobei } x_1,x_2,x_3 \in \R \und i \in \{1,2 \dots N\}
+\end{displaymath}
+Jedes Element wird dann durch jeweils vier Koordinatenindizes als $M \times 4, M\in\N$ Matrix gespeichert. Wobei die Koordinaten links herum durchgegangen werden.
+\begin{displaymath}
+ ELE[i,1:4] = E_i := (v_1,v_2,v_3,v_4) \text{ wobei } v_1,v_2,v_3,v_4 \in C \und i \in \{1,2 \dots M\}
+\end{displaymath}
+Aufgrund der Netzstabilität wollen wir maximal zwei Nachbarn an einer Kante zulassen, weshalb wir nur $2*4$ Einträge für die Nachbarschaftsrelationen pro Element benötigen. In der $M \times 8$ Matrix sind die Indizes der Nachbarn $n_1,n_2 \in E$ der Kante $(v_k, v_{k^{+1}}) | k\in \{1,2,3,4\}$ des Elements $E_i$ an den Stellen $NEI[i,[k,k+4]]$ gespeichert. Sollte eine Seite keinen Nachbarn haben markieren wir diese mit einer $0$. Außerdem wird an Kanten mit nur einem Nachbarn der erste Index auf den Nachbarn gesetzt und der zweite auf $0$.(Siehe Figur:\ref{exmpl3:nei:part})
+\\\noindent
+Ein ausführliches Beispiel ist in Figure \ref{exmpl3} dargestellt.
+\showMesh[Beispiel 3]{exmpl3}
+\begin{figure}[ht]
+\caption{Nachbarschaftsrelationen Element 4 aus \ref{exmpl3}}
+\label{exmpl3:nei:part}
+\centering
+ \subfloat[Lage]{\includegraphics{fig/Net_Neigh}}
+ \subfloat[Nachbarn]{\input{fig/exmpl3_nei_part}}
+\end{figure}
+
+
+
+
+
+\section{Verfeinern}
+Da wir im weiteren Verlauf sowohl adaptive also auch anisotrope Netzverfeinerung zulassen wollen, ist es sinnvoll eine Verfeinerungsfunktion zu implementieren, die alle Möglichen Teilungsprozesse auf einem Element unterstützt. Dabei sind vier nur wirklich relevant:
+\begin{enumerate}
+ \item keine Teilung
+ \item volle Teilung in vier gleich große Elemente
+ \item halbe Teilung in zwei gleichgroße horizontal liegende Elemente
+ \item halbe Teilung in zwei gleichgroße vertikal liegende Elemente
+\end{enumerate}
+Zusätzlich wurde auch Typ 5. belegt, welcher als Ergebnis eine volle Teilung 2. ausführt, diese aber schrittweise durch eine 3. Teilung und zwei 4. Teilungen.
+
+Damit jedem Element $E_i$ eine Teilungsart zugeordnet werden kann, haben wir einen Markierungsvektor $marked \in {1,2,3,4,5}^M$ eingeführt. Dabei entspricht $marked_i$ der Art der Teilung für das Element $E_i$. Um isotrope und auch uniforme Teilungen zu erleichtern kann statt dem Vektor $marked$ auch nur ein Skalar übergeben werden $marked \in {1,2\dots5}$, wodurch jedes Element mit der gewählten Art verfeinert wird.
+
+Relevant zum Verfeinern eines Netzes sind also die Koordinaten $COO$, Elemente $ELE$ sowie die Nachbarschaftsrelationen $NEI$ und der Markierungsvektor $marked$.
+
+Da wir später einen Fehlerschätzer berechnen wollen, ist es wichtig sich zu jedem Element seine Teilelemente zu merken. Dazu legen wir während der Teilung eine $M \times 4$ Matrix an, in der die maximal vier Elementindizes gespeichert sind. Wenn wir also ein Element in vier gleich große Teile verfeinern, so wird das neue Element links unten das erste sein und alle weiteren folgen im Kreis links herum. Teilen wir ein Element in zwei gleich große Elemente, so werden die doppelt belegten Quadranten auch doppelt eingetragen. Ein gar nicht geteiltes Element wird also vier mal den alten Indizes speichern. Dadurch wird sicher gestellt, dass das arithmetische Mittel über die Elemente immer gültig auszuführen ist.
+(Siehe Figur:\ref{exmpl3:f2s})
+
+$[COO_{fine}, ELE_{fine}, NEI_{fine}, F2S ] = refineQuad(COO, ELE, NEI, marked);$
+
+\begin{figure}[ht]
+\caption{VaterSohn}
+\centering
+\subfloat[VaterSohn aus \ref{exmpl2}]{\input{fig/exmpl2_f2s}}
+\subfloat[VaterSohn aus \ref{exmpl3}]{\input{fig/exmpl3_f2s}}
+\end{figure}
+
+\section{Fehlerschätzer}
+\begin{eqnarray*}
+\mu_h^2 & = & \norm{\varrho_h^{1/2}(\phi_{\frac h 2}-\Pi_h\phi_{\frac h 2})}_{\L^2(\T)}^2 \\
+& = &\sum_{T\in \T_h}\mu_h(T)^2\\
+\mu_h(T)^2 & = & \norm{\varrho_h^{1/2}(\phi_{\frac h 2}-\Pi_h\phi_{\frac h 2})}_{\L^2(\T)}^2 \\
+& = & h_{min}(T)\norm{\phi_{\frac h 2}-\Pi_h\phi_{\frac h 2}}_{\L^2(\T)}^2 \\
+&& T_j \in \tau_h, t_j^{(1)},\dots,t_j^{(4)} \in \tau_{\frac h 2} \\
+\phi_{\frac h 2}|_{T_j} & = &x_j^{(1)}\chi_{t_j^{(1)}}+\dots+x_j^{(4)}\chi_{t_j^{(4)}}\\
+\Pi_h\phi_{\frac h 2}|_{T_j}&=&\frac 1 {\abs{T_j}}\int_{T_j}\phi_{\frac h 2}d\Gamma\\
+& =&\frac 1 {\abs{T_j}}\sum_{k=1}^4x_j^{(k)}\int_{t^{(k)}}1d\Gamma\\
+& =&\frac 1 {\abs{T_j}}\sum_{k=1}^4x_j^{(k)}\abs{t_j^{(k)}}\\
+& =&\frac 1 {\abs{T_j}}\sum_{k=1}^4x_j^{(k)}\frac {\abs{T_j}} 4 \\
+& =&\frac 1 4 (x_j^{(1)} +\dots+x_j^{(4)} )=: m_j\\
+\norm{\phi_{\frac h 2} - \Pi_h\phi_{\frac h 2}}_{\L^2(T_j)}^2
+&=&\sum_{k=1}^4\norm{\phi_{\frac h 2} - \Pi_h\phi_{\frac h 2}}_{\L^2\left(t_j^{(k)}\right)}^2\\
+&=&\sum_{k=1}^4 \int_{t_j^{(k)}}(\phi_{\frac h 2} - \Pi_h\phi_{\frac h 2})^2d\Gamma\\
+&=&\sum_{k=1}^4 \int_{t_j^{(k)}}(x_j^{(k)}-m_j)^2d\Gamma\\
+&=&\sum_{k=1}^4 \abs{t_j^{(k)}}(x_j^{(k)}-m_j)^2\\
+&=&\frac {\abs{T_j}} 4\sum_{k=1}^4 (x_j^{(k)}-m_j)^2\\
+\mu_e(T_j)^2 & = & \frac{\why{ h_{min}(T_j)} \abs{T_j}}{4} \sum_{k=1}^4{(x_j^{(k)}-m_j)^2}
+\end{eqnarray*}
+
+$mu = computeEstSlpMuTilde(x_{fine}, COO, ELE, F2S);$
+
+\section{Markieren}
+Bestimme $M_e \subseteq T_e$ mit minimaler Kardinalität
+\begin{eqnarray*}
+\theta \sum_{T\in T_e} \mu_e(T)^2 & \leq & \sum_{T\in M_e} \mu_e(T)^2
+\end{eqnarray*}
+Zur anisotropen Verfeinerung wird weiterhin berechnet:
+\begin{eqnarray*}
+\begin{pmatrix}
+ C_j^{(1)}\\ C_j^{(2)}\\C_j^{(3)}\\ C_j^{(4)}
+\end{pmatrix}
+& = & \frac 1 4
+\begin{pmatrix}
+1 & 1 & 1 & 1\\
+1 & -1 & 1 & -1\\
+1 & 1 & -1 & -1\\
+1 & -1 & -1 & 1\\
+\end{pmatrix}
+\begin{pmatrix}
+ x_j^{(1)}\\ x_j^{(2)}\\x_j^{(3)}\\ x_j^{(4)}
+\end{pmatrix}
+\end{eqnarray*}
+
+\begin{eqnarray*}
+\eta \abs{ C_j^{(3)}} & \geq & \sqrt{\abs{ C_j^{(2)}}^2 + \abs{ C_j^{(4)}}^2}\\
+\eta \abs{ C_j^{(4)}} & \geq & \sqrt{\abs{ C_j^{(2)}}^2 + \abs{ C_j^{(3)}}^2}
+\end{eqnarray*}
+
+$xF2S := x_{fine}[F2S]$\\
+$marked = computeEstSlpMuTilde(xF2S, mu, theta, eta);$
+\end{document}
diff --git a/doc/fig/Net_Neigh.eps b/doc/fig/Net_Neigh.eps
new file mode 100644 (file)
index 0000000..0c21686
--- /dev/null
@@ -0,0 +1,272 @@
+%!PS-Adobe-3.0 EPSF-3.0
+%%Creator: cairo 1.10.2 (http://cairographics.org)
+%%CreationDate: Sun Mar  4 23:54:26 2012
+%%Pages: 1
+%%BoundingBox: 31 4 156 156
+%%DocumentData: Clean7Bit
+%%LanguageLevel: 2
+%%EndComments
+%%BeginProlog
+/cairo_eps_state save def
+/dict_count countdictstack def
+/op_count count 1 sub def
+userdict begin
+/q { gsave } bind def
+/Q { grestore } bind def
+/cm { 6 array astore concat } bind def
+/w { setlinewidth } bind def
+/J { setlinecap } bind def
+/j { setlinejoin } bind def
+/M { setmiterlimit } bind def
+/d { setdash } bind def
+/m { moveto } bind def
+/l { lineto } bind def
+/c { curveto } bind def
+/h { closepath } bind def
+/re { exch dup neg 3 1 roll 5 3 roll moveto 0 rlineto
+      0 exch rlineto 0 rlineto closepath } bind def
+/S { stroke } bind def
+/f { fill } bind def
+/f* { eofill } bind def
+/n { newpath } bind def
+/W { clip } bind def
+/W* { eoclip } bind def
+/BT { } bind def
+/ET { } bind def
+/pdfmark where { pop globaldict /?pdfmark /exec load put }
+    { globaldict begin /?pdfmark /pop load def /pdfmark
+    /cleartomark load def end } ifelse
+/BDC { mark 3 1 roll /BDC pdfmark } bind def
+/EMC { mark /EMC pdfmark } bind def
+/cairo_store_point { /cairo_point_y exch def /cairo_point_x exch def } def
+/Tj { show currentpoint cairo_store_point } bind def
+/TJ {
+  {
+    dup
+    type /stringtype eq
+    { show } { -0.001 mul 0 cairo_font_matrix dtransform rmoveto } ifelse
+  } forall
+  currentpoint cairo_store_point
+} bind def
+/cairo_selectfont { cairo_font_matrix aload pop pop pop 0 0 6 array astore
+    cairo_font exch selectfont cairo_point_x cairo_point_y moveto } bind def
+/Tf { pop /cairo_font exch def /cairo_font_matrix where
+      { pop cairo_selectfont } if } bind def
+/Td { matrix translate cairo_font_matrix matrix concatmatrix dup
+      /cairo_font_matrix exch def dup 4 get exch 5 get cairo_store_point
+      /cairo_font where { pop cairo_selectfont } if } bind def
+/Tm { 2 copy 8 2 roll 6 array astore /cairo_font_matrix exch def
+      cairo_store_point /cairo_font where { pop cairo_selectfont } if } bind def
+/g { setgray } bind def
+/rg { setrgbcolor } bind def
+/d1 { setcachedevice } bind def
+%%EndProlog
+11 dict begin
+/FontType 42 def
+/FontName /DejaVuSans def
+/PaintType 0 def
+/FontMatrix [ 1 0 0 1 0 0 ] def
+/FontBBox [ 0 0 0 0 ] def
+/Encoding 256 array def
+0 1 255 { Encoding exch /.notdef put } for
+Encoding 1 /uni0031 put
+Encoding 2 /uni0032 put
+Encoding 3 /uni0033 put
+Encoding 4 /uni0034 put
+Encoding 5 /uni0035 put
+Encoding 6 /uni0036 put
+Encoding 7 /uni0037 put
+Encoding 8 /uni0038 put
+/CharStrings 9 dict dup begin
+/.notdef 0 def
+/uni0031 1 def
+/uni0032 2 def
+/uni0033 3 def
+/uni0034 4 def
+/uni0035 5 def
+/uni0036 6 def
+/uni0037 7 def
+/uni0038 8 def
+end readonly def
+/sfnts [
+<00010000000a008000030020636d6170002bf0a7000006f8000000506376742000691d390000
+0748000001fe6670676d7134766a00000948000000ab676c79663a9b4078000000ac0000064c
+68656164f959a45f000009f400000036686865610cb8065a00000a2c00000024686d74782d85
+053d00000a50000000246c6f636100001d5800000a74000000286d6178700476067100000a9c
+00000020707265703b07f10000000abc0000056800020066fe96046605a400030007001a400c
+04fb0006fb0108057f0204002fc4d4ec310010d4ecd4ec301311211125211121660400fc7303
+1bfce5fe96070ef8f2720629000100e10000045a05d5000a004040154203a00402a005810700
+a009081f061c03001f010b10d44bb00f5458b9000100403859ecc4fcec31002fec32f4ecd4ec
+304b5358592201b40f030f04025d3721110535253311211521fe014afe990165ca014afca4aa
+047348b848fad5aa0000000100960000044a05f0001c009a4027191a1b03181c110504001105
+05044210a111940da014910400a00200100a02010a1c171003061d10fc4bb015544bb016545b
+4bb014545b58b90003ffc03859c4d4ecc0c011123931002fec32f4ecf4ec304b5358071005ed
+0705ed11173959220140325504560556077a047a05761b87190704000419041a041b051c7400
+7606751a731b741c82008619821a821b821ca800a81b115d005d25211521353600373e013534
+2623220607353e01333204151406070600018902c1fc4c73018d33614da7865fd3787ad458e8
+0114455b19fef4aaaaaa7701913a6d974977964243cc3132e8c25ca5701dfeeb00000001009c
+ffe3047305f000280070402e0015130a86091f862013a0150da00993061ca020932391068c15
+a329161c13000314191c2620101c03141f09062910fc4bb016544bb014545b58b90009ffc038
+59c4c4d4ecf4ec11173939310010ece4f4e4ec10e6ee10ee10ee10ee11123930014009641e61
+1f6120642104005d011e0115140421222627351e013332363534262b01353332363534262322
+0607353e01333204151406033f91a3fed0fee85ec76a54c86dbec7b9a5aeb6959ea39853be72
+73c959e6010c8e03251fc490ddf22525c33132968f8495a67770737b2426b42020d1b27cab00
+00020064000004a405d50002000d0081401d010d030d0003030d4200030b07a0050103810901
+0c0a001c0608040c0e10dc4bb00b544bb00d545b58b9000cffc03859d43cc4ec32113931002f
+e4d43cec321239304b5358071004c9071005c9592201402a0b002a0048005900690077008a00
+0716012b0026012b0336014e014f0c4f0d5601660175017a0385010d5d005d09012103331133
+152311231121350306fe0201fe35fed5d5c9fd5e0525fce303cdfc33a8fea00160c300000001
+009effe3046405d5001d005e4023041a071186101d1aa00714a010890d02a000810d8c07a41e
+171c010a031c000a10061e10fc014bb016544bb014545b58b90010ffc038594bb00f5458b900
+1000403859c4d4ec10c4ee310010e4e4f4ec10e6ee10fec410ee1112393013211521113e0133
+320015140021222627351e0133323635342623220607dd0319fda02c582cfa0124fed4feef5e
+c3685ac06badcacaad51a15405d5aafe920f0ffeeeeaf1fef52020cb3130b69c9cb624260000
+0002008fffe3049605f0000b0024005840241306000d860c00a01606a01c16a510a00c892291
+1c8c250c22091c191e131c03211f1b2510fcececf4ece4310010e4f4e4fce410ee10ee10ee11
+1239304014cb00cb01cd02cd03cd04cb05cb0607a41eb21e025d015d01220615141633323635
+342601152e01232202033e0133320015140023200011100021321602a4889f9f88889f9f0109
+4c9b4cc8d30f3bb26be10105fef0e2fefdfeee0150011b4c9b033bbaa2a1bbbba1a2ba0279b8
+2426fef2feef575dfeefebe6feea018d0179016201a51e000000000100a80000046805d50006
+00634018051102030203110405044205a0008103050301040100060710fcccc411393931002f
+f4ec304b5358071005ed071005ed5922014bb0165458bd00070040000100070007ffc0381137
+3859401258020106031a05390548056703b000b006075d005d13211501230121a803c0fde2d3
+01fefd3305d556fa81052b0000000003008bffe3048b05f0000b0023002f00434025180c00a0
+2706a01e2da012911e8c27a330180c242a1c15241c0f091c151b1e031c0f211b3010fcc4ecf4
+c4ec10ee10ee113939310010ece4f4ec10ee10ee393930012206151416333236353426252e01
+353424333216151406071e01151404232224353436131416333236353426232206028b90a5a5
+9090a6a5fea5829100ffdedffe918192a3fef7f7f7fef7a448918382939382839102c59a8787
+9a9b86879a5620b280b3d0d0b380b22022c68fd9e8e8d98fc601617482827474828200000000
+0002000300000000001400010000000000340004002000000004000400010000f008ffff0000
+f000ffff10000001000000000006001c00000000000900000001000200030004000500060007
+0008013500b800cb00cb00c100aa009c01a600b800660000007100cb00a002b20085007500b8
+00c301cb0189022d00cb00a600f000d300aa008700cb03aa0400014a003300cb000000d90502
+00f4015400b4009c01390114013907060400044e04b4045204b804e704cd0037047304cd0460
+0473013303a2055605a60556053903c5021200c9001f00b801df007300ba03e9033303bc0444
+040e00df03cd03aa00e503aa0404000000cb008f00a4007b00b80014016f007f027b0252008f
+00c705cd009a009a006f00cb00cd019e01d300f000ba018300d5009803040248009e01d500c1
+00cb00f600830354027f00000333026600d300c700a400cd008f009a0073040005d5010a00fe
+022b00a400b4009c00000062009c0000001d032d05d505d505d505f0007f007b005400a406b8
+0614072301d300b800cb00a601c301ec069300a000d3035c037103db0185042304a80448008f
+0139011401390360008f05d5019a0614072306660179046004600460047b009c000002770460
+01aa00e904600762007b00c5007f027b000000b4025205cd006600bc00660077061000cd013b
+01850389008f007b0000001d00cd074a042f009c009c0000077d006f0000006f0335006a006f
+007b00ae00b2002d0396008f027b00f600830354063705f6008f009c04e10266008f018d02f6
+00cd03440029006604ee00730000140000960000b707060504030201002c2010b002254964b0
+40515820c859212d2cb002254964b040515820c859212d2c20100720b00050b00d7920b8ffff
+5058041b0559b0051cb0032508b0042523e120b00050b00d7920b8ffff5058041b0559b0051c
+b0032508e12d2c4b505820b0fd454459212d2cb002254560442d2c4b5358b00225b002254544
+5921212d2c45442d2cb00225b0022549b00525b005254960b0206368208a108a233a8a10653a
+2d00000100000002547a5a7003545f0f3cf5001f080000000000ca6f827200000000ca6f8272
+f7d6fcae0d72095500000008000000010000000000010000076dfe1d00000de2f7d6fa510d72
+00010000000000000000000000000000000904cd0066051700e1051700960517009c05170064
+0517009e0517008f051700a80517008b0000000000000044000000b4000001b0000002980000
+035400000414000004ec000005780000064c0001000000090354002b0068000c000200100099
+000800000415021600080004b8028040fffbfe03fa1403f92503f83203f79603f60e03f5fe03
+f4fe03f32503f20e03f19603f02503ef8a4105effe03ee9603ed9603ecfa03ebfa03eafe03e9
+3a03e84203e7fe03e63203e5e45305e59603e48a4105e45303e3e22f05e3fa03e22f03e1fe03
+e0fe03df3203de1403dd9603dcfe03db1203da7d03d9bb03d8fe03d68a4105d67d03d5d44705
+d57d03d44703d3d21b05d3fe03d21b03d1fe03d0fe03cffe03cefe03cd9603cccb1e05ccfe03
+cb1e03ca3203c9fe03c6851105c61c03c51603c4fe03c3fe03c2fe03c1fe03c0fe03bffe03be
+fe03bdfe03bcfe03bbfe03ba1103b9862505b9fe03b8b7bb05b8fe03b7b65d05b7bb03b78004
+b6b52505b65d40ff03b64004b52503b4fe03b39603b2fe03b1fe03b0fe03affe03ae6403ad0e
+03acab2505ac6403abaa1205ab2503aa1203a98a4105a9fa03a8fe03a7fe03a6fe03a51203a4
+fe03a3a20e05a33203a20e03a16403a08a4105a096039ffe039e9d0c059efe039d0c039c9b19
+059c64039b9a10059b19039a1003990a0398fe0397960d0597fe03960d03958a410595960394
+930e05942803930e0392fa039190bb0591fe03908f5d0590bb039080048f8e25058f5d038f40
+048e25038dfe038c8b2e058cfe038b2e038a8625058a410389880b05891403880b0387862505
+8764038685110586250385110384fe038382110583fe0382110381fe0380fe037ffe0340ff7e
+7d7d057efe037d7d037c64037b5415057b25037afe0379fe03780e03770c03760a0375fe0374
+fa0373fa0372fa0371fa0370fe036ffe036efe036c21036bfe036a1142056a530369fe03687d
+036711420566fe0365fe0364fe0363fe0362fe03613a0360fa035e0c035dfe035bfe035afe03
+59580a0559fa03580a035716190557320356fe03555415055542035415035301100553180352
+1403514a130551fe03500b034ffe034e4d10054efe034d10034cfe034b4a13054bfe034a4910
+054a1303491d0d05491003480d0347fe0346960345960344fe0343022d0543fa0342bb03414b
+0340fe033ffe033e3d12053e14033d3c0f053d12033c3b0d053c40ff0f033b0d033afe0339fe
+033837140538fa033736100537140336350b05361003350b03341e03330d0332310b0532fe03
+310b03302f0b05300d032f0b032e2d09052e10032d09032c32032b2a25052b64032a2912052a
+25032912032827250528410327250326250b05260f03250b0324fe0323fe03220f0321011005
+2112032064031ffa031e1d0d051e64031d0d031c1142051cfe031bfa031a42031911420519fe
+031864031716190517fe031601100516190315fe0314fe0313fe031211420512fe0311022d05
+114203107d030f64030efe030d0c16050dfe030c0110050c16030bfe030a100309fe0308022d
+0508fe030714030664030401100504fe03401503022d0503fe0302011005022d0301100300fe
+0301b80164858d012b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b
+2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b
+2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b
+2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b
+2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b002b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b
+2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b
+2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b
+2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b
+2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b1d00>
+] def
+/f-0-0 currentdict end definefont pop
+%%Page: 1 1
+%%BeginPageSetup
+%%PageBoundingBox: 31 4 156 156
+%%EndPageSetup
+q 31 4 125 152 rectclip q
+0 g
+0.8 w
+0 J
+0 j
+[] 0.0 d
+4 M q 1 0 0 -1 0 160 cm
+32 29.145 96 94.855 re S Q
+q 1 0 0 -1 0 160 cm
+32 4.855 m 32.285 28.57 l 80 28.855 l 80 4.855 l 80 28.855 l 128.285 
+29.43 l 128.285 4.855 l S Q
+q 1 0 0 -1 0 160 cm
+128.57 29.145 m 155.715 28.855 l S Q
+q 1 0 0 -1 0 160 cm
+128.855 124 m 151.715 124 l S Q
+q 1 0 0 -1 0 160 cm
+32 124 m 32 155.43 l S Q
+q 1 0 0 -1 0 160 cm
+128 124 m 128 155.715 l S Q
+BT
+8 0 0 8 51.760217 39.0323 Tm
+/f-0-0 1 Tf
+<01>Tj
+ET
+BT
+8 0 0 8 120.208154 47.800415 Tm
+/f-0-0 1 Tf
+<02>Tj
+ET
+BT
+8 0 0 8 102.671899 119.925311 Tm
+/f-0-0 1 Tf
+<03>Tj
+ET
+BT
+8 0 0 8 34.789655 102.106219 Tm
+/f-0-0 1 Tf
+<04>Tj
+ET
+BT
+8 0 0 8 97.297894 39.0323 Tm
+/f-0-0 1 Tf
+<05>Tj
+ET
+BT
+8 0 0 8 120.77384 96.449365 Tm
+/f-0-0 1 Tf
+<06>Tj
+ET
+BT
+8 0 0 8 54.022961 121.905209 Tm
+/f-0-0 1 Tf
+<07>Tj
+ET
+BT
+8 0 0 8 35.638181 64.770984 Tm
+/f-0-0 1 Tf
+<08>Tj
+ET
+Q Q
+showpage
+%%Trailer
+count op_count sub {pop} repeat
+countdictstack dict_count sub {end} repeat
+cairo_eps_state restore
+%%EOF
diff --git a/doc/fig/Net_Neigh.pdf b/doc/fig/Net_Neigh.pdf
new file mode 100644 (file)
index 0000000..3de9aaa
Binary files /dev/null and b/doc/fig/Net_Neigh.pdf differ
diff --git a/doc/fig/Net_Neigh.svg b/doc/fig/Net_Neigh.svg
new file mode 100644 (file)
index 0000000..9bf6ac9
--- /dev/null
@@ -0,0 +1,378 @@
+<?xml version="1.0" encoding="UTF-8" standalone="no"?>
+<!-- Created with Inkscape (http://www.inkscape.org/) -->
+
+<svg
+   xmlns:dc="http://purl.org/dc/elements/1.1/"
+   xmlns:cc="http://creativecommons.org/ns#"
+   xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
+   xmlns:svg="http://www.w3.org/2000/svg"
+   xmlns="http://www.w3.org/2000/svg"
+   xmlns:sodipodi="http://sodipodi.sourceforge.net/DTD/sodipodi-0.dtd"
+   xmlns:inkscape="http://www.inkscape.org/namespaces/inkscape"
+   width="200"
+   height="200"
+   id="svg2"
+   version="1.1"
+   inkscape:version="0.48.2 r9819"
+   sodipodi:docname="Net_Split.svg">
+  <defs
+     id="defs4" />
+  <sodipodi:namedview
+     id="base"
+     pagecolor="#ffffff"
+     bordercolor="#666666"
+     borderopacity="1.0"
+     inkscape:pageopacity="0.0"
+     inkscape:pageshadow="2"
+     inkscape:zoom="2.8284271"
+     inkscape:cx="191.74258"
+     inkscape:cy="83.885664"
+     inkscape:document-units="px"
+     inkscape:current-layer="layer1"
+     showgrid="true"
+     showguides="true"
+     inkscape:guide-bbox="true"
+     inkscape:snap-global="true"
+     inkscape:window-width="1600"
+     inkscape:window-height="876"
+     inkscape:window-x="0"
+     inkscape:window-y="24"
+     inkscape:window-maximized="1"
+     inkscape:snap-grids="false"
+     inkscape:snap-to-guides="true"
+     inkscape:object-paths="false"
+     inkscape:object-nodes="true">
+    <inkscape:grid
+       type="xygrid"
+       id="grid3003"
+       empspacing="5"
+       visible="true"
+       enabled="true"
+       snapvisiblegridlinesonly="true" />
+    <sodipodi:guide
+       orientation="0,1"
+       position="100,100"
+       id="guide3015" />
+    <sodipodi:guide
+       orientation="1,0"
+       position="100,100"
+       id="guide3017" />
+    <sodipodi:guide
+       orientation="0,1"
+       position="80,80"
+       id="guide3021" />
+    <sodipodi:guide
+       orientation="0,1"
+       position="120,120"
+       id="guide3023" />
+    <sodipodi:guide
+       orientation="0,1"
+       position="100,155"
+       id="guide4138" />
+    <sodipodi:guide
+       orientation="0,1"
+       position="100,45"
+       id="guide4140" />
+    <sodipodi:guide
+       orientation="0,1"
+       position="100,175"
+       id="guide4142" />
+    <sodipodi:guide
+       orientation="0,1"
+       position="100,135"
+       id="guide4144" />
+    <sodipodi:guide
+       orientation="0,1"
+       position="100,65"
+       id="guide4146" />
+    <sodipodi:guide
+       orientation="0,1"
+       position="100,25"
+       id="guide4148" />
+  </sodipodi:namedview>
+  <metadata
+     id="metadata7">
+    <rdf:RDF>
+      <cc:Work
+         rdf:about="">
+        <dc:format>image/svg+xml</dc:format>
+        <dc:type
+           rdf:resource="http://purl.org/dc/dcmitype/StillImage" />
+        <dc:title></dc:title>
+      </cc:Work>
+    </rdf:RDF>
+  </metadata>
+  <g
+     inkscape:groupmode="layer"
+     id="layer2"
+     inkscape:label="RD"
+     style="display:none">
+    <path
+       style="fill:none;stroke:#000000;stroke-width:3;stroke-linecap:butt;stroke-linejoin:miter;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none"
+       d="m 190,10 -50,0 -10,70 50,0"
+       id="path3114"
+       inkscape:connector-curvature="0" />
+    <path
+       style="fill:none;stroke:#000000;stroke-width:3;stroke-linecap:butt;stroke-linejoin:miter;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none"
+       d="m 180,120 -50,0 -10,70 50,0"
+       id="path3116"
+       inkscape:connector-curvature="0" />
+    <path
+       sodipodi:type="arc"
+       style="fill:none;stroke:#000000;stroke-width:3;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none;display:inline"
+       id="path3915-9"
+       sodipodi:cx="-145"
+       sodipodi:cy="-35"
+       sodipodi:rx="5"
+       sodipodi:ry="5"
+       d="m -140,-35 a 5,5 0 1 1 -10,0 5,5 0 1 1 10,0 z"
+       transform="translate(274.89948,60)" />
+    <path
+       sodipodi:type="arc"
+       style="fill:none;stroke:#000000;stroke-width:3;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none;display:inline"
+       id="path3915-9-3"
+       sodipodi:cx="-145"
+       sodipodi:cy="-35"
+       sodipodi:rx="5"
+       sodipodi:ry="5"
+       d="m -140,-35 a 5,5 0 1 1 -10,0 5,5 0 1 1 10,0 z"
+       transform="translate(264.82935,170)" />
+  </g>
+  <g
+     inkscape:groupmode="layer"
+     id="layer3"
+     inkscape:label="LS"
+     style="display:inline" />
+  <g
+     inkscape:groupmode="layer"
+     id="layer4"
+     inkscape:label="LD"
+     style="display:none">
+    <path
+       style="fill:none;stroke:#000000;stroke-width:3;stroke-linecap:butt;stroke-linejoin:miter;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none"
+       d="m 10,190 50,0 10,-70 -50,0"
+       id="path3911"
+       inkscape:connector-curvature="0" />
+    <path
+       style="fill:none;stroke:#000000;stroke-width:3;stroke-linecap:butt;stroke-linejoin:miter;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none"
+       d="M 20,80 70,80 80,10 30,10"
+       id="path3913"
+       inkscape:connector-curvature="0" />
+    <path
+       sodipodi:type="arc"
+       style="fill:none;stroke:#000000;stroke-width:3;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none;display:inline"
+       id="path3915-9-32"
+       sodipodi:cx="-145"
+       sodipodi:cy="-35"
+       sodipodi:rx="5"
+       sodipodi:ry="5"
+       d="m -140,-35 a 5,5 0 1 1 -10,0 5,5 0 1 1 10,0 z"
+       transform="translate(230.89286,60)" />
+    <path
+       sodipodi:type="arc"
+       style="fill:none;stroke:#000000;stroke-width:3;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none;display:inline"
+       id="path3915-9-9"
+       sodipodi:cx="-145"
+       sodipodi:cy="-35"
+       sodipodi:rx="5"
+       sodipodi:ry="5"
+       d="m -140,-35 a 5,5 0 1 1 -10,0 5,5 0 1 1 10,0 z"
+       transform="translate(220.85072,170)" />
+  </g>
+  <g
+     inkscape:groupmode="layer"
+     id="layer6"
+     inkscape:label="KL"
+     style="display:none">
+    <path
+       sodipodi:type="arc"
+       style="fill:none;stroke:#000000;stroke-width:3;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none;display:inline"
+       id="path3915-9-8-0-9"
+       sodipodi:cx="-145"
+       sodipodi:cy="-35"
+       sodipodi:rx="5"
+       sodipodi:ry="5"
+       d="m -140,-35 a 5,5 0 1 1 -10,0 5,5 0 1 1 10,0 z"
+       transform="translate(216.94404,190)" />
+  </g>
+  <g
+     inkscape:groupmode="layer"
+     id="layer5"
+     inkscape:label="KR"
+     style="display:none">
+    <path
+       sodipodi:type="arc"
+       style="fill:none;stroke:#000000;stroke-width:3;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none;display:inline"
+       id="path3915-9-8-2"
+       sodipodi:cx="-145"
+       sodipodi:cy="-35"
+       sodipodi:rx="5"
+       sodipodi:ry="5"
+       d="m -140,-35 a 5,5 0 1 1 -10,0 5,5 0 1 1 10,0 z"
+       transform="translate(260.85187,190)" />
+  </g>
+  <g
+     inkscape:label="RS"
+     inkscape:groupmode="layer"
+     id="layer1"
+     transform="translate(0,-852.36218)"
+     style="display:inline">
+    <rect
+       style="fill:none;fill-rule:evenodd;stroke:#000000;stroke-width:1px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1"
+       id="rect3154"
+       width="120"
+       height="118.57143"
+       x="40"
+       y="36.42857"
+       transform="translate(0,852.36218)" />
+    <path
+       style="fill:none;stroke:#000000;stroke-width:1px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1"
+       d="M 40,6.0714286 40.357143,35.714286 100,36.071429 l 0,-30.0000004 0,30.0000004 60.35714,0.714285 0,-30.7142854"
+       id="path3924"
+       inkscape:connector-curvature="0"
+       transform="translate(0,852.36218)" />
+    <path
+       style="fill:none;stroke:#000000;stroke-width:1px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1"
+       d="m 160.71429,36.428571 33.92857,-0.357142"
+       id="path3926"
+       inkscape:connector-curvature="0"
+       transform="translate(0,852.36218)" />
+    <path
+       style="fill:none;stroke:#000000;stroke-width:1px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1"
+       d="m 161.07143,155 28.57143,0"
+       id="path3928"
+       inkscape:connector-curvature="0"
+       transform="translate(0,852.36218)" />
+    <path
+       style="fill:none;stroke:#000000;stroke-width:1px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1"
+       d="m 40,155 0,39.28571"
+       id="path3930"
+       inkscape:connector-curvature="0"
+       transform="translate(0,852.36218)" />
+    <path
+       style="fill:none;stroke:#000000;stroke-width:1px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1"
+       d="m 160,155 0,39.64286"
+       id="path3932"
+       inkscape:connector-curvature="0"
+       transform="translate(0,852.36218)" />
+    <text
+       xml:space="preserve"
+       style="font-size:40px;font-style:normal;font-weight:normal;line-height:125%;letter-spacing:0px;word-spacing:0px;fill:#000000;fill-opacity:1;stroke:none;font-family:Sans"
+       x="64.64286"
+       y="151.42857"
+       id="text3934"
+       sodipodi:linespacing="125%"
+       transform="translate(0,852.36218)"><tspan
+         sodipodi:role="line"
+         id="tspan3936"
+         x="64.64286"
+         y="151.42857"
+         style="font-size:10px;font-style:normal;font-variant:normal;font-weight:normal;font-stretch:normal;text-align:start;line-height:125%;writing-mode:lr-tb;text-anchor:start;font-family:Sans;-inkscape-font-specification:Sans" /></text>
+    <text
+       xml:space="preserve"
+       style="font-size:10px;font-style:normal;font-variant:normal;font-weight:normal;font-stretch:normal;text-align:start;line-height:125%;letter-spacing:0px;word-spacing:0px;writing-mode:lr-tb;text-anchor:start;fill:#000000;fill-opacity:1;stroke:none;font-family:Sans;-inkscape-font-specification:Sans"
+       x="64.700272"
+       y="151.20963"
+       id="text3942"
+       sodipodi:linespacing="125%"
+       transform="translate(0,852.36218)"><tspan
+         sodipodi:role="line"
+         id="tspan3944"
+         x="64.700272"
+         y="151.20963">1</tspan></text>
+    <text
+       xml:space="preserve"
+       style="font-size:10px;font-style:normal;font-variant:normal;font-weight:normal;font-stretch:normal;text-align:start;line-height:125%;letter-spacing:0px;word-spacing:0px;writing-mode:lr-tb;text-anchor:start;fill:#000000;fill-opacity:1;stroke:none;font-family:Sans;-inkscape-font-specification:Sans"
+       x="131.87541"
+       y="152.9774"
+       id="text3946"
+       sodipodi:linespacing="125%"
+       transform="translate(0,852.36218)"><tspan
+         sodipodi:role="line"
+         id="tspan3948"></tspan></text>
+    <text
+       xml:space="preserve"
+       style="font-size:10px;font-style:normal;font-variant:normal;font-weight:normal;font-stretch:normal;text-align:start;line-height:125%;letter-spacing:0px;word-spacing:0px;writing-mode:lr-tb;text-anchor:start;fill:#000000;fill-opacity:1;stroke:none;font-family:Sans;-inkscape-font-specification:Sans"
+       x="150.26019"
+       y="140.24948"
+       id="text3966"
+       sodipodi:linespacing="125%"
+       transform="translate(0,852.36218)"><tspan
+         sodipodi:role="line"
+         id="tspan3968"
+         x="150.26019"
+         y="140.24948">2</tspan></text>
+    <text
+       xml:space="preserve"
+       style="font-size:10px;font-style:normal;font-variant:normal;font-weight:normal;font-stretch:normal;text-align:start;line-height:125%;letter-spacing:0px;word-spacing:0px;writing-mode:lr-tb;text-anchor:start;fill:#000000;fill-opacity:1;stroke:none;font-family:Sans;-inkscape-font-specification:Sans"
+       x="128.33987"
+       y="50.093361"
+       id="text3970"
+       sodipodi:linespacing="125%"
+       transform="translate(0,852.36218)"><tspan
+         sodipodi:role="line"
+         id="tspan3972"
+         x="128.33987"
+         y="50.093361">3</tspan></text>
+    <text
+       xml:space="preserve"
+       style="font-size:10px;font-style:normal;font-variant:normal;font-weight:normal;font-stretch:normal;text-align:start;line-height:125%;letter-spacing:0px;word-spacing:0px;writing-mode:lr-tb;text-anchor:start;fill:#000000;fill-opacity:1;stroke:none;font-family:Sans;-inkscape-font-specification:Sans"
+       x="43.487068"
+       y="72.367226"
+       id="text3974"
+       sodipodi:linespacing="125%"
+       transform="translate(0,852.36218)"><tspan
+         sodipodi:role="line"
+         id="tspan3976"
+         x="43.487068"
+         y="72.367226">4</tspan></text>
+    <text
+       xml:space="preserve"
+       style="font-size:10px;font-style:normal;font-variant:normal;font-weight:normal;font-stretch:normal;text-align:start;line-height:125%;letter-spacing:0px;word-spacing:0px;writing-mode:lr-tb;text-anchor:start;fill:#000000;fill-opacity:1;stroke:none;font-family:Sans;-inkscape-font-specification:Sans"
+       x="121.62237"
+       y="151.20963"
+       id="text3978"
+       sodipodi:linespacing="125%"
+       transform="translate(0,852.36218)"><tspan
+         sodipodi:role="line"
+         id="tspan3980"
+         x="121.62237"
+         y="151.20963">5</tspan></text>
+    <text
+       xml:space="preserve"
+       style="font-size:10px;font-style:normal;font-variant:normal;font-weight:normal;font-stretch:normal;text-align:start;line-height:125%;letter-spacing:0px;word-spacing:0px;writing-mode:lr-tb;text-anchor:start;fill:#000000;fill-opacity:1;stroke:none;font-family:Sans;-inkscape-font-specification:Sans"
+       x="150.9673"
+       y="79.438293"
+       id="text3982"
+       sodipodi:linespacing="125%"
+       transform="translate(0,852.36218)"><tspan
+         sodipodi:role="line"
+         id="tspan3984"
+         x="150.9673"
+         y="79.438293">6</tspan></text>
+    <text
+       xml:space="preserve"
+       style="font-size:10px;font-style:normal;font-variant:normal;font-weight:normal;font-stretch:normal;text-align:start;line-height:125%;letter-spacing:0px;word-spacing:0px;writing-mode:lr-tb;text-anchor:start;fill:#000000;fill-opacity:1;stroke:none;font-family:Sans;-inkscape-font-specification:Sans"
+       x="67.528702"
+       y="47.618488"
+       id="text3986"
+       sodipodi:linespacing="125%"
+       transform="translate(0,852.36218)"><tspan
+         sodipodi:role="line"
+         id="tspan3988"
+         x="67.528702"
+         y="47.618488">7</tspan></text>
+    <text
+       xml:space="preserve"
+       style="font-size:10px;font-style:normal;font-variant:normal;font-weight:normal;font-stretch:normal;text-align:start;line-height:125%;letter-spacing:0px;word-spacing:0px;writing-mode:lr-tb;text-anchor:start;fill:#000000;fill-opacity:1;stroke:none;font-family:Sans;-inkscape-font-specification:Sans"
+       x="44.547726"
+       y="119.03627"
+       id="text3990"
+       sodipodi:linespacing="125%"
+       transform="translate(0,852.36218)"><tspan
+         sodipodi:role="line"
+         id="tspan3992"
+         x="44.547726"
+         y="119.03627">8</tspan></text>
+  </g>
+</svg>
diff --git a/doc/fig/exmpl1_coo.tex b/doc/fig/exmpl1_coo.tex
new file mode 100644 (file)
index 0000000..0f61052
--- /dev/null
@@ -0,0 +1,13 @@
+\begin{tabular}{>{\columncolor{gray}}rccc}
+\rowcolor{gray}
+ Index & x1 & x2 & x3\\
+   1 & 0.0000 & 0.0000 & 0.0000\\
+   2 & 0.0000 & 0.5000 & 0.0000\\
+   3 & 0.0000 & 1.0000 & 0.0000\\
+   4 & 0.5000 & 0.0000 & 0.0000\\
+   5 & 0.5000 & 0.5000 & 0.0000\\
+   6 & 0.5000 & 1.0000 & 0.0000\\
+   7 & 1.0000 & 0.0000 & 0.0000\\
+   8 & 1.0000 & 0.5000 & 0.0000\\
+   9 & 1.0000 & 1.0000 & 0.0000
+\end{tabular}
diff --git a/doc/fig/exmpl1_ele.tex b/doc/fig/exmpl1_ele.tex
new file mode 100644 (file)
index 0000000..4f55ff0
--- /dev/null
@@ -0,0 +1,9 @@
+\begin{tabular}{>{\columncolor{gray}}rcccc}
+\rowcolor{gray}
+ Index & c1 & c2 & c3 & c4\\
+   1 & 2 & 5 & 6 & 3\\
+   2 & 5 & 8 & 9 & 6\\
+   3 & 4 & 7 & 8 & 5\\
+   4 & 1 & 4 & 5 & 2
+\end{tabular}
+\label{exmpl1:ele}
\ No newline at end of file
diff --git a/doc/fig/exmpl1_f2s.tex b/doc/fig/exmpl1_f2s.tex
new file mode 100644 (file)
index 0000000..5037733
--- /dev/null
@@ -0,0 +1,16 @@
+\begin{tabular}{>{\columncolor{gray}}rcccc}
+\rowcolor{gray}
+ Index & e1 & e2 & e3 & e4\\
+   1 & 1 & 1 & 1 & 1\\
+   2 & 2 & 2 & 2 & 2\\
+   3 & 3 & 3 & 3 & 3\\
+   4 & 4 & 4 & 4 & 4\\
+   5 & 12 & 14 & 13 & 5\\
+   6 & 6 & 6 & 6 & 6\\
+   7 & 7 & 7 & 7 & 7\\
+   8 & 8 & 8 & 8 & 8\\
+   9 & 9 & 9 & 9 & 9\\
+   10 & 10 & 10 & 10 & 10\\
+   11 & 11 & 11 & 11 & 11
+\end{tabular}
+\label{exmpl1:f2s}
\ No newline at end of file
diff --git a/doc/fig/exmpl1_nei.tex b/doc/fig/exmpl1_nei.tex
new file mode 100644 (file)
index 0000000..8d7f3c9
--- /dev/null
@@ -0,0 +1,9 @@
+\begin{tabular}{>{\columncolor{gray}}rcccccccc}
+\rowcolor{gray}
+ Index & n1 & n2 & n3 & n4 & n5 & n6 & n7 & n8\\
+   1 & 4 & 2 & 0 & 0 & 0 & 0 & 0 & 0\\
+   2 & 3 & 0 & 0 & 1 & 0 & 0 & 0 & 0\\
+   3 & 0 & 0 & 2 & 4 & 0 & 0 & 0 & 0\\
+   4 & 0 & 3 & 1 & 0 & 0 & 0 & 0 & 0
+\end{tabular}
+\label{exmpl1:nei}
\ No newline at end of file
diff --git a/doc/fig/exmpl1_ref.eps b/doc/fig/exmpl1_ref.eps
new file mode 100644 (file)
index 0000000..f5c1c4a
--- /dev/null
@@ -0,0 +1,302 @@
+%!PS-Adobe-2.0 EPSF-1.2
+%%Creator: MATLAB, The MathWorks, Inc. Version 7.13.0.564 (R2011b). Operating System: Linux 3.0.0-16-generic #28-Ubuntu SMP Fri Jan 27 17:44:39 UTC 2012 x86_64.
+%%Title: ../doc/fig/exmpl1_ref.eps
+%%CreationDate: 03/13/2012  13:17:38
+%%DocumentNeededFonts: Helvetica
+%%DocumentProcessColors: Cyan Magenta Yellow Black
+%%Extensions: CMYK
+%%Pages: 1
+%%BoundingBox:    58   199   546   589
+%%EndComments
+
+%%BeginProlog
+% MathWorks dictionary
+/MathWorks 160 dict begin
+% definition operators
+/bdef {bind def} bind def
+/ldef {load def} bind def
+/xdef {exch def} bdef
+/xstore {exch store} bdef
+% operator abbreviations
+/c  /clip ldef
+/cc /concat ldef
+/cp /closepath ldef
+/gr /grestore ldef
+/gs /gsave ldef
+/mt /moveto ldef
+/np /newpath ldef
+/cm /currentmatrix ldef
+/sm /setmatrix ldef
+/rm /rmoveto ldef
+/rl /rlineto ldef
+/s {show newpath} bdef
+/sc {setcmykcolor} bdef
+/sr /setrgbcolor ldef
+/sg /setgray ldef
+/w /setlinewidth ldef
+/j /setlinejoin ldef
+/cap /setlinecap ldef
+/rc {rectclip} bdef
+/rf {rectfill} bdef
+% page state control
+/pgsv () def
+/bpage {/pgsv save def} bdef
+/epage {pgsv restore} bdef
+/bplot /gsave ldef
+/eplot {stroke grestore} bdef
+% orientation switch
+/portraitMode 0 def /landscapeMode 1 def /rotateMode 2 def
+% coordinate system mappings
+/dpi2point 0 def
+% font control
+/FontSize 0 def
+/FMS {/FontSize xstore findfont [FontSize 0 0 FontSize neg 0 0]
+  makefont setfont} bdef
+/reencode {exch dup where {pop load} {pop StandardEncoding} ifelse
+  exch dup 3 1 roll findfont dup length dict begin
+  { 1 index /FID ne {def}{pop pop} ifelse } forall
+  /Encoding exch def currentdict end definefont pop} bdef
+/isroman {findfont /CharStrings get /Agrave known} bdef
+/FMSR {3 1 roll 1 index dup isroman {reencode} {pop pop} ifelse
+  exch FMS} bdef
+/csm {1 dpi2point div -1 dpi2point div scale neg translate
+ dup landscapeMode eq {pop -90 rotate}
+  {rotateMode eq {90 rotate} if} ifelse} bdef
+% line types: solid, dotted, dashed, dotdash
+/SO { [] 0 setdash } bdef
+/DO { [.5 dpi2point mul 4 dpi2point mul] 0 setdash } bdef
+/DA { [6 dpi2point mul] 0 setdash } bdef
+/DD { [.5 dpi2point mul 4 dpi2point mul 6 dpi2point mul 4
+  dpi2point mul] 0 setdash } bdef
+% macros for lines and objects
+/L {lineto stroke} bdef
+/MP {3 1 roll moveto 1 sub {rlineto} repeat} bdef
+/AP {{rlineto} repeat} bdef
+/PDlw -1 def
+/W {/PDlw currentlinewidth def setlinewidth} def
+/PP {closepath eofill} bdef
+/DP {closepath stroke} bdef
+/MR {4 -2 roll moveto dup  0 exch rlineto exch 0 rlineto
+  neg 0 exch rlineto closepath} bdef
+/FR {MR stroke} bdef
+/PR {MR fill} bdef
+/L1i {{currentfile picstr readhexstring pop} image} bdef
+/tMatrix matrix def
+/MakeOval {newpath tMatrix currentmatrix pop translate scale
+0 0 1 0 360 arc tMatrix setmatrix} bdef
+/FO {MakeOval stroke} bdef
+/PO {MakeOval fill} bdef
+/PD {currentlinewidth 2 div 0 360 arc fill
+   PDlw -1 eq not {PDlw w /PDlw -1 def} if} def
+/FA {newpath tMatrix currentmatrix pop translate scale
+  0 0 1 5 -2 roll arc tMatrix setmatrix stroke} bdef
+/PA {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
+  0 0 1 5 -2 roll arc closepath tMatrix setmatrix fill} bdef
+/FAn {newpath tMatrix currentmatrix pop translate scale
+  0 0 1 5 -2 roll arcn tMatrix setmatrix stroke} bdef
+/PAn {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
+  0 0 1 5 -2 roll arcn closepath tMatrix setmatrix fill} bdef
+/vradius 0 def /hradius 0 def /lry 0 def
+/lrx 0 def /uly 0 def /ulx 0 def /rad 0 def
+/MRR {/vradius xdef /hradius xdef /lry xdef /lrx xdef /uly xdef
+  /ulx xdef newpath tMatrix currentmatrix pop ulx hradius add uly
+  vradius add translate hradius vradius scale 0 0 1 180 270 arc 
+  tMatrix setmatrix lrx hradius sub uly vradius add translate
+  hradius vradius scale 0 0 1 270 360 arc tMatrix setmatrix
+  lrx hradius sub lry vradius sub translate hradius vradius scale
+  0 0 1 0 90 arc tMatrix setmatrix ulx hradius add lry vradius sub
+  translate hradius vradius scale 0 0 1 90 180 arc tMatrix setmatrix
+  closepath} bdef
+/FRR {MRR stroke } bdef
+/PRR {MRR fill } bdef
+/MlrRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lry uly sub 2 div def
+  newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
+  rad rad scale 0 0 1 90 270 arc tMatrix setmatrix lrx rad sub lry rad
+  sub translate rad rad scale 0 0 1 270 90 arc tMatrix setmatrix
+  closepath} bdef
+/FlrRR {MlrRR stroke } bdef
+/PlrRR {MlrRR fill } bdef
+/MtbRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lrx ulx sub 2 div def
+  newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
+  rad rad scale 0 0 1 180 360 arc tMatrix setmatrix lrx rad sub lry rad
+  sub translate rad rad scale 0 0 1 0 180 arc tMatrix setmatrix
+  closepath} bdef
+/FtbRR {MtbRR stroke } bdef
+/PtbRR {MtbRR fill } bdef
+/stri 6 array def /dtri 6 array def
+/smat 6 array def /dmat 6 array def
+/tmat1 6 array def /tmat2 6 array def /dif 3 array def
+/asub {/ind2 exch def /ind1 exch def dup dup
+  ind1 get exch ind2 get sub exch } bdef
+/tri_to_matrix {
+  2 0 asub 3 1 asub 4 0 asub 5 1 asub
+  dup 0 get exch 1 get 7 -1 roll astore } bdef
+/compute_transform {
+  dmat dtri tri_to_matrix tmat1 invertmatrix 
+  smat stri tri_to_matrix tmat2 concatmatrix } bdef
+/ds {stri astore pop} bdef
+/dt {dtri astore pop} bdef
+/db {2 copy /cols xdef /rows xdef mul dup 3 mul string
+  currentfile exch readhexstring pop
+  dup 0 3 index getinterval /rbmap xdef
+  dup 2 index dup getinterval /gbmap xdef
+  1 index dup 2 mul exch getinterval /bbmap xdef pop pop}bdef
+/it {gs np dtri aload pop moveto lineto lineto cp c
+  cols rows 8 compute_transform 
+  rbmap gbmap bbmap true 3 colorimage gr}bdef
+/il {newpath moveto lineto stroke}bdef
+currentdict end def
+%%EndProlog
+
+%%BeginSetup
+MathWorks begin
+
+0 cap
+
+end
+%%EndSetup
+
+%%Page: 1 1
+%%BeginPageSetup
+%%PageBoundingBox:    58   199   546   589
+MathWorks begin
+bpage
+%%EndPageSetup
+
+%%BeginObject: obj1
+bplot
+
+/dpi2point 8.33333 def
+portraitMode 0150 5100 csm
+
+  334   190  4067  3251 MR c np
+85 dict begin %Colortable dictionary
+/c0 { 0.000000 0.000000 0.000000 sr} bdef
+/c1 { 1.000000 1.000000 1.000000 sr} bdef
+/c2 { 0.900000 0.000000 0.000000 sr} bdef
+/c3 { 0.000000 0.820000 0.000000 sr} bdef
+/c4 { 0.000000 0.000000 0.800000 sr} bdef
+/c5 { 0.910000 0.820000 0.320000 sr} bdef
+/c6 { 1.000000 0.260000 0.820000 sr} bdef
+/c7 { 0.000000 0.820000 0.820000 sr} bdef
+c0
+1 j
+1 sg
+   0    0 4801 3602 PR
+4.16667 w
+0 2935 3720 0 0 -2935 624 3205 4 MP
+PP
+-3720 0 0 2935 3720 0 0 -2935 624 3205 5 MP stroke
+2.77778 w
+DO
+SO
+4.16667 w
+0 sg
+ 624 3205 mt 4344 3205 L
+ 624 3205 mt  624  270 L
+ 624 3205 mt  624 3167 L
+%%IncludeResource: font Helvetica
+/Helvetica /ISOLatin1Encoding 83.3333 FMSR
+
+ 601 3306 mt 
+(0) s
+ 996 3205 mt  996 3167 L
+ 939 3306 mt 
+(0.1) s
+1368 3205 mt 1368 3167 L
+1311 3306 mt 
+(0.2) s
+1740 3205 mt 1740 3167 L
+1683 3306 mt 
+(0.3) s
+2112 3205 mt 2112 3167 L
+2055 3306 mt 
+(0.4) s
+2484 3205 mt 2484 3167 L
+2427 3306 mt 
+(0.5) s
+2856 3205 mt 2856 3167 L
+2799 3306 mt 
+(0.6) s
+3228 3205 mt 3228 3167 L
+3171 3306 mt 
+(0.7) s
+3600 3205 mt 3600 3167 L
+3543 3306 mt 
+(0.8) s
+3972 3205 mt 3972 3167 L
+3915 3306 mt 
+(0.9) s
+4344 3205 mt 4344 3167 L
+4321 3306 mt 
+(1) s
+ 624 3205 mt  661 3205 L
+ 554 3235 mt 
+(0) s
+ 624 2911 mt  661 2911 L
+ 485 2941 mt 
+(0.1) s
+ 624 2618 mt  661 2618 L
+ 485 2648 mt 
+(0.2) s
+ 624 2324 mt  661 2324 L
+ 485 2354 mt 
+(0.3) s
+ 624 2031 mt  661 2031 L
+ 485 2061 mt 
+(0.4) s
+ 624 1737 mt  661 1737 L
+ 485 1767 mt 
+(0.5) s
+ 624 1443 mt  661 1443 L
+ 485 1473 mt 
+(0.6) s
+ 624 1150 mt  661 1150 L
+ 485 1180 mt 
+(0.7) s
+ 624  857 mt  661  857 L
+ 485  887 mt 
+(0.8) s
+ 624  563 mt  661  563 L
+ 485  593 mt 
+(0.9) s
+ 624  270 mt  661  270 L
+ 554  300 mt 
+(1) s
+gs 624 270 3721 2936 MR c np
+/c8 { 0.000000 0.000000 1.000000 sr} bdef
+c8
+0 1467 -1860 0 0 -1467 1860 0 624 1737 5 MP stroke
+0 1467 -1860 0 0 -1467 1860 0 2484 1737 5 MP stroke
+0 1468 -1860 0 0 -1468 1860 0 2484 3205 5 MP stroke
+0 1468 -1860 0 0 -1468 1860 0 624 3205 5 MP stroke
+gr
+
+c8
+0 sg
+1554 1033 mt 
+(1) s
+3414 1033 mt 
+(2) s
+3414 2501 mt 
+(3) s
+1554 2501 mt 
+(4) s
+2462 3406 mt 
+(x) s
+ 428 1759 mt  -90 rotate
+(y) s
+90 rotate
+
+end %%Color Dict
+
+eplot
+%%EndObject
+
+epage
+end
+
+showpage
+
+%%Trailer
+%%EOF
diff --git a/doc/fig/exmpl1_ref.pdf b/doc/fig/exmpl1_ref.pdf
new file mode 100644 (file)
index 0000000..ff433d0
Binary files /dev/null and b/doc/fig/exmpl1_ref.pdf differ
diff --git a/doc/fig/exmpl2_coo.tex b/doc/fig/exmpl2_coo.tex
new file mode 100644 (file)
index 0000000..4597de8
--- /dev/null
@@ -0,0 +1,23 @@
+\begin{tabular}{>{\columncolor{gray}}rccc}
+\rowcolor{gray}
+ Index & x1 & x2 & x3\\
+   1 & 0.0000 & 0.0000 & 0.0000\\
+   2 & 0.0000 & 0.2500 & 0.0000\\
+   3 & 0.0000 & 0.5000 & 0.0000\\
+   4 & 0.0000 & 0.7500 & 0.0000\\
+   5 & 0.0000 & 1.0000 & 0.0000\\
+   6 & 0.2500 & 0.0000 & 0.0000\\
+   7 & 0.2500 & 0.2500 & 0.0000\\
+   8 & 0.2500 & 0.5000 & 0.0000\\
+   9 & 0.2500 & 0.7500 & 0.0000\\
+   10 & 0.2500 & 1.0000 & 0.0000\\
+   11 & 0.5000 & 0.0000 & 0.0000\\
+   12 & 0.5000 & 0.2500 & 0.0000\\
+   13 & 0.5000 & 0.5000 & 0.0000\\
+   14 & 0.5000 & 0.7500 & 0.0000\\
+   15 & 0.5000 & 1.0000 & 0.0000\\
+   16 & 1.0000 & 0.0000 & 0.0000\\
+   17 & 1.0000 & 0.5000 & 0.0000\\
+   18 & 1.0000 & 0.7500 & 0.0000\\
+   19 & 1.0000 & 1.0000 & 0.0000
+\end{tabular}
diff --git a/doc/fig/exmpl2_ele.tex b/doc/fig/exmpl2_ele.tex
new file mode 100644 (file)
index 0000000..8d7c7da
--- /dev/null
@@ -0,0 +1,16 @@
+\begin{tabular}{>{\columncolor{gray}}rcccc}
+\rowcolor{gray}
+ Index & c1 & c2 & c3 & c4\\
+   1 & 4 & 9 & 10 & 5\\
+   2 & 14 & 18 & 19 & 15\\
+   3 & 11 & 16 & 17 & 13\\
+   4 & 2 & 7 & 8 & 3\\
+   5 & 3 & 8 & 9 & 4\\
+   6 & 13 & 17 & 18 & 14\\
+   7 & 7 & 12 & 13 & 8\\
+   8 & 6 & 11 & 12 & 7\\
+   9 & 1 & 6 & 7 & 2\\
+   10 & 9 & 14 & 15 & 10\\
+   11 & 8 & 13 & 14 & 9
+\end{tabular}
+\label{exmpl2:ele}
\ No newline at end of file
diff --git a/doc/fig/exmpl2_f2s.tex b/doc/fig/exmpl2_f2s.tex
new file mode 100644 (file)
index 0000000..c7d66b4
--- /dev/null
@@ -0,0 +1,9 @@
+\begin{tabular}{>{\columncolor{gray}}rcccc}
+\rowcolor{gray}
+ Index & e1 & e2 & e3 & e4\\
+   1 & 5 & 11 & 10 & 1\\
+   2 & 6 & 6 & 2 & 2\\
+   3 & 3 & 3 & 3 & 3\\
+   4 & 9 & 8 & 7 & 4
+\end{tabular}
+\label{exmpl2:f2s}
\ No newline at end of file
diff --git a/doc/fig/exmpl2_nei.tex b/doc/fig/exmpl2_nei.tex
new file mode 100644 (file)
index 0000000..9da83b6
--- /dev/null
@@ -0,0 +1,16 @@
+\begin{tabular}{>{\columncolor{gray}}rcccccccc}
+\rowcolor{gray}
+ Index & n1 & n2 & n3 & n4 & n5 & n6 & n7 & n8\\
+   1 & 5 & 10 & 0 & 0 & 0 & 0 & 0 & 0\\
+   2 & 6 & 0 & 0 & 10 & 0 & 0 & 0 & 0\\
+   3 & 0 & 0 & 6 & 8 & 0 & 0 & 0 & 7\\
+   4 & 9 & 7 & 5 & 0 & 0 & 0 & 0 & 0\\
+   5 & 4 & 11 & 1 & 0 & 0 & 0 & 0 & 0\\
+   6 & 3 & 0 & 2 & 11 & 0 & 0 & 0 & 0\\
+   7 & 8 & 3 & 11 & 4 & 0 & 0 & 0 & 0\\
+   8 & 0 & 3 & 7 & 9 & 0 & 0 & 0 & 0\\
+   9 & 0 & 8 & 4 & 0 & 0 & 0 & 0 & 0\\
+   10 & 11 & 2 & 0 & 1 & 0 & 0 & 0 & 0\\
+   11 & 7 & 6 & 10 & 5 & 0 & 0 & 0 & 0
+\end{tabular}
+\label{exmpl2:nei}
\ No newline at end of file
diff --git a/doc/fig/exmpl2_ref.eps b/doc/fig/exmpl2_ref.eps
new file mode 100644 (file)
index 0000000..fe1ab01
--- /dev/null
@@ -0,0 +1,323 @@
+%!PS-Adobe-2.0 EPSF-1.2
+%%Creator: MATLAB, The MathWorks, Inc. Version 7.13.0.564 (R2011b). Operating System: Linux 3.0.0-16-generic #28-Ubuntu SMP Fri Jan 27 17:44:39 UTC 2012 x86_64.
+%%Title: ../doc/fig/exmpl2_ref.eps
+%%CreationDate: 03/13/2012  13:17:50
+%%DocumentNeededFonts: Helvetica
+%%DocumentProcessColors: Cyan Magenta Yellow Black
+%%Extensions: CMYK
+%%Pages: 1
+%%BoundingBox:    58   199   546   589
+%%EndComments
+
+%%BeginProlog
+% MathWorks dictionary
+/MathWorks 160 dict begin
+% definition operators
+/bdef {bind def} bind def
+/ldef {load def} bind def
+/xdef {exch def} bdef
+/xstore {exch store} bdef
+% operator abbreviations
+/c  /clip ldef
+/cc /concat ldef
+/cp /closepath ldef
+/gr /grestore ldef
+/gs /gsave ldef
+/mt /moveto ldef
+/np /newpath ldef
+/cm /currentmatrix ldef
+/sm /setmatrix ldef
+/rm /rmoveto ldef
+/rl /rlineto ldef
+/s {show newpath} bdef
+/sc {setcmykcolor} bdef
+/sr /setrgbcolor ldef
+/sg /setgray ldef
+/w /setlinewidth ldef
+/j /setlinejoin ldef
+/cap /setlinecap ldef
+/rc {rectclip} bdef
+/rf {rectfill} bdef
+% page state control
+/pgsv () def
+/bpage {/pgsv save def} bdef
+/epage {pgsv restore} bdef
+/bplot /gsave ldef
+/eplot {stroke grestore} bdef
+% orientation switch
+/portraitMode 0 def /landscapeMode 1 def /rotateMode 2 def
+% coordinate system mappings
+/dpi2point 0 def
+% font control
+/FontSize 0 def
+/FMS {/FontSize xstore findfont [FontSize 0 0 FontSize neg 0 0]
+  makefont setfont} bdef
+/reencode {exch dup where {pop load} {pop StandardEncoding} ifelse
+  exch dup 3 1 roll findfont dup length dict begin
+  { 1 index /FID ne {def}{pop pop} ifelse } forall
+  /Encoding exch def currentdict end definefont pop} bdef
+/isroman {findfont /CharStrings get /Agrave known} bdef
+/FMSR {3 1 roll 1 index dup isroman {reencode} {pop pop} ifelse
+  exch FMS} bdef
+/csm {1 dpi2point div -1 dpi2point div scale neg translate
+ dup landscapeMode eq {pop -90 rotate}
+  {rotateMode eq {90 rotate} if} ifelse} bdef
+% line types: solid, dotted, dashed, dotdash
+/SO { [] 0 setdash } bdef
+/DO { [.5 dpi2point mul 4 dpi2point mul] 0 setdash } bdef
+/DA { [6 dpi2point mul] 0 setdash } bdef
+/DD { [.5 dpi2point mul 4 dpi2point mul 6 dpi2point mul 4
+  dpi2point mul] 0 setdash } bdef
+% macros for lines and objects
+/L {lineto stroke} bdef
+/MP {3 1 roll moveto 1 sub {rlineto} repeat} bdef
+/AP {{rlineto} repeat} bdef
+/PDlw -1 def
+/W {/PDlw currentlinewidth def setlinewidth} def
+/PP {closepath eofill} bdef
+/DP {closepath stroke} bdef
+/MR {4 -2 roll moveto dup  0 exch rlineto exch 0 rlineto
+  neg 0 exch rlineto closepath} bdef
+/FR {MR stroke} bdef
+/PR {MR fill} bdef
+/L1i {{currentfile picstr readhexstring pop} image} bdef
+/tMatrix matrix def
+/MakeOval {newpath tMatrix currentmatrix pop translate scale
+0 0 1 0 360 arc tMatrix setmatrix} bdef
+/FO {MakeOval stroke} bdef
+/PO {MakeOval fill} bdef
+/PD {currentlinewidth 2 div 0 360 arc fill
+   PDlw -1 eq not {PDlw w /PDlw -1 def} if} def
+/FA {newpath tMatrix currentmatrix pop translate scale
+  0 0 1 5 -2 roll arc tMatrix setmatrix stroke} bdef
+/PA {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
+  0 0 1 5 -2 roll arc closepath tMatrix setmatrix fill} bdef
+/FAn {newpath tMatrix currentmatrix pop translate scale
+  0 0 1 5 -2 roll arcn tMatrix setmatrix stroke} bdef
+/PAn {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
+  0 0 1 5 -2 roll arcn closepath tMatrix setmatrix fill} bdef
+/vradius 0 def /hradius 0 def /lry 0 def
+/lrx 0 def /uly 0 def /ulx 0 def /rad 0 def
+/MRR {/vradius xdef /hradius xdef /lry xdef /lrx xdef /uly xdef
+  /ulx xdef newpath tMatrix currentmatrix pop ulx hradius add uly
+  vradius add translate hradius vradius scale 0 0 1 180 270 arc 
+  tMatrix setmatrix lrx hradius sub uly vradius add translate
+  hradius vradius scale 0 0 1 270 360 arc tMatrix setmatrix
+  lrx hradius sub lry vradius sub translate hradius vradius scale
+  0 0 1 0 90 arc tMatrix setmatrix ulx hradius add lry vradius sub
+  translate hradius vradius scale 0 0 1 90 180 arc tMatrix setmatrix
+  closepath} bdef
+/FRR {MRR stroke } bdef
+/PRR {MRR fill } bdef
+/MlrRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lry uly sub 2 div def
+  newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
+  rad rad scale 0 0 1 90 270 arc tMatrix setmatrix lrx rad sub lry rad
+  sub translate rad rad scale 0 0 1 270 90 arc tMatrix setmatrix
+  closepath} bdef
+/FlrRR {MlrRR stroke } bdef
+/PlrRR {MlrRR fill } bdef
+/MtbRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lrx ulx sub 2 div def
+  newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
+  rad rad scale 0 0 1 180 360 arc tMatrix setmatrix lrx rad sub lry rad
+  sub translate rad rad scale 0 0 1 0 180 arc tMatrix setmatrix
+  closepath} bdef
+/FtbRR {MtbRR stroke } bdef
+/PtbRR {MtbRR fill } bdef
+/stri 6 array def /dtri 6 array def
+/smat 6 array def /dmat 6 array def
+/tmat1 6 array def /tmat2 6 array def /dif 3 array def
+/asub {/ind2 exch def /ind1 exch def dup dup
+  ind1 get exch ind2 get sub exch } bdef
+/tri_to_matrix {
+  2 0 asub 3 1 asub 4 0 asub 5 1 asub
+  dup 0 get exch 1 get 7 -1 roll astore } bdef
+/compute_transform {
+  dmat dtri tri_to_matrix tmat1 invertmatrix 
+  smat stri tri_to_matrix tmat2 concatmatrix } bdef
+/ds {stri astore pop} bdef
+/dt {dtri astore pop} bdef
+/db {2 copy /cols xdef /rows xdef mul dup 3 mul string
+  currentfile exch readhexstring pop
+  dup 0 3 index getinterval /rbmap xdef
+  dup 2 index dup getinterval /gbmap xdef
+  1 index dup 2 mul exch getinterval /bbmap xdef pop pop}bdef
+/it {gs np dtri aload pop moveto lineto lineto cp c
+  cols rows 8 compute_transform 
+  rbmap gbmap bbmap true 3 colorimage gr}bdef
+/il {newpath moveto lineto stroke}bdef
+currentdict end def
+%%EndProlog
+
+%%BeginSetup
+MathWorks begin
+
+0 cap
+
+end
+%%EndSetup
+
+%%Page: 1 1
+%%BeginPageSetup
+%%PageBoundingBox:    58   199   546   589
+MathWorks begin
+bpage
+%%EndPageSetup
+
+%%BeginObject: obj1
+bplot
+
+/dpi2point 8.33333 def
+portraitMode 0150 5100 csm
+
+  334   190  4067  3251 MR c np
+85 dict begin %Colortable dictionary
+/c0 { 0.000000 0.000000 0.000000 sr} bdef
+/c1 { 1.000000 1.000000 1.000000 sr} bdef
+/c2 { 0.900000 0.000000 0.000000 sr} bdef
+/c3 { 0.000000 0.820000 0.000000 sr} bdef
+/c4 { 0.000000 0.000000 0.800000 sr} bdef
+/c5 { 0.910000 0.820000 0.320000 sr} bdef
+/c6 { 1.000000 0.260000 0.820000 sr} bdef
+/c7 { 0.000000 0.820000 0.820000 sr} bdef
+c0
+1 j
+1 sg
+   0    0 4801 3602 PR
+4.16667 w
+0 2935 3720 0 0 -2935 624 3205 4 MP
+PP
+-3720 0 0 2935 3720 0 0 -2935 624 3205 5 MP stroke
+2.77778 w
+DO
+SO
+4.16667 w
+0 sg
+ 624 3205 mt 4344 3205 L
+ 624 3205 mt  624  270 L
+ 624 3205 mt  624 3167 L
+%%IncludeResource: font Helvetica
+/Helvetica /ISOLatin1Encoding 83.3333 FMSR
+
+ 601 3306 mt 
+(0) s
+ 996 3205 mt  996 3167 L
+ 939 3306 mt 
+(0.1) s
+1368 3205 mt 1368 3167 L
+1311 3306 mt 
+(0.2) s
+1740 3205 mt 1740 3167 L
+1683 3306 mt 
+(0.3) s
+2112 3205 mt 2112 3167 L
+2055 3306 mt 
+(0.4) s
+2484 3205 mt 2484 3167 L
+2427 3306 mt 
+(0.5) s
+2856 3205 mt 2856 3167 L
+2799 3306 mt 
+(0.6) s
+3228 3205 mt 3228 3167 L
+3171 3306 mt 
+(0.7) s
+3600 3205 mt 3600 3167 L
+3543 3306 mt 
+(0.8) s
+3972 3205 mt 3972 3167 L
+3915 3306 mt 
+(0.9) s
+4344 3205 mt 4344 3167 L
+4321 3306 mt 
+(1) s
+ 624 3205 mt  661 3205 L
+ 554 3235 mt 
+(0) s
+ 624 2911 mt  661 2911 L
+ 485 2941 mt 
+(0.1) s
+ 624 2618 mt  661 2618 L
+ 485 2648 mt 
+(0.2) s
+ 624 2324 mt  661 2324 L
+ 485 2354 mt 
+(0.3) s
+ 624 2031 mt  661 2031 L
+ 485 2061 mt 
+(0.4) s
+ 624 1737 mt  661 1737 L
+ 485 1767 mt 
+(0.5) s
+ 624 1443 mt  661 1443 L
+ 485 1473 mt 
+(0.6) s
+ 624 1150 mt  661 1150 L
+ 485 1180 mt 
+(0.7) s
+ 624  857 mt  661  857 L
+ 485  887 mt 
+(0.8) s
+ 624  563 mt  661  563 L
+ 485  593 mt 
+(0.9) s
+ 624  270 mt  661  270 L
+ 554  300 mt 
+(1) s
+gs 624 270 3721 2936 MR c np
+/c8 { 0.000000 0.000000 1.000000 sr} bdef
+c8
+0 733 -930 0 0 -733 930 0 624 1003 5 MP stroke
+0 733 -1860 0 0 -733 1860 0 2484 1003 5 MP stroke
+0 1468 -1860 0 0 -1468 1860 0 2484 3205 5 MP stroke
+0 734 -930 0 0 -734 930 0 624 2471 5 MP stroke
+0 734 -930 0 0 -734 930 0 624 1737 5 MP stroke
+0 734 -1860 0 0 -734 1860 0 2484 1737 5 MP stroke
+0 734 -930 0 0 -734 930 0 1554 2471 5 MP stroke
+0 734 -930 0 0 -734 930 0 1554 3205 5 MP stroke
+0 734 -930 0 0 -734 930 0 624 3205 5 MP stroke
+0 733 -930 0 0 -733 930 0 1554 1003 5 MP stroke
+0 734 -930 0 0 -734 930 0 1554 1737 5 MP stroke
+gr
+
+c8
+0 sg
+1089  666 mt 
+(1) s
+3414  666 mt 
+(2) s
+3414 2501 mt 
+(3) s
+1089 2134 mt 
+(4) s
+1089 1400 mt 
+(5) s
+3414 1400 mt 
+(6) s
+2019 2134 mt 
+(7) s
+2019 2868 mt 
+(8) s
+1089 2868 mt 
+(9) s
+2019  666 mt 
+(10) s
+2019 1400 mt 
+(11) s
+2462 3406 mt 
+(x) s
+ 428 1759 mt  -90 rotate
+(y) s
+90 rotate
+
+end %%Color Dict
+
+eplot
+%%EndObject
+
+epage
+end
+
+showpage
+
+%%Trailer
+%%EOF
diff --git a/doc/fig/exmpl2_ref.pdf b/doc/fig/exmpl2_ref.pdf
new file mode 100644 (file)
index 0000000..2e36d3b
Binary files /dev/null and b/doc/fig/exmpl2_ref.pdf differ
diff --git a/doc/fig/exmpl3_coo.tex b/doc/fig/exmpl3_coo.tex
new file mode 100644 (file)
index 0000000..700d3b8
--- /dev/null
@@ -0,0 +1,28 @@
+\begin{tabular}{>{\columncolor{gray}}rccc}
+\rowcolor{gray}
+ Index & x1 & x2 & x3\\
+   1 & 0.0000 & 0.0000 & 0.0000\\
+   2 & 0.0000 & 0.2500 & 0.0000\\
+   3 & 0.0000 & 0.5000 & 0.0000\\
+   4 & 0.0000 & 0.6250 & 0.0000\\
+   5 & 0.0000 & 0.7500 & 0.0000\\
+   6 & 0.0000 & 1.0000 & 0.0000\\
+   7 & 0.1250 & 0.5000 & 0.0000\\
+   8 & 0.1250 & 0.6250 & 0.0000\\
+   9 & 0.1250 & 0.7500 & 0.0000\\
+   10 & 0.2500 & 0.0000 & 0.0000\\
+   11 & 0.2500 & 0.2500 & 0.0000\\
+   12 & 0.2500 & 0.5000 & 0.0000\\
+   13 & 0.2500 & 0.6250 & 0.0000\\
+   14 & 0.2500 & 0.7500 & 0.0000\\
+   15 & 0.2500 & 1.0000 & 0.0000\\
+   16 & 0.5000 & 0.0000 & 0.0000\\
+   17 & 0.5000 & 0.2500 & 0.0000\\
+   18 & 0.5000 & 0.5000 & 0.0000\\
+   19 & 0.5000 & 0.7500 & 0.0000\\
+   20 & 0.5000 & 1.0000 & 0.0000\\
+   21 & 1.0000 & 0.0000 & 0.0000\\
+   22 & 1.0000 & 0.5000 & 0.0000\\
+   23 & 1.0000 & 0.7500 & 0.0000\\
+   24 & 1.0000 & 1.0000 & 0.0000
+\end{tabular}
diff --git a/doc/fig/exmpl3_ele.tex b/doc/fig/exmpl3_ele.tex
new file mode 100644 (file)
index 0000000..e1f4e9d
--- /dev/null
@@ -0,0 +1,19 @@
+\begin{tabular}{>{\columncolor{gray}}rcccc}
+\rowcolor{gray}
+ Index & c1 & c2 & c3 & c4\\
+   1 & 5 & 14 & 15 & 6\\
+   2 & 19 & 23 & 24 & 20\\
+   3 & 16 & 21 & 22 & 18\\
+   4 & 2 & 11 & 12 & 3\\
+   5 & 4 & 8 & 9 & 5\\
+   6 & 18 & 22 & 23 & 19\\
+   7 & 11 & 17 & 18 & 12\\
+   8 & 10 & 16 & 17 & 11\\
+   9 & 1 & 10 & 11 & 2\\
+   10 & 14 & 19 & 20 & 15\\
+   11 & 12 & 18 & 19 & 14\\
+   12 & 3 & 7 & 8 & 4\\
+   13 & 8 & 13 & 14 & 9\\
+   14 & 7 & 12 & 13 & 8
+\end{tabular}
+\label{exmpl3:ele}
\ No newline at end of file
diff --git a/doc/fig/exmpl3_f2s.tex b/doc/fig/exmpl3_f2s.tex
new file mode 100644 (file)
index 0000000..3a0d666
--- /dev/null
@@ -0,0 +1,16 @@
+\begin{tabular}{>{\columncolor{gray}}rcccc}
+\rowcolor{gray}
+ Index & e1 & e2 & e3 & e4\\
+   1 & 1 & 1 & 1 & 1\\
+   2 & 2 & 2 & 2 & 2\\
+   3 & 3 & 3 & 3 & 3\\
+   4 & 4 & 4 & 4 & 4\\
+   5 & 12 & 14 & 13 & 5\\
+   6 & 6 & 6 & 6 & 6\\
+   7 & 7 & 7 & 7 & 7\\
+   8 & 8 & 8 & 8 & 8\\
+   9 & 9 & 9 & 9 & 9\\
+   10 & 10 & 10 & 10 & 10\\
+   11 & 11 & 11 & 11 & 11
+\end{tabular}
+\label{exmpl3:f2s}
\ No newline at end of file
diff --git a/doc/fig/exmpl3_nei.tex b/doc/fig/exmpl3_nei.tex
new file mode 100644 (file)
index 0000000..483a265
--- /dev/null
@@ -0,0 +1,19 @@
+\begin{tabular}{>{\columncolor{gray}}rcccccccc}
+\rowcolor{gray}
+ Index & n1 & n2 & n3 & n4 & n5 & n6 & n7 & n8\\
+   1 & 13 & 10 & 0 & 0 & 5 & 0 & 0 & 0\\
+   2 & 6 & 0 & 0 & 10 & 0 & 0 & 0 & 0\\
+   3 & 0 & 0 & 6 & 8 & 0 & 0 & 0 & 7\\
+   4 & 9 & 7 & 12 & 0 & 0 & 0 & 14 & 0\\
+   5 & 12 & 13 & 1 & 0 & 0 & 0 & 0 & 0\\
+   6 & 3 & 0 & 2 & 11 & 0 & 0 & 0 & 0\\
+   7 & 8 & 3 & 11 & 4 & 0 & 0 & 0 & 0\\
+   8 & 0 & 3 & 7 & 9 & 0 & 0 & 0 & 0\\
+   9 & 0 & 8 & 4 & 0 & 0 & 0 & 0 & 0\\
+   10 & 11 & 2 & 0 & 1 & 0 & 0 & 0 & 0\\
+   11 & 7 & 6 & 10 & 14 & 0 & 0 & 0 & 13\\
+   12 & 4 & 14 & 5 & 0 & 0 & 0 & 0 & 0\\
+   13 & 14 & 11 & 1 & 5 & 0 & 0 & 0 & 0\\
+   14 & 4 & 11 & 13 & 12 & 0 & 0 & 0 & 0
+\end{tabular}
+\label{exmpl3:nei}
\ No newline at end of file
diff --git a/doc/fig/exmpl3_nei_part.tex b/doc/fig/exmpl3_nei_part.tex
new file mode 100644 (file)
index 0000000..a1cc1a3
--- /dev/null
@@ -0,0 +1,6 @@
+\begin{tabular}{>{\columncolor{gray}}rcccccccc}
+\rowcolor{gray}
+ Index & n1 & n2 & n3 & n4 & n5 & n6 & n7 & n8\\
+   4 & 9 & 7 & 12 & 0 & 0 & 0 & 14 & 0
+\end{tabular}
+\label{exmpl3:nei}
diff --git a/doc/fig/exmpl3_ref.eps b/doc/fig/exmpl3_ref.eps
new file mode 100644 (file)
index 0000000..fc682fd
--- /dev/null
@@ -0,0 +1,332 @@
+%!PS-Adobe-2.0 EPSF-1.2
+%%Creator: MATLAB, The MathWorks, Inc. Version 7.13.0.564 (R2011b). Operating System: Linux 3.0.0-16-generic #28-Ubuntu SMP Fri Jan 27 17:44:39 UTC 2012 x86_64.
+%%Title: ../doc/fig/exmpl3_ref.eps
+%%CreationDate: 03/13/2012  13:17:56
+%%DocumentNeededFonts: Helvetica
+%%DocumentProcessColors: Cyan Magenta Yellow Black
+%%Extensions: CMYK
+%%Pages: 1
+%%BoundingBox:    58   199   546   589
+%%EndComments
+
+%%BeginProlog
+% MathWorks dictionary
+/MathWorks 160 dict begin
+% definition operators
+/bdef {bind def} bind def
+/ldef {load def} bind def
+/xdef {exch def} bdef
+/xstore {exch store} bdef
+% operator abbreviations
+/c  /clip ldef
+/cc /concat ldef
+/cp /closepath ldef
+/gr /grestore ldef
+/gs /gsave ldef
+/mt /moveto ldef
+/np /newpath ldef
+/cm /currentmatrix ldef
+/sm /setmatrix ldef
+/rm /rmoveto ldef
+/rl /rlineto ldef
+/s {show newpath} bdef
+/sc {setcmykcolor} bdef
+/sr /setrgbcolor ldef
+/sg /setgray ldef
+/w /setlinewidth ldef
+/j /setlinejoin ldef
+/cap /setlinecap ldef
+/rc {rectclip} bdef
+/rf {rectfill} bdef
+% page state control
+/pgsv () def
+/bpage {/pgsv save def} bdef
+/epage {pgsv restore} bdef
+/bplot /gsave ldef
+/eplot {stroke grestore} bdef
+% orientation switch
+/portraitMode 0 def /landscapeMode 1 def /rotateMode 2 def
+% coordinate system mappings
+/dpi2point 0 def
+% font control
+/FontSize 0 def
+/FMS {/FontSize xstore findfont [FontSize 0 0 FontSize neg 0 0]
+  makefont setfont} bdef
+/reencode {exch dup where {pop load} {pop StandardEncoding} ifelse
+  exch dup 3 1 roll findfont dup length dict begin
+  { 1 index /FID ne {def}{pop pop} ifelse } forall
+  /Encoding exch def currentdict end definefont pop} bdef
+/isroman {findfont /CharStrings get /Agrave known} bdef
+/FMSR {3 1 roll 1 index dup isroman {reencode} {pop pop} ifelse
+  exch FMS} bdef
+/csm {1 dpi2point div -1 dpi2point div scale neg translate
+ dup landscapeMode eq {pop -90 rotate}
+  {rotateMode eq {90 rotate} if} ifelse} bdef
+% line types: solid, dotted, dashed, dotdash
+/SO { [] 0 setdash } bdef
+/DO { [.5 dpi2point mul 4 dpi2point mul] 0 setdash } bdef
+/DA { [6 dpi2point mul] 0 setdash } bdef
+/DD { [.5 dpi2point mul 4 dpi2point mul 6 dpi2point mul 4
+  dpi2point mul] 0 setdash } bdef
+% macros for lines and objects
+/L {lineto stroke} bdef
+/MP {3 1 roll moveto 1 sub {rlineto} repeat} bdef
+/AP {{rlineto} repeat} bdef
+/PDlw -1 def
+/W {/PDlw currentlinewidth def setlinewidth} def
+/PP {closepath eofill} bdef
+/DP {closepath stroke} bdef
+/MR {4 -2 roll moveto dup  0 exch rlineto exch 0 rlineto
+  neg 0 exch rlineto closepath} bdef
+/FR {MR stroke} bdef
+/PR {MR fill} bdef
+/L1i {{currentfile picstr readhexstring pop} image} bdef
+/tMatrix matrix def
+/MakeOval {newpath tMatrix currentmatrix pop translate scale
+0 0 1 0 360 arc tMatrix setmatrix} bdef
+/FO {MakeOval stroke} bdef
+/PO {MakeOval fill} bdef
+/PD {currentlinewidth 2 div 0 360 arc fill
+   PDlw -1 eq not {PDlw w /PDlw -1 def} if} def
+/FA {newpath tMatrix currentmatrix pop translate scale
+  0 0 1 5 -2 roll arc tMatrix setmatrix stroke} bdef
+/PA {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
+  0 0 1 5 -2 roll arc closepath tMatrix setmatrix fill} bdef
+/FAn {newpath tMatrix currentmatrix pop translate scale
+  0 0 1 5 -2 roll arcn tMatrix setmatrix stroke} bdef
+/PAn {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
+  0 0 1 5 -2 roll arcn closepath tMatrix setmatrix fill} bdef
+/vradius 0 def /hradius 0 def /lry 0 def
+/lrx 0 def /uly 0 def /ulx 0 def /rad 0 def
+/MRR {/vradius xdef /hradius xdef /lry xdef /lrx xdef /uly xdef
+  /ulx xdef newpath tMatrix currentmatrix pop ulx hradius add uly
+  vradius add translate hradius vradius scale 0 0 1 180 270 arc 
+  tMatrix setmatrix lrx hradius sub uly vradius add translate
+  hradius vradius scale 0 0 1 270 360 arc tMatrix setmatrix
+  lrx hradius sub lry vradius sub translate hradius vradius scale
+  0 0 1 0 90 arc tMatrix setmatrix ulx hradius add lry vradius sub
+  translate hradius vradius scale 0 0 1 90 180 arc tMatrix setmatrix
+  closepath} bdef
+/FRR {MRR stroke } bdef
+/PRR {MRR fill } bdef
+/MlrRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lry uly sub 2 div def
+  newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
+  rad rad scale 0 0 1 90 270 arc tMatrix setmatrix lrx rad sub lry rad
+  sub translate rad rad scale 0 0 1 270 90 arc tMatrix setmatrix
+  closepath} bdef
+/FlrRR {MlrRR stroke } bdef
+/PlrRR {MlrRR fill } bdef
+/MtbRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lrx ulx sub 2 div def
+  newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
+  rad rad scale 0 0 1 180 360 arc tMatrix setmatrix lrx rad sub lry rad
+  sub translate rad rad scale 0 0 1 0 180 arc tMatrix setmatrix
+  closepath} bdef
+/FtbRR {MtbRR stroke } bdef
+/PtbRR {MtbRR fill } bdef
+/stri 6 array def /dtri 6 array def
+/smat 6 array def /dmat 6 array def
+/tmat1 6 array def /tmat2 6 array def /dif 3 array def
+/asub {/ind2 exch def /ind1 exch def dup dup
+  ind1 get exch ind2 get sub exch } bdef
+/tri_to_matrix {
+  2 0 asub 3 1 asub 4 0 asub 5 1 asub
+  dup 0 get exch 1 get 7 -1 roll astore } bdef
+/compute_transform {
+  dmat dtri tri_to_matrix tmat1 invertmatrix 
+  smat stri tri_to_matrix tmat2 concatmatrix } bdef
+/ds {stri astore pop} bdef
+/dt {dtri astore pop} bdef
+/db {2 copy /cols xdef /rows xdef mul dup 3 mul string
+  currentfile exch readhexstring pop
+  dup 0 3 index getinterval /rbmap xdef
+  dup 2 index dup getinterval /gbmap xdef
+  1 index dup 2 mul exch getinterval /bbmap xdef pop pop}bdef
+/it {gs np dtri aload pop moveto lineto lineto cp c
+  cols rows 8 compute_transform 
+  rbmap gbmap bbmap true 3 colorimage gr}bdef
+/il {newpath moveto lineto stroke}bdef
+currentdict end def
+%%EndProlog
+
+%%BeginSetup
+MathWorks begin
+
+0 cap
+
+end
+%%EndSetup
+
+%%Page: 1 1
+%%BeginPageSetup
+%%PageBoundingBox:    58   199   546   589
+MathWorks begin
+bpage
+%%EndPageSetup
+
+%%BeginObject: obj1
+bplot
+
+/dpi2point 8.33333 def
+portraitMode 0150 5100 csm
+
+  334   190  4067  3251 MR c np
+85 dict begin %Colortable dictionary
+/c0 { 0.000000 0.000000 0.000000 sr} bdef
+/c1 { 1.000000 1.000000 1.000000 sr} bdef
+/c2 { 0.900000 0.000000 0.000000 sr} bdef
+/c3 { 0.000000 0.820000 0.000000 sr} bdef
+/c4 { 0.000000 0.000000 0.800000 sr} bdef
+/c5 { 0.910000 0.820000 0.320000 sr} bdef
+/c6 { 1.000000 0.260000 0.820000 sr} bdef
+/c7 { 0.000000 0.820000 0.820000 sr} bdef
+c0
+1 j
+1 sg
+   0    0 4801 3602 PR
+4.16667 w
+0 2935 3720 0 0 -2935 624 3205 4 MP
+PP
+-3720 0 0 2935 3720 0 0 -2935 624 3205 5 MP stroke
+2.77778 w
+DO
+SO
+4.16667 w
+0 sg
+ 624 3205 mt 4344 3205 L
+ 624 3205 mt  624  270 L
+ 624 3205 mt  624 3167 L
+%%IncludeResource: font Helvetica
+/Helvetica /ISOLatin1Encoding 83.3333 FMSR
+
+ 601 3306 mt 
+(0) s
+ 996 3205 mt  996 3167 L
+ 939 3306 mt 
+(0.1) s
+1368 3205 mt 1368 3167 L
+1311 3306 mt 
+(0.2) s
+1740 3205 mt 1740 3167 L
+1683 3306 mt 
+(0.3) s
+2112 3205 mt 2112 3167 L
+2055 3306 mt 
+(0.4) s
+2484 3205 mt 2484 3167 L
+2427 3306 mt 
+(0.5) s
+2856 3205 mt 2856 3167 L
+2799 3306 mt 
+(0.6) s
+3228 3205 mt 3228 3167 L
+3171 3306 mt 
+(0.7) s
+3600 3205 mt 3600 3167 L
+3543 3306 mt 
+(0.8) s
+3972 3205 mt 3972 3167 L
+3915 3306 mt 
+(0.9) s
+4344 3205 mt 4344 3167 L
+4321 3306 mt 
+(1) s
+ 624 3205 mt  661 3205 L
+ 554 3235 mt 
+(0) s
+ 624 2911 mt  661 2911 L
+ 485 2941 mt 
+(0.1) s
+ 624 2618 mt  661 2618 L
+ 485 2648 mt 
+(0.2) s
+ 624 2324 mt  661 2324 L
+ 485 2354 mt 
+(0.3) s
+ 624 2031 mt  661 2031 L
+ 485 2061 mt 
+(0.4) s
+ 624 1737 mt  661 1737 L
+ 485 1767 mt 
+(0.5) s
+ 624 1443 mt  661 1443 L
+ 485 1473 mt 
+(0.6) s
+ 624 1150 mt  661 1150 L
+ 485 1180 mt 
+(0.7) s
+ 624  857 mt  661  857 L
+ 485  887 mt 
+(0.8) s
+ 624  563 mt  661  563 L
+ 485  593 mt 
+(0.9) s
+ 624  270 mt  661  270 L
+ 554  300 mt 
+(1) s
+gs 624 270 3721 2936 MR c np
+/c8 { 0.000000 0.000000 1.000000 sr} bdef
+c8
+0 733 -930 0 0 -733 930 0 624 1003 5 MP stroke
+0 733 -1860 0 0 -733 1860 0 2484 1003 5 MP stroke
+0 1468 -1860 0 0 -1468 1860 0 2484 3205 5 MP stroke
+0 734 -930 0 0 -734 930 0 624 2471 5 MP stroke
+0 367 -465 0 0 -367 465 0 624 1370 5 MP stroke
+0 734 -1860 0 0 -734 1860 0 2484 1737 5 MP stroke
+0 734 -930 0 0 -734 930 0 1554 2471 5 MP stroke
+0 734 -930 0 0 -734 930 0 1554 3205 5 MP stroke
+0 734 -930 0 0 -734 930 0 624 3205 5 MP stroke
+0 733 -930 0 0 -733 930 0 1554 1003 5 MP stroke
+0 734 -930 0 0 -734 930 0 1554 1737 5 MP stroke
+0 367 -465 0 0 -367 465 0 624 1737 5 MP stroke
+0 367 -465 0 0 -367 465 0 1089 1370 5 MP stroke
+0 367 -465 0 0 -367 465 0 1089 1737 5 MP stroke
+gr
+
+c8
+0 sg
+1089  666 mt 
+(1) s
+3414  666 mt 
+(2) s
+3414 2501 mt 
+(3) s
+1089 2134 mt 
+(4) s
+ 856 1217 mt 
+(5) s
+3414 1400 mt 
+(6) s
+2019 2134 mt 
+(7) s
+2019 2868 mt 
+(8) s
+1089 2868 mt 
+(9) s
+2019  666 mt 
+(10) s
+2019 1400 mt 
+(11) s
+ 856 1584 mt 
+(12) s
+1321 1217 mt 
+(13) s
+1321 1584 mt 
+(14) s
+2462 3406 mt 
+(x) s
+ 428 1759 mt  -90 rotate
+(y) s
+90 rotate
+
+end %%Color Dict
+
+eplot
+%%EndObject
+
+epage
+end
+
+showpage
+
+%%Trailer
+%%EOF
diff --git a/doc/fig/exmpl3_ref.pdf b/doc/fig/exmpl3_ref.pdf
new file mode 100644 (file)
index 0000000..27a0c96
Binary files /dev/null and b/doc/fig/exmpl3_ref.pdf differ