Kuramoto/OrderParameter.py at master · AndrewJSON/Kuramoto · 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

'''
 * OrderParameter.py
 *
 *   Created on:         07.02.2019
 *   Author:             Philippe Lehmann
 *
 * General description:
 *   xxx
'''

import numpy as np

class OrderParameter:

    def __init__(self):
        pass


    def SumUpOrderMatrix_Elements(self, _phases):

        OrderMatrix = self.SumUpOrderMatrixArray( _phases )
        OrderParam = np.sum(OrderMatrix)\
                         / (OrderMatrix.shape[0]*OrderMatrix.shape[0])
        return OrderParam


    def SumUpOrderMatrixArray(self, _phases):

        orderMatrixList = self.CreateOrderMatrixArray( _phases )
        orderMatrixShape = orderMatrixList[0].shape
        sumOrderMatrix = np.zeros(orderMatrixShape)

        for matrix in orderMatrixList:
            sumOrderMatrix += matrix
        #print(_omList)
        #print(_phases.shape[0])
        #print(sumOrderMatrix[49].shape)


        sumOrderMatrix = np.exp(1j*sumOrderMatrix)
        #sumOrderMatrix = np.abs(sumOrderMatrix)
        sumOrderMatrix = sumOrderMatrix / _phases.shape[0]

        #print(sumOrderMatrix)

        return sumOrderMatrix


    def CreateOrderMatrixArray(self, _phases):

        MatrixArray = []
        timeStepCount = _phases.shape[0]

        for ts in range(0,timeStepCount):
            singleTimeStepPhaseDiff = _phases[ts] - _phases[ts][:, np.newaxis]
            MatrixArray.append( singleTimeStepPhaseDiff )

        #print(timeStepCount)
        #print(OrderMatrix[999].shape)
        #print(OrderMatrix)

        return np.array(MatrixArray)

''' END '''