Discussion :

[kis93]

Bonjour,

J'ai un problème à la compilation de deux algorithmes écrits en C, sachant que:

mex -setup
Please choose your compiler for building external interface (MEX) files:

Would you like mex to locate installed compilers [y]/n? y

Select a compiler:
[1] Lcc C version 2.4 in C:\MATLAB7\sys\lcc

[0] None

Compiler: 1

Please verify your choices:

Compiler: Lcc C 2.4
Location: C:\MATLAB7\sys\lcc

Are these correct?([y]/n): y

Try to update options file: C:\Documents and Settings\GuyoM\Application Data\MathWorks\MATLAB\R14\mexopts.bat
From template: C:\MATLAB7\BIN\WIN32\mexopts\lccopts.bat

Done . . .

Apres la ligne de commande :

mex atgdmp_v2.c get_hdr_param.c

Le message d'erreur est :

atgdmp_v2.obj: multiple definition of _main
first definition is in file c:\matlab7\sys\lcc\lib\libcrt0.obj
Specified export _mexFunction is not defined
Missing exports. Aborting

C:\MATLAB7\BIN\WIN32\MEX.PL: Error: Link of 'atgdmp_v2.dll' failed.

??? Error using ==> mex
Unable to complete successfully

Je ne comprends pas car il n'y a qu'un seul 'main' dans les deux programmes utilisés

Je vous mets les codes en C au cas où quelqu'un aurait une idée

atgdmp_v2.c

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
 
