nano_exit

基礎的なことこそ、簡単な例が必要だと思うのです。

2019-05-01から1ヶ月間の記事一覧

高速フーリエ変換する時の端点の話。

高速フーリエ変換ことFFT(ファイナルファンタジータクティクスではありません)は、結局は離散フーリエ変換な訳で、ある範囲のデータが無限に周期的に繰り返しているとしてフーリエ変換してくれます。 この時に、端の点のどこで周期性を課すのか、ややこし…

FFTWで3次元高速フーリエ変換してみた。

前回の例を元に、3次元に拡張してみた。 koideforest.hatenadiary.com program fftw3_test_3d implicit none integer, parameter :: WP = selected_real_kind( 8 ) real( WP ), parameter :: PI = 3.141592653589793 complex( WP ), parameter :: I = ( 0, 1…

FFTWをfortranで使う際の自分用のメモ

FFTWをfortranで使う。本家の説明: FFTW 3.3.8: Fortran Examples参考にしたサイト: fftwの使い方:腰も砕けよ 膝も折れよ:So-net blog FortranのFFTWの使い方 - Qiita Hiroyuki R. Takahashi Fumio KUSUHARA -Fortran90で構造解析- program fftw3_test i…

Fortranの精度変更。

Fortranの精度変更の書式が自分の知らない感じだったので、自分用にまとめる。 Fortran Best Practices — Fortran90 1.0 documentation program test_presicion implicit none integer, parameter :: dp0 = kind( 0.d0 ) integer, parameter :: dp = selecte…

Makefileの自分用のメモ

FC = gfortran FFLAGS = INCLUDE = -I/opt/local/include LIBS = -L/opt/local/lib -lbrabra PROGRAM = main.out OBJS = sub.o SRCS = $(OBJS:%.o=%.f90) .Phony: all clean all: $(PROGRAM) $(PROGRAM): $(OBJS) $(FC) $(FFLAGS) -o $@ $^ $(INCLUDE) $(LIB…

pythonで三次元プロットする時の自分用メモ

三次元プロット用のx軸、y軸の配列。 import numpy as np N = 100 x_min, x_max = -1, 1 y_min, y_max = -1, 1 x1 = np.linspace( x_min, x_max, N ) y1 = np.linspace( y_min, y_max, N ) X, Y = np.meshgrid( x1, y1 ) 三次元プロット用のz軸データ。 def …