コンデンサの電位を計算

Fortran95

パリティ Vol.06 No.03 (1991) 62. のフォートラン入門より、平行平板コンデンサの作る二次元電場の電位を、ラプラス方程式のヤコビによる解法で解きます。forallを用い mask 配列を使ってコンデンサの電極は計算を飛ばすことにします。forall は並列代入操…

2012-04-24

Getting started in CM Fortran のサンプルプログラム

Fortran95

http://www.sthmuseum.org/downloads/CM5/GettingStartedinCMFortran.pdf 42ページ Primes:Second CM Fortran Versionこれは Fortran95 の文法にかなっているので無修正で動きます。今のコンピュータだとシリアル実行でも一瞬ｗ実行結果長いので見切れ画像…

2012-04-20

Fortranに於けるfilterとmap

Fortran95

Fortranは配列がいわゆる first class object になっており、配列の扱いに長けています。Fortran90 以降、ほとんどの組み込み関数は elemental 属性を持っており、定義上スカラーを引数として取る関数に配列を渡すと、結果を配列として返します。これは、LIS…

2012-01-12

partition function p(n)

Fortran95

整数の可能な分割数 p(n) は漸化式で求められます。p(n)はｎマスのヤング図（フェラーズ盤）の可能な数に対応しています。アンドリュース＆エリクソンの『整数の分割』という本に漸化式が導出されています（第五章四節）。結果だけ書くと、これより (p(0)=1…

2011-11-27

LU分解 Doolittle 法

Fortran95

基本的にガウス消去法と同じ方式になっています。実行結果 LU分解では精度がやや悪くなるようです。 Eliza’s father=Alfred P.DoolittleStation master=AlfeeYoung Freddy=Sherlock Holmes ソース・プログラム簡単に書けるのがうれしい。ほとんど実質３〜…

2011-11-06

オンライン Fortran 実行環境

Fortran95

sourcelair というオンライン実行環境が開発中のようで、Fortranも使えるようです。登録すると、ファイルの保存などもできます。今のところ設定ミスか Module が使えないようです。H23-11-08 メールしてみたところ、今のところ全言語で Module は使えないそ…

2011-10-31

挿入ソート改良

Fortran95

ここまで来ると更に簡単に。尻からではなく、頭から sort してゆく。 module m_isort2 implicit none contains subroutine isort3y( x ) real, intent(in out) :: x(:) integer :: j, k do k = 2, size(x) do j = 1, k - 1 if ( x(k) < x(j) ) exit end do x…

2011-10-31

挿入ソート

Fortran95

少し改良。配列のサイクリックシフトで。 module m_isort2 implicit none contains subroutine isort2y( x ) real, intent(in out) :: x(:) integer :: j, k, n n = size(x) if (n <= 1) return do k = n - 1, 1, -1 do j = k + 1, n if ( x(k) < x(j) ) exi…

2011-10-12

Modern Fortran Explained Chap.8 Exercises 3, 4

Fortran95

program exercise3 implicit none ! Variables real, target :: a(5:10) real, pointer :: p1(:), p2(:) ! Body of test5 p1 => a p2 => a(:) print *, associated(p1, a), lbound(p1) print *, associated(p2, a), lbound(p2) print *, associated(p1, p2),…

2011-10-11

メモ帳

Fortran95

program test4 implicit none ! Variables real, allocatable :: x(:, :) ! Body of test4 allocate( x(100, 100) ) x = 1.0 call desub(x) print *, allocated(x) call alsub(x) print *, allocated(x) print *, x stop contains subroutine desub(x) real,…