Discussion :

[ol9245]

Bonjour,

Je commence tout juste à interfacer MATLAB avec du codes C.
Pour l'instant, j'essaye de comprendre un peu l'API.

Je sais passer et récupérer des double. Voici le code de test qui me fait ça :

Code C :

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
#include <mex.h>
 
void mexFunction(int n_outputs, mxArray *p_outputs[], int n_inputs, const mxArray *p_inputs[]) {
 
    double a = mxGetScalar(p_inputs[0]);
    mexPrintf("a = %g\n", a) ;
 
//associate outputs 
    const mwSize di[2] = {1, 8} ;
    mxArray *c_out_m = p_outputs[0] = mxCreateNumericArray(2, di, mxDOUBLE_CLASS, mxREAL);
 
//associate pointers
    double *c = mxGetPr(c_out_m);
 
// do something
    for (int u=0; u!=8; u++) {
        c[u] = a ;
    }
 
    return;
}

Le joli résultat est celui-ci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
>> mex entrees_sorties.cpp
 
Warning: You are using gcc version "4.6.3-1ubuntu5)".  The version
         currently supported with MEX is "4.3.4".
         For a list of currently supported compilers see: 
         http://www.mathworks.com/support/compilers/current_release/
 
>> A = entrees_sorties(12.5)
a = 12.5
 
A =
 
   12.5000   12.5000   12.5000   12.5000   12.5000   12.5000   12.5000   12.5000

Ce que je n'arrive pas à faire : exactement la même chose mais sur des entiers (passer un tableau d'entiers en arguments, travailler en C sur ces entiers, retourner le résultat en entier, sans jamais dupliquer inutilement les données en double.

Merci de votre aide

[Invité]

Bonjour,

Qu'as-tu essayé pour les entiers ?
Qu'entends-tu par

sans jamais dupliquer inutilement les données en double.

?

[ol9245]

Merci pour ton passage.

Pour les entiers, je n'ai rien trouvé.
Je suis bloqué par l'instruction suivante (ligne 13 de mon premer post):

Code :

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

double *c = mxGetPr(c_out_m);

La fonction renvoie obligatoirement un double*. Je n'ai aps trouvé dans la doc d'équivalent qui renvoie un int* ou un "ce que je veux"*

J'ai donc fait ce petit script dans l'espoir que quelqu'un puisse me le convertir pour des entiers.

Je répète que j'en suis à mes tout-débuts dans la prise en main de cette API. Donc je ne crois pas que la solution soit ni loin ni compliquée. Simplement je ne l'ai pas vue dans la doc ni dans mes recherches.

[Invité]

doc mxGetData retournant un void * que tu peux caster par la suite.

[ol9245]

ça n'a pas marché

Code C :

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
#include <mex.h>
 
void mexFunction(int n_outputs, mxArray *p_outputs[], int n_inputs, const mxArray *p_inputs[]) {
 
    int a = mxGetScalar(p_inputs[0]);
    mexPrintf("a = %g\n", a) ;
 
//associate outputs 
    const mwSize di[2] = {1, 8} ;
    mxArray *c_out_m = p_outputs[0] = mxCreateNumericArray(2, di, mxDOUBLE_CLASS, mxREAL);
 
//associate pointers
    int *c = (int*) mxGetPr(c_out_m);
 
// do something
    for (int u=0; u!=8; u++) {
        c[u] = a ;
    }
 
    return;
}

Résultat :

Code :

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

1
2
3
4
5
6
7
8
9
10
>> a = int16(12) ; A = entrees_sorties(a)
a = 6.89886e-310
 
A =
 
  1.0e-312 *
 
    0.2546    0.2546    0.2546    0.2546         0         0         0         0
 
>>

[Jerome Briot]

Mais si :

Code C :

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
#include "mex.h"
 
void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]
        )
 
{
 
    unsigned char *pr, *pl;
    mwSize i, dims[2] = {0,0};
 
    if(!mxIsUint8)
        mexErrMsgTxt("Array must of class uint8");
 
    pr = (unsigned char *)mxGetData(prhs[0]);
 
    dims[0] = mxGetM(prhs[0]);
    dims[1] = mxGetN(prhs[0]);
 
    plhs[0] = mxCreateNumericArray(2, dims, mxUINT8_CLASS, mxREAL);
 
    pl = (unsigned char *)mxGetData(plhs[0]);
 
    for(i=0;i<dims[1]*dims[0];i++)
        pl[i] = pr[i]*2;
 
    mxSetData(plhs[0],pl);
 
 
}

Code :

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

1
2
3
4
5
6
7
8
9
10
11
>> A = randi(127,1,10,'uint8')
 
A =
 
  104  116   17  116   81   13   36   70  122  123
 
>> B = test(A)
 
B =
 
  208  232   34  232  162   26   72  140  244  246

[ol9245]

merci++ Dut

[Invité]

ça n'a pas marché

Code C :

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

double mxGetScalar(const mxArray *pm);

Code C :

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

double *mxGetPr(const mxArray *pm);

À part ça il est où mon mxGetData ?
Ensuite int16 est plutôt représenté par un short (enfin après cela peut dépendre du compilateur...).

