FORTRANとC言語プログラムで値の受け渡しがうまく行かないです。

前提・実現したいこと

cos(x)のxがa~bの面積を、矩形分積で求めるプログラムで、計算はFORTRAN側で行いたいです。
また、実行時にコマンドライン引数a,bを指定することで面積を求めるように実現したいです。
計算やコマンドラインからの受け取りは実現できているのですが、FORTRANへの受け渡しがうまくいかないです。。

ソースコード

FORTRAN77
1      DOUBLE PRECISION FUNCTION SEKIBUN(A, B) 
2      DOUBLE PRECISION X
3      DOUBLE PRECISION I
4      DOUBLE PRECISION AREA 
5      DOUBLE PRECISION DIV
6
7      X=B-A  
8      DIV=100000000.0
9      I=X/DIV
10
11      WRITE (*, *) A 
12      WRITE (*, *) B 
13
14      DO WHILE ( I .LT. B )     
15         AREA=AREA+((DCOS(I))*(X/DIV))
16         I=I+(X/DIV)
17      END DO
18
19      SEKIBUN=AREA
20
21      RETURN
22      END
23

C
1#include <stdio.h>
2#include <stdlib.h>
3
4extern double sekibun_(double *, double *);
5
6int main(int argc,char *argv[]){
7    double a_c,b_c,area;
8    a_c = atof(argv[1]);
9    b_c = atof(argv[2]);
10    printf("a_c==%lf,,,b_c==%lf\n",a_c,b_c);
11
12    area = sekibun_(&a_c, &b_c);
13    printf("面積%25.15lf\n",area);
14
15    return 0;
16}
17

出力結果(ターミナル）

 $ ./sekibun 0 1
a_c==0.000000,,,b_c==1.000000
   0.00000000    
   1.87500000    
面積        0.954085769656684

コマンドラインからの受け取りや積分の計算そのものはできているのですが、CからFORTRANの関数に値を渡す際に、1.0が1.875000のように変わってしまいます。

試したこと

printfで細々とデバッグしたのですが、関数に値を渡す直前ではしっかりと値も合っていて、渡した直後に値が変わってしまいます。。

行動規範の内容に同意します

回答1件

ベストアンサー

引数の型（バイト数）が一致していないことが原因だと思います。

値の引き渡す元（ C 側 )では倍精度( double )を宣言していますが、
受け側（ FORTRAN77 )での、引数 A,B に対しては、
倍精度( DOUBLE PRECISION )の宣言が見当たりません。

このため、FORTRAN77 側では、変数 A,B に対して暗黙の型宣言により
REAL 型(単精度）として扱われてしまい、引数を渡すとき、メモリ上で
バイト数が一致しないので、数値が変わったように見えているのだと思います。

C
1※ double で宣言して atof で 数値に変換
2    double a_c,b_c,area;
3    a_c = atof(argv[1]);
4    b_c = atof(argv[2]);
5
6    area = sekibun_(&a_c, &b_c);

FORTRAN77 側で DOUBLE の宣言を追加してみてはいかがですか？

FORTRAN77
1      DOUBLE PRECISION FUNCTION SEKIBUN(A, B) 
2      DOUBLE PRECISION X
3      DOUBLE PRECISION I
4      DOUBLE PRECISION AREA 
5      DOUBLE PRECISION DIV
6
7※引数を DOUBLEに宣言
8      DOUBLE PRECISION A
9      DOUBLE PRECISION B