SedimentHydraulics https://blog.hatena.ne.jp/SedimentHydraulics/ 趣味で計算流砂水理 Computational Sediment Hydraulics for Fun Learning https://computational-sediment-hyd.hatenablog.jp/ R 循環シフトとは 配列の最初と最後を循環してシフトするものです。 pythonのnumpy.rollで書くとこんな感じになります。 import numpy as np x = np.arange(5) x # array([0, 1, 2, 3, 4]) np.roll(x,2) # array([4, 0, 1, 2, 3]) np.roll(x,-1) # array([1, 2, 3, 4, 0]) いまいち何に使うのかわかりませんが(ビット演算とかでしょうか)、これを使えば周期境界の離散式がスッキリかけると思います。 ちなみにfortranだとcshiftです。 python nump… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fcomputational-sediment-hyd.hatenablog.jp%2Fentry%2F2020%2F10%2F06%2F210000" title="周期境界は配列の循環シフトを使えばソースコードがスッキリするかも - 趣味で計算流砂水理 Computational Sediment Hydraulics for Fun Learning " class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> https://chart.apis.google.com/chart?cht=tx&chl=%0A%5Cbegin%7Balign%7D%20%0A%20%20%20%20%5Cdfrac%7B%5Cpartial%20u%20%7D%7B%5Cpartial%20t%20%7D%20%26%3D%20-%20c%5Cdfrac%7B%5Cpartial%20u%20%7D%7B%5Cpartial%20x%20%7D%20%20%5C%5C%0A%20%20%20%20c%20%26%3D%20%5Crm%7Bconst%7D%0A%5Cend%7Balign%7D%0A Hatena Blog https://hatena.blog 2020-10-06 21:00:00 rich https://computational-sediment-hyd.hatenablog.jp/entry/2020/10/06/210000 1.0 100%