B.16.3 式

• 論理式の値
  • .FALSE. は 0
  • .TRUE. は 0 以外の任意の整数値
• 論理演算子
  • .not. は !
  • .and. は &&
  • .or. は ||
• 比較演算子
  • .eq. は ==
  • .ne. は !=
  • .gt. は >
  • .ge. は >=
  • .lt. は <
  • .le. は <=
• 配列について。
  • 括弧は (, ) から [, ] に

            real a(0:NN)
            ..
            x = a(3)
            ...
            ↓
            double a[NN+1];
            ..
            x = a[3];
    

  • FORTRAN では配列の添字は、任意の整数から任意の整数までと出来る。デフォー ルトでは 1 から。C では 0 から。

            real a(N)
            ...
            do i=1,N
              a(i)=0.0
            end do
            ...
            x=a(3)
            ↓
            double a[N];
            ...
            for (i=0; i<N; i++)
              a[i] = 0.0;
            ...
            x = a[2];
            または
            double a[N+1];
            ...
            for (i=1; i<=5; i++)
              a[i] = 0.0;
            ...
            x = a[3];
    

  • 多次元配列は FORTRAN では , で区切る。C では、[] を繰り返す。

            real mat(3,3)
            x = mat(i,j)
            ↓
            double mat[3][3];
            x = mat[i][j];
    

  • FORTRAN の多次元配列には整合配列の機能があるが、C にはない。 C で同じことをするには、次のいずれかの方法を用いる。
    • 十分大きな 2 次元配列を宣言して、その一部を使う。(「十分」これは関数 の可搬性が損なわれる。例えば、次の inv() はライブラリィに登録すること は出来ない。)

              real a(NDIM,NDIM)
              ...
              a(i,j)=むにゃ
              ...
              call inv(NDIM,3,a)
              ..
              subroutine inv(NDIM, N, a)
              integer NDIM,N
              real a(NDIM,*)
              ...
              ↓
              double a[NBIG][NBIG];
              ...
              a[i][j] = むにゃ
              ...
              inv(3,a);
              ...
              void inv(int n, double a[][NBIG])
      

    • 一次元配列を宣言して（あるいはポインターで同様のことをする）、添字 計算を自分でやる。（これはかなり面倒くさい。マクロを使って処理すべき？）

              double *a;
              ...
              a = malloc(sizeof(double) * n * n); if (a == NULL) ...
              ...
              a[i*3+j] = むにゃ
              ...
              inv(3,a);
              ...
              void inv(int n, double *a)
      

    • ポインター配列（あるいはポインターのポインター）を用いる。

              double *a[NDIM];
              for (i = 0; i < 3; i++) {
                 if ((a[i] = malloc(sizeof(double) * 3)) == NULL) {
                 }
              }
              ...
              a[i][j] = むにゃ
              ...
              inv(3,a);
              ...
              void inv(int n, double **a)
      

  • FORTRAN では、そのプログラム内で使う大きな配列は、主プログラムで宣言 するのが普通だが、C では動的にポインターで確保するのが普通。

         real a(0:NN,M)
         ...
         x = a(0,j)
         ↓
         double a[NN+1][M];
         ...
         x = a[0][j];
  

• FORTRAN の stop 文は C では exit() 関数の呼び出しに書き換える。 (exit() 関数の引数には小さな正整数を与えるが、正常終了の場合は 0, 異常 終了の場合は 0 以外の値を返すのが慣習である。)

          stop
         ↓
          exit(0);           正常終了の場合
  あるいは
          exit(1);           異常終了の場合
  

• FORTRAN の冪乗演算子は適当に

          L=j**2
          n=j**3
          m=j**4
          x=y**0.5
          z=y**(1.0/3.0)
  	↓
          L=j*j;
          n=j*j*j;
          tmp = j*j; m=tmp*tmp;
          x=sqrt(x);
          z=pow(x,1.0/3.0)
  

• FORTRAN の format 文の E, F, G 変換はそれぞれ printf() の書式文字列 内では %e, %f, %g に。 I 変換は printf() の書式文 字列内では %d に。
• FLOAT() のような型変換関数は何もしないでよい時もあるし、必要な場合は キャスト演算子で。

          integer n;
          real x,r
          ...
          x=float(n)
          r=float(1)/float(n)
          call mysub(a, float(n))
          ↓
          int n;
          double x,r
          ...
          x=n;                     自動的に型変換される。
          r=1.0/n;                 1.0 に合わせて n は double に型変換される
          mysub(a, (double)n);     もしも mysub() にプロトタイプ宣言があれば
                                   キャストは不要である。