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
|
#include <stdio.h>
#include <mlvalues.h>
#include <bigarray.h>
extern void filltab_(void);
extern void printtab_(float * data, int * dimx, int * dimy);
extern float ftab_[];
#define DIMX 6
#define DIMY 8
double ctab[DIMX][DIMY];
void filltab(void)
{
int x, y;
for (x = 0; x < DIMX; x++)
for (y = 0; y < DIMY; y++)
ctab[x][y] = x * 100 + y;
}
void printtab(double tab[DIMX][DIMY])
{
int x, y;
for (x = 0; x < DIMX; x++) {
printf("%3d", x);
for (y = 0; y < DIMY; y++)
printf(" %6.1f", tab[x][y]);
printf("\n");
}
}
value c_filltab(value unit)
{
filltab();
return alloc_bigarray_dims(BIGARRAY_FLOAT64 | BIGARRAY_C_LAYOUT,
2, ctab, DIMX, DIMY);
}
value c_printtab(value ba)
{
printtab(Data_bigarray_val(ba));
return Val_unit;
}
value fortran_filltab(value unit)
{
filltab_();
return alloc_bigarray_dims(BIGARRAY_FLOAT32 | BIGARRAY_FORTRAN_LAYOUT,
2, ftab_, 8, 6);
}
value fortran_printtab(value ba)
{
int dimx = Bigarray_val(ba)->dim[0];
int dimy = Bigarray_val(ba)->dim[1];
printtab_(Data_bigarray_val(ba), &dimx, &dimy);
return Val_unit;
}
|