#include <stdio.h>
#include <math.h>
main(int argc, char *argv[])
{
    FILE   *atg_fp, *reftrk_asc_fp, *reftrk_des_fp, *eqxl_fp;
    char   filename[100], mission[5];
    char   reftrk_asc_filename[100], reftrk_des_filename[100], eqxl_filename[100];
    int    cycle, pass, begin_rec, end_rec, start_of_data_offset;
    int    i, j, n, status, recsize, nchars, offset, nbytes;
    char   search_buf[200], search_param[200], *ptr;
    unsigned char header[3127], temp_buf[200];
    short  ssha[3127];
    int    whole, days, msecs;
    int    begin_cycle_days, begin_cycle_millisecs;
    double lat[3127], lon[3127], latA[3127], lonA[3127];
    double latD[3127], lonD[3127], eqxl[254], deltat;
    double reflat[3127], reflon[3127];
    double part, time, ti, integral, frac, z;
 
    recsize =  3127;
 
    if(argc != 7) {
        printf("Usage: %s <filename><mission><cycle><pass><begin_rec><end_rec>\n", argv[0]);
	printf("Example: atgdmp_v2 j1atgssha11.020 j1 20 1 1 50\n");
        exit(1);
    }
    strcpy(filename, argv[1]);
    strcpy(mission, argv[2]);
    cycle = atoi(argv[3]);
    pass = atoi(argv[4]);
    begin_rec = atoi(argv[5]);
    end_rec = atoi(argv[6]);
 
 
    if((strcmp(mission, "tp") == 0) && cycle >= 365 && cycle <= 368) {
        printf("Cycle must not include 365, 366, 367, or 368.\n");
        printf("These are drift cycles for which there are no\n");
        printf("along track gridded files available.\n");
        exit(1);
    }
    if((strcmp(mission, "tp") != 0 && strcmp(mission, "j1") != 0)) {
        printf("mission must be one of tp or j1\n");
        printf("where tp is TOPEX/Poseidon and j1 is Jason-1\n");
        exit(1);
    }
    if(begin_rec < 1 || end_rec > 3127) {
        printf("Input Error: start_rec and end_rec must be between 1 and 3127 and\n");
        exit(1);
    }
    if(begin_rec > end_rec) {
        printf("Input Error: start rec is greater than end_rec\n");
	exit(1);
    }
    printf("Dumping file %s for pass %03d, begin_rec = %04d, end_rec = %04d\n",
        filename, pass, begin_rec, end_rec);
    printf("\n");
    if((atg_fp = fopen(filename, "r")) == NULL) {
        printf("Can't open %s\n", filename);
        exit(1);
    }
    /****************/
    /* PRINT HEADER */
    /****************/
    if((nchars = fread(header, 1, 3127, atg_fp)) != 3127) {
        printf("Error reading header characters\n");
	printf("Number of characters read = %d not equal to 3127\n", nchars);
	exit(1);
    }
    /* Print Header */
    ptr = &header[0];
    for(i = 0; i < 15; i++) {
        j = 0;
        while((temp_buf[j] = *ptr++) != '\n') {
	    j++;
	}
	temp_buf[j] = '\0';
        printf("%s\n", temp_buf);
    }
    printf("\n");
 
    /* Get specific header values */
    strcpy(search_buf, "Begin_Cycle_Days");
    if((status = get_hdr_param(header, search_buf, search_param)) != 0) {
        printf("Bad status from get_j1gdr_hdr_param for %s\n", search_buf);
        exit(1);
    }
    begin_cycle_days = atoi(search_param);
    strcpy(search_buf, "Begin_Cycle_Millisecs");
    if((status = get_hdr_param(header, search_buf, search_param)) != 0) {
        printf("Bad status from get_j1gdr_hdr_param for %s\n", search_buf);
        exit(1);
    }
    begin_cycle_millisecs = atoi(search_param);
    strcpy(search_buf, "Deltat");
    if((status = get_hdr_param(header, search_buf, search_param)) != 0) {
        printf("Bad status from get_j1gdr_hdr_param for %s\n", search_buf);
        exit(1);
    }
    deltat = atof(search_param);
    printf("begin_cycle_days      = %d\n", begin_cycle_days);
    printf("begin_cycle_millisecs = %d\n", begin_cycle_millisecs);
    printf("deltat                = %6.4f seconds\n", deltat);
 
    /**************************************************************/
    /* Read reference latitudes and longitudes for ascending pass */
    /**************************************************************/
    strcpy(reftrk_asc_filename, "/data/pgs/pgsa/atg_v2/util/atgdmp/reftrk.ascending");
    if((reftrk_asc_fp = fopen(reftrk_asc_filename, "r")) == NULL) {
        printf("Can't open %s\n", reftrk_asc_filename);
        exit(1);
    }
    for(i = 0; i < 3127; i++) {
        if(fscanf(reftrk_asc_fp, "%lf %lf", &latA[i], &lonA[i]) == EOF) {
	    printf("EOF encountered while reading reference latitudes\n");
            printf("and longitudes for ascending passes from filename %s\n",
            reftrk_asc_filename);
	    exit(1);
	}
    }
    fclose(reftrk_asc_fp);
 
    /***************************************************************/
    /* Read reference latitudes and longitudes for descending pass */
    /***************************************************************/
    strcpy(reftrk_des_filename, "/data/pgs/pgsa/atg_v2/util/atgdmp/reftrk.descending");
    if((reftrk_des_fp = fopen(reftrk_des_filename, "r")) == NULL) {
        printf("Can't open %s\n", reftrk_des_filename);
        exit(1);
    }
    for(i = 0; i < 3127; i++) {
        if(fscanf(reftrk_des_fp, "%lf %lf", &latD[i], &lonD[i]) == EOF) {
	    printf("EOF encountered while reading reference latitudes\n");
            printf("and longitudes for descending passes from filename %s\n",
            reftrk_des_filename);
	    exit(1);
	}
    }
    fclose(reftrk_des_fp);
 
    /************************************/
    /* Read equator crossing longitudes */
    /************************************/
    if((strcmp(mission, "j1") == 0)) {
        strcpy(eqxl_filename, "/data/pgs/pgsa/atg_v2/util/atgdmp/eqxl_tp_then_j1.dat");
    }
    if((strcmp(mission, "tp") == 0)) {
        if(cycle <= 364)
	    strcpy(eqxl_filename, "/data/pgs/pgsa/atg_v2/util/atgdmp/eqxl_tp_then_j1.dat");
	else
	    strcpy(eqxl_filename, "/data/pgs/pgsa/atg_v2/util/atgdmp/eqxl_new_tp.dat");
    }
    if((eqxl_fp = fopen(eqxl_filename, "r")) == NULL) {
        printf("Can't open %s\n", eqxl_filename);
        exit(1);
    }
    for(i = 0; i < 254; i++) {
        if(fscanf(eqxl_fp, "%lf", &eqxl[i]) == EOF) {
	    printf("EOF encountered while reading equator crossing longitudes\n");
            printf("from filename %s\n", eqxl_filename);
	    exit(1);
	}
    }
    fclose(eqxl_fp);
 
    if((pass % 2) == 0) {  /* For even numbered passes */
        for(i = 0; i < 3127; i++) {
	  reflat[i] = latD[i];
	  reflon[i] = lonD[i];
        }
    }
    else {  /* For odd numbered passes */
        for(i = 0; i < 3127; i++) {
	  reflat[i] = latA[i];
	  reflon[i] = lonA[i];
        }
    }
    /**************************/
    /* Print table of results */
    /**************************/
    n = end_rec - begin_rec + 1;
 
    for(i = 0; i < 3127; i++) {
        /* Set longitude between -180 and 180 degrees */
        lon[i] = reflon[i] + eqxl[pass-1];
        if(lon[i] < -180.0) lon[i] += 360.0;
	if(lon[i] >= 180.0) lon[i] -= 360.0;
    }
    printf("Begin_days = %d\n", begin_cycle_days);
    printf("Begin_millisecs = %d\n", begin_cycle_millisecs);
    printf("eqxl file used = %s\n", eqxl_filename);
    printf("For pass %3d, equator crossing longitude = %10.6f deg\n", pass, eqxl[pass-1]);
    printf("\n");
    printf("==================================================\n");
    printf("REC    DAYS     MSECS      LAT        LON     SSHA\n");
    printf("==================================================\n");
 
    offset = 2*3127*pass;
    fseek(atg_fp, offset, SEEK_SET);
    if((nbytes = fread(ssha, 3127, 2, atg_fp)) == 0) {
        printf("Read error or end of file encountered\n");
        printf("file %s\n", filename);
        exit(1);
    }
    for(i = begin_rec; i <= end_rec; i++) {
	time = (double)(3127*(pass - 1)*deltat + (i-1) * deltat); /* time in sec */
	ti = begin_cycle_millisecs + time*1000.0;
	whole = (int)(ti/86400000.0);
	days = begin_cycle_days + whole;
	msecs = (long)ti;
	printf("%4d  %4d   %8d   %8.4f   %8.4f  %5d \n", 
	    i, days, msecs, reflat[i-1], lon[i-1], ssha[i-1]);
    }
    fclose(atg_fp);
    exit(0);
}

