67 lines
1.6 KiB
C
67 lines
1.6 KiB
C
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() */
|