RG/RG_dict.py at main · cchandre/RG · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
##################################################################################################
##               Dictionary of parameters: https://github.com/cchandre/RG                       ##
##################################################################################################

import numpy as xp

Method = 'region'
Iterates = 20

Nxy = 500
RelDist = 1e-7

## 2D -- golden mean
N = [[1, 1], [1, 0]]
omega0 = [(xp.sqrt(5)-1)/2, -1]
Omega = [1, 0]
K = ((0, 1, 0), (0, 1, 1))
AmpInf = [0, 0]
AmpSup = [0.12, 0.35]
CoordLine = [0.0, 0.028]
ModesLine = (1, 1)
DirLine = [1, 1]

## 3D -- spiral mean
# N = [[0, 1, 0], [0, 0, 1], [1, 0, -1]]
# sigma = 1.3247179572447460259
# omega0 = [sigma, sigma ** 2, 1.0]
# Omega = [1, 1, -1]
# K = ((0, 1, 0, 0), (0, 0, 1, 0), (0, 0, 0, 1))
# AmpInf = [0, 0, 0]
# AmpSup = [0.05, 0.25, 0.1]
# CoordLine = [0, 0.05]
# ModesLine = (1, 1, 0)
# DirLine = [1, 5, 0.1]

## 3D -- tau mean
# N = [[0, 1, -1],[1, -1, 1],[0, -1, 2]]
# tau = 0.445041867912629
# omega0 = [1.0, tau, 1.0 - tau - tau**2]
# Omega = [1, 1, -1]
# K = ((0, 0, -1, 1), (0, 1, -1, -1), (0, 0, 0, 1))
# AmpInf = [0.0, 0.0, 0.01]
# AmpSup = [7e-4, 6e-3, 0.01]
# CoordLine = [0.0, 0.002]
# ModesLine = (1, 1, 0)
# DirLine = [1, 5, 0.01]

L = 5
J = 5

ChoiceIm = 'AK2000'
Sigma = 0.6
Kappa = 0.1

CanonicalTransformation = 'expm_multiply'
MinStep = 0.05
AbsTol = 1e-2
RelTol = 1e-3

TolMax = 1e+4
TolMin = 1e-10
MaxLie = 5000

NormChoice = 'sum'
NormAnalytic = 1
Precision = 64

SaveData = True
PlotResults = False
Parallelization = (True, 100)

##################################################################################################
##                              DO NOT EDIT BELOW                                               ##
##################################################################################################
Precision = {32: xp.float32, 64: xp.float64, 128: xp.float128}.get(Precision, xp.float64)
param_dict = {'Method': 'compute_' + Method}
param_dict.update({
        'Iterates': Iterates,
        'Nxy': Nxy,
        'RelDist': RelDist,
        'N': xp.asarray(N, dtype=int),
        'omega0': xp.asarray(omega0, dtype=Precision),
        'Omega': xp.asarray(Omega, dtype=Precision),
        'K': K,
        'ModesK': [K[_] for _ in xp.nonzero(ModesLine)[0]],
        'AmpInf': AmpInf,
        'AmpSup': AmpSup,
        'CoordLine': CoordLine,
        'ModesLine': xp.asarray(ModesLine),
        'DirLine': xp.asarray(DirLine),
        'L': L,
        'J': J,
        'ChoiceIm': ChoiceIm,
        'Sigma': Sigma,
        'Kappa': Kappa,
        'CanonicalTransformation': 'self.' + CanonicalTransformation,
        'MinStep': MinStep,
        'AbsTol': AbsTol,
        'RelTol': RelTol,
        'TolMax': TolMax,
        'TolMin': TolMin,
        'MaxLie': MaxLie,
        'NormChoice': NormChoice,
        'NormAnalytic': NormAnalytic,
        'Precision': Precision,
        'SaveData': SaveData,
        'PlotResults': PlotResults,
        'Parallelization': Parallelization})
##################################################################################################