1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
|
/* Passing from fortran to C by value, using %VAL. */
#include <inttypes.h>
/* We used to #include <complex.h>, but this fails for some platforms
(like cygwin) who don't have it yet. */
#define complex __complex__
#define _Complex_I (1.0iF)
extern void f_to_f__ (float*, float, float*, float**);
extern void f_to_f8__ (double*, double, double*, double**);
extern void i_to_i__ (int*, int, int*, int**);
extern void i_to_i8__ (int64_t*, int64_t, int64_t*, int64_t**);
extern void c_to_c__ (complex float*, complex float, complex float*, complex float**);
extern void c_to_c8__ (complex double*, complex double, complex double*, complex double**);
extern void abort (void);
void
f_to_f__(float *retval, float a1, float *a2, float **a3)
{
if ( a1 != *a2 ) abort();
if ( a1 != **a3 ) abort();
a1 = 0.0;
*retval = *a2 * 2.0;
return;
}
void
f_to_f8__(double *retval, double a1, double *a2, double **a3)
{
if ( a1 != *a2 ) abort();
if ( a1 != **a3 ) abort();
a1 = 0.0;
*retval = *a2 * 2.0;
return;
}
void
i_to_i__(int *retval, int i1, int *i2, int **i3)
{
if ( i1 != *i2 ) abort();
if ( i1 != **i3 ) abort();
i1 = 0;
*retval = *i2 * 3;
return;
}
void
i_to_i8__(int64_t *retval, int64_t i1, int64_t *i2, int64_t **i3)
{
if ( i1 != *i2 ) abort();
if ( i1 != **i3 ) abort();
i1 = 0;
*retval = *i2 * 3;
return;
}
void
c_to_c__(complex float *retval, complex float c1, complex float *c2, complex float **c3)
{
if ( c1 != *c2 ) abort();
if ( c1 != *(*c3) ) abort();
c1 = 0.0 + 0.0 * _Complex_I;
*retval = (*c2) * 4.0;
return;
}
void
c_to_c8__(complex double *retval, complex double c1, complex double *c2, complex double **c3)
{
if ( c1 != *c2 ) abort();
if ( c1 != *(*c3) ) abort();
c1 = 0.0 + 0.0 * _Complex_I;;
*retval = (*c2) * 4.0;
return;
}
|
