リー代数の初歩として， -1 から 1 までを要素に持つ10000個のランダム行列についてヤコビの恒等式が成り立っていることを確かめるプログラムを書いた． ヤコビの恒等式： import numpy as np from sympy import * import time import random def com(A, B):…

0 準備 型の確認． type(None), type(0) (NoneType, int) 適切なパラメータを設定する． N = int(1e7) M = int(1e4) インポートも一応計測する． %%timeit import numpy as np #112 ns ± 0.579 ns per loop (mean ± std. dev. of 7 runs, 10000000 loops eac…

意外と詰まったし良いまとめがないので書く． 1 一番近い値の重複を考えない場合 1.1 最もシンプルなバージョン 探したい値が一つで，返したい index も一つの場合．これが最も簡単．探索されるデータは1次元配列を想定． import numpy as np def idx_of_the…

import %matplotlib inline import numpy as np import math from sympy import * from sympy.combinatorics import * init_printing(pretty_print=False) 関数 Python de 対称群 2 - Krypf’s Diary 参照． def plusone(x): add_one = [(i > 0) * (list(x)[i…

plusone() 関数の改良 どうもどこかで見たはずの，前回の tuple でやる方式は，スタンダードではないらしいので，list に変えて，更に拡張・改良する． 1 前回の関数 以前の関数．出力はコメントアウトで表現．名前は plusone_tuple() に変えた． import num…

Python 対称群 プログラムはデフォルトだと0から数え始めるので，普段使う1からの表示にする関数を作った． SymmetricGroup(3)[1] #Permutation(2)(0, 1) def plusone(x): X = tuple([0]) for i in range(n): a = tuple([tuple(x)[i] + 1]) X = X + a return…