-
Notifications
You must be signed in to change notification settings - Fork 1
/
arnold.py
37 lines (31 loc) · 1.23 KB
/
arnold.py
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
import matplotlib.pyplot as plt, copy
def arnold( inn, iterr ):
if (inn.ndim != 2):
error('Oly two dimensions allowed');
(m, n) = inn.shape;
if (m != n):
error('Arnold Transform is defined only for squares. Please complete empty rows or columns to make the square.');
out=copy.deepcopy(inn);
#print(type(out))
for c in range (1,iterr):
for i in range (0,n):
for j in range (0,n):
i_1=((i*1+j*1)%m);
j_1=((i*1+j*2)%m);
out[i_1, j_1] = inn[i, j];
return out;
def iarnold( inn, iterr ):
if (inn.ndim != 2):
error('Oly two dimensions allowed');
(m, n) = inn.shape;
if (m != n):
error('Arnold Transform is defined only for squares. Please complete empty rows or columns to make the square.');
out=copy.deepcopy(inn);
#print(type(out))
for c in range (1,iterr):
for i in range (0,n):
for j in range (0,n):
i_1=((i*2+j*(-1))%m);
j_1=((i*(-1)+j*1)%m);
out[i_1, j_1] = inn[i, j];
return out;