Discussion :

[Giansolo]

Bonjour à toutes et à tous,

Je suis en train de me mettre aux fichiers MEX pour transformer mes fonctions matlab en fonction C pour gagner de la vitesse.
Ca fait un bail que j'ai fait du C/C++ et j'avoue que je redécouvre ce langage.
J'ai écrit un petit programme pour essayer de faire fonctionner simplement les mécanismes du MEX, et il marche... une fois? il fait bien un carré, mais si je l'appelle 2 fois matlab crash (veuillez redémarrer matlab etc...).
J'ai donc du me planter quelque part dans la création d'espace mémoire, mais je ne sais pas trop ou je faute, pouvez vous m'aider?

Voici le code:

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
 
#include "mex.h"
 
void mexFunction(int nout, 
                mxArray *pout[], 
                int nin, 
                const mxArray *pin[] ) 
{
    double *x;
    double *y;
 
    *x = mxGetScalar(pin[0]);
 
    /* Preparation de la sortie de la fonction */
    pout[0] = mxCreateDoubleMatrix(1,1,mxREAL);
    y = mxGetPr(pout[0]);
 
    /* Calcul */
    y[0] = x[0] * x[0];
 
    /* Sortie de la fonction */
    mxSetPr(pout[0],y);
 
}

Je pense que l'erreur est dans la ligne: pout[0] = mxCreateDoubleMatrix(1,1,mxREAL); ne devrais-je pas utiliser un type mxDOUBLE?

Et si j'ai bien compris x et y sont des pointeurs vers des tableaux, ce qui sous entend que je peux les modifier directement en utilisant une fonction de type (calc(x,y))?

Merci par avance pour les informations que vous pourrez me communiquer...

Gian

[Jerome Briot]

Ton code ne contient pas vraiment d'erreur.

Il peut se simplifier comme ceci :

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
#include "mex.h"
 
void mexFunction(int nout, 
                mxArray *pout[], 
                int nin, 
                const mxArray *pin[] ) 
{
    double x;
    double *y;
 
    x = mxGetScalar(pin[0]);
 
    pout[0] = mxCreateDoubleMatrix(1,1,mxREAL);
    y = mxGetPr(pout[0]);
 
    y[0] = x * x;
 
}

ou encore :

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
#include "mex.h"
 
void mexFunction(int nout, 
                mxArray *pout[], 
                int nin, 
                const mxArray *pin[] ) 
{
    double x;
    double *y;
 
    x = mxGetScalar(pin[0]);
 
    /* Preparation de la sortie de la fonction */
    pout[0] = mxCreateDoubleScalar(x*x);
 
}

L'argument MXREAL de MXCREATEDOUBLEMATRIX sert juste à spécifier si le tableau contient des valeurs complexes ou non.

[Giansolo]

Merci pour ton retour,
Est-ce que par hasard, il y aurait une explication au fait que ce code plante? Je peux l'executer un fois par session de matlab, la 2ème fois matlab crash. Est-ce du à un problème de version de compilateur ou autre?

Merci par avance,

Gian

[Jerome Briot]

Je pense que l'erreur vient du fait que tu n'initialises pas le pointeur au début de ton code

Il est plus sûr de faire :

Code C :

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

1
2
3
4
5
double *pr, x;
 
pr = &x; 
 
*pr = mxGetScalar(pin[0]);

Sans ces lignes, le code ne fonctionne pas avec les compilateurs Lcc ou Open Watcom mais par contre, il fonctionne correctement avec le compilateur Microsoft Visual C++ Express. Serait-ce dû à une spécificité du compilateur Microsoft ?

De toute façon, pourquoi faire compliqué quand on peut faire simple :

Code C :

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

1
2
3
double x;
 
x = mxGetScalar(pin[0]);

[Giansolo]

ok, merci pour ton retour, je vais tester ca, comme je ne maitrise plus le C faut que je reregarde un peu tout ca...