kkatoのブログ

1日で使ったハードディスクの量の自己記録を更新しました。

あるプログラムをデバッグしていて、デバッグ出力をさせながらプログラムを動かしてそのまま昨日夜帰りました。
今日朝来てみると、ディスク容量不足で落ちてます。一晩で使った量は180Gバイト。

とほほ

C言語からfortranを呼び出すためのメモ。

普通の変数は参照渡しだからCではポインタになる、とかという点はいいのだが、２次元配列の渡しかたが難しい。

一般的な配列とparameterでは並び順が違うようなので気をつけなければならない。

以下サンプルコード：

a.f:

      module test
      integer i,j
      integer s
      integer,parameter:: param(2,5)= (/ 0,1,2,3,4,5,6,7,8,9 /)
      end module test
      use test
      integer a(2,5)
      integer b(2,5)
      s=0;
      do i=1,2
         do j=1,5
            a(i,j)=s
            s=s+10
         end do
      end do
      call f(a,b)
      do i=1,2
         do j=1,5
            write(*,*) 'b(',i,',',j,') = ',b(i,j)
         end do
      end do
      end

f.c:

#include <stdio.h>

extern test_mp_param_[][5];

void f_(int a[][2], int b[][2])
{
  int i,j;
  for (i=0; i<5; i++) {
    for (j=0; j<2; j++) {
      printf("a[%d][%d] = %d\n",i,j,a[i][j]);
      b[i][j] = a[i][j] + test_mp_param_[j][i];
    }
  }
  
}

a[i][j]とtest_mp_param_[j][i]で添字が逆になっている点に注意。
これのa.fとf.cをコンパイルしてリンクして実行する。ちなみにテストはgccとifortでやった。
以下実行結果：

a[0][0] = 0
a[0][1] = 50
a[1][0] = 10
a[1][1] = 60
a[2][0] = 20
a[2][1] = 70
a[3][0] = 30
a[3][1] = 80
a[4][0] = 40
a[4][1] = 90
 b(           1 ,           1 ) =            0
 b(           1 ,           2 ) =           11
 b(           1 ,           3 ) =           22
 b(           1 ,           4 ) =           33
 b(           1 ,           5 ) =           44
 b(           2 ,           1 ) =           55
 b(           2 ,           2 ) =           66
 b(           2 ,           3 ) =           77
 b(           2 ,           4 ) =           88
 b(           2 ,           5 ) =           99


この仕様、かなり気持ち悪いのだがなんとかならないものか。
おそらく歴史的事情やなにやらがあるんだろうが。

必要なライブラリでいろいろとうるさく言われたのでメモ。

sudo aptitude install libstdc++5 ia32-libs
とすればいいようです。

Tom Marioni

The Act of Drinking Beer with Friends Is the Highest Form of Art

作品名訳：「友達とビールを飲むという行為は芸術の最高の形態である」

そういえば私の友達でもこういうやついました。飲んだビールの缶を部屋に並べてたやつ。

彼もアーティストか。

食事のメニューを見て思ったこと。

やっぱり都市伝説だったんですね。安心しました。

メキシコのカンクンに仕事で来ています。 入国審査でbusiness!と答えるのが、はばかられるようなそんなところです。
目的は国際会議出席と、ことだまの開発ワークショップです。