test

Sun Studio 12：Fortran 编程指南

11.5.2 返回 COMPLEX 数据

COMPLEX 数据的互操作性情况在 SPARC 32 位和 64 位实现之间有所不同。

11.5.2.1 32 位平台

在 32 位平台上，返回 COMPLEX 或 DOUBLE COMPLEX 的 Fortran 函数等同于具有一个指向内存中返回值的附加第一参数的 C 函数。Fortran 函数及其相应的 C 函数的一般样式如下：

Fortran 函数 

C 函数 

 

COMPLEX FUNCTION CF( a1,a2,...,an)


cf_ (return, a1, a2, ..., an)
 struct { float r, i; } *return

表 11–13 返回 COMPLEX 数据的函数（32 位 SPARC）

Fortran 调用 C 

C 调用 Fortran 

 

COMPLEX U, V, RETCPX
EXTERNAL RETCPX
U = ( 7.0, -8.0)
V = RETCPX(U)
...

------------------------------

struct complex { float r, i; };
void retcpx_( temp, w )
struct complex *temp, *w;
{
   temp->r = w->r + 1.0;
   temp->i = w->i + 1.0;
   return;
}


struct complex { float r, i; };
struct complex c1, c2;<
struct complex *u=&c1, *v=&c2;
extern retfpx_();
u -> r = 7.0;
u -> i = -8.0;
retfpx_( v, u );
...

------------------------------

COMPLEX FUNCTION RETFPX(Z)
  COMPLEX Z
  RETFPX = Z + (1.0, 1.0)
  RETURN
END

11.5.2.2 64 位 SPARC 平台

在 64 位 SPARC 环境下，在浮点寄存器中返回 COMPLEX 值： 在 %f0 和 %f1 中返回 COMPLEX 和 DOUBLE COMPLEX，在 %f0、%f1、%f2 和 %f3 中返回 COMPLEX*32。对于 64 位 SPARC，返回结构（其字段均为浮点型）的 C 函数将在浮点寄存器中返回该结构，但条件是最多需要 4 个这样的寄存器进行此操作。在 64 位 SPARC 平台上，Fortran 函数及其相应的 C 函数的一般样式如下：

Fortran 函数 

C 函数 

COMPLEX FUNCTION CF(a1, a2, ..., an)

struct {float r,i;} cf_ (a1, a2, ..., an)

表 11–14 返回 COMPLEX 数据的函数（64 位 SPARC）

Fortran 调用 C 

 

COMPLEX U, V, RETCPX
EXTERNAL RETCPX
U = ( 7.0, -8.0)
V = RETCPX(U)
...

---------------------------------------------------------

struct complex {float r, i; };
struct complex retcpx_(struct complex *w )
{
    struct complex temp;
    temp.r = w->r + 1.0;
    temp.ii = w->i + 1.0;
    return (temp);
}

C 调用 Fortran 

 

struct complex { float r, i; };
struct complex c1, c2;
struct complex *u=&c1;
extern struct complex retfpx_(struct complex *);
  u -> r = 7.0;
  u -> i = -8.0;
  retfpx_( u );
...

---------------------------------------------------------

COMPLEX FUNCTION RETFPX(Z)
  COMPLEX Z
  RETFPX = Z + (1.0, 1.0)
  RETURN
END