ここまで来ると更に簡単に。尻からではなく、頭から 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…

少し改良。配列のサイクリックシフトで。 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…

配列生成子を使ったものは、配列はみ出しの場合分けをしなくて済むが、死ぬほど遅いｗ 素朴な方がいい。 module m_isort2 contains subroutine isort2( x ) real, intent(in out) :: x(:) real :: tmp integer :: i, k, n n = size(x) if (n <= 1) return do…

再帰を使ったので、すぐスタックあふれで死ぬｗ コンパイラオプション /assume:realloc_lhs が必要。 module m_isort implicit none contains recursive function isort(arr) result(res) real, intent(in) :: arr(:) real, allocatable :: res(:) if (size(…

今度は $15,000 に引き上げ！リーマンショックの時も $6,000 くらいに引き上げて学生とか貧民からド顰蹙くらってたのに。 大貧民・大富豪を地で行く展開ﾜﾛﾀ

ベストセラーになっただけあって読みやすく面白い。 華々しいところしか報じないジャーナリストへの不満が軍人側から何度か表明されるのも良い。そのためか縁の下の力持ちとして地味な仕事をしている兵隊さんの姿も描いている。また腐敗しきった支那人の有様…

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),…

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,…

Hamlet: >it has made me mad. I say we will have no >moe marriage. Moe Marriage まぁ moe は more の崩れなんですけど。