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