blob: bf076ce4a40650c81394bfed85230258d6a9833e (
plain)
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
|
double fabs (double);
/* interops with myftype_1 */
typedef struct {
int m, n;
float r;
} myctype_t;
/* interops with particle in f90 */
typedef struct particle
{
double x; /* x position */
double vx; /* velocity in x direction */
double y; /* y position */
double vy; /* velocity in y direction */
double z; /* z position */
double vz; /* velocity in z direction */
double m; /* mass */
}particle_t;
extern void abort(void);
void types_test(particle_t *my_particles, int num_particles);
/* declared in the fortran module bind_c_dts */
extern myctype_t myDerived;
int main(int argc, char **argv)
{
particle_t my_particles[100];
/* the fortran code will modify the middle particle */
my_particles[49].x = 1.0;
my_particles[49].vx = 1.0;
my_particles[49].y = 1.0;
my_particles[49].vy = 1.0;
my_particles[49].z = 1.0;
my_particles[49].vz = 1.0;
my_particles[49].m = 1.0;
myDerived.m = 1;
myDerived.n = 2;
myDerived.r = 3.0;
types_test(&(my_particles[0]), 100);
if(fabs(my_particles[49].x - 1.2) > 0.00000000)
abort();
if(fabs(my_particles[49].vx - 1.2) > 0.00000000)
abort();
if(fabs(my_particles[49].y - 1.2) > 0.00000000)
abort();
if(fabs(my_particles[49].vy - 1.2) > 0.00000000)
abort();
if(fabs(my_particles[49].z - 1.2) > 0.00000000)
abort();
if(fabs(my_particles[49].vz - 1.2) > 0.00000000)
abort();
if(fabs(my_particles[49].m - 1.2) > 0.00000000)
abort();
if(myDerived.m != 2)
abort();
if(myDerived.n != 3)
abort();
if(fabs(myDerived.r - 4.0) > 0.00000000)
abort();
return 0;
}/* end main() */
|
