Discussion :

[takout]

Bonjour à tous,

Après avoir écrit un mex-files qui fait la simple copie d'un tableau de double, je souhaiterai le rendre générique. C'est à dire quel que soit le type en entrée de mon tableau (uint16, int16,...) mon programme fonctionne et me retourne un résultat.
Pouvez-vous m'aider svp ?

voici mon mex-file :

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
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
#include "mex.h"
#include <stdio.h>
#include <stdlib.h>
 
#define res(i1,j1) res[(i1)+row*(j1)]
#define x(i1, j1) x[(i1)+row*(j1)]
 
 
    /*#######################************ Programme C hysteresis***********************####################*/
void copie(double *x,double *res, int row, int col)
{  
        int i1,j1;
 
       for(i1=0;i1<row;i1++)
       {
         for(j1=0;j1<col;j1++)
         {
         res(i1,j1)=x(i1,j1);
         }
       }
 
}
 
 
/* nrhs nombre de parametre d'entre
  nlhs nombre de parametre de sortie
  prhs tableau de pointeur d'un element en entre
  plhs tableau de pointeur de sortie
*/
void mexFunction(int nlhs,mxArray *plhs[],int nrhs,const mxArray *prhs[])
{
 
   double *a,*c; /* a  est le tableau en entre et c le tableau en sortie Matlab*/
    int row,col;
 
 
 
        /* Get inputs */ 
    row = mxGetM(prhs[0]); 
    col = mxGetN(prhs[0]);
 
   a = mxGetPr(prhs[0]); /* retourne la matrice pointée par prhs[0]*/
 
  /* Create a matrix for the return argument. */
  plhs[0] = mxCreateDoubleMatrix(row, col, mxREAL);
  c = mxGetPr(plhs[0]);
 
    /* main computation */
     copie(a,c,row,col);
 
}

[Jerome Briot]

Le plus simple consiste à tester la classe du tableau d'entrée avec mxIsClass et à écrire une sous-fonction "copie#class#" par classe supportée.

Tu utiliseras mxCreateNumericArray pour créer le tableau de sortie de la classe voulue.

[takout]

Ah ok en fait j'avais l'idée de faire une seul fonction et de faire un cast du tableau en entrée est ce possible (si le tableau par exemple n'est pas de type double) ?

[Jerome Briot]

C'est certainement possible, mais je commencerais quand même par la méthode que je préconise.

[gpcbitnik38]

Ah ok en fait j'avais l'idée de faire une seul fonction et de faire un cast du tableau en entrée est ce possible (si le tableau par exemple n'est pas de type double) ?

j'ai pas bien compris ce que tu veux faire,
si on te donne un t'ableau d'entier par exemple tu veux ressortir un tableau d'entier ? ou un tableau de double?

en utilisant a = mxGetPr(prhs[0]); tu retournera forcement un tableau de double
pour retouner un t'ableau d'un autre type tu doit passer par la fonction mxSetData

ensuite pour passer d'un type a un autre tu peux utiliser : reinterpret_cast < new_type > ( expression )

une fasson simple sans te casser la tete serai de caster toute ton entré en double de faire tes opération sur des double et a la fin de faire le cast dans l'autre sens
ça nécéssiterai un petit switch en entré pour savoir quel cast faire en sortie

[takout]

Ben en fait ce que je voulais faire c'est la dernière idée que tu as eu.
Écrire une fonction qui prend en entrée des doubles. Si lors de l'utilisation de cette fonction le tableau en entrée est un int par exemple j'opère un cast.

[Jerome Briot]

Cette solution a pour principal désavantage de consommer et d'avoir à gérer inutilement de la mémoire.

[takout]

Oui c'est vrai Dut à raison.

[gpcbitnik38]

Dut je parlais d'une solution simple à mettre en place pas optimale ^^

sinon je viens de me souvenir d'un truc
si tu veux juste faire une copie de la matrice donnée par ton .m pourquoi tu ne duplique pas l'entré avec mxDuplicateArray?

Code C :

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

1
2
3
4
5
6
7
8
9
10
#include "mex.h"
 
 
void mexFunction(int nlhs,mxArray *plhs[],int nrhs,const mxArray *prhs[])
{
 
 
  plhs[0] = mxDuplicateArray(prhs[0]);
 
}

[takout]

Je ne l'ai pas utilisé car c'étais juste un exemple pour illustrer mon problème.