4 次中央差分ルーチン ff_dx_sf のテスト計算をする. 下図に示す振幅の大きさが 1 の sin 関数の微分を計算して, その結果を解析的に微分した関数と比較する. 格子点数を 20, 50, 100 と変えて計算する.
|
|
| 微分される関数 |
| 計算領域 |
arare_sf-diff-ff_dx_sf_center4_R20.nc(微分される関数) arare_ff-diff-ff_dx_sf_center4_R20.nc(微分した関数, 誤差) |
| 全領域 |
arare_sf-diff-ff_dx_sf_center4_R20.nc(微分される関数) arare_ff-diff-ff_dx_sf_center4_R20.nc(微分した関数, 誤差) |
|
|
|
| 数値的に微分した関数 | 解析的に微分した関数 | 数値的に微分した関数の解析的に微分した関数に対する誤差 |
| 計算領域 |
arare_sf-diff-ff_dx_sf_center4_R50.nc(微分される関数) arare_ff-diff-ff_dx_sf_center4_R50.nc(微分した関数, 誤差) |
| 全領域 |
arare_sf-diff-ff_dx_sf_center4_R50.nc(微分される関数) arare_ff-diff-ff_dx_sf_center4_R50.nc(微分した関数, 誤差) |
|
|
|
| 数値的に微分した関数 | 解析的に微分した関数 | 数値的に微分した関数の解析的に微分した関数に対する誤差 |
| 計算領域 |
arare_sf-diff-ff_dx_sf_center4_R100.nc(微分される関数) arare_ff-diff-ff_dx_sf_center4_R100.nc(微分した関数, 誤差) |
| 全領域 |
arare_sf-diff-ff_dx_sf_center4_R100.nc(微分される関数) arare_ff-diff-ff_dx_sf_center4_R100.nc(微分した関数, 誤差) |
|
|
|
| 数値的に微分した関数 | 解析的に微分した関数 | 数値的に微分した関数の解析的に微分した関数に対する誤差 |