Beispiele ========= 3 Magnete --------- Die folgenden Abbildungen wurden aus der Datei `simple.xml` erstellt. Drei gleiche Magnete befinden sich unterhalb des Pendelkörpers in gleichem Abstand von der Ruhelage. Der Reibungskoeffizient beträgt :math:`\gamma = 0.1`, der Potentialexponent ist X = 2. In horizontaler Richtung ist :math:`\phi = [0, 2\pi]`, in vertikaler Richtung ist :math:`\vartheta = [\frac{\pi}{2}, \pi]` aufgetragen. Das entspricht gerade der unteren Halbkugel, in der sich der Pendelkörper bewegen kann. Einfluß der Auflösung auf die Rechenzeit ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ Die Abbildungen unterscheiden sich jeweils nur in der gewaehlten Aufloesung: von 400x400 (atol=rtol=1e-2) mit einer Rechenzeit von 91 Minuten über 400x400 (atol=rtol=1e-5) bei 126 Minuten und 800x800 bei 384 Minuten bis zu 1600x1600 mit 2663 Minuten! Die erste Abbildung zeigt außerdem, was für zu große Toleranzen bei der Integration passiert: es entsteht zu viel Rauschen. .. image:: ../plots/simple/400x400.png :width: 150pt .. image:: ../plots/simple/400x400-detail.png :width: 150pt .. image:: ../plots/simple/800x800.png :width: 150pt .. image:: ../plots/simple/1600x1600.png :width: 150pt Bei Projektion der Abbildungen auf die Kugeloberfläche entstehen anschaulichere Darstellungen: .. image:: ../plots/simple/400x400-sphere.png :width: 150pt .. image:: ../plots/simple/400x400-detail-sphere.png :width: 150pt .. image:: ../plots/simple/800x800-sphere.png :width: 150pt .. image:: ../plots/simple/1600x1600-sphere.png :width: 150pt Berechnet man nicht nur die untere Halbkugel, sondern auch die obere, so erhält man die folgende Karte und ihre Projektion auf die obere Halbkugel (Rechenzeit: 1024 Minuten): .. image:: ../plots/simple/simple-fullsphere.png :width: 300pt .. image:: ../plots/simple/simple-fullsphere-uppersphere.png :width: 150pt Einfluß der Reibung ^^^^^^^^^^^^^^^^^^^ Bei einem geringen Reibungskoeffizienten :math:`\gamma` braucht das Pendel eine vergleichsweise lange Zeit, um zur Ruhe zu kommen. Auf dem Weg dahin kann es sich mehrmals von Magnet zu Magnet umentscheiden. Insbesondere werden die Karten für größere Reibungskoeffizienten regelmäßiger. :math:`\gamma = 0.15` (800x800, 541 Minuten): .. image:: ../plots/simple/simple-fric0.15.png :width: 150pt .. image:: ../plots/simple/simple-fric0.15-sphere.png :width: 150pt :math:`\gamma = 0.20` (800x800, 612 Minuten): .. image:: ../plots/simple/simple-fric0.20.png :width: 150pt .. image:: ../plots/simple/simple-fric0.20-sphere.png :width: 150pt :math:`\gamma = 0.25` (800x800, 424 Minuten): .. image:: ../plots/simple/simple-fric0.25.png :width: 150pt .. image:: ../plots/simple/simple-fric0.25-sphere.png :width: 150pt :math:`\gamma = 0.30` (800x800, 395 Minuten): .. image:: ../plots/simple/simple-fric0.30.png :width: 150pt .. image:: ../plots/simple/simple-fric0.30-sphere.png :width: 150pt :math:`\gamma = 0.50` (1600x1600, 1300 Minuten): .. image:: ../plots/simple/simple-fric0.50-huge.png :width: 150pt .. image:: ../plots/simple/simple-fric0.50-huge-sphere.png :width: 150pt :math:`\gamma = 0.60` (800x800, 404 Minuten): .. image:: ../plots/simple/simple-fric0.60.png :width: 150pt .. image:: ../plots/simple/simple-fric0.60-sphere.png :width: 150pt :math:`\gamma = 1.00` (800x800, 268 Minuten): .. image:: ../plots/simple/simple-fric1.00.png :width: 150pt .. image:: ../plots/simple/simple-fric1.00-sphere.png :width: 150pt Einfluß der Magnetstärke ^^^^^^^^^^^^^^^^^^^^^^^^ Die folgende Karte unterscheidet sich von den vorherigen nur dadurch, daß die Magnetstärken wie folgt definiert sind: 0.9 (rot), 1.0 (grün), 1.1 (blau). .. image:: ../plots/strengths/strengths-800x800.png :width: 150pt .. image:: ../plots/strengths/strengths-800x800-sphere.png :width: 150pt (Rechenzeit: etwa 960 Minuten, Überlauf in der Zeitmessung ;-)) Einfluß des Abstands ^^^^^^^^^^^^^^^^^^^^ Nun wird der Abstand der Magnete von der Ruhelage variiert, so daß der rote, grüne und blaue Magnet in einem Abstand von 0.9, 1.0 und 1.1 angebracht ist. .. image:: ../plots/distances/distances-800x800.png :width: 150pt .. image:: ../plots/distances/distances-800x800-sphere.png :width: 150pt (Rechenzeit: etwa 800 Minuten, Überlauf der Zeitmessung mit `time.clock()`) Einfluß des Potentialexponenten ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ Karte mit einem Potentialexponenten von :math:`X = 3` (Rechenzeit etwa 840 Minuten): .. image:: ../plots/simple/simple-x3.png :width: 150pt .. image:: ../plots/simple/simple-x3-sphere.png :width: 150pt 4 Magnete --------- Nun sind vier Magnete unterhalb des Pendels aufgebaut. :math:`\gamma = 0.2` (800x800, 1699 Minuten): .. image:: ../plots/four/four-0.2.png :width: 150pt .. image:: ../plots/four/four-0.2-sphere.png :width: 150pt :math:`\gamma = 0.3` (800x400, 1278 Minuten): .. image:: ../plots/four/four-0.3-fullsphere.png :width: 300pt .. image:: ../plots/four/four-0.3-fullsphere-sphere.png :width: 150pt .. image:: ../plots/four/four-0.3-fullsphere-uppersphere.png :width: 150pt :math:`\gamma = 0.4` (800x400, parallel gerechnet, 11 Knoten, 48 Minuten): .. image:: ../plots/parallel/four-0.4-fullsphere.png :width: 300pt .. image:: ../plots/parallel/four-0.4-fullsphere-sphere.png :width: 150pt .. image:: ../plots/parallel/four-0.4-fullsphere-uppersphere.png :width: 150pt :math:`\gamma = 0.5` (800x400, parallel gerechnet, 10 Knoten, 43 Minuten): .. image:: ../plots/parallel/four-0.5-fullsphere.png :width: 300pt .. image:: ../plots/parallel/four-0.5-fullsphere-sphere.png :width: 150pt .. image:: ../plots/parallel/four-0.5-fullsphere-uppersphere.png :width: 150pt 6 Magnete --------- Sechs Magnete unterhalb des Pendels. :math:`\gamma = 0.2` (800x800, parallel gerechnet, 10 Knoten, 77 Minuten): .. image:: ../plots/parallel/six-0.20.png :width: 150pt .. image:: ../plots/parallel/six-0.20-sphere.png :width: 150pt :math:`\gamma = 0.3` (800x800, parallel gerechnet, 10 Knoten, 63 Minuten): .. image:: ../plots/parallel/six-0.3.png :width: 150pt .. image:: ../plots/parallel/six-0.3-sphere.png :width: 150pt :math:`\gamma = 0.4` (800x800, parallel gerechnet, 10 Knoten, 58 Minuten): .. image:: ../plots/parallel/six-0.4.png :width: 150pt .. image:: ../plots/parallel/six-0.4-sphere.png :width: 150pt :math:`\gamma = 0.5` (800x800, parallel gerechnet, 10 Knoten, 61 Minuten): .. image:: ../plots/parallel/six-0.5.png :width: 150pt .. image:: ../plots/parallel/six-0.5-sphere.png :width: 150pt :math:`\gamma = 0.6` (400x400, parallel gerechnet, 20 Knoten, 23 Minuten; Ausfall von 2 Knoten während der Berechnung!): .. image:: ../plots/parallel/six-0.6-small.png :width: 150pt .. image:: ../plots/parallel/six-0.6-small-sphere.png :width: 150pt