|
3 | 3 | Calculates the Pythagorean 3-body problem using different values for |
4 | 4 | the smoothing length in the n-body code. |
5 | 5 | """ |
6 | | -from __future__ import print_function |
| 6 | +# from __future__ import print_function |
7 | 7 | # import numpy |
8 | 8 | # import time |
9 | 9 |
|
@@ -57,11 +57,11 @@ def run_pyth(interface, tend=100, dt=0.125, parameters=[]): |
57 | 57 |
|
58 | 58 |
|
59 | 59 | if __name__ == "__main__": |
60 | | - codes_to_run = [('Hermite0, $\eta=0.03$', Hermite, [("dt_param", 0.03)]), |
61 | | - ('Hermite0, $\eta=0.01$', Hermite, [("dt_param", 0.01)]), |
62 | | - ('Hermite0, $\eta=0.003$', Hermite, [("dt_param", 0.003)]), |
63 | | - ('Hermite0, $\eta=0.001$', Hermite, [("dt_param", 0.001)])] |
64 | | - N = (len(codes_to_run) - 1) / 2 + 1 |
| 60 | + codes_to_run = [(r'Hermite0, $\eta=0.030$', Hermite, [("dt_param", 0.030)]), |
| 61 | + (r'Hermite0, $\eta=0.010$', Hermite, [("dt_param", 0.010)]), |
| 62 | + (r'Hermite0, $\eta=0.003$', Hermite, [("dt_param", 0.003)]), |
| 63 | + (r'Hermite0, $\eta=0.001$', Hermite, [("dt_param", 0.001)])] |
| 64 | + N = (len(codes_to_run) - 1) // 2 + 1 |
65 | 65 | f = pyplot.figure(figsize=(8, 4 * N)) |
66 | 66 |
|
67 | 67 | for i, (label, interface, parameters) in enumerate(codes_to_run): |
|
0 commit comments