@Dut : il me semble que mxSetData(plhs[0],pl); est inutile vu que l'on fait pl = (unsigned char *)mxGetData(plhs[0]); non ?
[EDIT] if(!mxIsUint8), il manque l'entrée prhs[0]

[Jerome Briot]

merci++ Dut

Winjerome a donné la solution, j'ai juste mis l’exemple

Pour plus d'informations sur les différentes classes dans les fichiers MEX, je conseille d'analyser le code C contenu dans le fichier explore.c dans le sous-dossier extern/examples/mex du dossier d'installation de MATLAB.

@Dut : il me semble que mxSetData(plhs[0],pl); est inutile vu que l'on fait pl = (unsigned char *)mxGetData(plhs[0]); non ?

Oui... je l'avais laissé passque je m'étais embrouillé entre pr et pl en écrivant le code vite fait

[Jerome Briot]

Ensuite int16 est plutôt représenté par un short (enfin après cela peut dépendre du compilateur...).

=> FAQ C : Les types et les variables

[ol9245]

Code C :

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

double mxGetScalar(const mxArray *pm);

Code C :

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

double *mxGetPr(const mxArray *pm);

À part ça il est où mon mxGetData ?
Ensuite int16 est plutôt représenté par un short (enfin après cela peut dépendre du compilateur...).

@Dut : il me semble que mxSetData(plhs[0],pl); est inutile vu que l'on fait pl = (unsigned char *)mxGetData(plhs[0]); non ?
[EDIT] if(!mxIsUint8), il manque l'entrée prhs[0]

Je comprends pas tout ce que tu dis

Là je débute à partir de zéro avec cette API donc le mieux pour moi est d'avoir un exemple simple, qui marche, que je comprend, et que je peux réutiliser pour mes besoins du moment. Après, je pourrai capitaliser peu à peu sur ces acquis. Pour l'instant, je ne sais pas la différence entre mxGetData, mxGetScalar etc. et ma manière d'apprendre ne consiste pas à se palucher des dizaines de pages de doc avant de cracher la première ligne de code : je cherche à récupérer des petits bouts de code qui marchent et à progresser step-by-step à partir de ça en fonction des besoins.

[Invité]

Je comprends pas tout ce que tu dis

Ce serait constructif de me dire quelles parties.

Se plonger dans la documentation est pourtant le meilleur moyen d'y arriver.
Et si tu veux des exemples, outre celui déjà fourni par Dut tu as justement ceux de la documentation. Rien que celle de mxGetData en contient 6 (dont explore.c mentionnée par Dut).

[ol9245]

Ce serait constructif de me dire quelles parties.

Se plonger dans la documentation est pourtant le meilleur moyen d'y arriver.
Et si tu veux des exemples, outre celui déjà fourni par Dut tu as justement ceux de la documentation. Rien que celle de mxGetData en contient 6 (dont explore.c mentionnée par Dut).

Bien sur, je lis les docs ! mais pas en linéaire en commençant par la page 1. Je ne crois pas avoir rien inventé. Tous ceux qui apprennent en travaillant doivent faire +/- la même chose. on a des échéances et des impératifs et au milieu de ça on cherche à en apprendre le plus possible sans trop impacter la productivité.

[Invité]

Lire la documentation en diagonale, qui ne l'a pas fait.

Pour te dire à quel point cela peut être rapide : tu va dans la documentation générale des MEX-Files (Google "MATLAB MEX-Files" saura t'y conduire rapidement ) puis dans la partie C/C++ Examples tu peux facilement voir Handle 8-, 16-, and 32-Bit Data où tu peux voir les deux premiers paragraphes :

You can create and manipulate signed and unsigned 8-, 16-, and 32-bit data from within your MEX-files. The MATLAB API provides a set of functions that support these data types. The API function mxCreateNumericArray constructs an unpopulated N-dimensional numeric array with a specified data size. Refer to the entry for mxClassID in the online reference pages for a discussion of how the MATLAB API represents these data types.

Once you have created an unpopulated MATLAB array of a specified data type, you can access the data using mxGetData and mxGetImagData. These two functions return pointers to the real and imaginary data. You can perform arithmetic on data of 8-, 16-, or 32-bit precision in MEX-files and return the result to MATLAB, which recognizes the correct data class.

avec un petit fichier exemple à consulter.

Voilà en même pas 1 minute, tu as tout sous la main

[Jerome Briot]

Ce qu'il faut bien comprendre, c'est que MATLAB ne connaît que les mxArray.

Donc en entrée et en sortie d'un fichier MEX, il n'y a que des (pointeurs sur) mxArray

Ensuite, pour y accéder à partir du langage C, on utilise des fonctions toutes faites (comme mxGetPr ou MxGetM par exemple) qui vont "piocher" directement dans la structure du mxArray les informations utiles.

[ol9245]

merci infiniment à vous deux pour l'aide.
c'est vraiment ce dont j'avais besoin.

et merci à winjerome pour les liens sur les pages exactes dont j'avais besoin (et que je n'avais pas trouvées tout seul, shame on me)

je mets "résolu". Je reviendrai poster éventuellement le code final pour documenter le fil en vue de futurs internautes qui tomberaient dessus.