Makefile: Das Bauen der Bilder korrigiert / vereinfacht.
[diplomarbeit.git] / diplomarbeit.tex
1 \documentclass[a4paper,11pt]{article}
2 \usepackage[utf8]{inputenc}
3 \usepackage{ngerman}
4 \usepackage{fancyhdr}
5 \usepackage{geometry}
6 \usepackage{amsmath, bbm}
7 \usepackage{amsfonts}
8 \usepackage{amssymb}
9 \usepackage{listings}
10 \usepackage{graphicx}
11 \usepackage{url}
12 \usepackage[table]{xcolor}
13 %\usepackage{longtable}
14 \usepackage{subfigure}
15 \usepackage{icomma}
16
17 \usepackage{tikz}
18 \usetikzlibrary{arrows,shapes}
19
20 % Fuer mathtoolsset
21 \usepackage{mathtools}
22
23 \geometry{paper=a4paper,margin=30mm}
24
25 \pagestyle{fancy}
26 %\fancyhf{}
27 %\fancyhead[LO,LE]{"Ubung zu Computational Intelligence}
28 %\fancyhead[CO,CE]{2006-05-15}
29 %\fancyhead[RO,RE]{Florian Forster (2099894)}
30
31 \title{Evolutionäre Optimierung von Sortiernetzwerken}
32 \author{Florian Forster}
33 \date{\today}
34
35 \newcommand{\false}{\textsc{False}}
36 \newcommand{\true}{\textsc{True}}
37 \newcommand{\todo}[1]{{\bf TODO:} #1}
38 \newcommand{\qed}{\hfill $\Box$ \par \bigskip}
39
40 \newcommand{\oes}[1]{\ensuremath{\operatorname{OES}\left(#1\right)}}
41 \newcommand{\bs}[1]{\ensuremath{\operatorname{BS}\left(#1\right)}}
42 \newcommand{\ps}[1]{\ensuremath{\operatorname{PS}\left(#1\right)}}
43 \newcommand{\oem}[1]{\ensuremath{\operatorname{OEM}\left(#1\right)}}
44 \newcommand{\bm}[1]{\ensuremath{\operatorname{BM}\left(#1\right)}}
45 \newcommand{\oet}[1]{\ensuremath{\operatorname{OET}\left(#1\right)}}
46
47 \newcommand{\gcell}{\cellcolor{green!10}}
48 \newcommand{\Gcell}{\cellcolor{green!10!white!95!black}}
49
50 \newtheorem{definition}{Definition}
51 \newtheorem{satz}{Satz}
52
53 % Zeige Nummern nur bei referenzierten Gleichungen an.
54 \mathtoolsset{showonlyrefs=true}
55
56 \begin{document}
57
58 \tikzstyle{vertex}   = [circle,draw,thick,fill=black,minimum size=5,inner sep=0pt]
59 \tikzstyle{comp}     = [draw,thick,-]
60 \tikzstyle{compup}   = [draw,thick,->]
61 \tikzstyle{compdown} = [draw,thick,<-]
62 \tikzstyle{edge}     = [draw,thick,-]
63 \tikzstyle{diredge}  = [draw,thick,->]
64 \tikzstyle{prob}     = [font=\tiny]
65
66 \tikzstyle{edge minimum} = [edge,color=blue!20]
67 \tikzstyle{edge maximum} = [edge,color=red!20]
68 \tikzstyle{vertex active minimum} = [vertex,color=blue!50, fill=blue!50]
69 \tikzstyle{vertex active maximum} = [vertex,color=red!50, fill=red!50]
70 \tikzstyle{vertex active minimum maximum} = [vertex,color=violet!50, fill=violet!50]
71 \tikzstyle{vertex inactive minimum} = [vertex,color=blue!20, fill=blue!20]
72 \tikzstyle{vertex inactive maximum} = [vertex,color=red!20, fill=red!20]
73 \tikzstyle{vertex inactive minimum maximum} = [vertex,color=black!20, fill=black!20]
74 \tikzstyle{comp active minimum} = [comp]
75 \tikzstyle{comp active maximum} = [comp]
76 \tikzstyle{comp active minimum maximum} = [comp,color=black!20]
77 \tikzstyle{comp inactive minimum} = [comp,color=blue!20]
78 \tikzstyle{comp inactive maximum} = [comp,color=red!20]
79 \tikzstyle{comp inactive minimum maximum} = [comp,color=black!20]
80
81 \tikzstyle{red box}   = [draw,-,color=red, top color=red!2,bottom color=red!10]
82 \tikzstyle{blue box}  = [draw,-,color=blue,top color=blue!2,bottom color=blue!10]
83 \tikzstyle{green box} = [draw,-,color=teal,top color=teal!2,bottom color=teal!10]
84 \tikzstyle{gray box}  = [draw,-,color=black, top color=black!2,bottom color=black!10]
85
86 \maketitle
87 \begin{abstract}
88 Sortiernetzwerke erweisen sich als sehr schwieriges Optimierungsproblem. Zwar
89 ist das Konzept leicht verständlich und schnell erklärt, effiziente und
90 schnelle Sortiernetzwerke zu finden oder zu konstruieren bleibt aber eine
91 Herausforderung.
92
93 Diese Arbeit verwendet „Schnitte“ oder „Leitungselimination“ und
94 Mischer-Netz\-werke, um evolutionäre Algorithmen anzugeben, deren Individuen
95 die Menge der gültigen Sortiernetzwerke nie verlassen. Bisherige Ansätze
96 bewegten sich in der Regel in der Menge aller Komparatornetzwerke und suchten
97 dort nach Sortiernetzwerken. Nach dem Vorstellen der zwei Algorithmen
98 \textsc{SN-Evolution} und \textsc{SN-Evolution-Cut} sowie einiger Ergebnisse,
99 die mit diesen Algorithmen erzielt werden konnten, wird -- basierend auf dem
100 evolutionären Algorithmus \textsc{SN-Evolution} -- eine Markov-Kette für
101 Sortiernetzwerke angegeben.
102 \end{abstract}
103 \newpage
104
105 \section*{Eidesstattliche Erklärung}
106
107 Ich versichere, dass ich die vorliegende wissenschaftliche Arbeit
108 selbstständig verfasst und keine anderen als die angegebenen Hilfsmittel
109 verwendet habe. Die Stellen der Arbeit, die anderen Werken dem Wortlaut oder
110 dem Sinn nach entnommen sind, wurden unter Angabe der Quelle als Entlehnung
111 deutlich gemacht. Das Gleiche gilt auch für beigegebene Skizzen und
112 Darstellungen. Diese Arbeit hat in gleicher oder ähnlicher Form meines Wissene
113 nach noch keiner Prüfungsbehörde vorgelegen.
114 \\[3cm]
115 München, den 20.~März 2011,\\
116 \\[1cm]
117 Florian Forster
118 \newpage
119
120 \tableofcontents
121
122 \newpage
123 \section{Motivation und Einleitung}
124
125 \subsection{Motivation}\label{sect:motivation}
126
127 \emph{Sortiernetzwerke} sind ein theoretisches Konstrukt, dass auch von
128 Personen ohne Zugang zum Thema, beziehungsweise der theoretischen Informatik,
129 schnell verstanden werden kann. Eine Einführung wird in
130 Abschnitt~\ref{sect:einleitung_sortiernetzwerke} gegeben. Nichtsdestotrotz ist
131 das Finden von Sortiernetzwerken sowie das Beweisen, dass ein
132 Komparatornetzwerk jede beliebige Eingabe sortiert, im Allgemeinen sehr
133 schwierig und nach heutigem Kenntnisstand nur mit exponentiellem Aufwand zu
134 bewältigen.
135
136 Einfacher ist der Korrektheitsbeweis bei konstruktiven Verfahren, da hier die
137 Konstruktionsvorschrift genutzt werden kann, um die Korrektheit für beliebige
138 Eingabegrößen zu beweisen. Im Abschnitt~\ref{sect:konstruktive_netzwerke}
139 geschieht dies beispielsweise für zwei von \emph{Kenneth~E. Batcher} 1968
140 gefundene Konstruktionsvorschriften.
141
142 Um effiziente und schnelle Sortiernetzwerke zu finden, wurden schon wiederholt
143 Computer und automatische Suchverfahren eingesetzt. Bisherige Ansätze
144 versuchen meist in der Menge aller Komparatornetzwerke jene zu finden, die
145 die Sortiereigenschaft besitzen und aus wenigen Komparatoren bestehen. Die
146 Eigenschaft, jede Eingabepermutation zu sortieren, ist also ein
147 Optimierungsziel und nicht durch das Vorgehen gewährleistet. Dafür können
148 theoretisch alle Sortiernetzwerke durch diese Algorithmen gefunden werden --
149 genügend Laufzeit vorausgesetzt.
150
151 In dieser Arbeit werden Methoden verwendet, die die Menge der Sortiernetzwerke
152 nie verlassen, dafür aber auch nicht alle existierenden Sortiernetzwerke
153 erzeugen können. So muss für ein gefundenes Komparatornetzwerk nicht mehr
154 nachgewiesen werden, dass es jede beliebige Eingabe sortiert, weil diese
155 Eigenschaft durch das Verfahren sichergestellt ist.
156
157 \subsection{Einleitung}\label{sect:einleitung}
158
159 \subsubsection{Sortiernetzwerke}\label{sect:einleitung_sortiernetzwerke}
160
161 \emph{Komparatoren} sind die Bausteine, die \emph{Komparatornetzwerken}
162 zugrunde liegen. Sie haben zwei Eingänge über die sie zwei Zahlen erhalten
163 können und zwei Ausgänge, auf denen die Zahlen wieder ausgegeben werden. Dabei
164 sind die Ausgänge im Gegensatz zu den Eingängen unterscheidbar, da die größere
165 der beiden Zahlen immer auf dem einen, die kleinere der beiden Zahlen
166 immer auf dem anderen Ausgang ausgegeben wird.
167
168 Kombiniert man mehrere \emph{Komparatoren} in der Form miteinander, dass die
169 Ausgänge eines Komparators mit Eingängen weiterer Komparatoren verbunden sind,
170 erhält man ein {\em Komparatornetzwerk}.
171
172 \begin{figure}
173   \begin{center}
174     \input{images/einfaches_komparatornetzwerk.tex}
175   \end{center}
176   \caption{Einfaches Komparatornetzwerk mit 4~Ein- beziehungsweise Ausgängen,
177     bestehend aus 5~Komparatoren.}
178   \label{fig:einfaches_komparatornetzwerk}
179 \end{figure}
180
181 Abbildung~\ref{fig:einfaches_komparatornetzwerk} zeigt ein einfaches
182 \emph{Komparatornetzwerk} aus fünf Komparatoren. Insgesamt gibt es vier
183 verschiedene Eingänge und vier Ausgänge. Die Ein- und Ausgänge werden durch
184 eine horizontale Linie dargestellt und als \emph{Leitung} bezeichnet. Die
185 \emph{Komparatoren} sind durch vertikale Pfeile dargestellt und verbinden je
186 zwei verschiedene \emph{Leitungen} miteinander. Die Verbindungsstellen von
187 \emph{Leitungen} und \emph{Komparatoren} sind zur besseren Übersichtlichkeit
188 durch schwarze Punkte symbolisiert.
189
190 Auf der linken Seite befinden sich die Eingänge. Hier wird eine Zahlenfolge in
191 das Netzwerk hinein gegeben. Jeder Komparator vergleicht die Zahlen „auf“ den
192 beiden Leitungen, die er verbindet. Nach einem Komparator befindet sich die
193 kleinere Zahl immer auf der Leitung, auf die der Pfeil zeigt, die größere Zahl
194 befindet sich auf der Leitung, auf der der Pfeil seinen Ursprung hat. Wenn in
195 einem Komparatornetzwerk alle Komparatoren in die gleiche Richtung zeigen --
196 die Konvention in dieser Arbeit ist, dass das Minimum auf der unteren Leitung
197 ausgegeben wird -- werden die Pfeile durch einfache Linien ersetzt. Zu diesen
198 sogenannten \emph{Standard-Netzwerken} siehe auch
199 Abschnitt~\ref{sect:normalisieren}.
200
201 Komparatoren, die \emph{unterschiedliche} Leitungen miteinander vergleichen,
202 können gleichzeitig angewendet werden. Das Beispiel in
203 Abbildung~\ref{fig:einfaches_komparatornetzwerk} verwendet diesen Umstand und
204 vergleicht die zwei oberen und die zwei unteren Leitungen gleichzeitig. Eine
205 Gruppe von Komparatoren, die gleichzeitig angewendet werden können, nennt man
206 eine \emph{Schicht} des Komparatornetzwerks. Die \emph{Geschwindigkeit} eines
207 Komparatornetzwerks ist gleichbedeutend mit der Anzahl der Schichten, in die
208 sich die Komparatoren mindestens gruppieren lassen, da sie die Anzahl der
209 benötigten parallelen Schritte darstellt.
210
211 \emph{Komparatornetzwerke}, die für \emph{jede} Eingabefolge eine Ausgabe
212 erzeugen, die der Sortierung der Eingabe entspricht, heißen
213 \emph{Sortiernetzwerke}. Das in
214 Abbildung~\ref{fig:einfaches_komparatornetzwerk} gezeigte Komparatornetzwerk
215 ist \emph{kein} Sortiernetzwerk: Die Eingabefolge ${(1, 2, 3, 4)}$ führt zur
216 Ausgabe ${(2, 1, 3, 4)}$ -- die bestehenden Sortierung wird also sogar
217 zerstört.
218
219 \begin{figure}
220   \begin{center}
221     \input{images/09-e2-c24-allbut1.tex}
222   \end{center}
223   \caption{Ein \emph{Komparatornetzwerk} mit 9~Eingängen und 24~Komparatoren,
224   die in 8~Schichten angeordnet sind. Das Netzwerk sortiert alle Eingaben, bei
225   denen das Minimum nicht auf dem mittleren Eingang liegt.}
226   \label{fig:09-e2-c24-allbut1}
227 \end{figure}
228 Zu beweisen, dass ein gegebenes Komparatornetzwerk die Sortiereigenschaft {\em
229 nicht} hat, ist mit einem gegebenen Gegenbeispiel einfach möglich. Das
230 Komparatornetzwerk wird auf das Gegenbeispiel angewendet und anschließend wird
231 überprüft, ob die Ausgabe sortiert ist. Ist sie es nicht, heißt das, dass es
232 mindestens eine Eingabefolge gibt, die nicht sortiert wird. Entsprechend der
233 Definition handelt es sich bei dem \emph{Komparatornetzwerk} folglich
234 \emph{nicht} um ein \emph{Sortiernetzwerk}. Ein solches Gegenbeispiel für ein
235 gegebenes Komparatornetzwerk zu finden ist nach heutigem Kenntnisstand jedoch
236 nicht \emph{effizient}\footnote{In diesem Zusammenhang heißt \emph{effizient},
237 dass keine Algorithmen bekannt sind, die eine polynomielle Laufzeit (in
238 Abhängigkeit von der Eingabelänge) haben.} möglich.
239
240 Beispielsweise sortiert das im Rahmen dieser Arbeit entdeckte
241 Komparatornetzwerk in Abbildung~\ref{fig:09-e2-c24-allbut1} viele der 362.880
242 möglichen Eingabepermutationen. Mit dem Gegenbeispiel $(3, 5, 2, 1, 0, 7, 4,
243 8, 6)$ lässt sich jedoch leicht beweisen, dass das Komparatornetzwerk die
244 Sortiereigenschaft \emph{nicht} besitzt, da es in diesem Fall die Folge $(1,
245 0, 2, 3, 4, 5, 6, 7, 8)$ ausgibt.
246
247 Insgesamt gibt es $n!$~Permutationen von $n$~Elementen. Wenn ein
248 Komparatornetzwerk die Sortiereigenschaft besitzt, bildet es alle diese
249 Permutationen auf die sortierte Reihenfolge ab. Allerdings wächst $n!$
250 so schnell, dass ein Ausprobieren aller möglichen Permutationen schon bei
251 16~Leitungen praktisch nicht mehr zu bewerkstelligen
252 ist.\footnote{1.307.674.368.000 Permutationen}
253
254 \label{sect:0-1-prinzip}
255 Glücklicherweise reicht es aus, alle möglichen 0-1-Folgen zu überprüfen, wie
256 \textit{Donald~E. Knuth} in \cite{KNUTH} zeigt. Die Beweisidee ist folgende:
257 Angenommen ein Komparatornetzwerk sortiert alle 0-1-Folgen und es gibt eine
258 Permutation $E = (e_0, \dots, e_{n-1})$ beliebiger Zahlen, die nicht sortiert
259 wird. Die Ausgabefolge sei $A = (a_0, \dots, a_{n-1})$. Sei $i$ eine Position
260 in der Ausgabe, die die Sortierbedingung verletzt:
261 \begin{displaymath}
262   a_0 \leqq a_1 \leqq \dots \leqq a_{i-1} > a_i \dots
263 \end{displaymath}
264 Die Eingabe kann mittels
265 \begin{displaymath}
266   \hat{e}_j = \left\{
267     \begin{array}{cl}
268       0 & e_j \leqq a_i \\
269       1 & e_j > a_i
270     \end{array} \right.
271 \end{displaymath}
272 auf eine 0-1-Folge abgebildet werden, die entsprechend der Annahme vom
273 Komparatornetzwerk sortiert wird. Allerdings verändert diese Abbildung das
274 Verhalten jedes einzelnen Komparators nicht: Wenn bei der Permutation eine
275 Zahl größer als $a_i$ und eine Zahl kleiner oder gleich $a_i$ verglichen
276 wurden, liegen jetzt entsprechend eine Null und eine Eins an, die genauso
277 vertauscht werden oder nicht, wie das bei der Permutation der Fall war. Liegen
278 zwei Nullen oder zwei Einsen an, entsprechen sie zwei Zahlen kleiner als
279 $a_i$, beziehungsweise zwei Zahlen größer oder gleich $a_i$. Da im Fall der
280 0-1-Folge zwei gleiche Zahlen am Komparator anliegen, dürfen wir davon
281 ausgehen, dass sich der Komparator so verhält, wie er sich bei der Permutation
282 verhalten hat -- ohne das Ergebnis zu beeinflussen. Entsprechend müssen an den
283 Ausgängen $i-1$ und $i$ eine Null und eine Eins in der falschen Reihenfolge
284 ankommen. Das steht im Widerspruch zu der Annahme, dass alle 0-1-Folgen
285 sortiert werden.
286
287 Im Gegensatz zum Überprüfen aller möglichen Permutationen, was mit dem Aufwand
288 $\Theta\left(\sqrt{n}\left(\frac{n}{e}\right)^n\right)$ verbunden ist, besitzt
289 das Überprüfen aller 0-1-Folgen „nur“ den Aufwand $\Theta(2^n)$. Entsprechend
290 ist dieses Verfahren nicht \emph{effizient} -- ein schnelleres Verfahren ist
291 bisher allerdings nicht bekannt.
292
293 Um zu überprüfen, ob ein Komparatornetzwerk mit 16~Leitungen die
294 Sortiereigenschaft besitzt, sind mit dieser Methode nur 65.536 Tests notwendig
295 -- eine Zahl, die für aktuelle Prozessoren keine Herausforderung darstellt.
296 Für die Überprüfung eines Komparatornetzwerks mit 32~Leitungen sind jedoch
297 bereits etwa 4,3~Milliarden Tests notwendig, die einen Rechner durchaus
298 mehrere Stunden beschäftigen. Das ist deshalb problematisch, weil die im
299 Folgenden vorgestellten \emph{Evolutionären Algorithmen} eine entsprechende
300 Überprüfung in jeder Iteration durchführen müssten. Wenn die Überprüfung eines
301 Zwischenergebnisses drei Stunden in Anspruch nimmt, sind für eine Million
302 Iterationen etwa 340~Jahre Rechenzeit notwendig. Selbst wenn die Berechnung
303 auf 1000~Computern mit je 4~Prozessoren verteilt wird, wird etwa ein Monat für
304 einen Lauf benötigt.
305
306 \subsubsection{Evolutionäre Algorithmen}
307
308 Viele {\em kombinatorische Optimierungsprobleme} sind schwer zu lösen -- die
309 entsprechenden Entscheidungsprobleme liegen oft in der Komplexitätsklasse
310 $\mathcal{NP}$-vollständig. Das heißt, dass keine Verfahren bekannt sind, die
311 diese Probleme effizient exakt lösen. Sollte sich herausstellen, dass diese
312 Probleme außerhalb der Komplexitätsklasse~$\mathcal{P}$ liegen, wäre eine
313 Konsequenz, dass es für diese Probleme keine effizienten exakten Algorithmen
314 gibt. Stellt sich hingegen heraus, dass diese Probleme neben
315 $\mathcal{NP}$-vollständig auch in der Komplexitätsklasse~\textit{P} liegen,
316 gibt es effiziente Algorithmen. Es ist jedoch wahrscheinlich, dass die
317 Zeitkonstanten solcher Algorithmen sehr groß wären, so dass der praktische
318 Nutzen fraglich bleibt.
319
320 Aus diesem Grund besteht die Notwendigkeit, einen Kompromiss einzugehen: Statt
321 die \emph{optimale Lösung}, beziehungsweise eine der \emph{optimalen Lösungen}
322 als einzige Ausgabe des Algorithmus zuzulassen, wird eine "`möglichst gute"'
323 Lösung ausgegeben. Dafür verringert sich die Laufzeit des Algorithmus. Viele
324 dieser Optimierungsalgorithmen orientieren sich an Vorgängen in der Natur.
325 Beispielsweise imitieren die „Ameisenalgorithmen“ das Verhalten von Ameisen
326 auf der Futtersuche, um kurze Rundreisen auf Graphen zu berechnen.
327
328 Bei {\em Evolutionären Algorithmen} stand die Evolution Pate. Die Grundidee
329 ist, bekannte Lösungen zu neuen -- unter Umständen besseren -- Lösungen zu
330 kombinieren. Dabei bedient man sich der in der Evolutionstheorie etablierten
331 Nomenklatur, beispielsweise werden konkrete Lösungen für ein Problem als {\em
332 Individuen} bezeichnet.
333
334 Die Vorgehensweise lässt sich abstrakt wie folgt beschreiben: Aus einer
335 bestehenden Lösungsmenge, der {\em Population}, werden zufällig Lösungen
336 ausgesucht {\em (Selektion)} und zu einer neuen Lösung kombiniert ({\em
337 Rekombination}). Unter Umständen wird die neue Lösung noch zufällig
338 verändert {\em (Mutation)}, bevor sie in die bestehende Lösungsmenge
339 eingefügt wird. Die verwendeten Wahrscheinlichkeiten, beispielsweise bei der
340 {\em Selektion}, sind dabei nicht zwangsläufig gleichverteilt -- üblicherweise
341 werden bessere Lösungen bevorzugt. Zur Bewertung dient die sogenannte {\em
342 Gütefunktion}.
343
344 Nicht alle Probleme eignen sich für diese Strategie. Zum einen muss es möglich
345 sein, eine initiale Population zur Verfügung zu stellen, da diese als Basis
346 aller weiteren Operationen dient. Das ist häufig keine große Einschränkung, da
347 es oft einfach ist, {\em irgendeine} Lösung anzugeben. Die angegebenen
348 Algorithmen verwenden als einfache initiale Lösung häufig das
349 \emph{Odd-Even-Transpositionsort}-Netzwerk, das in
350 Abschnitt~\ref{sect:odd_even_transpositionsort} beschrieben wird. Zum anderen
351 muss eine Methode für die Rekombination existieren. Das ist insbesondere dann
352 problematisch, wenn {\em Nebenbedingungen} eingehalten werden müssen. Die in
353 dieser Arbeit verwendeten Rekombinationsmethoden sind so gewählt, dass die
354 Nebenbedingungen nicht verletzt werden.
355
356 Beim Aussuchen von zufälligen Lösungen aus der Population, der
357 \emph{Selektion}, werden gute Lösungen bevorzugt. Wie sehr diese Lösungen
358 bevorzugt werden, hat einen starken Einfluss auf das Verhalten des
359 Algorithmus. Werden gute Lösungen stark bevorzugt, konvergiert der Algorithmus
360 schnell gegen ein (lokales) Optimum. Dieses \textit{Exploitation} (Englisch
361 für „Ausnutzung“) genannte Verhalten sorgt dafür, dass sich der Algorithmus
362 schnell auf eine Lösung festlegt und andere, möglicherweise bessere lokale
363 Optima nicht mehr findet. Werden gute Lösungen hingegen nur wenig bevorzugt,
364 erforscht der Algorithmus den Lösungsraum in viele Richtungen. Dieses
365 \textit{Exploration} (Englisch für „Erforschung“) genannte Verhalten sorgt
366 zwar dafür, dass der Algorithmus langsamer auf ein Optimum zusteuert, dafür
367 findet er aber in der Regel bessere Lösungen. Die Rolle, die
368 \textit{Exploitation} und \textit{Exploration} bei evolutionären
369 Optimierungsalgorithmen spielen, wird von \textit{Eiben} und
370 \textit{Schippers} in~\cite{ES1998} untersucht.
371
372 Die Parameter evolutionärer Algorithmen so einzustellen, dass sich ein guter
373 Mittelweg zwischen den beiden Extremen einstellt, ist eine Aufgabe, die sich
374 nur experimentell lösen lässt. Die genauen Parameter hängen nicht nur vom
375 eigentlichen Algorithmus, sondern auch vom konkreten Problem ab, so dass sich
376 beispielsweise bei der Optimierung von Sortiernetzwerken die Parameter
377 zwischen verschiedenen Leitungszahlen stark unterscheiden. Einen Überblick
378 geben \textit{Kalyanmoy Deb} und \textit{Samir Agrawal} in~\cite{DA1998}.
379
380 Die Erforschung (\textit{Exploration}) kann von einem weiteren Mechanismus
381 unterstützt werden, der ebenfalls der Evolutionslehre entliehen ist, der
382 \emph{Mutation}. Dabei werden Lösungen zufällig verändert, so dass auch andere
383 Lösungen „in der Nähe“ von direkten Nachfolgern erreicht werden können. Das
384 hilft insbesondere bei der intensiven Suche in der Nähe eines lokalen Optimums
385 aber auch beim „Ausbrechen“ aus lokalen Optima und finden noch besserer
386 Lösungen.
387
388 Bei \emph{Sortiernetzwerken} ist eine \emph{Mutation} leider immer damit
389 verbunden, dass anschließend die Sortiereigenschaft des resultierenden
390 \emph{Komparatornetzwerks} wieder überprüft werden muss, da selbst das
391 Hinzufügen eines zufälligen Komparators diese Eigenschaft zer\-stö\-ren kann.
392 Beim Suchen möglichst effizienter Netzwerke ist das zufällige Entfernen von
393 Komparatoren interessanter, was die Sortiereigenschaft fast immer aufhebt.
394
395 Die im Folgenden beschriebenen Algorithmen mutieren (verändern) daher nicht
396 die \emph{Sortiernetzwerke} selbst, sondern verzichten entweder ganz auf
397 Mutation oder mutieren lediglich Transformationen von Sortiernetzwerken, die
398 die Sortiereigenschaft erhalten. Transformationen von Sortiernetzwerken werden
399 in Abschnitt~\ref{sect:tranformation} beschrieben, ein Algorithmus, der
400 Mutation einsetzt, wird in Abschnitt~\ref{sect:sn-evolution-cut} vorgestellt.
401
402 \begin{figure}
403   \begin{center}
404     \input{images/16-hillis.tex}
405   \end{center}
406   \caption{Das 16-Sortiernetzwerk, das \textit{Hillis} in~\cite{H1990} angibt.
407   Es besteht aus 61~Komparatoren in 11~Schichten.}
408   \label{fig:16-hillis}
409 \end{figure}
410 Evolutionäre Algorithmen wurden bereits mehrfach eingesetzt, um
411 Sortiernetzwerke zu untersuchen. \textit{W.~Daniel Hillis} verwendete
412 \emph{Co-Evolution} um neben Komparatornetzwerken auch „schwierige Eingaben“
413 zu optimieren~\cite{H1990}. Diese \emph{Parasiten} genannten Eingaben wurden
414 daran gemessen, bei wie vielen Komparatornetzwerken sie beweisen konnten, dass
415 sie keine Sortiernetzwerke sind. So mussten bei neuen Individuen
416 (Komparatornetzwerken) nicht alle 0-1-Folgen, sondern nur erfolgreiche
417 Parasiten (schwierige Eingaben) überprüft werden. Auf diese Art und Weise
418 gelang es \textit{Hillis} ein 16-Sortiernetzwerk mit 61~Komparatoren
419 anzugeben, das in Abbildung~\ref{fig:16-hillis} zu sehen ist.
420
421 \textit{Michael~L. Harrison} und \textit{James~A. Foster} nutzten ebenfalls
422 \emph{Co-Evolution} in~\cite{HF2004}, um die Stabilität von Sortiernetzwerken
423 zu erhöhen. Die zweite Population bestand aus \emph{Fehlern} -- der
424 Information, welche Komparatoren defekt, beziehungsweise inaktiv sind. So
425 generierte Sortiernetzwerke können auch dann noch für viele Eingaben eine
426 korrekt sortierte Ausgabe erzeugen, wenn ein oder mehrere Komparatoren
427 (zufällig) entfernt werden.
428
429 \begin{figure}
430   \centering
431   \subfigure{\input{images/13-juille-0.tex}}
432   \subfigure{\input{images/13-juille-1.tex}}
433   \caption{13-Sortiernetzwerke, die von \textit{Hugues Juillé} mithilfe des
434   END-Algorithmus gefunden wurden. Sie bestehen jeweils aus 45~Komparatoren in
435   10~Schichten.}
436   \label{fig:13-juille}
437 \end{figure}
438 \textit{Hugues Juillé} entwickelte ein Verfahren, das er \emph{Evolving
439 Non-Determinism} (END) nannte~\cite{J1995}. Dabei handelt es sich nicht um
440 einen der \emph{Evolutionären Algorithmen}, wie sie hier vorgestellt wurden,
441 sondern um eine verteilte, probabilistische Breitensuche, die an die
442 \emph{Strahlsuche} (englisch: \textit{beam search}), ein Verfahren der
443 Künstlichen Intelligenz, angelehnt ist. Die aufwendigste Operation bei diesem
444 Ansatz ist die Bewertungsfunktion, die abschätzt, wie viele Komparatoren zu
445 einem Komparatornetzwerk hinzugefügt werden müssen, um ein Sortiernetzwerk zu
446 erhalten. Mit diesem Ansatz gelang es \textit{Juillé} zwei 13-Sortiernetzwerke
447 anzugeben, die mit 45~Komparatoren effizienter sind als alle bis dahin
448 bekannten (Abbildung~\ref{fig:13-juille}).
449
450 \newpage
451 \section[Konstruktionsverfahren]{Konstruktionsverfahren für Sortiernetzwerke}
452 \label{sect:konstruktive_netzwerke}
453
454 Die bekannten Konstruktionsverfahren für Sortiernetzwerke, insbesondere ein
455 häufig verwendeter Baustein, sogenannte \emph{Mischer}\footnote{Eine
456 Fehlübersetzung aus dem Englischen, von \textit{to~merge} (Deutsch:
457 zusammenfügen). Da der Begriff des "`mischens"' beziehungsweise der
458 "`Mischer"' in der Literatur sehr weit verbreitet ist, werden diese Begriffe
459 in dieser Arbeit trotzdem verwendet.}, bilden die Grundlage für die
460 beschriebenen evolutionären Algorithmen beziehungsweise dienen als initiale
461 Eingabe. Im Folgenden werden daher drei Konstruktionsverfahren vorgestellt.
462
463 \subsection{Das Odd-Even-Transpositionsort-Netzwerk}
464 \label{sect:odd_even_transpositionsort}
465
466 Das Sortiernetzwerk {\em Odd-Even-Transpositionsort} (OET) ist eines der
467 einfachsten Sortiernetzwerke. Es besteht aus $n$~{\em Schichten}, die jede
468 "`Leitung"' abwechselnd mit den benachbarten Leitungen verbindet.
469 Abbildung~\ref{fig:odd-even-transposition-08} zeigt das OET-Netzwerk für
470 ${n = 8}$ Leitungen.
471
472 \begin{figure}
473   \begin{center}
474     \input{images/oe-transposition-8.tex}
475   \end{center}
476   \caption{Das \emph{Odd-Even-Transpositionsort}-Netzwerk mit 8~Eingängen.}
477   \label{fig:odd-even-transposition-08}
478 \end{figure}
479
480 Dass das \emph{Odd-Even-Transpositionsort}-Netzwerk tatsächlich jede beliebige
481 Eingabe sortiert, ist nicht offensichtlich. Leicht zu sehen ist jedoch, dass
482 sowohl das Minimum als auch das Maximum durch das im Netzwerk enthaltene
483 Treppenmuster auf die unterste beziehungsweise oberste Leitung gelangt.
484
485 Die Sortiereigenschaft größerer OET-Netzwerke lässt sich rekursiv beweisen,
486 indem man $\operatorname{OET}(n)$ auf $\operatorname{OET}(n-1)$ durch
487 Herausschneiden einer Leitung reduziert. In
488 Abschnitt~\ref{sect:leitungen_entfernen} wird das Vorgehen im Detail
489 beschrieben, Abbildung~\ref{fig:oe-transposition-cut} zeigt das
490 Herausschneiden einer Leitung aus $\operatorname{OET}(8)$. Die Rekursion endet
491 beim \emph{Odd-Even-Transpositionsort}-Netzwerk mit drei Eingängen, bei dem
492 das Minimum auf der untersten, das Maximum auf der obersten und der mittlere
493 Wert auf der mittleren Leitung landet. Folglich ist die Ausgabe bei
494 $\operatorname{OET}(3)$ sortiert.
495
496 Das \emph{Odd-Even-Transpositionsort}-Netzwerk ist weder in Bezug auf die
497 Anzahl der Komparatoren noch in Bezug auf die Anzahl der Schichten, in denen
498 sich die Komparatoren anordnen lassen, effizient. Es benötigt ${\frac12 n
499 (n-1)} = \Theta(n^2)$~Komparatoren, die in $n$~Schichten angeordnet sind.
500 Die im Folgenden vorgestellten Sortiernetzwerke benötigen deutlich weniger Komparatoren,
501 ($\Theta(n \log (n)^2)$), die in weniger Schichten,
502 ($\Theta(\log (n)^2)$), angeordnet sind.
503
504 Das Interessante am OET-Netzwerk ist seine einfache Konstruktion. Einige der
505 folgenden Algorithmen benötigen ein möglichst einfaches Sortiernetzwerk als
506 Starteingabe, auf dessen Basis sie versuchen optimierte Sortiernetzwerke zu
507 finden. Häufig dient $\operatorname{OET}(n)$ als Eingabe für diese
508 Algorithmen.
509
510 Außerdem bedienen sich die Algorithmen der Technik des Herausschneidens einer,
511 beziehungsweise mehrerer Leitungen, um die Anzahl der Leitungen eines
512 Sortiernetzwerks zu reduzieren. Die Technik wird in Detail im
513 Abschnitt~\ref{sect:leitungen_entfernen} beschrieben.
514
515 \subsection{Das bitone Mergesort-Netzwerk}
516
517 Das \emph{bitone Mergesort}-Netzwerk ($\operatorname{BS}(n)$) ist ein
518 Sortiernetzwerk, das 1968 von \emph{Kenneth~E. Batcher} in~\cite{B1968}
519 veröffentlicht wurde. Es ist deutlich effizienter als das
520 Odd-Even-Transposi\-tionsort-Netzwerk -- sowohl in Bezug auf die Anzahl der
521 Komparatoren als auch der benötigten Zeit, also der Anzahl der Schichten.
522
523 Das Sortiernetzwerk basiert auf einem Komparatornetzwerk, welches zwei
524 sortierte Listen zusammenfügen (Englisch: \textit{to~merge}) kann. Dieser
525 \emph{„bitone Mischer“} (Englisch: \textit{bitonic merger}) genannte Baustein
526 verleiht dem Sortiernetzwerk seinen Namen.
527
528 Da das Sortiernetzwerk rekursiv definiert ist, betrachten wir hier nur die
529 Instanzen des Netzwerks, deren Leitungszahl $n = 2^d$ eine Zweierpotenz ist.
530 Es ist jedoch möglich, das Sortiernetzwerk für beliebige~$n$ zu erzeugen.
531
532 \subsubsection{Der bitone Mischer}\label{sect:der_bitone_mischer}
533
534 Das \emph{bitone Mergesort}-Netzwerk basiert auf dem sogenannten \emph{bitonen
535 Mischer} $\operatorname{BM}(n)$, einem Kom\-parator-Netzwerk, das eine
536 beliebige \emph{bitone Folge} in eine sortierte Liste umordnen kann. Eine
537 \emph{bitone Folge} ist eine monoton steigende Folge, gefolgt von einer
538 monoton absteigenden Folge, oder ein zyklischer Shift davon.
539 Abbildung~\ref{fig:beispiel-biton} zeigt die vier prinzipiellen Möglichkeiten,
540 die durch zyklische Shifts entstehen können. Die wichtigsten Varianten für das
541 \emph{bitone Mergesort}-Netzwerk zeigen die
542 Abbildungen~\ref{fig:beispiel-biton-0} und~\ref{fig:beispiel-biton-1}. Sie
543 erhält man, wenn man eine aufsteigend und eine absteigend sortierte Liste
544 aneinanderhängt. Bei den anderen beiden Formen ist wichtig zu beachten, dass
545 das letzte Element nicht größer (Abbildung~\ref{fig:beispiel-biton-2}) beziehungsweise
546 kleiner (Abbildung~\ref{fig:beispiel-biton-3}) als das erste Element der Folge
547 sein darf.
548
549 \begin{figure}
550   \centering
551   \subfigure[aufsteigend, absteigend]{\input{images/beispiel-biton-0.tex}\label{fig:beispiel-biton-0}}
552   \subfigure[absteigend, aufsteigend]{\input{images/beispiel-biton-1.tex}\label{fig:beispiel-biton-1}}
553   \subfigure[aufsteigend, absteigend, aufsteigend]{\input{images/beispiel-biton-2.tex}\label{fig:beispiel-biton-2}}
554   \subfigure[absteigend, aufsteigend, absteigend]{\input{images/beispiel-biton-3.tex}\label{fig:beispiel-biton-3}}
555   \caption{Beispiele bitoner Folgen.}
556   \label{fig:beispiel-biton}
557 \end{figure}
558
559 \begin{figure}
560   \centering
561   \subfigure[normal]{\input{images/bitonic-merge.tex}\label{fig:bitonic-merge-normal}}
562   \qquad
563   \subfigure[trichter]{\input{images/bitonic-merge-trichter.tex}\label{fig:bitonic-merge-tricheter}}
564   \caption{Schematischer Aufbau des bitonen Mischers: Jedes Element der
565   aufsteigenden Folge $u_0, u_1, \ldots$ wird mit dem entsprechenden Element
566   der absteigend sortierten Folge $v_0, v_1, \ldots$ verglichen. Die beiden
567   resultierenden Teilfolgen sind wiederum biton.}
568   \label{fig:bitonic-merge-schema}
569 \end{figure}
570
571 Der Mischer funktioniert folgendermaßen: Gegeben sind zwei Folgen mit je
572 ${m = \frac{n}{2}}$ Elementen, $U = \left(u_0, u_1, \ldots, u_{m-1}\right)$ und
573 $V = \left(v_0, v_1, \ldots, v_{m-1}\right)$. Die Folge $U$ sei aufsteigend
574 sortiert, die Folge $V$ sei absteigend sortiert:
575 \begin{eqnarray}
576  u_0 \leqq u_1 \leqq &\ldots& \leqq u_{m-1} \\
577  v_0 \geqq v_1 \geqq &\ldots& \geqq v_{m-1}
578 \end{eqnarray}
579 Im ersten Schritt werden nun jeweils die Elemente an den gleichen relativen
580 Positionen verglichen und ggf. vertauscht:
581 \begin{equation}
582 u_i \longleftrightarrow v_i, \quad 0 \leqq i < m
583 \end{equation}
584 Sei $j \in \{0 \ldots m\}$ der Index der ersten Elemente $u_j$ und $v_j$, die
585 durch den gemeinsamen Komparator vertauscht werden. Unter der Annahme, dass
586 Elemente nur vertauscht werden wenn, sie ungleich sind, muss ${u_j > v_j}$
587 gelten. Mit $u_j \leqq u_{j+1}$ und $v_j \geqq v_{j+1}$ folgt daraus $u_{j+1}
588 > v_{j+1}$. Es werden also alle Elemente $u_k$ und $v_k$ mit $k \geqq j$
589 vertauscht. $j = m$ bezeichnet den Fall, in dem das größte Element der
590 "`linken"' Folge $u_{m-1}$ kleiner ist als das kleinste Element der
591 "`rechten"' Folge $v_{m-1}$. Daraus folgt, dass sich das Resultat in zwei
592 bitone Folgen aufteilen lässt: Eine aufsteigende~/ absteigende Folge und eine
593 absteigende~/ aufsteigende Folge. Abbildung~\ref{fig:bitonic-merge-normal}
594 zeigt die Situationen vor und nach diesem Schritt des Mischers schematisch.
595
596 Um die Folge vollständig zu sortieren, müssen anschließend die beiden
597 resultierenden bitonen Folgen sortiert werden. Die geschieht ebenfalls
598 mithilfe des bitonen Mischers, mit zwei Instanzen von
599 $\operatorname{BM}(\frac{n}{2})$. Diese rekursive Definition endet mit dem
600 bitonen Mischer mit zwei Leitungen, $\operatorname{BM}(2)$, der als
601 Komparator-Netzwerk mit einem Komparator zwischen den beiden Leitungen
602 definiert ist.
603
604 Der bitone Mischer kann auch zwei aufsteigende Folgen sortieren. Dazu ist
605 lediglich eine etwas modifizierte Vergleichs-Kaskade im ersten Schritt
606 notwendig. Die folgenden, kleineren Mischer erhalten als Eingabe wieder eine
607 „echte“ bitone Folge. Abbildung~\ref{fig:bitonic-merge-tricheter} zeigt das
608 Schema des bitonen Mischers für zwei aufsteigend sortierte Folgen. Durch das
609 Umdrehen einer Folge verändert sich das Muster der Komparatoren ein wenig:
610 Statt an eine Treppe erinnert das Muster nun an einen Trichter.
611
612 Da sich die Anzahl der Leitungen in jedem Rekursionsschritt halbiert, endet
613 die Rekursion nach $\log(n)$~Schritten. In jedem Rekursionsschritt werden
614 $\frac{n}{2}$~Komparatoren eingefügt, so dass der gesamte Mischer aus
615 $\frac{1}{2} n \log(n) = \Theta\left(n \log(n)\right)$~Komparatoren
616 besteht, die in $\log(n)$~Schichten angeordnet werden können.
617
618 \subsubsection{Das bitone Mergesort-Netzwerk}
619
620 Ebenso wie der bitone Mischer $\operatorname{BM}(n)$ ist auch das \emph{bitone
621 Mergesort-Netzwerk} $\operatorname{BS}(n)$ rekursiv definiert. Es setzt sich
622 zusammen aus zwei Instanzen des bitonen Mergesort-Netzwerks halber Größe
623 $\bs{\frac{n}{2}}$ für je die Hälfte der Eingänge, sowie dem bitonen Mischer
624 für $n$~Leitungen $\operatorname{BM}(n)$. Das Rekursionsende ist das bitone
625 Mergesort-Netzwerk mit nur einer Leitung $\operatorname{BS}(1)$, welches als
626 leeres Komparatornetzwerk definiert ist. Entsprechend sind die
627 Komparatornetzwerke $\operatorname{BM}(2)$ und $\operatorname{BS}(2)$
628 identisch.
629
630 Bei der Konstruktion kommt die trichterförmige Anordnung der Komparatoren
631 (Abbildung~\ref{fig:bitonic-merge-tricheter}) gelegen, weil so die beiden
632 rekursiven Sortiernetzwerke in die gleiche Richtung sortieren können und so
633 alle Komparatoren in die gleiche Richtung zeigen.
634
635 \begin{figure}
636   \begin{center}
637   \input{images/batcher-8.tex}
638   \end{center}
639   \caption{\bs{8}, Batchers \emph{bitones Mergesort}-Netzwerk für 8~Eingänge.
640     Markiert sind die beiden Instanzen von \bs{4} (rot), die beiden bitonen
641     Mischer~\bm{4} (blau) und die Komparatoren, die im letzten rekursiven
642     Schritt hinzugefügt wurden (grün).}
643   \label{fig:bitonic-08}
644 \end{figure}
645
646 Das Sortiernetzwerk~$\operatorname{BS}(8)$ ist in
647 Abbildung~\ref{fig:bitonic-08} zu sehen. Eingezeichnet sind ebenfalls die
648 beiden Instanzen des Netzwerks~$\operatorname{BS}(4)$ (rot) sowie der bitone
649 Mischer~$\operatorname{BM}(8)$ (blau). Die trichterförmige Komparator-Kaskade,
650 die die bitone Eingabefolge in zwei bitone Ausgabefolgen transformiert, ist
651 grün hinterlegt.
652
653 Das \emph{bitone Mergesort}-Netzwerk mit einer Leitungszahl $n = 2^d$, die
654 eine Zweierpotenz ist, besteht aus $\frac{1}{4} n \log(n) \log(n+1) =
655 \Theta\left(n (log (n))^2\right)$ Komparatoren, die in $\frac{1}{2}
656 \log(n) \log(n+1) = \Theta(\log(n)^2)$ Schichten angeordnet sind.
657
658 \subsection{Das Odd-Even-Mergesort-Netzwerk}
659
660 Obwohl der Name ähnlich klingt, haben das \emph{Odd-Even-Mergesort}-Netzwerk
661 (OES) und das \emph{Odd-Even-Transpositionsort}-Netzwerk (siehe
662 Abschnitt~\ref{sect:odd_even_transpositionsort}) wenig gemein. Vielmehr ist
663 OES dem \emph{bitonen Mergesort}-Netzwerk, das im vorherigen Abschnitt
664 vorgestellt wurde, ähnlich: Auch dieses Sortiernetzwerk ist von
665 \textit{Kenneth~E. Batcher} gefunden worden und ist ebenfalls in~\cite{B1968}
666 beschrieben und initial analysiert worden. Eine weitere Gemeinsamkeit besteht
667 darin, dass es ebenfalls rekursiv durch einen Mischer definiert ist.
668
669 \subsubsection{Der \emph{Odd-Even}-Mischer}\label{sect:der_odd_even_mischer}
670
671 Der \emph{Odd-Even}-Mischer $\operatorname{OEM}(n,m)$ ist ein
672 Komparatornetzwerk, das zwei sortierte Folgen mit $n$ beziehungsweise $m$
673 Elementen zu einer sortierten Ausgabefolge mit $N = n+m$~Elementen
674 zusammenfügen kann. Dabei kommt es mit weniger Vergleichen aus als der
675 \emph{bitone Mischer}, der im Abschnitt~\ref{sect:der_bitone_mischer}
676 vorgestellt wurde. Im allgemeinen Fall, wenn die Anzahl der Leitungen keine
677 Zweierpotenz ist, kann das \emph{bitonic-Merge}-Netzwerk schneller sein 
678 als das \emph{Odd-Even-Merge}-Netzwerk.~\cite{KNUTH}
679
680 Der \emph{Odd-Even}-Mischer selbst ist ebenfalls rekursiv aufgebaut: Die
681 Eingabe für den Mischer mit $N = n + m$ Leitungen besteht aus den beiden
682 sortierten Folgen $U = \left(u_0, u_1, \ldots, u_{n-1}\right)$ und
683 $V = \left(v_0, v_1, \ldots, v_{m-1}\right)$. Die gesamte Eingabe sei
684 $W = \left(w_0, w_1, \ldots, w_{N-1}\right)$ mit:
685 \begin{equation}
686 w_i = \left\{ \begin{array}{ll}
687         u_i,     & i < n \\
688         v_{i-n}, & i \geqq n
689       \end{array} \right.,
690       \quad 0 \leqq i < N
691 \end{equation}
692
693 \begin{figure}
694   \begin{center}
695   \input{images/oe-merge.tex}
696   \end{center}
697   \caption{Schematischer Aufbau des \emph{Odd-Even-Merge}-Netzwerks. Die
698     beiden Dreiecke symbolisieren die beiden sortierten Folgen $U$ und $V$,
699     die Blöcke darunter die rekursiven Mischer mit etwa der Hälfte der
700     Leitungen. Im Vergleich zum \emph{bitonen Mischer} für 8~Leitungen kommt
701     dieses Schema mit einem Komparator weniger aus. Der Effekt wird durch den
702     rekursiven Aufbau verstärkt.}
703   \label{fig:oe-merge}
704 \end{figure}
705
706 Diese werden in insgesamt vier sortierte Folgen aufgeteilt, je eine Liste der
707 geraden Indizes und je eine Liste der ungeraden Indizes.
708 \begin{eqnarray}
709   U_{\textrm{gerade}}   &=& \left(u_0, u_2, u_4, \ldots\right) \\
710   U_{\textrm{ungerade}} &=& \left(u_1, u_3, u_5, \ldots\right) \\
711   V_{\textrm{gerade}}   &=& \left(v_0, v_2, u_4, \ldots\right) \\
712   V_{\textrm{ungerade}} &=& \left(v_1, v_3, u_5, \ldots\right)
713 \end{eqnarray}
714
715 Die geraden Folgen $U_{\textrm{gerade}}$ und $V_{\textrm{gerade}}$,
716 beziehungsweise die ungeraden Folgen $U_{\textrm{ungerade}}$ und
717 $V_{\textrm{ungerade}}$ werden rekursiv von kleineren \emph{Odd-Even}-Mischern
718 zusammengefügt, so dass sich am Ausgang der Mischer die Folgen
719 \begin{eqnarray}
720   W_{\textrm{gerade}}   &=& \left(w_0, w_2, w_4, \ldots\right) \\
721   W_{\textrm{ungerade}} &=& \left(w_1, w_3, w_5, \ldots\right)
722 \end{eqnarray}
723 ergeben.
724
725 Anschließend werden die Komparatoren zwischen benachbarten Leitungen
726 hinzugefügt,
727 \begin{equation}
728   w_{2i-1} \longleftrightarrow w_{2i}, \quad 1 \leqq i < \frac{N}{2}
729 \end{equation}
730 die die Folge~$W$ sortieren. Den schematischen Aufbau des
731 \emph{Odd-Even}-Mischers zeigt Abbildung~\ref{fig:oe-merge}.
732
733 Leider bricht die Rekursion nicht so schön ab, wie das beim {\em bitonen
734 Mischer} der Fall gewesen ist. Insbesondere für ${n = m = 1}$ würde --
735 entsprechend der Konstruktionsvorschrift -- ein leeres Netzwerk entstehen, was
736 offensichtlich nicht korrekt wäre. Die Abbruchbedingungen für den rekursiven
737 Aufbau lauten:
738 \begin{itemize}
739   \item Falls ${n = 0}$ oder ${m = 0}$: Das Netzwerk ist leer.
740   \item Falls ${n = 1}$ und ${m = 1}$: Das Netzwerk besteht aus einem
741   einzelnen Komparator.
742 \end{itemize}
743
744 Mit dem {\em 0-1-Prinzip} lässt sich zeigen, dass die resultierende Folge
745 sortiert ist. Da $U$ und $V$ sortiert sind, ist die Anzahl der Nullen in den
746 geraden Teilfolgen $U_{\textrm{gerade}}$, beziehungsweise
747 $V_{\textrm{gerade}}$ größer oder gleich der Anzahl der Nullen in den
748 ungeraden Teilfolgen $U_{\textrm{ungerade}}$ beziehungsweise
749 $V_{\textrm{ungerade}}$ --~die Einsen verhalten sich entsprechend umgekehrt.
750 Das trifft demnach auch auf die Folgen $W_{\textrm{gerade}}$ und
751 $W_{\textrm{ungerade}}$ entsprechend zu:
752 \begin{eqnarray}
753   \left|W_{\textrm{gerade}}\right|_0
754   &=& \left|U_{\textrm{gerade}}\right|_0
755     + \left|V_{\textrm{gerade}}\right|_0
756    =  \left\lceil \frac{1}{2} \left|U\right|_0 \right\rceil
757    +  \left\lceil \frac{1}{2} \left|V\right|_0 \right\rceil \\
758   \left|W_{\textrm{ungerade}}\right|_0
759   &=& \left|U_{\textrm{ungerade}}\right|_0
760     + \left|V_{\textrm{ungerade}}\right|_0
761    =  \left\lfloor \frac{1}{2} \left|U\right|_0 \right\rfloor
762    +  \left\lfloor \frac{1}{2} \left|V\right|_0 \right\rfloor
763 \end{eqnarray}
764 Daraus folgt, dass $W_{\textrm{gerade}}$ $0$, $1$ oder $2$ Nullen mehr enthält
765 als $W_{\textrm{ungerade}}$. In den ersten beiden Fällen ist die "`verzahnte"'
766 Ausgabe der beiden kleineren Mischer bereits sortiert. Nur im letzten Fall,
767 wenn $W_{\textrm{gerade}}$ zwei Nullen mehr enthält als
768 $W_{\textrm{ungerade}}$, muss genau eine Vertauschung stattfinden, um die
769 Ausgabe zu sortieren. Diese wird von den Komparatoren ausgeführt, die
770 benachbarte Leitungen miteinander vergleichen. Die jeweiligen Situationen sind
771 in Abbildung~\ref{fig:oe-post-recursive} dargestellt.
772
773 \begin{figure}
774   \centering
775   \subfigure[$\left|W_{\textrm{gerade}}\right|_0 - \left|W_{\textrm{ungerade}}\right|_0 = 0$]{\input{images/oe-post-recursive-diff0.tex}}
776   \qquad
777   \subfigure[$\left|W_{\textrm{gerade}}\right|_0 - \left|W_{\textrm{ungerade}}\right|_0 = 1$]{\input{images/oe-post-recursive-diff1.tex}}
778   \qquad
779   \subfigure[$\left|W_{\textrm{gerade}}\right|_0 - \left|W_{\textrm{ungerade}}\right|_0 = 2$]{\input{images/oe-post-recursive-diff2.tex}}
780   \caption{Die drei Situationen, die nach dem Verzahnen der Ausgaben der
781   kleineren \emph{Odd-Even}-Mischer entstehen können. Ist die Differenz der
782   Anzahl der Nullen gleich $0$ oder $1$, ist die Folge bereits sortiert. Im
783   letzten Fall stellt einer der Komparatoren sicher, dass das Ergebnis
784   sortiert ist.}
785   \label{fig:oe-post-recursive}
786 \end{figure}
787
788 Da die Teilfolgen $U$ und $V$ in jedem Rekursionsschritt etwa halbiert werden,
789 bricht die Rekursion nach $\Theta\left(\log (n) + \log (m)\right)$
790 Schritten ab. Die exakte Anzahl der benötigten Rekursionsschritte (und damit
791 Schichten im Mischer-Netzwerk), hängt von der längeren der beiden
792 Eingabefolgen ab und beträgt $1 + \lceil \log\left(\max(n, m)\right) \rceil$.
793
794 Die Anzahl der Komparatoren $K(n,m)$, die $\operatorname{OEM}(n,m)$ im
795 allgemeinen Fall verwendet, hängt gemäß der rekursiven Definition von der
796 Länge der Eingabefolgen, $n$ und $m$ ab:
797 \begin{displaymath}
798   K(n,m) = \left\{ \begin{array}{ll}
799     nm, & \mathrm{falls} \quad nm \leqq 1 \\
800     K\left(\left\lceil \frac{n}{2} \right\rceil, \left\lceil \frac{m}{2} \right\rceil\right)
801     + K\left(\left\lfloor \frac{n}{2} \right\rfloor, \left\lfloor \frac{m}{2} \right\rfloor\right)
802     + \left\lfloor \frac{1}{2} (m + n - 1) \right\rfloor & \mathrm{falls} \quad nm > 1
803   \end{array} \right.
804 \end{displaymath}
805 Leider ist es schwierig, diese allgemeine Formel in einer geschlossenen Form
806 anzugeben. Aus der Anzahl der Rekursionsschritte ist jedoch leicht erkennbar,
807 dass $K(n,m)$ in $\Theta(N \log (N))$ enthalten ist.
808
809 Für den wichtigen Spezialfall, dass $n = m = 2^{d-1}$ beträgt, lässt sich die
810 Anzahl der Komparatoren im Vergleich zum \emph{bitonen Mischer} angeben: Der
811 erste Rekursionsschritt der OEM-Konstruktion fügt
812 $\left\lfloor \frac{1}{2} (m + n - 1) \right\rfloor = \frac{N}{2} - 1$
813 Komparatoren ein -- einen Komparator weniger als der \emph{bitone Mischer} in
814 diesem Schritt. Das selbe gilt für die rekursiv verwendeten kleineren Mischer
815 $\operatorname{OEM}(\frac{n}{2}, \frac{n}{2})$ und so weiter bis
816 einschließlich $\operatorname{OEM}(2, 2)$, von denen es $2, 4, \dots,
817 \frac{N}{4} = 2^{\log(N)-2}$ Instanzen gibt. Insgesamt werden
818 \begin{displaymath}
819   \sum_{i=0}^{\log(N)-2} 2^i = 2^{\log(N) - 1} - 1 = \frac{N}{2} - 1 = n - 1
820 \end{displaymath}
821 Komparatoren eingespart. Damit ergibt sich
822 \begin{displaymath}
823   K\left(n = 2^{d-1}, n = 2^{d-1}\right) = \frac{1}{2} N \log(N) - \frac{N}{2} + 1
824 \end{displaymath}
825 für die Anzahl der Komparatoren, die von $\operatorname{OEM}(N = 2^d)$
826 benötigt werden.
827
828 \subsubsection{Das Odd-Even-Mergesort-Netzwerk}
829
830 Das \emph{Odd-Even-Mergesort}-Netzwerk $\operatorname{OES}(n)$ besteht, wie
831 das \emph{bitone Mergesort}-Netzwerk, rekursiv aus kleineren Varianten von
832 sich selbst und einem abschließenden \emph{Odd-Even}-Mischer. Die
833 effizientesten Sortiernetzwerke in Bezug auf Komparator- und Schichtzahl
834 entstehen, wenn die Anzahl der Leitungen jeweils halbiert wird. Somit besteht
835 \oes{n} aus $\oes{\left\lceil\frac{n}{2}\right\rceil}$,
836 $\oes{\left\lfloor\frac{n}{2}\right\rfloor}$ und
837 $\oem{\left\lceil\frac{n}{2}\right\rceil, \left\lfloor\frac{n}{2}\right\rfloor}$.
838 Die Rekursion endet mit $\operatorname{OES}(1)$ und $\operatorname{OES}(0)$,
839 die als leere Komparatornetzwerke definiert sind.
840
841 \begin{figure}
842   \begin{center}
843   \input{images/oe-mergesort-8.tex}
844   \end{center}
845   \caption{Das {\em Odd-Even-Mergesort-Netzwerk} für 8~Eingänge. Markiert
846   sind die Instanzen von $\operatorname{OES}(4)$ (rot), die beiden
847   \emph{Odd-Even}-Mischer $\operatorname{OEM}(4)$ für gerade und ungerade
848   Leitungen (blau) und die im ersten Rekursionsschritt hinzugefügten
849   Komparatoren zwischen benachbarten Leitungen (grün).}
850   \label{fig:odd-even-mergesort-08}
851 \end{figure}
852
853 In Abbildung~\ref{fig:odd-even-mergesort-08} ist das \oes{8}-Sortiernetzwerk
854 zu sehen. Rot markiert sind die beiden rekursiven Instanzen
855 $\operatorname{OES}(4)$. Die anderen Blöcke stellen den
856 \emph{Odd-Even}-Mischer für 8~Leitungen dar: die beiden blauen Blöcke sind
857 die rekursiven Instanzen von $\operatorname{OEM}(4)$, der grüne Block markiert
858 die Komparatoren, die im ersten Rekursionsschritt hinzugefügt werden.
859
860 Im Allgemeinen ist die Anzahl der Komparatoren, die vom
861 \emph{Odd-Even-Mergesort-Netz\-werk} verwendet wird, $k(n)$, direkt aus der
862 Definition, beziehungsweise der Konstruktionsanleitung abzulesen:
863 \begin{displaymath}
864   k(n) = k\left(\left\lceil\frac{n}{2}\right\rceil\right)
865        + k\left(\left\lfloor\frac{n}{2}\right\rfloor\right)
866        + K\left(\left\lceil\frac{n}{2}\right\rceil, \left\lfloor\frac{n}{2}\right\rfloor\right)
867 \end{displaymath}
868 Da es schwierig ist für $K(n,m)$ eine geschlossene Form anzugeben, ist eine
869 geschlossene Darstellung von $k(n)$ ebenfalls nicht ohne weiteres möglich. Es
870 ist allerdings bekannt, dass $k(n)$ in $\Theta\left(n \left(\log
871 (n)\right)^2\right)$ enthalten ist.
872
873 Für den wichtigen Spezialfall, dass $n = 2^d$ eine Zweierpotenz ist, kann die
874 Anzahl der Komparatoren wieder explizit angegeben werden. \textit{Kenneth
875 Batcher} zeigt in~\cite{B1968}, dass in diesem Fall
876 \begin{displaymath}
877   k(n = 2^d) = \frac{1}{4} n \left(\log (n)\right)^2 - \frac{1}{4}n\log(n) + n - 1
878 \end{displaymath}
879 gilt.
880
881 %\subsection{Das Pairwise-Sorting-Netzwerk}
882 %
883 %Das \emph{Pairwise-Sorting}-Netzwerk \ps{n} ist eine Konstruktionsvorschrift
884 %für Sortiernetzwerke, die 1992 von \textit{Ian Parberry} in seiner Arbeit „The
885 %Pairwise Sorting Network“ \cite{P1992} definiert wurde. Wenn die Anzahl der
886 %Leitungen $n = 2^d$ eine Zweierpotenz ist, hat das
887 %\emph{Pairwise-Sorting}-Netzwerk die selbe Effizienz und Geschwindigkeit wie
888 %das \emph{Odd-Even-Mergesort}-Netzwerk.
889
890 \newpage
891 \section{Transformation von Sortiernetzwerken}
892 \label{sect:tranformation}
893
894 \subsection{Komprimieren}
895
896 Komparatoren, die unterschiedliche Leitungen miteinander vergleichen, können
897 gleichzeitig ausgewertet werden, wie bereits in
898 Abschnitt~\ref{sect:einleitung_sortiernetzwerke} beschrieben. Durch manche
899 Transformationen, insbesondere das Entfernen einer Leitung wie in
900 Abschnitt~\ref{sect:leitungen_entfernen} beschrieben, kann es vorkommen, dass
901 die Komparatoren eines Sortiernetzwerks nicht mehr in der kleinstmöglichen
902 Anzahl von \emph{Schichten} angeordnet sind. Unter \emph{Komprimierung} wird
903 eine (Neu-)Gruppierung der Komparatoren verstanden, die jeden Komparator so
904 früh wie möglich ausführt. So entsteht die kleinstmögliche Anzahl von
905 \emph{Schichten}, in die sich ein Sortiernetzwerk unterteilen lässt.
906
907 Diese Anzahl ist insbesondere beim automatisierten Bewerten von
908 Komparatornetzwerken interessant, wie in Abschnitt~\ref{sect:sn-evolution:bewertung}
909 beschrieben. Die Anzahl der Schichten kann künstlich vergrößert werden, indem
910 Komparatoren später angewendet werden. Deshalb sollte vor einer Bewertung, die
911 die Anzahl der Schichten als Bewertungskriterium verwendet, immer eine
912 Komprimierung durchgeführt werden.
913
914 \subsection{Normalisieren}
915 \label{sect:normalisieren}
916
917 \begin{figure}
918   \centering
919   \subfigure[$S(8)$ (nach Konstruktion)]{\input{images/batcher-8-nonstd.tex}\label{fig:bitonic-nonstd}}
920   \subfigure[$S(8)$ (normalisiert)]{\input{images/batcher-8-std.tex}\label{fig:bitonic-std}}
921   \caption{Jedes Sortiernetzwerk kann in ein Standard-Sortiernetzwerk
922   transformiert werden. Gezeigt ist das bitone Sortiernetzwerk nach der
923   intuitiven Konstruktion und die normalisierte Variante.}
924   \label{fig:beispiel_normalisieren}
925 \end{figure}
926
927 Ein \emph{Standard-Sortiernetzwerk} oder \emph{normalisiertes Sortiernetzwerk}
928 ist ein Sortiernetzwerk, dessen Komparatoren alle in die selbe Richtung
929 zeigen.\footnote{Die Konvention in dieser Arbeit ist, dass in diesem Fall alle
930 Pfeile nach unten zeigen. Das Minimum wird auf der untersten, das Maximum auf
931 der obersten Leitung ausgegeben.} Jedes Sortiernetzwerk kann in eine
932 normaliesierte Variante transformiert werden. Dazu gibt beispielsweise
933 \emph{Donald~E. Knuth} in~\cite{KNUTH} einen Algorithmus an.
934
935 Abbildung~\ref{fig:beispiel_normalisieren} stellt das \emph{bitone
936 Mergesort}-Netzwerk in zwei Varianten dar. Abbildung~\ref{fig:bitonic-nonstd}
937 zeigt das Netzwerk nach der Konstruktionsvorschrift, siehe auch
938 Abbildung~\ref{fig:bitonic-merge-normal}: In den ersten drei Schichten werden
939 die untere und die obere Hälfte gegenläufig sortiert. Das heißt, dass nach
940 drei Schritten die eine Hälfte auf- und die andere Hälfte absteigend sortiert
941 ist. In den Schichten~4 bis~6 folgt der bitone Mischer entsprechend der
942 rekursiven Definition.
943
944 In Abbildung~\ref{fig:bitonic-std} ist die normalisierte Version des bitonen
945 Mergesort-Netzwerks zu sehen. Alle Komparatoren zeigen hier in die selbe
946 Richtung. Statt dem typischen „Treppenmuster“ sind abwechselnd das Treppen-
947 und das Trichtermuster zu sehen.
948
949 \subsection{Zwei Netzwerke kombinieren}
950
951 Um Sortiernetzwerke als \emph{Individuen} evolutionärer Algorithmen verwenden
952 zu können, muss es möglich sein, zwei Sortiernetzwerke zu einem neuen
953 Sortiernetzwerk zusammenzufassen.
954
955 Diese Technik wurde in den vorangegangen Abschnitten bereits verwendet,
956 beispielsweise um zwei \emph{bitone Mergesort}-Netzwerke mit jeweils der
957 halben Leitungszahl, $\operatorname{BS}\left(\frac{n}{2}\right)$, zu einem
958 einzigen Sortiernetzwerk $\operatorname{BS}(n)$ zu kombinieren. Auch das
959 \emph{Odd-Even-Mergesort}-Netzwerk $\operatorname{OES}(n)$ wurde auf diese Art
960 und Weise rekursiv aufgebaut.
961
962 Die vorgestellten \emph{Mischer} erwarten als Eingabe zwei bereits sortierte
963 Folgen. \emph{Wie} diese Folgen sortiert wurden ist unerheblich. Entsprechend
964 können wir beliebige Sortiernetzwerke einsetzen, um die beiden Eingabefolgen
965 zu sortieren und die Ausgaben mit einem der beschriebenen Mischer
966 zusammenfügen.
967
968 Beispielsweise kann die Ausgabe von zwei \emph{bitonen Mergesort-Netzwerken}
969 $\operatorname{BS}(8)$ mit je 8~Leitungen mit dem
970 \emph{Odd-Even-Merge}-Netzwerk $\operatorname{OEM(8,8)}$ zu einer sortierten
971 Gesamtfolge zusammengefügt werden. Das resultierende Sortiernetzwerk besitzt
972 73~Komparatoren (zum Vergleich: $\operatorname{BS}(16)$ benötigt
973 80~Komparatoren, $\operatorname{OES}(16)$ nur 63).
974
975 Verbesserungen der Effizienz (die Anzahl der benötigten Komparatoren),
976 beziehungsweise der Geschwindigkeit (die Anzahl der Schichten) eines „kleinen“
977 Sortiernetzwerks, übertragen sich direkt auf das resultierende Gesamtnetzwerk.
978 Das \emph{Odd-Even-Mergesort}-Netzwerk $\operatorname{OES}(9)$ benötigt
979 beispielsweise 26~Komparatoren, die in 9~Schichten angeordnet sind. Es sind
980 allerdings Sortiernetzwerke mit 9~Eingängen bekannt, die lediglich
981 25~Komparatoren in 7~Schichten benötigen. Wenn zwei dieser Netzwerke mit dem
982 \emph{Odd-Even}-Mischer kombiniert werden, entsteht ein 18-Sortiernetzwerk,
983 das aus 80~Komparatoren in 11~Schichten besteht. Damit ist das resultierende
984 Netzwerk genauso schnell wie das Sortiernetzwerk mit 18~Eingängen, das
985 \textit{Sherenaz~W. Al-Haj Baddar} und \textit{Kenneth~E. Batcher} in ihrer
986 Arbeit „An 11-Step Sorting Network for 18~Elements“~\cite{BB2009} vorstellen,
987 benötigt aber 6~Komparatoren weniger. $\operatorname{OES}(18)$ benötigt
988 82~Komparatoren in 13~Schichten.
989
990 Das Zusammenfassen von zwei Sortiernetzwerken durch Hintereinanderausführung
991 ist nicht sinnvoll: Da die Ausgabe des ersten Sortiernetzwerks bereits
992 sortiert ist, ist das zweite Sortiernetzwerk überflüssig. Eine
993 Aneinanderreihung der Art „die ersten $x$~Schichten des einen, dann die
994 letzten $y$~Schichten des anderen Sortiernetzwerks“ zerstören im Allgemeinen
995 die Sortiereigenschaft. Die Sortiereigenschaft des resultierenden
996 Komparatornetzwerks müsste überprüft werden, was nach heutigem Wissensstand
997 nur mit exponentiellem Aufwand möglich ist.
998
999 \subsection{Leitungen entfernen}
1000 \label{sect:leitungen_entfernen}
1001
1002 Im vorherigen Abschnitt wurde gezeigt, dass es mithilfe von \emph{Mischern}
1003 möglich ist, aus zwei Sortiernetzwerken mit je $n$~Eingängen
1004 ein neues Sortiernetzwerk mit $2n$~Eingängen zu erzeugen. Für einen
1005 beabsichtigen \emph{evolutionären Algorithmus} ist es jedoch notwendig, dass
1006 sich die Anzahl der Eingänge nicht verändert. Es soll wieder ein
1007 Sortiernetzwerk mit $n$~Eingängen entstehen.
1008
1009 Man kann ein gegebenes Sortiernetzwerk mit $n$~Eingängen auf ein
1010 Sortiernetzwerk mit ${n-1}$~Leitungen verkleinern, indem man eine Leitung
1011 „eliminiert“. Dazu wird angenommen, dass das Minimum oder das Maximum an einem
1012 bestimmten Eingang anliegt. Der Weg, den das Minimum beziehungsweise das
1013 Maximum durch das Sortiernetzwerk nimmt, ist eindeutig bestimmt und endet an
1014 einem der „Ränder“, also auf der Leitung mit dem höchsten oder dem niedrigsten
1015 Index. Insbesondere ist bekannt, welche Komparatoren „berührt“ werden und
1016 welche dafür sorgen, dass der Wert die Leitung wechselt, da das Minimum jeden
1017 Vergleich „verliert“ und das Maximum jeden Vergleich „gewinnt“. Die
1018 Abbildung~\ref{fig:oe-transposition-cut0} zeigt den Weg eines Maximums durch
1019 das \emph{Odd-Even-Transpositionsort}-Netzwerk.
1020
1021 Im ersten Schritt wird eine Leitung ausgewählt und Maximum oder Minimum auf
1022 dieser Leitung angenommen. Dadurch ist der Weg durch das Sortiernetzwerk
1023 eindeutig festgelegt.
1024
1025 \begin{figure}
1026   \centering
1027   \subfigure[Auf der Leitung~4 wird $-\infty$ angelegt. Dadurch ist der Pfad
1028   durch das Sortiernetzwerk eindeutig festgelegt.]{\input{images/oe-transposition-cut0.tex}\label{fig:oe-transposition-cut0}}
1029   \subfigure[Komparatoren, die einen Wechsel der Leitungen bewirken, werden
1030   durch sich kreuzende Leitungen ersetzt.]{\input{images/oe-transposition-cut1.tex}\label{fig:oe-transposition-cut1}}
1031   \subfigure[Leitung~4 wurde entfernt. Übrig bleibt ein Sortiernetzwerk mit
1032   7~Leitungen.]{\input{images/oe-transposition-cut2.tex}\label{fig:oe-transposition-cut2}}
1033   \subfigure[Die Leitungen wurden wieder gerade eingezeichnet und die
1034   Komparatoren regelmäßig angeordnet. Blau eingezeichnet ist \oet{7}.]{\input{images/oe-transposition-cut3.tex}\label{fig:oe-transposition-cut3}}
1035   \caption{Eine Leitung wird aus dem
1036   \emph{Odd-Even-Transpositionsort}-Netzwerk \oet{8} entfernt: Auf der rot
1037   markierten Leitung wird $-\infty$ angelegt. Da der Wert bei jedem Komparator
1038   nach unten weiter gegeben wird, ist der Pfad fest vorgegeben. Da die
1039   restlichen Werte trotzdem noch richtig sortiert werden müssen, kann dieser
1040   Pfad heraus getrennt werden. In der letzten Abbildung ist \oet{7} markiert.}
1041   \label{fig:oe-transposition-cut}
1042 \end{figure}
1043
1044 Im nächsten Schritt werden alle beteiligten Komparatoren gelöscht,
1045 beziehungsweise ersetzt: Komparatoren, die {\em nicht} zu einem Wechsel der
1046 Leitung geführt haben, werden ersatzlos gelöscht. Diese Komparatoren sind in
1047 Abbildung~\ref{fig:oe-transposition-cut0} grün markiert. Die Komparatoren, die
1048 zum Wechsel der Leitung geführt haben, werden durch sich kreuzende Leitungen
1049 ersetzt. Das Resultat ist eine Leitung, auf der das Minimum beziehungsweise
1050 das Maximum angenommen wird, die an unterster oder oberster Stelle endet und
1051 die keine Komparatoren mehr berührt
1052 (Abbildung~\ref{fig:oe-transposition-cut1}).
1053
1054 Die Werte auf den verbleibenden $(n-1)$~Leitungen müssen vom restlichen
1055 Komparatornetzwerk immer noch sortiert werden: Es wurde lediglich die
1056 \emph{Position}  des Minimums oder des Maximums in der Eingabe angenommen. Ein
1057 Sortiernetzwerk muss die Eingabe sortieren, unabhängig davon auf welcher
1058 Leitung das Minimum oder das Maximum liegt. Das Sortiernetzwerk unter diese
1059 Annahme auszuwerten -- über die verbleibenden Eingänge wurde keine Aussage
1060 getroffen. Entsprechend müssen die verbleibenden Ausgänge eine sortierte Liste
1061 mit $(n-1)$~Elementen darstellen.
1062
1063 Wird die Minimum- beziehungsweise Maximum-Leitung entfernt, wie in
1064 Abbildung~\ref{fig:oe-transposition-cut2} dargestellt, bleibt das
1065 Sortiernetzwerk für $(n-1)$~Leitungen übrig. Je nachdem, ob auf einer Leitung
1066 ein Minimum oder ein Maximum angenommen wird, wird das eliminieren einer
1067 Leitung auf diese Art und Weise als \emph{Minimum-Schnitt}, beziehungsweise
1068 \emph{Maximum-Schnitt} bezeichnet.
1069
1070 Die letzte Abbildung, \ref{fig:oe-transposition-cut3}, zeigt das
1071 Sortiernetzwerk wieder mit den üblichen geraden Leitungen und die rot
1072 markierten Komparatoren sind verschoben, so dass sich eine kompaktere
1073 Darstellung ergibt. Außerdem ist das
1074 \emph{Odd-Even-Transpositionsort}-Netzwerk für sieben Werte markiert. Der
1075 zusätzliche Komparator vor dem \oet{7} hat keinen Einfluss auf die Ausgabe und
1076 kann entfernt werden.
1077
1078 Durch das Ersetzen von Komparatoren durch gekreuzte Leitungen werden häufig
1079 \emph{Nicht-Standard-Sortiernetzwerke} erzeugt. Im Anschluss an einen
1080 \emph{Schnitt} empfiehlt es sich deshalb, das Sortiernetzwerk zu
1081 \emph{normalisieren}, wie in Abschnitt~\ref{sect:normalisieren} beschrieben.
1082
1083 \subsubsection{Anzahl möglicher und unterschiedlicher Schnittmuster}
1084 \label{sect:anzahl_schnittmuster}
1085
1086 Der Eliminierungsschritt kann iterativ angewendet werden, um aus einem
1087 Sortiernetzwerk mit $n$~Ein\-gängen Sortiernetzwerke mit $n-1$, $n-2$,
1088 $n-3$,~\dots Eingängen zu erzeugen. Insbesondere können auf diese Art und
1089 Weise Sortiernetzwerke mit $2n$~Eingängen auf Sortiernetzwerke mit
1090 $n$~Eingängen reduziert werden. Als \emph{$k$-Schnittmuster} bezeichnet man
1091 die $k$~Minimum- und Maximum-Schnitte, die nacheinander angewendet ein
1092 $n$-Sortiernetzwerk auf ein ${(n-k)}$-Sortiernetz\-werk reduzieren.
1093
1094 Zwei Schnittmuster heißen \emph{äquivalent} bezüglich~$S$, wenn ihre Anwendung
1095 auf das Sortiernetzwerk~$S$ das selbe Ergebnis liefert. Ansonsten heißen die
1096 Schnittmuster \emph{unterschiedlich} bezüglich~$S$. 
1097
1098 Bei einem Sortiernetzwerk mit $n$~Eingängen gibt es $2n$~Möglichkeiten eine
1099 Leitung zu entfernen: Auf jeder der $n$~Leitungen kann sowohl das Minimum als
1100 auch das Maximum angenommen werden. Wendet man das Verfahren iterativ an, um
1101 ein $n$-Sortiernetzwerk auf ein ${(n-k)}$-Sortiernetzwerk zu reduzieren,
1102 ergeben sich insgesamt
1103 \begin{displaymath}
1104   \prod_{i=n}^{1+n-k} 2i = 2^k \cdot \frac{n!}{(n-k)!}
1105   \quad (n > m)
1106 \end{displaymath}
1107 \emph{mögliche} Schnittmuster. Diese Schnittmuster sind nicht alle
1108 unterschiedlich. Wird beispielsweise das Minimum auf der untersten Leitung
1109 und das Maximum auf der obersten Leitung eines Standard-Sortiernetzwerks
1110 angenommen, führen beide möglichen Schnitt-Reihenfolgen zum selben Ergebnis.
1111
1112 \textit{Moritz Mühlenthaler} zeigt in seiner Arbeit~\cite{M2009}, dass es
1113 möglich ist, mehrere Eingänge gleichzeitig mit Minimum beziehungsweise Maximum
1114 vorzubelegen, ohne die Menge der erreichbaren Sortiernetzwerke einzuschränken.
1115 Dadurch wird die Anzahl der möglichen Schnittmuster reduziert, die Menge der
1116 so erzeugbaren Sortiernetzwerke bleibt aber unverändert. Die Anzahl der
1117 möglichen Schnittmuster setzt sich zusammen aus der Anzahl von Möglichkeiten,
1118 $k$~Leitungen aus $n$~Leitungen auszuwählen, und die möglichen Minimum-~/
1119 Maximum-Muster. Damit ergibt sich folgende Formel für die Anzahl der möglichen
1120 Schnittmuster:
1121 \begin{equation}\label{eqn:anzahl_schnittmuster}
1122   2^k \cdot \left( \begin{array}{c} n \\ k \end{array} \right)
1123   = 2^{k} \cdot \frac{n!}{k! (n-k)!}
1124   = 2^{k} \cdot \frac{n!}{(n-k)!} \cdot \frac{1}{k!}
1125   \quad (1 \leqq k < n)
1126 \end{equation}
1127
1128 Die Anzahl der möglichen Schnittmuster wird mit der Anzahl der zu entfernenden
1129 Leitungen sehr schnell sehr groß. Um ein Sortiernetzwerk mit 32~Eingängen auf
1130 ein Sortiernetzwerk mit 16~Eingängen zu reduzieren, ist ein Schnittmuster mit
1131 16~Schnitten notwendig, für das es bereits etwa ${3,939 \cdot 10^{13}}$
1132 Möglichkeiten gibt. Ein Ausprobieren aller Möglichkeiten ist für große
1133 Netzwerke nicht oder nur unter erheblichem Ressourcenaufwand möglich.
1134
1135 Die Anzahl der \emph{unterschiedlichen} Schnittmuster ist allerdings kleiner
1136 als die Anzahl der \emph{möglichen} Schnittmuster. Für jeden Komparator auf
1137 der ersten Stufe gibt es neun verschiedene Eingangskonfigurationen: Für beide
1138 Eingänge gibt es drei mögliche Eingangswerte, Minimum, Maximum und
1139 unspezifiziert. Es gibt drei Konfigurationen, bei denen an beiden Eingängen
1140 der gleiche Wert angelegt wird, und sechs Konfigurationen, bei denen sich die
1141 Werte unterscheiden.
1142
1143 Bei diesen letzten sechs Konfigurationen werden je zwei auf das selbe
1144 Ausgangsmuster abgebildet, weil die Position des Minimums beziehungsweise des
1145 Maximums durch den Komparator vorgegeben wird. Das heißt, dass die neun
1146 unterschiedlichen Eingangsmuster nur sechs unterschiedliche Ausgangsmuster
1147 erzeugen. In der zweiten und allen folgenden Schichten kann man diesen
1148 Unterschied nicht mehr erkennen. In allen sechs Fällen, in denen sich die
1149 Eingänge unterscheiden, wird anschließend der Komparator entfernt, so dass
1150 sich die Resultate auch in der ersten Schicht nicht unterscheiden.
1151
1152 \begin{figure}
1153   \begin{center}
1154     \includegraphics[viewport=0 0 425 262,width=15cm]{images/count-cuts-16.pdf}
1155   \end{center}
1156   \caption{Anzahl der \emph{unterschiedlichen} Sortiernetzwerke, die durch
1157   8-Schnittmuster aus $\operatorname{OES}(16)$, $\operatorname{BS}(16)$ und
1158   $\operatorname{PS}(16)$ hervorgegangen sind. Die Anzahl der
1159   unterschiedlichen Netzwerke nach $10^6$~Iterationen ist 3519 für das
1160   \emph{Odd-Even-Mergesort}-Netzwerk, 4973 für das \emph{bitone
1161   Mergesort}-Netzwerk und 18764 für das \emph{Pairwise-Sorting}-Netzwerk.}
1162   \label{fig:count-cuts-16}
1163 \end{figure}
1164
1165 Alleine durch Betrachten der ersten Schicht von Komparatoren konnte die Anzahl
1166 der \emph{unterschiedlichen} Schnittmuster auf höchstens $\frac{2}{3}$ der
1167 \emph{möglichen} Schnittmuster reduziert werden. Um die Anzahl der
1168 \emph{unterschiedlichen} Schnittmuster experimentell zu ermitteln, wurden je
1169 eine Million zufällige 8-Schnittmuster auf die 16-Sortiernetzwerke \oes{16},
1170 \bs{16} und \ps{16} angewandt. Anschließend wurde mithilfe einer Hashtabelle
1171 überprüft, ob das resultierende Sortiernetzwerk schon von einem
1172 \emph{äquivalenten} Schnittmuster erzeugt wurde. Falls das Sortiernetzwerk
1173 noch nicht in der Hashtabelle enthalten war, wurde der Zähler für
1174 unterschiedliche Schnittmuster erhöht und das Sortiernetzwerk eingefügt.
1175
1176 Abbildung~\ref{fig:count-cuts-16} trägt die Anzahl der
1177 \emph{unterschiedlichen} Schnittmuster gegen die Anzahl der zufälligen
1178 Schnittmuster auf. Klar zu sehen ist, dass sich die Anzahl der erzeugten
1179 Sortiernetzwerke nach $500.000$~Iterationen nur noch gering verändert und der
1180 Wert nach $1.000.000$~Iterationen allem Anschein nach dem Endwert schon sehr
1181 nahe ist.
1182
1183 Die Anzahl der möglichen 8-Schnittmuster ist entsprechend der
1184 Formel~\eqref{eqn:anzahl_schnittmuster} 3.294.720. Diese möglichen
1185 Schnittmuster führen aber nur zu wenigen \emph{unterschiedlichen}
1186 Sortiernetzwerken: 3519 ($\approx 0,1\%$) im Fall des
1187 \emph{Odd-Even-Mergesort}-Netzwerks, 4973 ($\approx 0,15\%$) beim
1188 \emph{bitonen Mergesort}-Netzwerk und 18764 ($\approx 0,57\%$) beim
1189 \emph{Pairwise-Sorting}-Netzwerk. Zwar ist es möglich, dass mehr Iterationen
1190 die Anzahl der unterschiedlichen Schnittmuster noch wachsen lässt. Die Graphen
1191 in Abbildung~\ref{fig:count-cuts-16} geben jedoch Grund zu der Annahme, dass
1192 die Anzahl dieser zusätzlichen, unterschiedlichen Schnittmuster
1193 vernachlässigbar klein ist.
1194
1195 Bedingt durch die sehr große Anzahl möglicher Schnittmuster ist dieses
1196 Experiment für größere Sortiernetzwerke nicht sinnvoll durchführbar. Die
1197 Hashtabelle würde mehr Arbeitsspeicher benötigen als in derzeitigen Rechnern
1198 vorhanden ist, bevor ein entsprechender Graph den linearen Bereich für
1199 „kleine“ x-Werte verlässt.
1200
1201 Um die Anzahl der unterschiedlichen Schnittmuster trotzdem abschätzen zu
1202 können, kann man sich einer stochastischen Methode bedienen, der sogenannten
1203 \emph{Monte-Carlo-Methode}, die \textit{Rolf Wanka} in~\cite{W2006} für
1204 schwierige Zählprobleme vorstellt. Zunächst generiert man eine Menge~$S$ von
1205 $k$~unterschiedlichen Schnittmustern. Anschließend werden $n$~Schnittmuster
1206 zufällig erzeugt und überprüft, ob sie in der Menge~$S$ enthalten sind. Unter
1207 der Annahme, dass auf diese Art und Weise Sortiernetzwerke zufällig und
1208 gleichverteilt erzeugt werden, entspricht das Verhältnis der zufälligen
1209 Schnittmuster, die in $S$ enthalten sind, und $n$ gleich dem Verhältnis von
1210 $k$ und der Anzahl der unterschiedlichen Schnittmuster insgesamt. Damit kann
1211 die Anzahl der unterschiedlichen Schnittmuster abgeschätzt werden.
1212
1213 \begin{figure}
1214   \begin{center}
1215     \includegraphics[viewport=0 0 425 262,width=15cm]{images/collisions-10000-1000000-32.pdf}
1216   \end{center}
1217   \caption{Abschätzung der unterschiedlichen Schnittmuster mit der
1218   \emph{Monte-Carlo-Methode} für $\operatorname{OES}(32)$ und
1219   $\operatorname{BS}(32)$.}
1220   \label{fig:collisions-10000-1000000-32}
1221 \end{figure}
1222
1223 In Abbildung~\ref{fig:collisions-10000-1000000-32} ist das Ergebnis des
1224 Monte-Carlo-Algorithmus für 16-Schnittmuster zu sehen, die auf
1225 $\operatorname{OES}(32)$ und $\operatorname{BS}(32)$ angewandt wurden: Von
1226 jedem Sortiernetzwerk wurden zunächst eine Menge~$S$ von 10.000
1227 \emph{unterschiedlichen} Schnittmustern erzeugt. Anschließend wurden 1.000.000
1228 zufällige Schnittmuster erzeugt und der Anteil der zufälligen Schnittmuster,
1229 die \emph{äquivalent} zu einem in~$S$ enthalten Schnittmuster sind, berechnet.
1230 Für $\operatorname{OES}(32)$ war dieser Anteil etwa $0,19 \%$, für
1231 $\operatorname{BS}(32)$ etwa $0,29 \%$. Das ergibt eine Abschätzung von $5,2
1232 \cdot 10^6$ unterschiedlichen 16-Schnittmustern für $\operatorname{OES}(32)$
1233 und $3,4 \cdot 10^6$ für $\operatorname{BS}(32)$.
1234
1235 \begin{figure}
1236   \begin{center}
1237     \includegraphics[viewport=0 0 425 262,width=15cm]{images/collisions-100000-1000000-32-ps.pdf}
1238   \end{center}
1239   \caption{Abschätzung der unterschiedlichen Schnittmuster mit der
1240   \emph{Monte-Carlo-Methode} für $\operatorname{PS}(32)$. 385 von 1.000.000
1241   zufälligen Schnittmustern waren äquivalent zu einem Schnittmuster in einer
1242   Menge von 100.000. Daraus ergibt sich eine Schätzung von $2,6 \cdot 10^8$
1243   unterschiedlichen Schnittmustern.}
1244   \label{fig:collisions-100000-1000000-32-ps}
1245 \end{figure}
1246
1247 Im vorherigen Abschnitt wurde das \emph{Pairwise-Sorting}-Netzwerk
1248 $\operatorname{PS}(32)$ nicht betrachtet, da es für dieses Netzwerk viel mehr
1249 unterschiedliche 16-Schnittmuster gibt als für $\operatorname{OES}(32)$ und
1250 $\operatorname{BS}(32)$. In Anbetracht der Tatsache, dass die Anzahl der
1251 unterschiedlichen 8-Schnittmuster für $\operatorname{PS}(16)$ in
1252 Abbildung~\ref{fig:count-cuts-16} bereits mehr als dreimal größer war als die
1253 Anzahl für $\operatorname{OES}(16)$ beziehungsweise $\operatorname{BS}(16)$,
1254 ist dieser Umstand wenig verwunderlich. Entsprechend hätte man in einem
1255 kombinierten Graphen keine Details mehr erkennen können. Aufgrund der hohen
1256 Anzahl unterschiedlicher Schnittmuster, wurde für das gleiche Experiment mit
1257 $\operatorname{PS}(32)$ eine initiale Menge von 100.000 unterschiedlichen
1258 Schnittmustern erzeugt. Trotzdem wurden nach 1.000.000 Iterationen nur 385
1259 Schnittmuster gefunden, die zu einem Schnittmuster in der Menge äquivalent
1260 waren. Daraus ergibt sich eine Abschätzung von $2,6 \cdot 10^8$
1261 unterschiedlichen Schnittmustern -- zwei Zehnerpotenzen mehr als bei den
1262 vorherigen Sortiernetzwerken, aber immer noch fünf Zehnerpotenzen kleiner als
1263 die Anzahl der \emph{möglichen} Schnittmuster.
1264
1265 \newpage
1266 \section{Der \textsc{SN-Evolution}-Algorithmus}
1267 \label{sect:sn-evolution}
1268
1269 Der \textsc{SN-Evolution}-Algorithmus ist ein \emph{evolutionärer
1270 Algorithmus}, der die in den vorherigen Abschnitten beschriebenen Mischer
1271 (Abschnitt~\ref{sect:konstruktive_netzwerke}) und Schnittmuster
1272 (Abschnitt~\ref{sect:leitungen_entfernen}) verwendet, um „möglichst gute“
1273 Sortiernetzwerke zu erzeugen. Was ein „gutes“ Sortiernetzwerk ausmacht, wird
1274 in Abschnitt~\ref{sect:sn-evolution:bewertung} behandelt. Informationen zur Implementierung
1275 von \textsc{SN-Evolution} befinden sich in
1276 Abschnitt~\ref{sect:implementierung}.
1277
1278 \subsection{Bewertungsfunktion}\label{sect:sn-evolution:bewertung}
1279
1280 Um Sortiernetzwerke überhaupt optimieren zu können, muss zunächst die
1281 {\em Güte} eines Netzwerks definiert werden. Prinzipiell gibt es zwei Ziele,
1282 die bei Sortiernetzwerken verfolgt werden können:
1283 \begin{itemize}
1284   \item Möglichst wenige Komparatoren („effizient“)
1285   \item Möglichst wenige Schichten („schnell“)
1286 \end{itemize}
1287
1288 \begin{figure}
1289   \centering
1290   \subfigure[16-Sortiernetzwerk aus 60~Komparatoren in 10~Schichten. Das
1291   Netzwerk wurde von \textit{M.~W. Green} konstruiert und 1969 in~\cite{G1972}
1292   veröffentlicht.]{\input{images/16-green.tex}\label{fig:16-green}}
1293   \subfigure[16-Sortiernetzwerk aus 61~Komparatoren in 9~Schichten. Das
1294   Netzwerk wurde von \textit{D. Van~Voorhis} 1974 in~\cite{V1974}
1295   veröffentlicht.]{\input{images/16-voorhis.tex}\label{fig:16-voorhis}}
1296   \caption{Das effizienteste und das schnellste Sortiernetzwerk für
1297   16~Leitungen, das derzeit bekannt ist.}
1298   \label{fig:16-best-known}
1299 \end{figure}
1300 Diese Ziele führen im Allgemeinen zu unterschiedlichen Netzwerken.
1301 Beispielsweise besteht das \emph{effizienteste} bekannte Sortiernetzwerk für
1302 16~Eingänge aus 60~Komparatoren in 10~Schichten. Es ist in
1303 Abbildung~\ref{fig:16-green} dargestellt. Das \emph{schnellste} bekannte
1304 16-Sortiernetzwerk besteht aus 61~Komparatoren in nur 9~Schichten und ist in
1305 Abbildung~\ref{fig:16-voorhis} zu sehen.
1306
1307 \textsc{SN-Evolution} verwendet eine Gütefunktion, die die beiden Ziele
1308 "`effizient"' und "`schnell"' berücksichtigen kann. Sie hat die folgende
1309 generelle Form:
1310 \begin{equation}
1311   \operatorname{Guete}(S) = w_{\mathrm{Basis}}
1312                     + w_{\mathrm{Komparatoren}} \cdot \left|S\right|_\mathrm{Komparatoren}
1313                     + w_{\mathrm{Schichten}} \cdot \left|S\right|_\mathrm{Schichten}
1314 \end{equation}
1315 Die Parameter $w_{\mathrm{Komparatoren}}$ und $w_{\mathrm{Schichten}}$ dienen
1316 dabei der Festlegung des Optimierungsziels. Wenn einer der beiden Parameter
1317 gleich Null ist, wird nur das jeweils andere Ziel verfolgt. Sind beide
1318 Parameter gleich Null, werden alle Netzwerke mit der gleich Güte bewertet --
1319 jegliche Ergebnisse sind dann rein zufälliger Natur.\footnote{Dass dies nicht
1320 so schlecht ist wie man intuitiv vermuten könnte, zeigt der
1321 \textsc{SN-Markov}-Algorithmus in Abschnitt~\ref{sect:markov}.}
1322
1323 Da möglichst effiziente und schnelle Sortiernetzwerke gefunden werden sollen,
1324 ist ein kleiner Wert von $\operatorname{Guete}(S)$ besser als ein großer Wert.
1325 Das heißt, dass das Ziel von \textsc{SN-Evolution} ist,
1326 $\operatorname{Guete}(S)$ zu \emph{minimieren}.
1327
1328 Mit dem Parameter $w_{\mathrm{Basis}}$ kann auf die Selektion Einfluss
1329 genommen werden. Ist er groß, wird der relative Unterschied der Güten
1330 verschiedener Netzwerke kleiner, was die {\em Exploration}, das Absuchen des
1331 gesamten Lösungsraums, begünstigt. Wählt man $w_{\mathrm{Basis}}$ hingegen
1332 klein -- in Abhängigkeit von den anderen beiden Parametern sind auch negative
1333 Werte möglich -- werden die relativen Unterschiede groß. Dadurch wird die {\em
1334 Exploitation}, das Streben zu (lokalen) Optima, verstärkt. In~\cite{WW2002}
1335 geben \textit{Karsten und Nicole Weicker} einen Überblick über
1336 Selektionsmethoden und Rekombinationsmöglichkeiten.
1337
1338 Diese Parameter haben einen großen Einfluss auf die Geschwindigkeit, mit der
1339 der \textsc{SN-Evolution}-Algorithmus konvergiert und ob er tatsächlich gute
1340 Lösungen findet oder sich in \emph{lokalen} Optima "`verfängt"'. Leider gibt
1341 es kein Patentrezept für die Wahl der Parameter, so dass für verschiedene
1342 Leitungszahlen und Mischer-Typen experimentiert werden muss.
1343
1344 Als guter Standardansatz für \textsc{SN-Evolution} haben sich die folgenden
1345 Werte herausgestellt:
1346 \begin{eqnarray*}
1347   w_{\mathrm{Basis}}        &=& 0 \\
1348   w_{\mathrm{Komparatoren}} &=& 1 \\
1349   w_{\mathrm{Schichten}}    &=& \left|S\right|_\mathrm{Leitungen}
1350 \end{eqnarray*}
1351 Sofern nicht anders angegeben, werden diese Werte im Folgenden zur Bewertung
1352 von Sortiernetzwerken verwendet. Die Bewertungsfunktion bevorzugt mit diesen
1353 Konstanten \emph{schnelle} Sortiernetzwerke, da das Einsparen einer Schicht
1354 ein höheres Gewicht als das Einsparen von Komparatoren hat.
1355
1356 Wenn der \textsc{SN-Evolution}-Algorithmus nach \emph{effizienten}
1357 Sortiernetzwerken suchen soll, werden alternative Werte für die Konstanten der
1358 Bewertungsfunktion verwendet. Die Werte
1359 \begin{eqnarray*}
1360   w_{\mathrm{Basis}}        &=& 0 \\
1361   w_{\mathrm{Komparatoren}} &=& 2 \\
1362   w_{\mathrm{Schichten}}    &=& 1
1363 \end{eqnarray*}
1364 geben dem Einsparen eines Komparators ein höheres Gewicht als dem Einsparen
1365 einer Schicht.
1366
1367 \subsection{Selektion}
1368
1369 Als \emph{Selektion} wird der Vorgang bezeichnet, der zwei Individuen zufällig
1370 aus der Population auswählt. Sie werden im folgenden Schritt miteinander
1371 rekombiniert. Die Auswahl der Individuen erfolgt zufällig, aber nicht
1372 gleichverteilt. So sorgt die \emph{Selektion} dafür, dass bessere Individuen
1373 eine größere Wahrscheinlichkeit haben zur nächsten Generation beizutragen.
1374 Diese Ungleichbehandlung von Individuen verschiedener Güte ist der Grund für
1375 das Streben des Algorithmus nach besseren Lösungen.
1376
1377 Obwohl dieser Vorteil für gute Individuen intuitiv als sehr gering erscheint,
1378 passiert es häufig, dass die Ausnutzung \emph{(Exploitation)} überhand gewinnt
1379 und der Algorithmus vorschnell in Richtung eines lokalen Optimums optimiert.
1380
1381 Die in \textsc{SN-Evolution} implementierte Selektion eines Individuums lässt
1382 sich mit Pseudocode wie folgt beschreiben:
1383 \begin{verbatim}
1384   Gütesumme := 0
1385   Auswahl := (leer)
1386   
1387   für jedes Individuum in Population
1388   {
1389     reziproke Güte := 1.0 / Guete(Individuum)
1390     Wahrscheinlichkeit P := reziproke Güte / (Gütesumme + reziproke Güte)
1391     Gütesumme := Gütesumme + reziproke Güte
1392   
1393     mit Wahrscheinlichkeit P
1394     {
1395       Auswahl := Individuum
1396     }
1397   }
1398   gib Auswahl zurück
1399 \end{verbatim}
1400
1401 Diese Auswahl wird zweimal ausgeführt, um zwei Individuen für die
1402 Rekombination zu erhalten. Das heißt, dass die Individuen bei
1403 \textsc{SN-Evolution} stochastisch unabhängig voneinander ausgewählt werden.
1404
1405 \subsection{Rekombination}
1406 \label{sect:sn-evolution:rekombination}
1407
1408 Bei der Rekombination werden zwei Individuen --~hier Sortiernetzwerke~-- zu
1409 einer neuen Lösung kombiniert. Geeignete Mischer, um die beiden Netzwerke zu
1410 einem Netzwerk mit $2n$~Leitungen zusammenzufügen, sind zum Beispiel der {\em
1411 bitonen Mischer} (Abschnitt~\ref{sect:der_bitone_mischer}) und der
1412 \emph{Odd-Even}-Mischer (Abschnitt~\ref{sect:der_odd_even_mischer}),
1413 Anschließend werden $n$~Leitungen mit einem zufälligen $n$-Schnittmuster wie
1414 in Abschnitt~\ref{sect:leitungen_entfernen} beschrieben entfernt.
1415
1416 Dieses Verfahren hat den großen Vorteil, dass es die Sortiereigenschaft
1417 erhält. Entsprechend muss nicht aufwendig überprüft werden, ob das
1418 Komparatornetzwerk die Sortiereigenschaft besitzt. Der Nachteil ist, dass
1419 nicht alle Sortiernetzwerke auf diese Art und Weise erzeugt werden können.
1420
1421 \subsection{Mutation}
1422
1423 Zu einem vollständigen evolutionären Algorithmus gehört außerdem die Mutation
1424 --~eine zufällige Veränderung einer Lösung. Leider ist es nicht möglich ein
1425 Sortiernetzwerk zufällig zu verändern und dabei die Sortiereigenschaft zu
1426 erhalten. Selbst das \emph{Hinzufügen} eines zufälligen Komparators kann diese
1427 Eigenschaft zerstören.
1428
1429 Nach einer Mutation müsste man überprüfen, ob das neue Komparatornetzwerk die
1430 Sortiereigenschaft noch besitzt. Nach heutigem Wissenstand ist diese
1431 Überprüfung nur mit exponentiellem Aufwand möglich, etwa durch das
1432 Ausprobieren aller $2^n$~Bitmuster, wie in Abschnitt~\ref{sect:0-1-prinzip}
1433 beschrieben.
1434
1435 Um das Potenzial einer Mutation abzuschätzen wurde in \textsc{SN-Evolution}
1436 eine Überprüfung eingebaut: Unmittelbar vor dem Einfügen in die Population
1437 überprüft eine Funktion die Notwendigkeit jedes einzelnen Komparators. Dazu
1438 wird nacheinander jeder Komparator entfernt und überprüft, ob das verbleibende
1439 Netzwerk die Sortiereigenschaft noch besitzt. Trotz des hohen Rechenaufwands
1440 -- bei 16-Sortiernetzwerken sind gut 4~Millionen Tests notwendig, um alle
1441 Komparatoren zu überprüfen -- waren die Ergebnisse ernüchternd: Nach circa
1442 1~Million Iterationen mit 16-Sortiernetzwerken fand der so modifizierte
1443 Algorithmus keinen einzigen Komparator, den er hätte entfernen können. Daher
1444 wurde beim \textsc{SN-Evolution}-Algorithmus auf eine Mutation verzichtet.
1445
1446 \subsection[Bitoner Mischer]{Versuche mit dem bitonen Mischer}
1447
1448 Wenn \textsc{SN-Evolution} mit dem \emph{Odd-Even-Transpositionsort}-Netzwerk
1449 als Eingabe gestartet wird und in der Rekombinationsphase den \emph{bitonen
1450 Mischer} verwendet, gibt der Algorithmus \emph{effiziente} und in einigen
1451 Fällen \emph{schnelle} Sortiernetzwerke aus. Die Ergebnisse des
1452 \textsc{SN-Evolution}-Algorithmus unter Verwendung des \emph{bitonen Mischers}
1453 sind in Tabelle~\ref{tbl:sn-ev-bm-fast} zusammengefasst.
1454
1455 \begin{table}\label{tbl:sn-ev-bm-fast}
1456 \begin{center}
1457 \rowcolors{4}{black!5}{}
1458 \begin{tabular}{|r|r|r|r|r|}
1459 \hline
1460 Leitungen & \multicolumn{2}{l|}{\textsc{SN-EV} mit \bm{n}} & \multicolumn{2}{|l|}{\bs{n}} \\
1461 \cline{2-5}
1462     ($n$) & Komp. & Schichten & Komp. & Schichten \\
1463 \hline
1464         8 & \gcell  20 &         6 &  24 &  6 \\
1465         9 & \Gcell  26 &         8 &  28 &  8 \\
1466        10 & \gcell  31 & \gcell  8 &  33 &  9 \\
1467        11 & \Gcell  37 & \Gcell  9 &  39 & 10 \\
1468        12 & \gcell  42 & \gcell  9 &  46 & 10 \\
1469        13 & \Gcell  48 &        10 &  53 & 10 \\
1470        14 & \gcell  54 &        10 &  61 & 10 \\
1471        15 & \Gcell  61 &        10 &  70 & 10 \\
1472        16 & \gcell  67 &        10 &  80 & 10 \\
1473        17 & \Gcell  76 &        12 &  85 & 12 \\
1474        18 & \gcell  87 & \gcell 12 &  91 & 13 \\
1475        19 & \Gcell  93 & \Gcell 13 &  98 & 14 \\
1476        20 & \gcell 104 & \gcell 13 & 106 & 14 \\
1477        21 & \Gcell 109 & \Gcell 14 & 114 & 15 \\
1478        22 & \gcell 118 & \gcell 14 & 123 & 15 \\
1479        23 & \Gcell 129 & \Gcell 14 & 133 & 15 \\
1480        24 & \gcell 133 &        15 & 144 & 15 \\
1481 \hline
1482 \end{tabular}
1483 \caption{Übersicht über die Ergebnisse des \textsc{SN-Evolution}-Algorithmus
1484   unter Verwendung des \emph{bitonen Merge}-Netzwerks \bm{n}. Der Algorithmus
1485   wurde mit dem \emph{Odd-Even-Transpositionsort}-Netzwerk \oet{n} gestartet
1486   und nach 2.500.000 Iterationen beendet. Die Bewertungsfunktion nutzte die
1487   Konstanten $w_{\mathrm{Basis}} = 0$, $w_{\mathrm{Komparatoren}} = 1$,
1488   $w_{\mathrm{Schichten}} = n$.}
1489 \end{center}
1490 \end{table}
1491
1492 Alle Sortiernetzwerke, die von \textsc{SN-Evolution} in dieser Konfiguration
1493 gefunden wurden, waren \emph{effizienter} als das \emph{bitone
1494 Mergesort}-Netzwerk \bs{n}, das ebenfalls auf dem \emph{bitonen
1495 Merge}-Netzwerk \bm{n} beruht. Zum Beispiel benötigt das in
1496 Abbildung~\ref{fig:16-e1-bitonic-1296542566} dargestellte 16-Sortiernetzwerk
1497 67~Komparatoren, 13~Komparatoren weniger als \bs{n}.
1498
1499 \begin{figure}
1500   \begin{center}
1501     \input{images/16-e1-bitonic-1296542566.tex}
1502   \end{center}
1503   \caption{Sortiernetzwerk mit 16~Leitungen und 67~Komparatoren in
1504     10~Schichten. Das Netzwerk wurde von dem Algorithmus
1505     \textsc{SN-Evolution} unter Verwendung des \emph{bitonen Mischers}
1506     erzeugt.}
1507   \label{fig:16-e1-bitonic-1296542566}
1508 \end{figure}
1509
1510 Wenn die Gütefunktion so gewählt ist, dass sie schnelle Sortiernetzwerke
1511 bevorzugt, werden in einigen Fällen Netzwerke zurückgegeben, die
1512 \emph{schneller} und \emph{effizienter} als \bs{n} sind. Das
1513 19-Sortiernetzwerk in Abbildung~\ref{fig:19-e1-bm-fast} besitzt beispielsweise
1514 nur 13~Schichten und benötigt damit einen parallelen Schritt weniger als
1515 \bs{19}.
1516
1517 \begin{figure}
1518   \begin{center}
1519     \input{images/19-e1-bm-fast.tex}
1520   \end{center}
1521   \caption{Sortiernetzwerk mit 19~Leitungen und 93~Komparatoren in
1522     13~Schichten. Das Netzwerk wurde von dem Algorithmus \textsc{SN-Evolution}
1523     unter Verwendung des \emph{bitonen Mischers} erzeugt.}
1524   \label{fig:19-e1-bm-fast}
1525 \end{figure}
1526
1527 \subsection[Odd-Even-Mischer]{Versuche mit dem Odd-Even-Mischer}
1528
1529 Die folgenden Ergebnisse wurden erzielt, indem \textsc{SN-Evolution} mit dem
1530 \emph{Odd-Even-Transpositionsort}-Netzwerk als Eingabe gestartet wurde und in
1531 der Rekombinationsphase das \emph{Odd-Even-Merge}-Netzwerk verwendete. So
1532 erzeugt der Algorithmus entweder Sortiernetzwerke, die genauso schnell und
1533 effizient wie das \oes{n}-Netzwerk, oder Sortiernetzwerke, die schneller aber
1534 weniger effizient als das \oes{n}-Netzwerk sind. Die Ergebnisse von
1535 \textsc{SN-Evolution} mit dem \emph{Odd-Even}-Mischer sind in
1536 Tabelle~\ref{tbl:sn-ev-oem-fast} zusammengefasst.
1537
1538 \begin{table}\label{tbl:sn-ev-oem-fast}
1539 \begin{center}
1540 \rowcolors{4}{black!5}{}
1541 \begin{tabular}{|r|r|r|r|r|}
1542 \hline
1543 Leitungen & \multicolumn{2}{l|}{\textsc{SN-EV} mit \oem{n}} & \multicolumn{2}{|l|}{\oes{n}} \\
1544 \cline{2-5}
1545           & Komp. & Schichten & Komp. & Schichten \\
1546 \hline
1547         8 &   19 &         6 &         19 &         6 \\
1548         9 &   26 &         8 &         26 &         8 \\
1549        10 &   31 &         9 &         31 &         9 \\
1550        11 &   38 & \Gcell  9 & \Gcell  37 &        10 \\
1551        12 &   43 & \gcell  9 & \gcell  41 &        10 \\
1552        13 &   48 &        10 &         48 &        10 \\
1553        14 &   53 &        10 &         53 &        10 \\
1554        15 &   59 &        10 &         59 &        10 \\
1555        16 &   63 &        10 &         63 &        10 \\
1556        17 &   74 &        12 &         74 &        12 \\
1557        18 &   82 &        13 &         82 &        13 \\
1558        19 &   93 & \Gcell 13 & \Gcell  91 &        14 \\
1559        20 &   97 &        14 &         97 &        14 \\
1560        21 &  108 & \Gcell 14 & \Gcell 107 &        15 \\
1561        22 &  117 & \gcell 14 & \gcell 114 &        15 \\
1562        23 &  127 & \Gcell 14 & \Gcell 122 &        15 \\
1563        24 &  128 &        15 & \gcell 127 &        15 \\
1564 \hline
1565 \end{tabular}
1566 \caption{Übersicht über die Ergebnisse des \textsc{SN-Evolution}-Algorithmus
1567   unter Verwendung des \emph{Odd-Even-Merge}-Netzwerks \oem{n}. Der
1568   Algorithmus wurde mit dem \emph{Odd-Even-Transpositionsort}-Netzwerk \oet{n}
1569   gestartet und nach 2.500.000 Iterationen beendet. Die Bewertungsfunktion
1570   nutzte die Konstanten $w_{\mathrm{Basis}} = 0$, $w_{\mathrm{Komparatoren}} =
1571   1$, $w_{\mathrm{Schichten}} = n$.}
1572 \end{center}
1573 \end{table}
1574
1575 Im vorherigen Abschnitt wurde gezeigt, dass der
1576 \textsc{SN-Evolution}-Algorithmus unter Verwendung des \emph{bitonen Mischers}
1577 Sortiernetzwerke erzeugen kann, die effizienter als das rekursiv aus dem
1578 \emph{bitonen Mischer} aufgebaute \emph{bitone Mergesort}-Netzwerk sind.
1579 Dieses Ergebnis lies sich mit dem \emph{Odd-Even-Merge}-Netzwerk nicht
1580 erzielen. Die Sortiernetzwerke, die \textsc{SN-Evolution} unter Verwendung des
1581 \emph{Odd-Even-Merge}-Netzwerks findet, erreichen das
1582 \emph{Odd-Even-Mergesort}-Netzwerk bezüglich Effizienz, übertreffen es aber
1583 nicht. Ein Beispiel für ein entsprechendes Sortiernetzwerk ist in
1584 Abbildung~\ref{fig:16-e1-oem-fast} dargestellt.
1585
1586 \begin{figure}
1587   \begin{center}
1588     \input{images/16-e1-oem-fast.tex}
1589   \end{center}
1590   \caption{Sortiernetzwerk mit 16~Leitungen und 63~Komparatoren in
1591     10~Schichten. Das Netzwerk wurde von dem Algorithmus
1592     \textsc{SN-Evolution} unter Verwendung des \emph{Odd-Even}-Mischers
1593     erzeugt.}
1594   \label{fig:16-e1-oem-fast}
1595 \end{figure}
1596
1597 Mit einer Gütefunktion, die schnelle Sortiernetzwerke bevorzugt, ist es auch
1598 mit dem \emph{Odd-Even}-Mischer möglich, dass \textsc{SN-Evolution}
1599 Sortiernetzwerke zurück gibt, die schneller als \oes{n} sind. Beispielsweise
1600 benötigt das 19-Sortiernetzwerk, das in Abbildung~\ref{fig:19-e1-oem-fast}
1601 dargestellt ist, nur 13~Schichten, während \oes{19} 14~Schichten benötigt.
1602
1603 \begin{figure}
1604   \begin{center}
1605     \input{images/19-e1-oem-fast.tex}
1606   \end{center}
1607   \caption{Sortiernetzwerk mit 19~Leitungen und 93~Komparatoren in
1608     13~Schichten. Das Netzwerk wurde von dem Algorithmus \textsc{SN-Evolution}
1609     unter Verwendung des \emph{Odd-Even}-Mischers erzeugt.}
1610   \label{fig:19-e1-oem-fast}
1611 \end{figure}
1612
1613 %\begin{figure}
1614 %\begin{center}
1615 %\input{images/08-e2-1237993371.tex}
1616 %\end{center}
1617 %\caption{{\tt images/08-e2-1237993371.tex}: 19~Komparatoren in 6~Schichten}
1618 %\label{fig:08-e2-1237993371}
1619 %\end{figure}
1620 %
1621 %\begin{figure}
1622 %\begin{center}
1623 %\input{images/09-e2-1237997073.tex}
1624 %\end{center}
1625 %\caption{{\tt images/09-e2-1237997073.tex}: 25~Komparatoren in 8~Schichten}
1626 %\label{fig:09-e2-1237997073}
1627 %\end{figure}
1628 %
1629 %\begin{figure}
1630 %\begin{center}
1631 %\input{images/09-e2-1237999719.tex}
1632 %\end{center}
1633 %\caption{{\tt images/09-e2-1237999719.tex}: 25~Komparatoren in 7~Schichten}
1634 %\label{fig:09-e2-1237999719}
1635 %\end{figure}
1636 %
1637 %\begin{figure}
1638 %\begin{center}
1639 %\input{images/10-e2-1239014566.tex}
1640 %\end{center}
1641 %\caption{{\tt images/10-e2-1239014566.tex}: 29~Komparatoren in 8~Schichten}
1642 %\label{fig:10-e2-1239014566}
1643 %\end{figure}
1644
1645 \subsection{Zufälliger Mischer}
1646
1647 Die Ergebnisse der beiden vorhergehenden Abschnitte zeigen, dass für einige
1648 Leitungszahlen der \emph{bitone Mischer} und für andere Leitungszahlen der
1649 \emph{Odd-Even}-Mischer bessere Ergebnisse liefert. Beispielsweise hat das
1650 Netzwerk für $n = 18$ bei Verwendung des \emph{bitone Mischers} nur
1651 12~Schichten, bei Verwendung des \emph{Odd-Even}-Mischers hingegen nur
1652 82~Komparatoren. Daher liegt die Idee nahe, beide Mischer-Netzwerke zu nutzen,
1653 um das beste Ergebnis beider Konstruktionen zu erreichen.
1654 \textsc{SN-Evolution} kann zu diesem Zweck beim Zusammenfügen zweier
1655 Individuen zufällig zwischen dem \emph{bitonen Mischer} und dem
1656 \emph{Odd-Even}-Mischer wählen. Die Ergebnisse von \textsc{SN-Evolution} bei
1657 einer zufälligen Wahl des Mischers in der Rekombinationsphase sind in
1658 Tabelle~\ref{tbl:sn-ev-rnd-fast} zusammengefasst.
1659
1660 \begin{table}\label{tbl:sn-ev-rnd-fast}
1661 \begin{center}
1662 \rowcolors{4}{black!5}{}
1663 \begin{tabular}{|r|r|r|r|r|r|r|}
1664 \hline
1665 Leitungen & \multicolumn{2}{l|}{\textsc{SN-EV} mit \bm{n}}
1666           & \multicolumn{2}{l|}{\textsc{SN-EV} mit \oem{n}}
1667           & \multicolumn{2}{l|}{\textsc{SN-EV} mit Zufall} \\
1668 \cline{2-7}
1669     ($n$) & Komp. & Schichten & Komp. & Schichten & Komp. & Schichten \\
1670 \hline
1671         8 &         20 &         6 & \gcell  19 &         6 & \gcell  19 &         6 \\
1672         9 &         26 &         8 &         26 &         8 &         26 &         8 \\
1673        10 &         31 & \gcell  8 &         31 &         9 &         31 & \gcell  8 \\
1674        11 & \Gcell  37 &         9 &         38 &         9 & \Gcell  37 &         9 \\
1675        12 &         42 &         9 &         43 &         9 & \gcell  41 &         9 \\
1676        13 &         48 &        10 &         48 &        10 &         48 &        10 \\
1677        14 &         54 &        10 & \gcell  53 &        10 & \gcell  53 &        10 \\
1678        15 &         61 &        10 & \Gcell  59 &        10 & \Gcell  59 &        10 \\
1679        16 &         67 &        10 & \gcell  63 &        10 &         64 &        10 \\
1680        17 &         76 &        12 & \Gcell  74 &        12 & \Gcell  74 &        12 \\
1681        18 &         87 & \gcell 12 & \gcell  82 &        13 &         83 & \gcell 12 \\
1682        19 &         93 &        13 &         93 &        13 & \Gcell  92 &        13 \\
1683        20 &        104 & \gcell 13 & \gcell  97 &        14 &        101 & \gcell 13 \\
1684        21 &        109 &        14 &        108 &        14 & \Gcell 107 &        14 \\
1685        22 &        118 &        14 &        117 &        14 & \gcell 116 &        14 \\
1686        23 &        129 &        14 & \Gcell 127 &        14 &        128 &        14 \\
1687        24 &        133 &        15 & \gcell 128 &        15 &        130 &        15 \\
1688 \hline
1689 \end{tabular}
1690 \caption{Übersicht über die Ergebnisse des \textsc{SN-Evolution}-Algorithmus
1691   unter Verwendung der beiden Mischer-Netzwerke. Der Algorithmus wurde mit dem
1692   \emph{Odd-Even-Transpositionsort}-Netzwerk \oet{n} gestartet und nach
1693   2.500.000 Iterationen beendet. Die Bewertungsfunktion nutzte die Konstanten
1694   $w_{\mathrm{Basis}} = 0$, $w_{\mathrm{Komparatoren}} = 1$ und
1695   $w_{\mathrm{Schichten}} = n$.}
1696 \end{center}
1697 \end{table}
1698
1699 Bei einigen Leitungszahlen kann der Algorithmus durch die Verfügbarkeit beider
1700 Mi\-scher-Netzwerke Sortiernetzwerke zurückgeben, die effizienter als die
1701 vorherigen Ergebnisse sind. Beispielsweise ist das 19-Sortiernetzwerk in
1702 Abbildung~\ref{fig:19-e1-rnd-fast} mit 92~Komparatoren effizienter als die
1703 19-Sortiernetzwerke, die mit nur einem der beiden Mischer-Netzwerke erreicht
1704 wurden (Abbildungen~\ref{fig:19-e1-bm-fast} und~\ref{fig:19-e1-oem-fast}).
1705
1706 \begin{figure}
1707   \begin{center}
1708     \input{images/19-e1-rnd-fast.tex}
1709   \end{center}
1710   \caption{Sortiernetzwerk mit 19~Leitungen und 92~Komparatoren in
1711     13~Schichten. Das Netzwerk wurde von dem Algorithmus \textsc{SN-Evolution}
1712     unter Verwendung des \emph{bitonen Mischers} und des
1713     \emph{Odd-Even}-Mischers erzeugt.}
1714   \label{fig:19-e1-rnd-fast}
1715 \end{figure}
1716
1717 Die Ergebnisse anderer Leitungszahlen erreichen die Geschwindigkeit der
1718 Ergebnisse, die mit dem \emph{bitonen Mischer} erzielt wurden. Die Effizienz
1719 liegt zwischen den Ergebnissen, die mit dem \emph{bitonen Mischer} erzielt
1720 wurden, und den Ergebnissen, die mit dem \emph{Odd-Even}-Mischer erzielt
1721 wurden. Beispielsweise ist das 18-Sortiernetzwerk in
1722 Abbildung~\ref{fig:18-e1-rnd-fast} so schnell wie das Ergebnis, das mit dem
1723 \emph{bitonen Mischer} ausgegeben wurde. Mit 83~Komparatoren liegt die
1724 Effizienz des Sortiernetzwerks zwischen den Ergebnissen, die mit dem
1725 \emph{bitonen Mischer} (87~Komparatoren), beziehungsweise dem
1726 \emph{Odd-Even}-Mischer (82~Komparatoren) erreicht werden konnten.
1727
1728 \begin{figure}
1729   \begin{center}
1730     \input{images/18-e1-rnd-fast.tex}
1731   \end{center}
1732   \caption{Sortiernetzwerk mit 18~Leitungen und 83~Komparatoren in
1733     12~Schichten. Das Netzwerk wurde von dem Algorithmus \textsc{SN-Evolution}
1734     unter Verwendung des \emph{bitonen Mischers} und des
1735     \emph{Odd-Even}-Mischers erzeugt.}
1736   \label{fig:18-e1-rnd-fast}
1737 \end{figure}
1738
1739 In einigen Fällen hat \textsc{SN-Evolution} in dieser Konfiguration
1740 Sortiernetzwerke ausgegeben, die weniger effizient und genauso schnell wie die
1741 bisherigen Ergebnisse unter Verwendung des \emph{Odd-Even}-Mischers sind.
1742 Prinzipiell könnte der Algorithmus in jeder Iteration zufällig den
1743 \emph{Odd-Even}-Mischers auswählen, um die selektierten Individuen zu
1744 rekombinieren. Das heißt, das die Ergebnisse auch bei einer zufälligen Wahl
1745 des Mischer-Netzwerks theoretisch erreicht werden können. Allerdings sind
1746 unter Umständen mehr Iterationen notwendig, bis die gleiche Effizienz erreicht
1747 wird.
1748
1749 %\input{shmoo-aequivalenz.tex}
1750
1751 \newpage
1752 \section{Der \textsc{SN-Markov}-Algorithmus}
1753 \label{sect:markov}
1754
1755 Der evolutionäre \textsc{SN-Evolution}-Algorithmus aus dem vorherigen
1756 Abschnitt verwendet immer zwei zufällige Sortiernetzwerke („Individuen“) aus
1757 einer Population. Da die beiden „Eltern“ zufällig und unabhängig voneinander
1758 ausgewählt werden, kann es vorkommen, dass das selbe Sortiernetzwerk zweimal
1759 verwendet und mit sich selbst kombiniert wird.
1760
1761 Macht man diesen Spezialfall zum Regelfall, kombiniert das aktuelle Netzwerk
1762 \emph{immer} mit sich selbst und eliminiert anschließend die Hälfte aller
1763 Leitungen, lassen sich einige interessante Beobachtungen anstellen. Netzwerke,
1764 die aus einem Netzwerk $S_0$ durch die beschriebene Kombination von $S_0$ mit
1765 sich selbst und anschließendem Eliminieren der Hälfte der Leitungen hervorgehen
1766 können, heißen \emph{Nachfolger} von $S_0$.
1767
1768 Beim beschriebenen Vorgehen kann man die Sortiernetzwerke als Knoten in einem
1769 (gerichteten) Graphen betrachten. Zwei Knoten $V_0$ und $V_1$, die zwei
1770 Sortiernetzwerke $S_0$ und $S_1$ repräsentieren, sind genau dann mit einer
1771 Kante ${E_{0,1} = (V_0, V_1)}$ verbunden, wenn $S_1$ ein \emph{Nachfolger} von
1772 $S_0$ ist, das heißt, dass $S_1$ durch die Rekombination von $S_0$ mit sich
1773 selbst erzeugt werden kann.
1774
1775 Wie in Abschnitt~\ref{sect:anzahl_schnittmuster} beschrieben, ist die Anzahl
1776 der \emph{unterschiedlichen} Schnittmuster und damit die Anzahl der Nachfolger
1777 sehr groß. Bei den untersuchten 16-Sortiernetzwerken lag die Anzahl der
1778 Nachfolger zwar noch unter 20.000, bei den untersuchten
1779 32-Sortier\-netz\-werken wurden jedoch bereits bis zu $2,6 \cdot 10^8$
1780 unterschiedliche Schnittmuster geschätzt.
1781
1782 Der Algorithmus {\sc SN-Markov} legt auf diesem Nachfolger-Graph einen
1783 zufälligen Weg (englisch: \textit{random walk}) zurück. Er startet auf einem
1784 gegebenen Sortiernetzwerk. Um von einem Sortiernetzwerk zum Nächsten zu
1785 gelangen, rekombiniert der Algorithmus das aktuelle Sortiernetzwerk mit sich
1786 selbst und erhält so einen zufälligen Nachfolger. In Pseudocode lässt sich der
1787 Algorithmus wie folgt beschreiben:
1788
1789 \begin{verbatim}
1790   Netzwerk := Eingabe
1791   
1792   für n Iterationen
1793   {
1794     Nachfolger := kombiniere (Netzwerk, Netzwerk)
1795     Netzwerk   := Nachfolger
1796   }
1797   
1798   gib Netzwerk zurück
1799 \end{verbatim}
1800
1801 Die Graphen in Abbildung~\ref{fig:markov-comparators} zeigen die Anzahl der
1802 Komparatoren der Sortiernetzwerke, die \textsc{SN-Markov} auf seinem
1803 zufälligen Pfad durchläuft (rot). Für jeden Graphen wurde der
1804 \textsc{SN-Markov}-Algorithmus auf einem entsprechenden
1805 \emph{Odd-Even-Transpositionsort}-Netzwerk gestartet und hat mindestens
1806 1.000.000~Iterationen durchlaufen. In jedem Schritt wurde die Anzahl der
1807 Komparatoren des Sortiernetzwerks bestimmt und ein entsprechender Zähler
1808 erhöht. In Abbildung~\ref{fig:markov-comparators} ist die resultierende
1809 prozentuale Verteilung zu sehen.
1810
1811 Ebenfalls in die Graphen der Abbildung~\ref{fig:markov-comparators}
1812 eingezeichnet ist eine \emph{Gamma-Verteilung} (grün), die die gemessenen
1813 Daten gut annähert. Die Gamma-Verteilung verwendet einen Offset~$\delta$, der
1814 um Eins kleiner als die kleinste erreichte Komparatorzahl gewählt wurde.
1815 Beispielsweise war die kleinste erreichte Komparatorzahl bei
1816 16-Sortiernetzwerken~63, entsprechend wurde der Offset $\delta = 63 - 1$
1817 gesetzt und die Gamma-Verteilung $g(x - 62)$ eingezeichnet. Die Parameter $k$
1818 und $\theta$, die eine Gamma-Verteilung charakterisieren, wurden mit einem
1819 Fitting-Algorithmus bestimmt. Der konkrete Offset ist als Parameter~$\delta$
1820 unter den Graphen angegeben.
1821
1822 \begin{figure}
1823   \centering
1824   \subfigure[12 Leitungen, $k = 8,267$, $\theta = 0,962$, $\delta = 40$]{\includegraphics[viewport=0 0 425 262,width=7cm]{images/markov-comparators-12-pct.pdf}}
1825   \subfigure[14 Leitungen, $k = 9,522$, $\theta = 0,867$, $\delta = 52$]{\includegraphics[viewport=0 0 425 262,width=7cm]{images/markov-comparators-14-pct.pdf}}
1826   \subfigure[16 Leitungen, $k = 17,939$, $\theta = 1,091$, $\delta = 62$]{\includegraphics[viewport=0 0 425 262,width=7cm]{images/markov-comparators-16-pct.pdf}}
1827   \subfigure[18 Leitungen, $k = 10,724$, $\theta = 0,766$, $\delta = 81$]{\includegraphics[viewport=0 0 425 262,width=7cm]{images/markov-comparators-18-pct.pdf}}
1828   \caption{Anzahl der Komparatoren von Sortiernetzwerken,
1829   die von {\sc SN-Markov} durchlaufen wurden (rot). Ebenfalls eingezeichnet
1830   ist jeweils eine \emph{Gamma-Verteilung} (grün), die eine gute Näherung der
1831   gemessenen Daten darstellt.}
1832   \label{fig:markov-comparators}
1833 \end{figure}
1834
1835 \begin{figure}
1836   \begin{center}
1837     \includegraphics[viewport=0 0 425 262,width=15cm]{images/comparison-comparators-16.pdf}
1838   \end{center}
1839   \caption{Anzahl der Komparatoren, die 16-Sortiernetzwerke von
1840   \textsc{SN-Markov} und \textsc{SN-Evolution} (mit dem
1841   \emph{Odd-Even}-Mischer und dem \emph{bitonen Mischer}) besaßen.}
1842   \label{fig:comparison-comparators}
1843 \end{figure}
1844
1845 Der Graph in Abbildung~\ref{fig:comparison-comparators} zeigt, dass der
1846 \textsc{SN-Markov}-Algorithmus nicht schlechter ist als der
1847 \textsc{SN-Evolution}-Algo\-rith\-mus. Analog zu dem Versuch mit
1848 \textsc{SN-Markov}, wurde beim \textsc{SN-Evolution}-Algorithmus die Anzahl
1849 der Komparatoren jedes neuen Individuums ermittelt und gespeichert. Als
1850 Startnetzwerk diente bei beiden Algorithmen das
1851 \emph{Odd-Even-Transpositionsort}-Netzwerk \oet{16}. Der Graph zeigt auf der
1852 x-Achse die Anzahl der Komparatoren, auf der y-Achse die Häufigkeit, mit der
1853 ein Sortiernetzwerk mit dieser Komparatorzahl durch die Rekombination erzeugt
1854 wurde. Die Ergebnisse von \textsc{SN-Evolution} unterscheiden sich außerdem je
1855 nach verwendetem Mischer-Netzwerk -- \oem{32}, beziehungsweise \bm{32}.
1856
1857 Sowohl der \textsc{SN-Markov}-Algorithmus, der das
1858 \emph{Odd-Even-Merge}-Netzwerk verwendet, als auch \textsc{SN-Evolution} mit
1859 \oem{32} erreichen eine Komparatorzahl von~63 und finden Sortiernetzwerke, die
1860 bezüglich Effizienz und Geschwindigkeit identisch zu \oes{16} sind.
1861 Interessanterweise erzeugt \textsc{SN-Markov} derartige Netzwerke häufiger:
1862 Während nur $0,000017 \%$ der Individuen von \textsc{SN-Evolution} mit
1863 63~Komparatoren auskamen, ist die Rate bei \textsc{SN-Markov} mit $0,000335
1864 \%$ rund 20~mal höher.
1865
1866 Erwartungsgemäß sind die besten Netzwerke, die \textsc{SN-Evolution} mit dem
1867 \emph{bitonen Mischer} findet, aus 67~Komparatoren aufgebaut. Überraschend ist
1868 jedoch, dass in dieser Konfiguration Sortiernetzwerke auftreten können, die
1869 mehr Komparatoren besitzen als \emph{Odd-Even-Transpositionsort}. \oet{16}
1870 ist aus 120~Komparatoren aufgebaut. Bei dem Lauf, der die Daten für
1871 Abbildung~\ref{fig:comparison-comparators} lieferte, trat auch jeweils ein
1872 Sortiernetzwerk mit 121 und 124~Komparatoren auf. Dass Sortiernetzwerke mit so
1873 vielen Komparatoren im Verlauf des Experiments selbst nach über 100~Millionen
1874 Iterationen nicht noch einmal erzeugt wurden, ist vermutlich ein Phänomen, das
1875 mit der Initialisierung durch das \emph{Odd-Even-Transpositionsort}-Netzwerk
1876 zusammenhängt.
1877
1878 %\begin{figure}
1879 %  \begin{center}
1880 %  \includegraphics[viewport=0 0 425 262,width=15cm]{images/markov-comparators-14-pct.pdf}
1881 %  \end{center}
1882 %  \caption{Anzahl der Komparatoren von Sortiernetzwerken (mit 14~Leitungen),
1883 %  die von {\sc SN-Markov} durchlaufen wurden. Grün eingezeichnet ist die
1884 %  \emph{Gamma-Verteilung} $f(x - 52)$ mit $k = 9,522$ und $\theta = 0,867$.}
1885 %  \label{fig:markov-comparators-14}
1886 %\end{figure}
1887 %
1888 %\begin{figure}
1889 %  \begin{center}
1890 %  \includegraphics[viewport=0 0 425 262,width=15cm]{images/markov-comparators-16-pct.pdf}
1891 %  \end{center}
1892 %  \caption{Anzahl der Komparatoren von Sortiernetzwerken (mit 16~Leitungen),
1893 %  die von {\sc SN-Markov} durchlaufen wurden. Grün eingezeichnet ist die
1894 %  \emph{Gamma-Verteilung} $f(x - 62)$ mit $k = 17,939$ und $\theta = 1,091$.}
1895 %  \label{fig:markov-comparators-16}
1896 %\end{figure}
1897 %
1898 %\begin{figure}
1899 %  \begin{center}
1900 %  \includegraphics[viewport=0 0 425 262,width=15cm]{images/markov-comparators-18-pct.pdf}
1901 %  \end{center}
1902 %  \caption{Anzahl der Komparatoren von Sortiernetzwerken (mit 18~Leitungen),
1903 %  die von {\sc SN-Markov} durchlaufen wurden. Grün eingezeichnet ist die
1904 %  \emph{Gamma-Verteilung} $f(x - 81)$ mit $k = 10,724$ und $\theta = 0,766$.}
1905 %  \label{fig:markov-comparators-18}
1906 %\end{figure}
1907
1908 %\begin{figure}
1909 %  \begin{center}
1910 %  \includegraphics[viewport=0 0 425 262,width=15cm]{images/markov-cycles-16.pdf}
1911 %  \end{center}
1912 %  \caption{Zyklen, die beim \textit{Random Walk} des
1913 %  \textsc{SN-Markov}-Algorithmus detektiert wurden. Auf der x-Achse sind die
1914 %  Anzahl der Schritte, die \textsc{SN-Markov} zurückgelegt hat, auf der
1915 %  y-Achse die Längen der gefundenen Zyklen aufgetragen. Das initiale
1916 %  Start-Sortiernetzwerk war $\operatorname{OET}(16)$.}
1917 %  \label{fig:markov-cycles-16}
1918 %\end{figure}
1919
1920 \newpage
1921 \section[\textsc{SN-Evolution-Cut}]{Der \textsc{SN-Evolution-Cut}-Algorithmus}
1922 \label{sect:sn-evolution-cut}
1923
1924 Das Programm \textsc{SN-Evolution-Cut} implementiert einen evolutionären
1925 Algorithmus, der zu einem gegebenen Sortiernetzwerk und einer gewünschten
1926 Leitungszahl ein Schnittmuster sucht, dass ein Sortiernetzwerk mit einer
1927 möglichst geringen Anzahl von Komparatoren und Schichten ergibt. Zur Bewertung
1928 von Sortiernetzwerken siehe auch Abschnitt~\ref{sect:sn-evolution:bewertung}.
1929
1930 Der \textsc{SN-Evolution-Cut}-Algorithmus verwendet \emph{Schnittmuster}, die
1931 in Abschnitt~\ref{sect:anzahl_schnittmuster} definiert wurden, als Individuen.
1932 Ein Individuum besteht aus einer Liste von $n$~Zahlen, die entweder 1, $-1$
1933 oder 0 sind. Dieser Werte entsprechen Maximum, Minimum und unbelegt. Bei einem
1934 $k$-Schnittmuster sind genau $k$ Zahlen ungleich Null.
1935
1936 Um zwei Individuen zu rekombinieren werden die ersten $r$~Werte des einen
1937 Schnittmusters und die letzten ${n-r}$~Schnitte des zweiten Schnittmusters
1938 verwendet. $r$ ist eine Zufallsvariable mit $0 \leqq r \leqq n$. Anschließend
1939 werden zufällig Werte auf Null beziehungsweise 1 oder $-1$ gesetzt, um die
1940 Anzahl der Schnitte zu korrigieren.
1941
1942 Die Mutation vertauscht entweder die Werte von zwei zufälligen Positionen oder
1943 multipliziert den Wert einer Leitung mit $-1$, um die Schnittrichtung zu
1944 invertieren.
1945
1946 Die Eingabe für \textsc{SN-Evolution-Cut} ist ein $n$-Sortiernetzwerk und eine
1947 Zahl $k$, $1 \leqq k < n$, die angibt wie viele Leitungen entfernt werden
1948 sollen. Der Rückgabewert des \textsc{SN-Evolution-Cut}-Algorithmus ist ein
1949 \emph{$k$-Schnittmuster}. Wird das Schnittmuster auf das Sortiernetzwerk, mit
1950 dem der Algorithmus gestartet wurde, angewendet, entsteht ein möglichst
1951 schnelles und effizientes Sortiernetzwerk mit $m = n - k$ Leitungen. Da mit
1952 dem Eingabe-Netzwerk und dem zurückgegebenen $k$-Schnittmuster das
1953 $m$-Sortiernetzwerk eindeutig bestimmt ist, werden im Folgenden sowohl das
1954 $k$-Schnittmuster als auch das $m$-Sortiernetzwerk als Ausgabe von
1955 \textsc{SN-Evolution-Cut} bezeichnet.
1956
1957 \subsection[Bitones Mergesort-Netzwerk]{Versuche mit dem bitonen Mergesort-Netzwerk}
1958 \label{sect:sn-evolution-cut:bs}
1959
1960 % Effizienz
1961
1962 Wenn der \textsc{SN-Evolution-Cut}-Algorithmus mit dem \emph{bitonen
1963 Mergesort}-Netzwerk \bs{n} gestartet wird und $k$~Leitungen entfernen soll,
1964 ergeben die gefundenen Schnittmuster in vielen Fällen effizientere Netzwerke
1965 als \bs{n-k}. Wird \textsc{SN-Evolution-Cut} beispielsweise mit \bs{22} und $k
1966 = 6$ gestartet, resultiert das gefundene Schnittmuster in einem
1967 Sortiernetzwerk mit 67~Komparatoren, 13~Komparatoren weniger als \bs{16}
1968 benötigt. Eines der Sortiernetzwerke, die auf diese Art und Weise generiert
1969 wurde, ist in Abbildung~\ref{fig:16-ec-from-bs22} zu sehen.
1970
1971 % Beispiel Effizienz
1972
1973 \begin{figure}
1974   \begin{center}
1975     \input{images/16-ec-from-bs22.tex}
1976   \end{center}
1977   \caption{Sortiernetzwerk mit 16~Leitungen und 67~Komparatoren in
1978     10~Schichten. Das Netzwerk wurde von dem Algorithmus
1979     \textsc{SN-Evolution-Cut} aus dem \emph{bitonen Mergesort}-Netzwerk
1980     $\operatorname{BS}(22)$ durch das 6-Schnittmuster $\operatorname{MIN}(4,
1981     10, 17)$, $\operatorname{MAX}(7, 15, 20)$ erzeugt.}
1982   \label{fig:16-ec-from-bs22}
1983 \end{figure}
1984
1985 Eine Übersicht über die Effizienz der Ergebnisse, die mit dem \emph{bitonen
1986 Mergesort}-Netzwerk als Eingabe für \textsc{SN-Evolution-Cut} erzielt wurden,
1987 gibt Tabelle~\ref{tbl:ec-bs-efficiency}. \textsc{SN-E\-vo\-lu\-tion-Cut} wurde
1988 mit \bs{n}, $n = 9 \dots 24$ und $k = 1 \dots (n-8)$ gestartet. Die Konstanten
1989 der Bewertungsfunktion waren $w_{\mathrm{Basis}} = 0$,
1990 $w_{\mathrm{Komparatoren}} = 1$ und $w_{\mathrm{Schichten}} = n$. In jeder
1991 Zeile befinden sich die Ergebnisse für ein Eingabenetzwerk, in den Spalten
1992 befinden sich die Ergebnisse für eine Leitungszahl $m=n-k$ des
1993 Ausgabenetzwerks. In den Zellen stehen jeweils die Anzahl der Komparatoren des
1994 resultierenden Netzwerks. Die letzte Zeile enthält die Anzahl der
1995 Komparatoren, die \bs{m} benötigt, um die Ergebnisse besser einordnen zu
1996 können.
1997
1998 \begin{table}
1999   \begin{center}
2000     \rowcolors{2}{black!5}{}
2001     \begin{tabular}{|r|rrrrrrrrrrrrrrrr|}
2002     \hline
2003        &  8 &  9 & 10 & 11 & 12 & 13 & 14 & 15 & 16 & 17 & 18 & 19 &  20 &  21 &  22 &  23 \\
2004     \hline
2005     9  & 21 &    &    &    &    &    &    &    &    &    &    &    &     &     &     &     \\
2006     10 & 20 & 27 &    &    &    &    &    &    &    &    &    &    &     &     &     &     \\
2007     11 & 20 & 27 & 32 &    &    &    &    &    &    &    &    &    &     &     &     &     \\
2008     12 & 20 & 26 & 32 & 39 &    &    &    &    &    &    &    &    &     &     &     &     \\
2009     13 & 20 & 26 & 32 & 39 & 45 &    &    &    &    &    &    &    &     &     &     &     \\
2010     14 & 20 & 26 & 32 & 39 & 45 & 53 &    &    &    &    &    &    &     &     &     &     \\
2011     15 & 20 & 26 & 32 & 39 & 45 & 53 & 61 &    &    &    &    &    &     &     &     &     \\
2012     16 & 20 & 26 & 32 & 39 & 45 & 53 & 61 & 70 &    &    &    &    &     &     &     &     \\
2013     17 & 20 & 26 & 32 & 37 & 43 & 50 & 57 & 65 & 74 &    &    &    &     &     &     &     \\
2014     18 & 20 & 26 & 31 & 37 & 43 & 49 & 56 & 63 & 71 & 82 &    &    &     &     &     &     \\
2015     19 & 20 & 26 & 31 & 37 & 43 & 48 & 55 & 62 & 70 & 79 & 88 &    &     &     &     &     \\
2016     20 & 20 & 26 & 32 & 37 & 44 & 48 & 55 & 61 & 68 & 77 & 86 & 95 &     &     &     &     \\
2017     21 & 20 & 26 & 32 & 37 & 44 & 48 & 55 & 61 & 68 & 77 & 85 & 94 & 103 &     &     &     \\
2018     22 & 20 & 26 & 31 & 37 & 42 & 48 & 54 & 61 & 67 & 77 & 84 & 93 & 102 & 112 &     &     \\
2019     23 & 20 & 26 & 31 & 37 & 42 & 48 & 54 & 61 & 68 & 76 & 84 & 93 & 102 & 112 & 122 &     \\
2020     24 & 20 & 26 & 32 & 37 & 42 & 48 & 54 & 61 & 68 & 76 & 84 & 93 & 102 & 112 & 122 & 133 \\
2021     \hline
2022 \bs{m} & 24 & 28 & 33 & 39 & 46 & 53 & 61 & 70 & 80 & 85 & 91 & 98 & 106 & 114 & 123 & 133 \\
2023     \hline
2024     \end{tabular}
2025   \end{center}
2026   \caption{Anzahl der Komparatoren der Ergebnisse von
2027     \textsc{SN-Evolution-Cut} mit verschiedenen Größen des \emph{bitonen
2028     Mergesort}-Netzwerks und unterschiedlichen Werten für~$k$. Jede Zeile gibt
2029     die Ergebnisse für ein Eingabenetzwerk \bs{n} an, jede Spalte enthält die
2030     Ergebnisse für $m=n-k$, die Anzahl der Leitungen des Ausgabenetzwerks.}
2031   \label{tbl:ec-bs-efficiency}
2032 \end{table}
2033
2034 Zu sehen ist, dass jedes einzelne Ergebnis von \textsc{SN-Evolution-Cut}
2035 mindestens so effizient wie das \emph{bitone Mergesort}-Netzwerk mit der
2036 gleichen Leitungszahl ist. Außerdem enthält jede Spalte (mit Ausnahme von
2037 $m=23$) ein Ergebnis, das effizienter als \bs{m} ist.
2038
2039 In zahlreichen Fällen reicht das Entfernen einer einzigen Leitung aus, um ein
2040 effizientes Ergebnis zu erzielen. Das Ergebnis, das \textsc{SN-Evolution-Cut}
2041 gestartet mit \bs{20} und $k = 1$ erreicht, benötigt mit 95~Komparatoren
2042 3~weniger als \bs{19}.
2043
2044 Bei anderen Größen ergeben erst größere~$k$ effiziente Sortiernetzwerke,
2045 beispielsweise bei $m = 10$: erst für $n = 18$, $k = 8$ wird ein
2046 Sortiernetzwerk mit 31~Komparatoren gefunden.
2047
2048 % Geschwindigkeit
2049
2050 Bei einigen Werten für die Ziel-Leitungsanzahl $m$ kann der
2051 \textsc{SN-Evolution-Cut}-Algorithmus Ergebnisse erzielen, die schneller als
2052 das entsprechende \emph{bitone Mergesort}-Netzwerk \bs{m} sind. In
2053 Tabelle~\ref{tbl:ec-bs-speed} ist die Anzahl der Schichten, die die Ergebnisse
2054 von \textsc{SN-Evolution-Cut} benötigen, um die Eingabe zu sortieren,
2055 aufgelistet. Jede Zeile enthält die Ergebnisse für ein Eingabenetzwerk \bs{n},
2056 jede Spalte enthält die Ergebnisse für eine Ziel-Leitungszahl $m = n-k$. Die
2057 Zellen enthalten die Anzahl der Schichten des jeweiligen Ergebnis-Netzwerks.
2058
2059 \begin{table}
2060   \begin{center}
2061     \rowcolors{2}{black!5}{}
2062 \begin{tabular}{|r|rrrrrrrrrrrrrrrr|}
2063 \hline
2064     &   8 &   9 &  10 &  11 &  12 &  13 &  14 &  15 &  16 &  17 &  18 &  19 &  20 &  21 &  22 &  23 \\
2065 \hline
2066   9 &   6 &     &     &     &     &     &     &     &     &     &     &     &     &     &     &     \\
2067  10 &   6 &   8 &     &     &     &     &     &     &     &     &     &     &     &     &     &     \\
2068  11 &   6 &   8 &   9 &     &     &     &     &     &     &     &     &     &     &     &     &     \\
2069  12 &   6 &   8 &   9 &  10 &     &     &     &     &     &     &     &     &     &     &     &     \\
2070  13 &   6 &   8 &   9 &  10 &  10 &     &     &     &     &     &     &     &     &     &     &     \\
2071  14 &   6 &   8 &   9 &  10 &  10 &  10 &     &     &     &     &     &     &     &     &     &     \\
2072  15 &   6 &   8 &   9 &  10 &  10 &  10 &  10 &     &     &     &     &     &     &     &     &     \\
2073  16 &   6 &   8 &   9 &  10 &  10 &  10 &  10 &  10 &     &     &     &     &     &     &     &     \\
2074  17 &   6 &   8 &   8 &   9 &  10 &  10 &  10 &  10 &  10 &     &     &     &     &     &     &     \\
2075  18 &   6 &   8 &   8 &   9 &   9 &  10 &  10 &  10 &  10 &  12 &     &     &     &     &     &     \\
2076  19 &   6 &   8 &   8 &   9 &   9 &  10 &  10 &  10 &  10 &  12 &  13 &     &     &     &     &     \\
2077  20 &   6 &   8 &   8 &   9 &   9 &  10 &  10 &  10 &  10 &  12 &  13 &  14 &     &     &     &     \\
2078  21 &   6 &   8 &   8 &   9 &   9 &  10 &  10 &  10 &  10 &  12 &  13 &  14 &  14 &     &     &     \\
2079  22 &   6 &   8 &   8 &   9 &   9 &  10 &  10 &  10 &  10 &  12 &  13 &  14 &  14 &  15 &     &     \\
2080  23 &   6 &   8 &   8 &   9 &   9 &  10 &  10 &  10 &  10 &  12 &  13 &  14 &  14 &  15 &  15 &     \\
2081  24 &   6 &   8 &   8 &   9 &   9 &  10 &  10 &  10 &  10 &  12 &  13 &  14 &  14 &  15 &  15 &  15 \\
2082 \hline
2083 \bs{m}& 6 &   8 &   9 &  10 &  10 &  10 &  10 &  10 &  10 &  12 &  13 &  14 &  14 &  15 &  15 &  15 \\
2084 \hline
2085 \end{tabular}
2086   \end{center}
2087   \caption{Anzahl der Schichten der Ergebnisse von
2088     \textsc{SN-Evolution-Cut} mit verschiedenen Größen des \emph{bitonen
2089     Mergesort}-Netzwerks und unterschiedlichen Werten für~$k$. Jede Zeile gibt
2090     die Ergebnisse für ein Eingabenetzwerk \bs{n} an, jede Spalte enthält die
2091     Ergebnisse für $m=n-k$, die Anzahl der Leitungen des Ausgabenetzwerks.}
2092   \label{tbl:ec-bs-speed}
2093 \end{table}
2094
2095 Für die Ziel-Leitungszahlen 9, 10 und 11 wurden Schnittmuster gefunden, die
2096 schnelle Sortiernetzwerke erzeugen. Beispiele für schnelle Sortiernetzwerke,
2097 die mit den von \textsc{SN-Evolution-Cut} ausgegebenen Schnittmustern erzeugt
2098 werden können, sind in Abbildung~\ref{fig:ec-bs-fast_networks} dargestellt.
2099
2100 % Beispiel Geschwindigkeit
2101
2102 \begin{figure}
2103   \centering
2104   \subfigure[10-Sortiernetzwerk aus 31~Komparatoren in 8~Schichten. Das
2105   Netzwerk wurde von \textsc{SN-Evolution-Cut} aus \bs{19} erzeugt.]{\input{images/10-ec-from-bs19-fast.tex}\label{fig:10-ec-from-bs19-fast}}
2106   \subfigure[11-Sortiernetzwerk aus 37~Komparatoren in 9~Schichten. Das
2107   Netzwerk wurde von \textsc{SN-Evolution-Cut} aus \bs{18} erzeugt.]{\input{images/11-ec-from-bs18-fast.tex}\label{fig:11-ec-from-bs18-fast}}
2108   \subfigure[12-Sortiernetzwerk aus 42~Komparatoren in 9~Schichten. Das
2109   Netzwerk wurde von \textsc{SN-Evolution-Cut} aus \bs{22} erzeugt.]{\input{images/12-ec-from-bs22-fast.tex}\label{fig:12-ec-from-bs22-fast}}
2110   \subfigure[19-Sortiernetzwerk aus 92~Komparatoren in 13~Schichten. Das
2111   Netzwerk wurde von \textsc{SN-Evolution-Cut} aus \bs{37} erzeugt.]{\input{images/19-ec-from-bs37-fast.tex}\label{fig:19-ec-from-bs37-fast}}
2112   \caption{Für einige Ziel-Leitungszahlen, unter anderem $m \in \{10, 11,
2113   12, 19\}$, kann der \textsc{SN-Evolution-Cut}-Algorithmus Sortiernetzwerke
2114   erzeugen, die \emph{schneller} und \emph{effizienter} als \bs{m} sind.}
2115   \label{fig:ec-bs-fast_networks}
2116 \end{figure}
2117
2118 Bei der Betrachtung der Effizienz wurde festgestellt, dass oft schon das
2119 Entfernen einer einzigen Leitung zu eines effizienteren Ergebnis als \bs{m}
2120 führt. Bei der Geschwindigkeit ist die Anzahl der Leitungen, die entfernt
2121 werden müssen, um schnellere Netzwerke zu erzielen, größer. Um eine Schicht
2122 einzusparen waren bei $m = 10$ und $m = 11$ $k = 6$ Schnitte notwendig. Bei $m
2123 = 9$ war sogar ein 7-Schnittmuster notwendig, um die Anzahl der Schichten zu
2124 reduzieren. Für schnelle \emph{und} effiziente Netzwerke musste $k$ teilweise
2125 noch größer gewählt werden.
2126
2127 % Detaillierte Betrachtung fuer m = 19
2128
2129 Die Effizienz und Geschwindigkeit der Sortiernetzwerke, die von
2130 \textsc{SN-Evolution-Cut} aus dem \emph{bitonen Mergesort}-Netzwerk erzeugten
2131 werden, ist für $m = 19$ und $n = 20 \dots 38$ ($k = 1 \dots 19$) in
2132 Tabelle~\ref{tbl:ec-bs-19} aufgelistet. Erst, wenn $k \geqq 6$ ist, wird im
2133 Vergleich zu \bs{19} eine Schicht eingespart. Für $n = 36$ ($k = 17$) und $n =
2134 37$ ($k = 18$) werden Sortiernetzwerke ausgegeben, die schneller als \bs{19}
2135 und \oes{19} sind und nur einen Komparator mehr als \oes{19} benötigen. Ein
2136 Beispiel für ein solches Netzwerk ist in
2137 Abbildung~\ref{fig:19-ec-from-bs37-fast} zu sehen.
2138
2139 \begin{table}
2140   \begin{center}
2141     \rowcolors{2}{black!5}{}
2142     \begin{tabular}{|r|r|r|}
2143     \hline
2144     $n$ & Komp. & Schichten \\
2145     \hline
2146           20 & 95 & 14 \\
2147           21 & 94 & 14 \\
2148           22 & 93 & 14 \\
2149           23 & 93 & 14 \\
2150           24 & 93 & 14 \\
2151           25 & 96 & 13 \\
2152           26 & 96 & 13 \\
2153           27 & 96 & 13 \\
2154           28 & 96 & 13 \\
2155           29 & 95 & 13 \\
2156           30 & 96 & 13 \\
2157           31 & 95 & 13 \\
2158           32 & 96 & 13 \\
2159           33 & 93 & 13 \\
2160           34 & 94 & 13 \\
2161           35 & 93 & 13 \\
2162           \rowcolor{green!10}
2163           36 & 92 & 13 \\
2164           \rowcolor{green!10!white!95!black}
2165           37 & 92 & 13 \\
2166           38 & 93 & 13 \\
2167     \hline
2168     \bs{19}  & 98 & 14 \\
2169     \oes{19} & 91 & 14 \\
2170     \hline
2171     \end{tabular}
2172   \end{center}
2173   \caption{Anzahl der Komparatoren und Schichten von 19-Sortiernetzwerken, die
2174     von \textsc{SN-Evolution-Cut} aus \bs{n}, $n = 20, \dots, 38$ erzeugt
2175     wurden. Für $k \geqq 6$ ergeben sich Sortiernetzwerke, die schneller als
2176     \bs{19} sind. Mit $k \in \{14, 16, 19\}$ erreichen die Ergebnisse mit
2177     13~Schichten die Effizienz der vorherigen
2178     Ergebnisse mit 14~Schichten, mit $k = 17$ und $k = 18$ wird diese
2179     Effizienz noch übertroffen. Ein 19-Sortiernetzwerk, das aus \bs{37}
2180     auf diese Art erzeugt wurde, ist in
2181     Abbildung~\ref{fig:19-ec-from-bs37-fast} dargestellt.}
2182   \label{tbl:ec-bs-19}
2183 \end{table}
2184
2185 % 2-er Potenz
2186
2187 \textit{Moritz Mühlenthaler} und \textit{Rolf Wanka} zeigen in~\cite{MW2010},
2188 wie ein \emph{bitoner Mischer} $\bm{n = 2^d}$, der nach Batchers Methode
2189 konstruiert wurde, durch systematisches Entfernen von Leitungen in einen
2190 ebenfalls bitonen Mischer mit der Hälfte der Leitungen transformiert werden
2191 kann, so dass dieser alternative Mischer im Vergleich zu $\bm{\frac{n}{2} =
2192 2^{d-1}}$ Komparatoren einspart.
2193
2194 Basierend auf diesen alternativen Mischern geben \textit{Mühlenthaler} und
2195 \textit{Wanka} eine Konstruktionsvorschrift für Sortiernetzwerke an, die
2196 gegenüber \bs{n} ${\frac{1}{4}n(\log n - 1)}$ Komparatoren einspart.
2197 Beispielsweise wird ein 16-Sortiernetzwerk angegeben, das nur 68~Komparatoren
2198 benötigt. Dieses Netzwerk ist in Abbildung~\ref{fig:16-muehlenthaler}
2199 dargestellt.
2200
2201 \begin{figure}
2202   \begin{center}
2203     \input{images/16-muehlenthaler.tex}
2204   \end{center}
2205   \caption{Sortiernetzwerk mit 16~Leitungen und 68~Komparatoren in
2206     10~Schichten. Das Netzwerk wurde 2010 von \textit{Mühlenthaler} und
2207     \textit{Wanka} aus optimierten bitonen Mischern konstruiert und
2208     in~\cite{MW2010} veröffentlicht.}
2209   \label{fig:16-muehlenthaler}
2210 \end{figure}
2211
2212 \begin{figure}
2213   \begin{center}
2214     \input{images/16-ec-from-bs32.tex}
2215   \end{center}
2216   \caption{Visualisierung eines 16-Schnittmusters, das von
2217     \textsc{SN-Evolution-Cut} für das \emph{bitone Mergesort}-Netzwerk \bs{32}
2218     berechnet wurde. Das resultierende Sortiernetzwerk besteht aus
2219     68~Komparatoren in 10~Schichten und ist in
2220     Abbildung~\ref{fig:16-ec-from-bs32-normalized} als
2221     Standard-Sortiernetzwerk dargestellt.}
2222   \label{fig:16-ec-from-bs32}
2223 \end{figure}
2224
2225 \begin{figure}
2226   \begin{center}
2227     \input{images/16-ec-from-bs32-normalized.tex}
2228   \end{center}
2229   \caption{Sortiernetzwerk mit 16~Leitungen und 68~Komparatoren in
2230     10~Schichten. Das Netzwerk wurde mit einem 16-Schnittmuster, das von
2231     \textsc{SN-Evolution-Cut} berechnet wurde, aus dem \emph{bitone
2232     Mergesort}-Netzwerk \bs{32} erzeugt. Das Schnittmuster ist in
2233     Abbildung~\ref{fig:16-ec-from-bs32} dargestellt.}
2234   \label{fig:16-ec-from-bs32-normalized}
2235 \end{figure}
2236
2237 Startet man {\sc SN-Evolution-Cut} mit dem \emph{bitonen Mergesort}-Netzwerk
2238 $\operatorname{BS}(32)$ und der Vorgabe 16~Leitungen zu entfernen, liefert der
2239 Algorithmus Sortiernetzwerke, die ebenfalls aus 68~Komparatoren bestehen. Ein
2240 16-Sortiernetzwerk, das auf diese Weise generiert wurde, ist in den
2241 Abbildungen~\ref{fig:16-ec-from-bs32} und~\ref{fig:16-ec-from-bs32-normalized}
2242 zu sehen. Abbildung~\ref{fig:16-ec-from-bs32} zeigt $\operatorname{BS}(32)$
2243 und das Schnittmuster ${\operatorname{MIN}(0, 5, 9, 11, 15, 17, 20, 22, 26,
2244 29, 30)}$, ${\operatorname{MAX}(2, 4, 13, 19, 24)}$, das durch
2245 \textsc{SN-Evolution-Cut} gefunden wurde.
2246 Abbildung~\ref{fig:16-ec-from-bs32-normalized} zeigt das 16-Sortiernetzwerk
2247 nachdem das Schnittmuster angewendet und das Netzwerk normalisiert wurde.
2248 % Eine Ähnlichkeit zu $\operatorname{BS}(32)$ oder $\operatorname{BS}(16)$ ist
2249 % in diesem Netzwerk nicht mehr erkennbar -- insbesondere die ersten Schichten
2250 % des Netzwerks scheinen rein zufällig zu sein.
2251
2252 \begin{figure}
2253   % 0:MAX 1:MAX 4:MIN 6:MAX 9:MAX 11:MAX 14:MIN 15:MAX 18:MAX 19:MAX 21:MAX
2254   % 23:MIN 24:MAX 25:MAX 30:MIN 31:MIN 32:MAX 34:MAX 36:MIN 37:MAX 40:MAX
2255   % 43:MAX 46:MIN 47:MAX 48:MAX 49:MAX 54:MIN 55:MAX 56:MAX 58:MIN 60:MAX
2256   % 63:MAX
2257   \begin{center}
2258     \input{images/32-ec-from-bs64.tex}
2259   \end{center}
2260   \caption{Sortiernetzwerk mit 32~Leitungen und 206~Komparatoren in
2261     15~Schichten. Das Netzwerk wurde von dem Algorithmus
2262     \textsc{SN-Evolution-Cut} aus dem bitonen Mergesort-Netzwerk
2263     $\operatorname{BS}(64)$ durch 32~Schnitte erzeugt. Das zugehörige
2264     Schnittmuster ist
2265     $\operatorname{MIN}(4, 14, 23, 30, 31, 36, 46, 54, 58)$,
2266     $\operatorname{MAX}(0, 1, 6, 9, 11, 15, 18, 19, 21, 24, 25, 32, 34, 37,
2267     40, 43, 47, 48, 49, 55, 56, 60, 63)$.}
2268   \label{fig:32-ec-from-bs64}
2269 \end{figure}
2270
2271 Wenn \textsc{SN-Evolution-Cut} mit dem \bs{64}-Netzwerk und $k = 32$ gestartet
2272 wird, findet der Algorithmus 32-Sortiernetzwerke, die effizienter sind als
2273 32-Sortiernetzwerke, die nach \textit{Mühlenthalers} und \textit{Wankas}
2274 Methode konstruiert werden. Ein von \textsc{SN-Evolution-Cut} aus \bs{64}
2275 generiertes 32-Sortiernetzwerk ist in Abbildung~\ref{fig:32-ec-from-bs64}
2276 dargestellt. Das \emph{bitone Mergesort}-Netzwerk \bs{32} benötigt
2277 240~Komparatoren, ein aus den optimierten Mischern aufgebautes Netzwerk
2278 verbessert die Effizienz auf 208~Komparatoren. Das Ergebnis von
2279 \textsc{SN-Evolution-Cut} kommt mit nur 206~Komparatoren aus. Die
2280 Geschwindigkeit aller genannten Sortiernetzwerke ist mit 15 parallelen
2281 Schritten identisch.
2282
2283 Wenn die Leitungszahl des Eingabenetzwerks keine Zweierpotenz ist, kann
2284 \textsc{SN-Evo\-lu\-tion-Cut} auch 16-Sortiernetzwerke erzeugen, die diese
2285 Effizienz unterbieten. Das geht aus den Daten in
2286 Tabelle~\ref{tbl:ec-bs-efficiency} hervor. Ein 16-Sortiernetzwerk mit
2287 67~Komparatoren, das von \textsc{SN-Evolution-Cut} generiert wurde, ist in
2288 Abbildung~\ref{fig:16-ec-from-bs22} dargestellt.
2289
2290 Leider sind die Schnittmuster, die \textsc{SN-Evolution-Cut} ausgibt, sehr
2291 unregelmäßig. Bisher ist es nicht gelungen eine Konstruktionsanweisung für
2292 gute Schnittmuster anzugeben.
2293
2294 Entscheidend für das Ergebnis eines Schnittmusters scheint beim \emph{bitonen
2295 Mergesort}-Netzwerk die Aufteilung der Minimum- und Maximumschnitte zu sein.
2296 Von Hundert 16-Schnittmustern für $\operatorname{BS}(32)$, die in
2297 Sortiernetzwerken mit 68~Komparatoren in 10~Schichten resultieren, hatten 73
2298 ein Verhältnis von $5/11$, 13 hatten ein Verhältnis von $4/12$ und 14 hatten
2299 ein Verhältnis von $3/13$ Minimum- beziehungsweise Maximumschnitten. Da sich
2300 die Schnittmuster aufgrund der Symmetrie des \emph{bitonen
2301 Mergesort}-Netzwerks leicht invertieren lassen, ist eine Fallunterscheidung --
2302 mehr Minimum- oder mehr Maximumschnitte -- nicht notwendig.
2303
2304 Dass die Sortiernetzwerke, die mit den Schnittmustern von
2305 \textsc{SN-Evolution-Cut} entstehen, keine erkennbare Struktur haben, ist
2306 jedoch kein Eigenschaft des Algorithmus, sondern hängt insbesondere von der
2307 Eingabe ab. Wird \textsc{SN-Evolution-Cut} beispielsweise mit dem
2308 \emph{Odd-Even-Transpositionsort}-Netzwerk $\operatorname{OET}(n)$ und
2309 $k$~Schnitten gestartet, so ist das beste Ergebnis immer das
2310 $\operatorname{OET}(n-k)$-Netzwerk. 
2311
2312 \subsection[Odd-Even-Mergesort-Netzwerk]{Versuche mit dem Odd-Even-Mergesort-Netzwerk}
2313 \label{sect:sn-evolution-cut:oes}
2314
2315 Wird \textsc{SN-Evolution-Cut} mit dem \emph{Odd-Even-Mergesort}-Netzwerk
2316 \oes{n} gestartet, gibt der Algorithmus meist Sortiernetzwerke zurück, die
2317 genauso effizient und schnell wie das entsprechende
2318 \emph{Odd-Even-Mergesort}-Netzwerk \oes{m} sind. Die Effizienz der
2319 Sortiernetzwerke, die mit Schnittmustern von \textsc{SN-Evolution-Cut} aus
2320 \oes{n} entstehen können, zeigt Tabelle~\ref{tbl:ec-oes-efficiency}
2321 tabellarisch.
2322
2323 \begin{table}
2324   \begin{center}
2325     \rowcolors{2}{black!5}{}
2326 \begin{tabular}{|r|rrrrrrrrrrrrrrrr|}
2327 \hline
2328     &   8 &   9 &  10 &  11 &  12 &  13 &  14 &  15 &  16 &  17 &  18 &  19 &  20 &  21 &  22 &  23 \\
2329 \hline
2330   9 &  19 &     &     &     &     &     &     &     &     &     &     &     &     &     &     &     \\
2331  10 &  19 &  26 &     &     &     &     &     &     &     &     &     &     &     &     &     &     \\
2332  11 &  19 &  26 &  31 &     &     &     &     &     &     &     &     &     &     &     &     &     \\
2333  12 &  19 &  26 &  31 &  37 &     &     &     &     &     &     &     &     &     &     &     &     \\
2334  13 &  19 &  26 &  31 &  37 &  41 &     &     &     &     &     &     &     &     &     &     &     \\
2335  14 &  19 &  26 &  31 &  37 &  41 &  48 &     &     &     &     &     &     &     &     &     &     \\
2336  15 &  19 &  26 &  31 &  37 &  41 &  48 &  53 &     &     &     &     &     &     &     &     &     \\
2337  16 &  19 &  26 &  31 &  37 &  41 &  48 &  53 &  59 &     &     &     &     &     &     &     &     \\
2338  17 &  19 &  26 &  31 &  38 &  41 &  48 &  53 &  59 &  63 &     &     &     &     &     &     &     \\
2339  18 &  19 &  26 &  31 &  38 &  43 &  48 &  53 &  59 &  63 &  74 &     &     &     &     &     &     \\
2340  19 &  19 &  26 &  31 &  38 &  43 &  48 &  53 &  59 &  63 &  74 &  82 &     &     &     &     &     \\
2341  20 &  19 &  26 &  31 &  38 &  43 &  48 &  53 &  59 &  63 &  74 &  82 &  91 &     &     &     &     \\
2342  21 &  19 &  26 &  31 &  38 &  43 &  48 &  53 &  59 &  63 &  74 &  82 &  91 &  97 &     &     &     \\
2343  22 &  19 &  26 &  31 &  38 &  43 &  48 &  53 &  59 &  63 &  74 &  82 &  91 &  97 & 107 &     &     \\
2344  23 &  19 &  26 &  31 &  38 &  43 &  48 &  53 &  59 &  63 &  74 &  82 &  91 &  97 & 107 & 114 &     \\
2345  24 &  19 &  26 &  31 &  38 &  43 &  48 &  53 &  59 &  63 &  74 &  82 &  91 &  97 & 107 & 114 & 122 \\
2346 \hline
2347 \oes{m}&19&  26 &  31 &  37 &  41 &  48 &  53 &  59 &  63 &  74 &  82 &  91 &  97 & 107 & 114 & 122 \\
2348 \hline
2349 \end{tabular}
2350   \end{center}
2351   \caption{Anzahl der Komparatoren der Ergebnisse von
2352     \textsc{SN-Evolution-Cut} mit verschiedenen Größen des
2353     \emph{Odd-Even-Mergesort}-Netzwerks und unterschiedlichen Werten für~$k$.
2354     Jede Zeile gibt die Ergebnisse für ein Eingabenetzwerk \oes{n} an, jede
2355     Spalte enthält die Ergebnisse für $m=n-k$, die Anzahl der Leitungen des
2356     Ausgabenetzwerks.}
2357   \label{tbl:ec-oes-efficiency}
2358 \end{table}
2359
2360 \begin{figure}
2361   \centering
2362   \subfigure[11-Sortiernetzwerk aus 38~Komparatoren in 9~Schichten. Das
2363   Netzwerk wurde von \textsc{SN-Evolution-Cut} aus \oes{17} erzeugt.]{\input{images/11-ec-from-oes17-fast.tex}\label{fig:11-ec-from-oes17-fast}}
2364   \subfigure[12-Sortiernetzwerk aus 43~Komparatoren in 9~Schichten. Das
2365   Netzwerk wurde von \textsc{SN-Evolution-Cut} aus \oes{18} erzeugt.]{\input{images/12-ec-from-oes18-fast.tex}\label{fig:12-ec-from-oes18-fast}}
2366   \caption{Für einige Ziel-Leitungszahlen, unter anderem $m = 10$ und $m =
2367   11$, kann der \textsc{SN-Evolution-Cut}-Algorithmus Sortiernetzwerke
2368   erzeugen, die \emph{schneller} aber weniger \emph{effizient} als \oes{m}
2369   sind.}
2370   \label{fig:ec-oes-fast_networks}
2371 \end{figure}
2372
2373 Die Bewertungsfunktion, die \textsc{SN-Evolution-Cut} verwendet, bevorzugt
2374 schnelle Sortiernetzwerke. Dadurch kann es vorkommen, dass ein
2375 $m$-Sortiernetzwerk, das durch ein von \textsc{SN-Evolution-Cut} ausgegebenes
2376 Schnittmuster entsteht, schneller als \oes{m} ist. Diese Geschwindigkeit
2377 war allerdings in allen beobachteten Fällen nur dann möglich, wenn
2378 zusätzliche Komparatoren in Kauf genommen wurden. In den
2379 Tabellen~\ref{tbl:ec-oes-efficiency} und~\ref{tbl:ec-oes-speed} ist dieser
2380 Fall für $m = 11$ und $k \geqq 6$, beziehungsweise $m = 12$ und $k \geqq 6$ zu
2381 beobachten. Die entsprechenden schnellen Sortiernetzwerke sind in
2382 Abbildung~\ref{fig:ec-oes-fast_networks} dargestellt.
2383
2384 Wie beim \emph{bitonen Mergesort}-Netzwerk reicht auch beim
2385 \emph{Odd-Even-Mergesort}-Netzwerk ein einziger Schnitt nicht aus, um die
2386 Geschwindigkeit gegenüber \oes{m} zu verbessern. Bei $m = 11$ und $m = 12$ war
2387 jeweils mindestens ein 6-Schnittmuster notwendig, um eine höhere
2388 Geschwindigkeit zu erreichen.
2389
2390 In Tabelle~\ref{tbl:ec-oes-19} sind die Ergebnisse von
2391 \textsc{SN-Evolution-Cut} für \oes{n}, $n = 20$ und $m = 19$ ($k = 1 \dots
2392 19$) aufgelistet. Mit $k = 10$ wird das erste mal ein schnelles
2393 19-Sortiernetzwerk mit 13~Schichten ausgegeben. Mit $k \geqq 11$ sind die
2394 resultierenden Netzwerke mit 93~Komparatoren effizienter als das Ergebnis mit
2395 $k = 10$, das 95~Komparatoren benötigt. Das Ergebnis, das auf Basis des
2396 \emph{bitonen Mergesort}-Netzwerks erreicht wurde (92~Komparatoren in
2397 13~Schichten, siehe Tabelle~\ref{tbl:ec-bs-19}), wird nicht erreicht.
2398
2399 \begin{table}
2400   \begin{center}
2401     \rowcolors{2}{black!5}{}
2402 \begin{tabular}{|r|rrrrrrrrrrrrrrrr|}
2403 \hline
2404     &   8 &   9 &  10 &  11 &  12 &  13 &  14 &  15 &  16 &  17 &  18 &  19 &  20 &  21 &  22 &  23 \\
2405 \hline
2406   9 &   6 &     &     &     &     &     &     &     &     &     &     &     &     &     &     &     \\
2407  10 &   6 &   8 &     &     &     &     &     &     &     &     &     &     &     &     &     &     \\
2408  11 &   6 &   8 &   9 &     &     &     &     &     &     &     &     &     &     &     &     &     \\
2409  12 &   6 &   8 &   9 &  10 &     &     &     &     &     &     &     &     &     &     &     &     \\
2410  13 &   6 &   8 &   9 &  10 &  10 &     &     &     &     &     &     &     &     &     &     &     \\
2411  14 &   6 &   8 &   9 &  10 &  10 &  10 &     &     &     &     &     &     &     &     &     &     \\
2412  15 &   6 &   8 &   9 &  10 &  10 &  10 &  10 &     &     &     &     &     &     &     &     &     \\
2413  16 &   6 &   8 &   9 &  10 &  10 &  10 &  10 &  10 &     &     &     &     &     &     &     &     \\
2414  17 &   6 &   8 &   9 &   9 &  10 &  10 &  10 &  10 &  10 &     &     &     &     &     &     &     \\
2415  18 &   6 &   8 &   9 &   9 &   9 &  10 &  10 &  10 &  10 &  12 &     &     &     &     &     &     \\
2416  19 &   6 &   8 &   9 &   9 &   9 &  10 &  10 &  10 &  10 &  12 &  13 &     &     &     &     &     \\
2417  20 &   6 &   8 &   9 &   9 &   9 &  10 &  10 &  10 &  10 &  12 &  13 &  14 &     &     &     &     \\
2418  21 &   6 &   8 &   9 &   9 &   9 &  10 &  10 &  10 &  10 &  12 &  13 &  14 &  14 &     &     &     \\
2419  22 &   6 &   8 &   9 &   9 &   9 &  10 &  10 &  10 &  10 &  12 &  13 &  14 &  14 &  15 &     &     \\
2420  23 &   6 &   8 &   9 &   9 &   9 &  10 &  10 &  10 &  10 &  12 &  13 &  14 &  14 &  15 &  15 &     \\
2421  24 &   6 &   8 &   9 &   9 &   9 &  10 &  10 &  10 &  10 &  12 &  13 &  14 &  14 &  15 &  15 &  15 \\
2422 \hline
2423 \oes{m}& 6 &  8 &   9 &  10 &  10 &  10 &  10 &  10 &  10 &  12 &  13 &  14 &  14 &  15 &  15 &  15 \\
2424 \hline
2425 \end{tabular}
2426   \end{center}
2427   \caption{Anzahl der Schichten der Ergebnisse von
2428     \textsc{SN-Evolution-Cut} mit verschiedenen Größen des
2429     \emph{Odd-Even-Mergesort}-Netzwerks und unterschiedlichen Werten für~$k$.
2430     Jede Zeile gibt die Ergebnisse für ein Eingabenetzwerk \oes{n} an, jede
2431     Spalte enthält die Ergebnisse für $m=n-k$, die Anzahl der Leitungen des
2432     Ausgabenetzwerks.}
2433   \label{tbl:ec-oes-speed}
2434 \end{table}
2435
2436 \begin{table}
2437   \begin{center}
2438     \rowcolors{2}{black!5}{}
2439     \begin{tabular}{|r|r|r|}
2440       \hline
2441       $n$ & Komp. & Schichten \\
2442       \hline
2443       20  &  91 & 14 \\
2444       21  &  91 & 14 \\
2445       22  &  91 & 14 \\
2446       23  &  91 & 14 \\
2447       24  &  91 & 14 \\
2448       25  &  91 & 14 \\
2449       26  &  91 & 14 \\
2450       27  &  91 & 14 \\
2451       28  &  91 & 14 \\
2452       29  &  95 & 13 \\
2453       \rowcolor{green!10}
2454       30  &  93 & 13 \\
2455       \rowcolor{green!10!white!95!black}
2456       31  &  93 & 13 \\
2457       \rowcolor{green!10}
2458       32  &  93 & 13 \\
2459       \rowcolor{green!10!white!95!black}
2460       33  &  93 & 13 \\
2461       \rowcolor{green!10}
2462       34  &  93 & 13 \\
2463       \rowcolor{green!10!white!95!black}
2464       35  &  93 & 13 \\
2465       \rowcolor{green!10}
2466       36  &  93 & 13 \\
2467       \rowcolor{green!10!white!95!black}
2468       37  &  93 & 13 \\
2469       \rowcolor{green!10}
2470       38  &  93 & 13 \\
2471       \hline
2472  \bs{19}  &  98 & 14 \\
2473  \oes{19} &  91 & 14 \\
2474       \hline
2475     \end{tabular}
2476   \end{center}
2477   \caption{Komparatoren und Schichten von Sortiernetzwerken, die von
2478     \textsc{SN-Evolution-Cut} mit \oes{n} und $k = n - 19$ ermittelt wurden. Erst mit $k = 10$
2479     ist es möglich gegenüber \oes{19} eine Schicht einzusparen. Dafür ist die
2480     Effizienz von 91~Komparatoren nicht mehr erreichbar.}
2481   \label{tbl:ec-oes-19}
2482 \end{table}
2483
2484 % 2-er Potenzen
2485
2486 In Abschnitt~\ref{sect:anzahl_schnittmuster} wurde bereits untersucht, wie
2487 viele \emph{unterschiedliche} 16-Schnittmuster die konstruierten
2488 Sortiernetzwerke $\operatorname{OES}(32)$, $\operatorname{BS}(32)$ und
2489 $\operatorname{PS}(32)$ besitzen. Eines der Ergebnisse war, dass von diesen
2490 Sortiernetzwerken das \emph{Odd-Even-Mergesort}-Netzwerk die wenigsten
2491 unterschiedlichen 16-Schnittmuster besitzt -- nur etwa $5,2$~Millionen.
2492 Entsprechend ist es wenig verwunderlich, dass \textsc{SN-Evolution-Cut}
2493 gestartet mit $\operatorname{OES}(32)$ sehr schnell\footnote{Ein
2494 entsprechendes Ergebnis wird meist nach 20.000 bis 100.000 Iterationen
2495 geliefert. Bei dieser Problemgröße erreicht die Implementierung (siehe
2496 Abschnitt~\ref{sect:implementierung}) etwa 20.000 Iterationen pro Sekunde auf
2497 derzeitigen Computern.} ein gutes 16-Schnittmuster findet.
2498
2499 Eines der 16-Schnittmuster für \oes{32}, die ein Sortiernetzwerk erzeugen, das
2500 bezüglich Effizienz und Geschwindigkeit identisch zu \oes{16} ist, ist
2501 $\operatorname{MIN}(1, 6, 11, 14, 17, 23, 26, 29)$, $\operatorname{MAX}(2, 7,
2502 8,$ $13, 18, 21, 27, 31)$. Das Schnittmuster ist in
2503 Abbildung~\ref{fig:16-ec-from-oes32-cut} veranschaulicht, das resultierende
2504 Netzwerk ist in Abbildung~\ref{fig:16-ec-from-oes32} zu sehen.
2505
2506 \begin{figure}
2507   \begin{center}
2508     \input{images/16-ec-from-oes32-cut.tex}
2509   \end{center}
2510   \caption{Visualisierung eines 16-Schnittmusters, das auf
2511   $\operatorname{OES}(32)$ angewendet ein Sortiernetzwerk ergibt, das
2512   bezüglich Geschwindigkeit und Effizienz identisch zu \oes{16} ist. Das
2513   resultierende Sortiernetzwerk ist in Abbildung~\ref{fig:16-ec-from-oes32}
2514   dargestellt.}
2515   \label{fig:16-ec-from-oes32-cut}
2516 \end{figure}
2517
2518 \begin{figure}
2519   \begin{center}
2520     \input{images/16-ec-from-oes32.tex}
2521   \end{center}
2522   \caption{16-Sortiernetzwerk mit 63~Komparatoren in 10~Schichten. 
2523     Das Netzwerk wurde aus dem \emph{Odd-Even-Mergesort}-Netzwerk \oes{32} mit
2524     einem 16-Schnittmuster erzeugt, das von \textsc{SN-Evolution-Cut}
2525     berechnet wurde. Das Schnittmuster ist in
2526     Abbildung~\ref{fig:16-ec-from-oes32-cut} dargestellt.}
2527   \label{fig:16-ec-from-oes32}
2528 \end{figure}
2529
2530 % Regelmaessiges Schnittmuster fuer n = 2^d
2531
2532 Bei diesem Schnittmuster fällt auf, dass es für jeweils vier Eingänge (0--3,
2533 4--7, \dots, 28--31) einen Minimum- und einen Maximumschnitt gibt. Aus dieser
2534 Beobachtung kann das regelmäßige Schnittmuster
2535 \begin{displaymath}
2536 \textit{Eingang}_i = \left\{ \begin{array}{rl}
2537    \infty & \quad \textrm{falls } i \bmod 4 = 0 \\
2538   -\infty & \quad \textrm{falls } i \bmod 4 = 3 \\
2539         ? & \quad \mathrm{sonst}
2540   \end{array} \right.
2541 \end{displaymath}
2542 abgeleitet werden. Es entfernt die Hälfte der Leitungen, vorausgesetzt die
2543 Anzahl der Leitungen ist durch Vier teilbar. Das Schnittmuster erzeugt
2544 effiziente Netzwerke, wenn die Anzahl der Leitungen $n = 2^d$ eine
2545 Zweierpotenz ist. Ein 32-Sortiernetzwerk, das mit diesem Schnittmuster aus
2546 \oes{64} erzeugt wurde, ist in Abbildung~\ref{fig:32-ec-from-oes64} zu sehen.
2547
2548 \begin{figure}
2549   \begin{center}
2550     \input{images/32-ec-from-oes64.tex}
2551   \end{center}
2552   \caption{32-Sortiernetzwerk mit 191~Komparatoren in 15~Schichten. 
2553     Das Netzwerk wurde mit einem regelmäßigen Schnittmuster aus dem
2554     \emph{Odd-Even-Mergesort}-Netzwerk \oes{64} erzeugt.}
2555   \label{fig:32-ec-from-oes64}
2556 \end{figure}
2557
2558 Wenn die Anzahl der Leitungen keine Zweierpotenz ist, erreichen die durch
2559 dieses regelmäßige Schnittmuster erzeugten Sortiernetzwerke die Effizienz des
2560 \emph{Odd-Even-Mergesort}-Netzwerks nicht. Wendet man das Schnittmuster
2561 beispielsweise auf \oes{24} an, so erhält man ein Sortiernetzwerk mit
2562 43~Komparatoren -- \oes{12} kommt mit 41~Komparatoren aus. Die Geschwindigkeit
2563 beider Sortiernetzwerke ist mit 10~Schichten identisch.
2564
2565 % SN-Evolution-Cut vs. regelmaessiges Schnittmuster
2566
2567 Wird der \textsc{SN-Evolution-Cut}-Algorithmus mit \oes{n}, $n = 22, 23, 24$
2568 und $k = 10, 11, 12$ gestartet, hängt die Ausgabe von der verwendeten
2569 Gütefunktion ab. Werden effiziente Netzwerke bevorzugt, findet der Algorithmus
2570 Schnittmuster wie das 12-Schnittmus\-ter $\operatorname{MIN}(6, 7, 8, 9, 16, 17,
2571 20, 22)$, $\operatorname{MAX}(2, 4, 12, 14)$ für \oes{24}, dessen Ergebnis in
2572 Abbildung~\ref{fig:12-ec-from-oes24-efficient} zu sehen ist. Das resultierende
2573 Sortiernetzwerk besteht aus 41~Komparatoren, die in 10~Schichten angeordnet
2574 werden können. Damit ist das Netzwerk bezüglich Effizienz und Geschwindigkeit
2575 gleichauf mit \oes{12}. Werden hingegen schnelle Sortiernetzwerke bevorzugt,
2576 werden stattdessen Schnittmuster wie $\operatorname{MIN}(6, 7, 11, 12, 15,
2577 16)$, $\operatorname{MAX}(1, 3, 10, 17, 20, 23)$ ausgegeben. Das Ergebnis
2578 dieses Schnittmusters ist in Abbildung~\ref{fig:12-ec-from-oes24-fast} zu
2579 sehen, weitere Ergebnisse für diese Gütefunktion sind in den
2580 Tabellen~\ref{tbl:ec-oes-efficiency} und~\ref{tbl:ec-oes-speed} zusammengefasst.
2581 Das resultierende 12-Sortiernetzwerk besteht aus 43~Komparatoren, die
2582 in 9~Schichten angeordnet sind. Das resultierende Netzwerk zwar nicht so
2583 effizient wie \oes{12}, dafür aber schneller als \oes{12} und \bs{12}.
2584
2585 \begin{figure}
2586   \centering
2587   \subfigure[Effizientes 12-Sortiernetzwerk aus 41~Komparatoren in
2588   10~Schichten, das von \textsc{SN-Evolution-Cut} aus dem
2589   \emph{Odd-Even-Mergesort}-Netzwerk generiert
2590   wurde.]{\input{images/12-ec-from-oes24-efficient.tex}\label{fig:12-ec-from-oes24-efficient}}
2591   \subfigure[Schnelles 12-Sortiernetzwerk aus 43~Komparatoren in 9~Schichten,
2592   das von \textsc{SN-Evolution-Cut} aus dem \emph{Odd-Even-Mergesort}-Netzwerk
2593   generiert
2594   wurde.]{\input{images/12-ec-from-oes24-fast.tex}\label{fig:12-ec-from-oes24-fast}}
2595   \caption{Startet man \textsc{SN-Evolution-Cut} mit \oes{24}, hängt das
2596   Ergebnis von der Bewertungsfunktion ab.}
2597   \label{fig:12-ec-from-oes24}
2598 \end{figure}
2599
2600 \subsection[Pairwise-Sorting-Netzwerk]{Versuche mit dem Pairwise-Sorting-Netzwerk}
2601
2602 Eine weitere interessante Eingabe für \textsc{SN-Evolution-Cut} ist das
2603 \emph{Pairwise-Sorting}-Netzwerk \ps{n}, das \textit{Ian
2604 Parberry} in seiner Arbeit „The Pairwise Sorting Network“ \cite{P1992}
2605 definiert. Einerseits wurde in Abschnitt~\ref{sect:anzahl_schnittmuster}
2606 gezeigt, dass es für \ps{n} sehr viele \emph{unterschiedliche} Schnittmuster
2607 gibt. Andererseits sind die Sortiernetzwerke, die nach \textit{Parberrys}
2608 Methode erzeugt werden, genauso schnell und effizient wie das
2609 \emph{Odd-Even-Mergesort}-Netzwerk, wenn die Leitungszahl $n = 2^d$ eine
2610 Zweierpotenz ist.
2611
2612 % Effizienz
2613
2614 Für viele Kombinationen von \ps{n} und $k$ sind die Sortiernetzwerke, die
2615 \textsc{SN-Evolution-Cut} ausgibt, weniger effizient als das entsprechende
2616 \ps{m}-Netzwerk. Für einige Kombinationen werden jedoch auch effizientere
2617 Netzwerke generiert, beispielsweise für $m = 11$ und $m = 20$. Die Effizienz
2618 der Sortiernetzwerke, die \textsc{SN-Evolution-Cut} auf Basis des
2619 \emph{Pairwise-Sorting}-Netzwerks berechnet, ist tabellarisch in
2620 Tabelle~\ref{tbl:ec-ps-efficiency} dargestellt. 
2621
2622 \begin{table}
2623   \begin{center}
2624     \rowcolors{2}{black!5}{}
2625 \begin{tabular}{|r|rrrrrrrrrrrrrrrr|}
2626 \hline
2627     &   8 &   9 &  10 &  11 &  12 &  13 &  14 &  15 &  16 &  17 &  18 &  19 &  20 &  21 &  22 &  23 \\
2628 \hline
2629   9 &  20 &     &     &     &     &     &     &     &     &     &     &     &     &     &     &     \\
2630  10 &  20 &  27 &     &     &     &     &     &     &     &     &     &     &     &     &     &     \\
2631  11 &  20 &  28 &  32 &     &     &     &     &     &     &     &     &     &     &     &     &     \\
2632  12 &  20 &  28 &  32 &  38 &     &     &     &     &     &     &     &     &     &     &     &     \\
2633  13 &  19 &  27 &  31 &  37 &  41 &     &     &     &     &     &     &     &     &     &     &     \\
2634  14 &  19 &  27 &  31 &  37 &  41 &  48 &     &     &     &     &     &     &     &     &     &     \\
2635  15 &  19 &  27 &  31 &  37 &  41 &  48 &  53 &     &     &     &     &     &     &     &     &     \\
2636  16 &  19 &  27 &  31 &  37 &  41 &  48 &  53 &  59 &     &     &     &     &     &     &     &     \\
2637  17 &  21 &  29 &  32 &  39 &  43 &  51 &  57 &  64 &  68 &     &     &     &     &     &     &     \\
2638  18 &  22 &  29 &  32 &  39 &  43 &  52 &  58 &  65 &  69 &  80 &     &     &     &     &     &     \\
2639  19 &  23 &  29 &  32 &  39 &  43 &  52 &  58 &  65 &  69 &  80 &  88 &     &     &     &     &     \\
2640  20 &  23 &  29 &  32 &  39 &  43 &  52 &  58 &  65 &  69 &  80 &  88 &  97 &     &     &     &     \\
2641  21 &  20 &  30 &  34 &  38 &  44 &  51 &  57 &  64 &  74 &  82 &  87 &  96 & 102 &     &     &     \\
2642  22 &  20 &  30 &  34 &  38 &  46 &  51 &  57 &  64 &  72 &  82 &  89 &  96 & 102 & 112 &     &     \\
2643  23 &  20 &  27 &  34 &  38 &  42 &  51 &  57 &  66 &  72 &  83 &  89 &  97 & 102 & 112 & 119 &     \\
2644  24 &  20 &  27 &  34 &  38 &  42 &  51 &  57 &  64 &  72 &  82 &  89 &  96 & 102 & 112 & 119 & 127 \\
2645 \hline
2646 \ps{m}&19 &  27 &  32 &  38 &  42 &  48 &  53 &  59 &  63 &  79 &  88 &  97 & 103 & 112 & 119 & 127 \\
2647 \hline
2648 \end{tabular}
2649   \end{center}
2650   \caption{Anzahl der Komparatoren der Ergebnisse von
2651     \textsc{SN-Evolution-Cut} mit verschiedenen Größen des
2652     \emph{Pairwise-Sorting}-Netzwerks und unterschiedlichen Werten für~$k$.
2653     Jede Zeile gibt die Ergebnisse für ein Eingabenetzwerk \ps{n} an, jede
2654     Spalte enthält die Ergebnisse für $m=n-k$, die Anzahl der Leitungen des
2655     Ausgabenetzwerks.}
2656   \label{tbl:ec-ps-efficiency}
2657 \end{table}
2658
2659 % Beispiel Effizienz
2660
2661 Zwei Ergebnisse, die effizienter als die entsprechenden
2662 \emph{Pairwise-Sorting}-Netzwerke sind, zeigt
2663 Abbildung~\ref{fig:ec-ps-efficient_networks}. Sie erreichen die
2664 Geschwindigkeit und Effizienz des \emph{Odd-Even-Mergesort}-Netzwerks mit der
2665 entsprechenden Leitungszahl. Bei größeren Netzwerken, beispielsweise $m = 19$,
2666 ist dies mit der in Tabelle~\ref{tbl:ec-ps-efficiency} dargestellten Größe der
2667 Schnittmuster noch nicht zu beobachten.
2668
2669 \begin{figure}
2670   \centering
2671   \subfigure[11-Sortiernetzwerk aus 37~Komparatoren in 10~Schichten. Das
2672     Netzwerk wurde von \textsc{SN-Evolution-Cut} aus dem
2673     \emph{Pairwise-Sorting}-Netzwerk \ps{13}
2674     erzeugt.]{\input{images/11-ec-from-ps13.tex}}
2675   \subfigure[12-Sortiernetzwerk aus 41~Komparatoren in 10~Schichten. Das
2676     Netzwerk wurde von \textsc{SN-Evolution-Cut} aus dem
2677     \emph{Pairwise-Sorting}-Netzwerk \ps{16}
2678     erzeugt.]{\input{images/12-ec-from-ps16.tex}}
2679   \caption{Zwei effiziente Sortiernetzwerke, die durch Schnittmuster, die von
2680     \emph{SN-Evolution-Cut} berechnet wurden, aus dem
2681     \emph{Pairwise-Sorting}-Netzwerk \ps{n} erzeugt wurden.}
2682   \label{fig:ec-ps-efficient_networks}
2683 \end{figure}
2684
2685 % Wie viele Schnitte man braucht.
2686
2687 Bei welchen Parametern \textsc{SN-Evolution-Cut} effiziente
2688 19-Sortiernetzwerke findet, ist Tabelle~\ref{tbl:ec-ps-19} zu entnehmen. Für
2689 $n = 31$, $k = 12$ und $n = 32$, $k = 13$ werden 19-Sortiernetzwerke mit der
2690 selben Effizienz und Geschwindigkeit wie \oes{19} erzeugt. Das
2691 19-Sortiernetzwerk, das auf diese Art und Weise aus \ps{32} erzeugt wurde, ist
2692 in Abbildung~\ref{fig:19-ec-from-ps32} dargestellt.
2693
2694 \begin{table}
2695   \begin{center}
2696     \rowcolors{2}{black!5}{}
2697     \begin{tabular}{|r|r|r|}
2698     \hline
2699     $n$ & Komp. & Schichten \\
2700     \hline
2701           20 &  97 & 15 \\
2702           21 &  96 & 15 \\
2703           22 &  96 & 15 \\
2704           23 &  97 & 14 \\
2705           24 &  96 & 14 \\
2706           25 &  93 & 14 \\
2707           26 &  92 & 14 \\
2708           27 &  94 & 14 \\
2709           28 &  94 & 14 \\
2710           29 &  92 & 14 \\
2711           30 &  92 & 14 \\
2712           \rowcolor{green!10!white!95!black}
2713           31 &  91 & 14 \\
2714           \rowcolor{green!10}
2715           32 &  91 & 14 \\
2716           33 & 101 & 15 \\
2717           34 & 104 & 15 \\
2718           35 & 106 & 15 \\
2719           36 & 107 & 15 \\
2720           37 & 106 & 15 \\
2721           38 & 102 & 15 \\
2722     \hline
2723      \ps{19} &  97 & 15 \\
2724     \oes{19} &  91 & 14 \\
2725     \hline
2726     \end{tabular}
2727   \end{center}
2728   \caption{Anzahl der Komparatoren und Schichten von 19-Sortiernetzwerken, die
2729     von \textsc{SN-Evolution-Cut} aus \ps{n}, $n = 20, \dots, 38$ erzeugt
2730     wurden.}
2731   \label{tbl:ec-ps-19}
2732 \end{table}
2733
2734 \begin{figure}
2735   \begin{center}
2736     \input{images/19-ec-from-ps32.tex}
2737   \end{center}
2738   \caption{Sortiernetzwerk mit 19~Leitungen und 91~Komparatoren in
2739     14~Schichten. Das Netzwerk wurde von dem Algorithmus
2740     \textsc{SN-Evolution-Cut} aus dem \emph{Pairwise-Sorting}-Netzwerk
2741     $\operatorname{PS}(32)$ erzeugt.}
2742   \label{fig:19-ec-from-ps32}
2743 \end{figure}
2744
2745 An den Daten in Tabelle~\ref{tbl:ec-ps-19} fällt auf, dass die Effizienz und
2746 Geschwindigkeit der Ergebnisse für $n > 32$ schlechter werden. Das
2747 \emph{Pairwise-Sorting}-Netzwerk ist Leitungszahlen, die Zweierpotenzen sind,
2748 besonders effizient und schnell. Um der Vermutung nachzugehen, dass der
2749 \textsc{SN-Evolution-Cut}-Algorithmus für $\ps{n = 2^d}$ besonders effiziente
2750 Schnittmuster findet, wurde \textsc{SN-Evolution-Cut} mit \ps{32} und $k = 1
2751 \dots 16$ gestartet. Die Ergebnisse sind in Tabelle~\ref{tbl:ec-ps-32}
2752 zusammengefasst.
2753
2754 \begin{table}
2755   \begin{center}
2756     \rowcolors{2}{black!5}{}
2757     \begin{tabular}{|r|r|r|}
2758     \hline
2759     $m$ & Komp. & Schichten \\
2760     \hline
2761      16 &    69 &        11 \\
2762      17 &    77 &        13 \\
2763      18 &    89 &        13 \\
2764      19 &    91 &        14 \\
2765      20 &    97 &        14 \\
2766      21 &   107 &        15 \\
2767      22 &   114 &        15 \\
2768      23 &   122 &        15 \\
2769      24 &   127 &        15 \\
2770      25 &   138 &        15 \\
2771      26 &   146 &        15 \\
2772      27 &   155 &        15 \\
2773      28 &   161 &        15 \\
2774      29 &   171 &        15 \\
2775      30 &   178 &        15 \\
2776      31 &   186 &        15 \\
2777     \hline
2778     \end{tabular}
2779   \end{center}
2780   \caption{Anzahl der Komparatoren und Schichten von $m$-Sortiernetzwerken,
2781     die von \textsc{SN-Evolution-Cut} nach 5.000.000 Iterationen aus \ps{32}
2782     erzeugt wurden.}
2783   \label{tbl:ec-ps-32}
2784 \end{table}
2785
2786 % Geschwindigkeit
2787
2788 Die Schnittmuster, die \textsc{SN-Evolution-Cut} für das
2789 \emph{Pairwise-Sorting}-Netzwerk berechnet, können zu schnelleren
2790 Sortiernetzwerken als \ps{m} führen. Beispielsweise konnte aus \ps{24} ein
2791 18-Sortiernetzwerk erzeugt werden, das mit 13~Schichten zwei parallele
2792 Schritte im Vergleich zu \ps{18} einspart. Eine Darstellung dieses
2793 Sortiernetzwerks befindet sich in Abbildung~\ref{fig:18-ec-from-ps24}. Für
2794 andere $m$ wurde die Geschwindigkeit des \emph{Pairwise-Sorting}-Netzwerks
2795 nicht übertroffen und im Fall von $m = 16$ wurde die Geschwindigkeit nicht
2796 einmal erreicht.
2797
2798 \begin{table}
2799   \begin{center}
2800     \rowcolors{2}{black!5}{}
2801     \begin{tabular}{|r|rrrrrrrrrrrrrrrr|}
2802     \hline
2803         &   8 &   9 &  10 &  11 &  12 &  13 &  14 &  15 &  16 &  17 &  18 &  19 &  20 &  21 &  22 &  23 \\
2804     \hline
2805       9 &   6 &     &     &     &     &     &     &     &     &     &     &     &     &     &     &     \\
2806      10 &   6 &  10 &     &     &     &     &     &     &     &     &     &     &     &     &     &     \\
2807      11 &   6 &   9 &  10 &     &     &     &     &     &     &     &     &     &     &     &     &     \\
2808      12 &   6 &   8 &   9 &  10 &     &     &     &     &     &     &     &     &     &     &     &     \\
2809      13 &   6 &   8 &   9 &  10 &  10 &     &     &     &     &     &     &     &     &     &     &     \\
2810      14 &   6 &   8 &   9 &  10 &  10 &  10 &     &     &     &     &     &     &     &     &     &     \\
2811      15 &   6 &   8 &   9 &  10 &  10 &  10 &  10 &     &     &     &     &     &     &     &     &     \\
2812      16 &   6 &   8 &   9 &  10 &  10 &  10 &  10 &  10 &     &     &     &     &     &     &     &     \\
2813      17 &   7 &   8 &   9 &  10 &  10 &  11 &  11 &  11 &  11 &     &     &     &     &     &     &     \\
2814      18 &   7 &   8 &   9 &  10 &  10 &  11 &  11 &  11 &  11 &  15 &     &     &     &     &     &     \\
2815      19 &   7 &   8 &   9 &  10 &  10 &  11 &  11 &  11 &  11 &  15 &  15 &     &     &     &     &     \\
2816      20 &   7 &   8 &   9 &  10 &  10 &  11 &  11 &  11 &  11 &  15 &  15 &  15 &     &     &     &     \\
2817      21 &   6 &   8 &   9 &  10 &  10 &  11 &  11 &  11 &  12 &  14 &  15 &  15 &  15 &     &     &     \\
2818      22 &   6 &   8 &   9 &  10 &  10 &  11 &  11 &  11 &  12 &  14 &  14 &  15 &  15 &  15 &     &     \\
2819      23 &   6 &   9 &   9 &  10 &  10 &  10 &  11 &  11 &  11 &  13 &  14 &  14 &  15 &  15 &  15 &     \\
2820      24 &   6 &   9 &   9 &  10 &  10 &  10 &  11 &  11 &  11 &  13 &  13 &  14 &  14 &  15 &  15 &  15 \\
2821      \hline
2822  \ps{m} &   6 &  10 &  10 &  10 &  10 &  10 &  10 &  10 &  10 &  15 &  15 &  15 &  15 &  15 &  15 &  15 \\
2823     \hline
2824     \end{tabular}
2825   \end{center}
2826   \caption{Anzahl der Schichten der Ergebnisse von
2827     \textsc{SN-Evolution-Cut} mit verschiedenen Größen des
2828     \emph{Pairwise-Sorting}-Netzwerks und unterschiedlichen Werten für~$k$.
2829     Jede Zeile gibt die Ergebnisse für ein Eingabenetzwerk \ps{n} an, jede
2830     Spalte enthält die Ergebnisse für $m=n-k$, die Anzahl der Leitungen des
2831     Ausgabenetzwerks.}
2832   \label{tbl:ec-ps-speed}
2833 \end{table}
2834
2835 % Beispiel Geschwindigkeit
2836
2837 \begin{figure}
2838   \begin{center}
2839     \input{images/18-ec-from-ps24.tex}
2840   \end{center}
2841   \caption{Sortiernetzwerk mit 18~Leitungen und 89~Komparatoren in
2842     13~Schichten. Das Netzwerk wurde von dem Algorithmus
2843     \textsc{SN-Evolution-Cut} aus dem \emph{Pairwise-Sorting}-Netzwerk
2844     $\operatorname{PS}(24)$ erzeugt.}
2845   \label{fig:18-ec-from-ps24}
2846 \end{figure}
2847
2848 % 2-er Potenz
2849
2850 Die Ergebnisse, die \textsc{SN-Evolution-Cut} erzielte, wenn das gegebene
2851 Sortiernetzwerk das \emph{bitone Mergesort}-Netzwerk war
2852 (Abschnitt~\ref{sect:sn-evolution-cut:bs}), waren sehr wirr. Beispielsweise
2853 ist bei dem Netzwerk in Abbildung~\ref{fig:32-ec-from-bs64} nicht ersichtlich,
2854 wie und warum es jede beliebige Eingabe sortiert.
2855
2856 Bei dem \emph{Pairwise-Sorting}-Netzwerk $\ps{n=2^d}$ ist das anders. Startet
2857 man \textsc{SN-Evolution-Cut} mit $\operatorname{PS}(32)$ und der Vorgabe,
2858 16~Leitungen zu entfernen, kann der Algorithmus ein Sortiernetzwerk
2859 zurückgeben, das die gleiche Anzahl Komparatoren und Schichten wie
2860 $\operatorname{PS}(16)$ und $\operatorname{OES}(16)$ hat. Eines dieser
2861 Sortiernetzwerke ist in Abbildung~\ref{fig:16-ec-from-ps32} dargestellt.
2862 Dieses Ergebnis demonstriert, dass sich die Ergebnisse in den gezeigten
2863 Tabellen oft durch zusätzliche Iterationen verbessern lassen.
2864
2865 \begin{figure}
2866   \begin{center}
2867     \input{images/16-ec-from-ps32.tex}
2868   \end{center}
2869   \caption{Sortiernetzwerk mit 16~Leitungen und 63~Komparatoren in
2870     10~Schichten. Das Netzwerk wurde von dem Algorithmus
2871     \textsc{SN-Evolution-Cut} aus dem \emph{Pairwise-Sorting-Netzwerk}
2872     $\operatorname{PS}(32)$ durch 16~Schnitte erzeugt.}
2873   \label{fig:16-ec-from-ps32}
2874 \end{figure}
2875
2876 Obwohl das \emph{Pairwise-Sorting}-Netzwerk den \emph{Odd-Even}-Mischer nicht
2877 einsetzt und auch nicht auf einem Mischer basiert, ist das
2878 \emph{Odd-Even-Merge}-Netzwerk $\operatorname{OEM}(8,8)$ im Sortiernetzwerk in
2879 Abbildung~\ref{fig:16-ec-from-ps32} eindeutig erkennbar (Schichten~7--10). In
2880 den Schichten~1--6 erkennt man zwei unabhängige Sortiernetzwerke, die
2881 strukturell sehr ähnlich zu $\operatorname{PS}(8)$ sind -- lediglich die
2882 Schichten~1 und~2 sowie 4~und~5 sind vertauscht.
2883
2884 \begin{figure}
2885   \begin{center}
2886     \input{images/32-pairwise-cut-16-pairwise.tex}
2887   \end{center}
2888   \caption{Das \ps{32}-Netzwerk mit 8~Maximum- und 8~Minimumschnitten. Gut zu
2889     sehen sind die verbleibenden Komparatoren, die das \ps{16}-Netzwerk
2890     bilden.}
2891   \label{fig:ps16-from-ps32}
2892 \end{figure}
2893
2894 Für das \emph{Pairwise-Sorting-Netzwerk} ist es vergleichsweise einfach
2895 regelmäßige Schnittmuster anzugeben, die aus dem Netzwerk ein kleineres
2896 schnelles und effizientes Sortiernetzwerk erzeugen. Beispielsweise führt das
2897 einfache Schnittmuster
2898 \begin{displaymath}
2899 \textit{Eingang}_i = \left\{ \begin{array}{rl}
2900   -\infty & \quad \textrm{falls } i < \frac{1}{4} n \\
2901    \infty & \quad \textrm{falls } i \geqq \frac{3}{4} n \\
2902         ? & \quad \mathrm{sonst}
2903   \end{array} \right.
2904 \end{displaymath}
2905 für $\operatorname{PS}\left(n = 2^d\right)$ zum Sortiernetzwerk
2906 $\operatorname{PS}\left(\frac{1}{2}n\right)$. Die Art und Weise, mit der
2907 dieses Schnittmuster Komparatoren eliminiert und welche Komparatoren das
2908 verbleibende Netzwerk ausmachen, ist in Abbildung~\ref{fig:ps16-from-ps32}
2909 dargestellt. Die matt blauen und roten Leitungen und Komparatoren sind
2910 diejenigen, die Aufgrund eines Minimums oder eines Maximums im resultierenden
2911 Netzwerk nicht mehr enthalten sind. Da die Minima und Maxima bereits auf den
2912 „richtigen“ Leitungen angelegt werden, müssen keine Leitungen vertauscht
2913 werden und das Ergebnis ist bereits normalisiert. Daher ist das resultierende
2914 Netzwerk in schwarz gut zu erkennen.
2915
2916 \begin{figure}
2917   \begin{center}
2918     \input{images/16-pairwise.tex}
2919   \end{center}
2920   \caption{Das $\operatorname{PS}(16)$-Sortiernetzwerk mit 8~Schnitten
2921     ($\operatorname{MIN}(0, 2, 4, 6), \operatorname{MAX}(9, 11, 13, 15)$). Das
2922     resultierende 8-Sortiernetzwerk ist $\operatorname{OES}(8)$.}
2923   \label{fig:16-pairwise}
2924 \end{figure}
2925
2926 Ein Spezialfall ergibt sich, wenn \textsc{SN-Evolution-Cut} auf
2927 $\operatorname{PS}(16)$ angewendet wird: In diesem Fall kann ein
2928 8-Schnittmuster ausgegeben werden, das \emph{Odd-Even-Mergesort}-Netzwerk
2929 \oes{8} aus \ps{16} erzeugt.. Für größere Sortiernetzwerke ist dies hingegen
2930 nicht mehr möglich, beispielsweise kann $\operatorname{PS}(32)$ nicht durch
2931 ein 16-Schnittmuster in \oes{16} konvertiert werden. Die Verwandtschaft von
2932 $\operatorname{PS}(n)$ und \oes{n} untersucht \textit{Moritz Mühlenthaler}
2933 ausführlich in~\cite{M2009}.
2934
2935 \newpage
2936 \section{Fazit und Ausblick}
2937
2938 Mit dem Entfernen von Leitungen aus bekannten Sortiernetzwerken lassen sich
2939 interessante Ergebnisse erzielen. Dies zeigte \textit{Moritz Mühlenthaler}
2940 bereits in~\cite{M2009}. Die in dieser Arbeit vorgestellten Methoden und
2941 Resultate machen deutlich, dass sich mit diesem Verfahren noch weitere
2942 interessante Beobachtungen machen lassen.
2943
2944 Das \emph{Odd-Even-Mergesort}-Netzwerk wird sowohl von \textsc{SN-Evolution},
2945 \textsc{SN-Evolution-Cut} und \textsc{SN-Markov} erreicht. Wenn die Anzahl der
2946 Leitungen keine Zweierpotenz ist, kann gegebenenfalls ein schnelleres
2947 Sortiernetzwerk erzeugt werden. Einige Beispiele hierfür wurden in
2948 Abschnitt~\ref{sect:sn-evolution-cut:oes} aufgezeigt.
2949
2950 Das \emph{bitone Mergesort}-Netzwerk kann in Bezug auf Effizienz von den
2951 vorgestellten Algorithmen übertroffen werden. Der Algorithmus
2952 \textsc{SN-Evolution-Cut} kann das Ergebnis von \textit{Mühlenthaler} und
2953 \textit{Wanka} (\cite{MW2010}) für ein 16-Sortiernetzwerk reproduzieren und
2954 für ein 32-Sortiernetzwerk sogar noch übertreffen. Der
2955 \textsc{SN-Evolution}-Algorithmus fand 16-Sortiernetzwerke, die gegenüber dem
2956 Ergebnis von \textsc{SN-Evolution-Cut} beziehungsweise~\cite{MW2010} einen
2957 weiteren Komparator einsparen.
2958
2959 Leider weisen die Sortiernetzwerke, die von den angegebenen Algorithmen
2960 zurückgegeben werden, keine Struktur auf, die sich zur Angabe einer
2961 Konstruktionsanweisung eigenen würde. Für das \emph{Pairwise-Sorting}- und das
2962 \emph{Odd-Even-Mergesort}-Netzwerk mit Zweierpotenzen als Leitungszahl wurden
2963 regelmäßige Schnittmuster angegeben. Diese ergeben Sortiernetzwerke, die so
2964 schnell und effizient sind wie die vergleichbaren \oes{n} und \ps{n}
2965 Netzwerke.
2966
2967 Die Anzahl der \emph{unterschiedlichen} Schnitte von verschiedenen
2968 Sortiernetzwerken wurde experimentell bestimmt und gezeigt, dass es deutlich
2969 weniger \emph{unterschiedliche} als \emph{mögliche} Schnittmuster gibt. Das
2970 bedeutet im Umkehrschluss, dass die gewonnenen Sortiernetzwerke mit mehreren
2971 Schnittmustern erreicht werden können.
2972
2973 Die Möglichkeiten der Optimierung von Sortiernetzwerken mit
2974 \emph{Evolutionären Algorithmen} sind durch die in dieser Arbeit vorgestellten
2975 Herangehensweisen bei weitem nicht erschöpft. Im Folgenden werden Ansätze
2976 umrissen, mit denen an die Untersuchungen in dieser Arbeit nahtlos angeknüpft
2977 werden könnte.
2978
2979 \subsection{Ausblick: Das \textit{Pairwise-Sorting}-Netzwerk und \textsc{SN-Evolution}}
2980
2981 Die aktuelle Implementierung von \textsc{SN-Evolution}
2982 (Abschnitte~\ref{sect:sn-evolution}
2983 beziehungsweise~\ref{sect:implementierung}) kann sowohl den \emph{bitonen
2984 Mischer} als auch den \emph{Odd-Even}-Mischer verwenden, um zwei Individuen zu
2985 rekombinieren. Das \emph{Pairwise-Sorting}-Netzwerk verwendet zwar keinen
2986 Mischer, es ist aber ebenfalls rekursiv über kleinere Versionen von sich
2987 selbst definiert. Das heißt, dass \ps{n} aus zwei Instanzen von
2988 $\ps{\frac{n}{2}}$ und zusätzlichen Komparatoren besteht, die die Eingabe für
2989 die kleineren Sortiernetzwerke vorbereiten und anschließend für eine sortierte
2990 Ausgabe sorgen. Anstelle von $\ps{\frac{n}{2}}$ können beliebige
2991 Sortiernetzwerke mit $\frac{n}{2}$~Leitungen verwendet werden.
2992
2993 Dies ließe sich für \textsc{SN-Evolution} nutzen, um zwei Individuen zu
2994 rekombinieren. Da es für das \emph{Pairwise-Sorting}-Netzwerk sehr viele
2995 \emph{unterschiedliche} Schnittmuster gibt
2996 (Abschnitt~\ref{sect:anzahl_schnittmuster}), ist es möglich, dass die
2997 Verwendung dieser Rekombinationsmethode neue Ergebnisse ermöglicht. Leider
2998 wird die Aussicht auf Erfolg durch die Tatsache geschmälert, dass keine
2999 $n$-Schnittmuster für \ps{2n} gefunden werden konnten, die zu besseren
3000 $n$-Sortiernetzwerken als \ps{n} führen.
3001
3002 \subsection{Ausblick: Kooperation von \textsc{SN-Evolution} und \textsc{SN-Evolution-Cut}}
3003
3004 Ähnlich zu der parasitären \emph{Co-Evolution}, die \textit{W.~Daniel Hillis}
3005 in~\cite{H1990} beschreibt, könnte man versuchen, die Algorithmen
3006 \textsc{SN-Evolution} und \textsc{SN-Evolution-Cut} zu kombinieren. Nach dem
3007 Zusammenfügen von zwei $n$-Sortiernetzwerken könnte der Algorithmus
3008 \textsc{SN-Evolution-Cut} beispielsweise einen möglichst guten Schnitt für
3009 \emph{dieses} Netzwerk ermitteln. Da sich die Lösungen, die Evolutionäre
3010 Algorithmen in ihre Population aufnehmen, in den ersten Schritten rasch
3011 verbessern, könnten selbst weniger Iterationen von \textsc{SN-Evolution-Cut}
3012 die Zwischenlösungen von \textsc{SN-Evolution} deutlich verbessern.
3013
3014 Alternativ könnte man -- analog zur Herangehensweise von \textit{Hillis} --
3015 eine zweite Population von Schnittmustern evolvieren, die für die
3016 Sortiernetzwerke in der Population von \textsc{SN-Evolution} besonders gut
3017 funktionieren. In jeder Iteration wendet man alle oder eine zufällige Menge
3018 Schnittmuster auf das zusammengeführte Netzwerk an und gibt dem besten
3019 Ergebnis den Zuschlag. Anschließend erfährt das entsprechende Schnittmuster
3020 eine Aufwertung, so dass es wahrscheinlicher wird, dass \emph{dieses}
3021 Schnittmuster zur nächsten Generation beiträgt. Im Gegensatz zum Ansatz der
3022 parasitären Eingaben entsteht eine \emph{Synergie} zweier Populationen, die
3023 das Gesamtergebnis oder zumindest die Konvergenzgeschwindigkeit verbessern
3024 könnte.
3025
3026 \newpage
3027 \section{Implementierung}
3028 \label{sect:implementierung}
3029
3030 Alle in dieser Arbeit beschriebenen Versuche wurden mit einer eigens
3031 entwickelten C-Bibliothek, \textit{libsortnetwork}, und zugehörigen
3032 Kommandozeilen-Programmen durchgeführt. Die Bibliothek wurde unter der
3033 \textit{GNU Lesser General Public License} (LGPL) in der Version~2.1
3034 veröffentlicht; die Kommandozeilen-Programme, die in vielen Fällen lediglich
3035 Funktionalität der Bibliothek auf der Kommandozeile zur Verfügung stellen,
3036 stehen unter der \textit{GNU General Public License}, Version~2. Diese
3037 Lizenzen räumen einem Benutzer weitreichende Rechte ein, unter anderem das
3038 Programm beliebig zu verwenden, zu studieren, zu verändern sowie veränderte
3039 und unveränderte Kopien zu veröffentlichen.
3040
3041 Die Programmierschnittstelle (API) der Bibliothek orientiert sich an
3042 Paradigmen der \textit{objektorientierten Programmierung}. Beispielsweise kann
3043 mit der Funktion \texttt{sn\_network\_ create()} ein neues Zustands-Objekt
3044 erzeugt werden, für das mehrere Manipulations-Methoden, zum Beispiel
3045 \texttt{sn\_network\_comparator\_add()}, zur Verfügung stehen. Auf diese Art
3046 und Weise kann die Bibliothek leicht erweitert werden, ohne dass bestehende
3047 Programme angepasst werden müssen.
3048
3049 Die meisten Kommandozeilen-Programmen lesen ein Komparatornetzwerk von der
3050 Standard-Eingabe und schreiben ihr Ergebnis auf die Standard-Ausgabe. Um
3051 Beispielsweise eine \emph{normalisierte} Variante des \emph{bitonen
3052 Mergesort}-Netzwerks \bs{42} zu erzeugen, kann folgendes Kommando verwendet
3053 werden:
3054 \begin{verbatim}
3055   $ sn-bitonicsort 42 | sn-normalize >sn-42
3056 \end{verbatim}
3057 Dieses Prinzip, kleine Programme \emph{eine} Aufgabe erledigen zu lassen und
3058 es einfach zu ermöglichen, Programme zu verketten, ist eines der
3059 Grundprinzipien des UNIX-Be\-triebs\-sys\-tems. Es hat sich in den letzten
3060 Jahrzehnten und beim Verfassen dieser Arbeit als sehr flexibel und mächtig
3061 erwiesen.
3062
3063 Funktionen, die von Kommandozeilen-Programmen zur Verfügung gestellt werden,
3064 sind unter anderem das Erzeugen des \emph{Odd-Even-Mergesort}-, \emph{bitonen
3065 Mergesort}- und \emph{Pairwise-Sorting}-Netzwerks, das Normalisieren von
3066 Sortiernetzwerken, Anwendung von Schnittmustern auf Sortiernetzwerke und
3067 Anwendung eines Komparatornetzwerks auf eine Eingabepermutation. Das
3068 Darstellen von Sortiernetzwerken wird ebenfalls angeboten, beispielsweise
3069 wurden die Sortiernetzwerke in dieser Arbeit mit dem Kommando \texttt{sn-tex}
3070 visualisiert.
3071
3072 \textit{libsortnetwork} kann unter der Web-Adresse
3073 \url{http://octo.it/libsortnetwork/} unentgeltlich heruntergeladen werden.
3074
3075 \newpage
3076 \bibliography{references}
3077 \bibliographystyle{plain}
3078
3079 %\listoffigures
3080
3081 \end{document}
3082
3083 % vim: set shiftwidth=2 softtabstop=2 tabstop=8 fdm=marker tw=78 spelllang=de :