三浦ノート

自分の経験したことを検索可能にしていくブログ.誰かの役に立ってくれれば嬉しいです.

python

剰余の関数プロット

atcoder の D - Remainder Reminder のような,剰余 a%b を調べる問題で a と b のどちらを固定するかで議論の難しさがとても違ったので(自分は最初 a を固定して考えていて分からなくなった),この際,剰余の関数グラフをプロットして見てみる. 特別な結…

2次元イジングモデルのオンサーガー解の数値計算

無限体積の2次元イジングモデルの厳密解であるオンサーガー解を数値的に計算する. 1スピン当たりの自由エネルギー \begin{align} -\beta f &= \log(2\cosh 2K) + \frac{1}{2\pi ^ 2} \int ^ \pi _ 0 d\omega _ 1 \int ^ \pi _ 0 d\omega _ 2 \log(1 - k ^ 2…

numpy.gradientで配列の微分をするとき,両端の微分の精度を上げるには

numpy.gradientで配列を微分するとき,配列の境界部分(始端終端)の値だけ目的の正しい微分値が出ないことがあった. すぐに解決したのでメモです. 境界部分での微分はedge_orderというオプションで精度を選べる. numpy.gradient — NumPy v1.16 Manual 何…

python3 のループ文のelse節についてのチュートリアル

python3.6のチュートリアルでループ制御でのelse節について知ったのでメモです. docs.python.jp ループのelse節の動作についてはチュートリアルに書いてある通りでたぶん理解できたと思います. チュートリアルにあったサンプルプログラムについて考えてみ…