def normal_form_game_payoff(): sum_1 = np.exp(lamda * (level_k_probabilities_1(k - 1, (player + 1)%2, lamda2)[0] * normal_form_game_player[0][0][0] +level_k ...
See demo.cpp, where you can input the matrix A (whose Smith Normal Form you wish to compute) and pass it into the calculator, giving its Smith Normal Form Decomposition.