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