get_hdr_param.c

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
 
#include <stdio.h>
#include <string.h>
int get_hdr_param(char *buf, char *search_param, char *hdr_value)
{
    FILE *fp;
    int  i, nchars;
    char *ptr;
    unsigned char temp_buf[200];
 
    /* For example, suppose that the input, search_param, contains
       the string "Product_File_Name".  The actual header record is:
 
       "Product_File_Name = Ja1_GDR_2PaP001_053.NASA"
 
       and we want to extract the value "Ja1_GDR_2PaP001_053.NASA" from it.
 
       buf contains the whole header of the j1gdr file.
       Then buf is searched for the string containing 
       "Product_File_Name", where a pointer is returned pointing to the 
       first character, "P", of "Product_File_Name".
 
       Then the whole header record "Product_File_Name = Ja1_GDR_2PaP001_053.NASA"
       is copied to parameter temp_buf.  From this parameter, temp_buf, the 
       value of the parameter, "Ja1_GDR_2PaP001_053.NASA" is copied to the
       output buffer, hdr_value.
    */
    /* Search for Product_File_Name */
    ptr = strstr(buf, search_param);
    if(ptr == NULL) {
        printf("GET_J1GDR_HDR_PARAM: No match was found for search of %s\n", search_param);
	return -1;
    }
    else {
	/* Now copy header record to temp_buf */
	i = 0;
	while((temp_buf[i] = *ptr++) != ';')
            i++;
        temp_buf[i] = '\0';
 
	/* Now extract the value out of temp_param */
	ptr = strstr(temp_buf, "=");
        ptr++; /* point to next character after equal sign */
 
        /* Point to first non_blank char after equal sign */
 
        while(*ptr == ' ') {
            ptr++;
	    i++;
        }
	i = 0;
        while((hdr_value[i] = *ptr++) != '\0') {
            i++;
        }
        hdr_value[i] = '\0';   
    }
    return 0;
}

[kis93]

Finalement, j'ai ré écrit les algorithmes sous matlab. Charité bien ordonnée.....Pour la peine, je vais travailler demain !

Si quelqu'un a tout de même une idée concernant le problème de compilation de mes deux codes en C je suis toujours preneur.

[paradize3]

Finalement, j'ai ré écrit les algorithmes sous matlab. Charité bien ordonnée.....Pour la peine, je vais travailler demain !

Si quelqu'un a tout de même une idée concernant le problème de compilation de mes deux codes en C je suis toujours preneur.

Le point d'entree d'un fichier mex est
extern "C" __declspec(dllexport) void mexFunction( int nlhs,mxArray *plhs[],int nrhs, mxArray *prhs[])
et non main().

plus d'info dans la doc sur les mex files!


salutations,

Greg

[kis93]

Merci Greg

Par contre, j'ai du mal comprendre quelquechose mais je ne vois pas pourquoi je dois changer la ligne de commande "main". je voulais juste pouvoir compiler les deux algorithme ecrit en C avec matlab et non utiliser des fonctions C dans un programme matlab. Je ne sais pas si c'est tres clair je m'y perd aussi