IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

MATLAB Discussion :

remplissage matrice sparse


Sujet :

MATLAB

  1. #1
    Membre averti
    Homme Profil pro
    Étudiant
    Inscrit en
    Décembre 2011
    Messages
    265
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Rhône (Rhône Alpes)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Décembre 2011
    Messages : 265
    Points : 352
    Points
    352
    Par défaut remplissage matrice sparse
    Bonjour,
    je cherche a créer une matrice sparse depuis un .mex probleme la matrice retournée est vide alors qu'elle devrai contenir size[0] element non nul.

    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
    52
    53
    54
    55
    56
    57
    58
     
    void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
    {
        double*ind,*t1,*t2,*size,*L,*sr;
        mxArray *inda,*t1a,*t2a,*sizea,*La,*sparse;
        mwIndex *irs,*jcs;
        inda=mxDuplicateArray(prhs[0]);
        t1a=mxDuplicateArray(prhs[1]);
        t2a=mxDuplicateArray(prhs[2]);
        La=mxDuplicateArray(prhs[3]);
        sizea=mxDuplicateArray(prhs[4]);
        ind=mxGetPr(inda);
        t1=mxGetPr(t1a);
        t2=mxGetPr(t2a);
        L=mxGetPr(La);
        size=mxGetPr(sizea);
        mxComplexity cmplx=mxREAL;
     
        sparse=plhs[0]=mxCreateSparse(size[0],size[1]*size[2],size[0],cmplx);
        sr  = mxGetPr(sparse);
        irs = mxGetIr(sparse);
        jcs = mxGetJc(sparse);
        int k;
        k=0;
        for(int j=0;j<size[0];j++)
        {
            jcs[j] = k;
            if(1-t2[j]>0)
            {
                printf("\n 1-t2[%d]=%f",j,1-t2[j]);
                sr[k]=1-t2[j];
                irs[k]=ind[j];
                k++;
            }
            if(t1[j]>0)
            {
                printf("\n t1[%d]=%f",j,t1[j]);
                sr[k]=t1[j];
                irs[k]=ind[j]+size[2]+1;
                k++;
            }
            if(t2[j]-t1[j]>0)
            {
                printf("\n t2[%d]-t1[%d]=%f",j,j,t2[j]-t1[j]);
                sr[k]=t2[j]-t1[j];
                irs[k]=ind[j]+L[j];
                k++;
            }
     
     
        }
        jcs[(int)size[0]]=k;
     
        mexCallMATLAB(0,NULL,0,NULL,"pause");
     
     
     
    }

    les printf m'indique que le code tourne bien, l'algo rempli bien le champ sr jcs et irs
    je me suis inspiré de l'exemple fournis par matlab (fulltosparse.c) et je ne vois pas d'ou vien le probleme

  2. #2
    Membre averti
    Homme Profil pro
    Étudiant
    Inscrit en
    Décembre 2011
    Messages
    265
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Rhône (Rhône Alpes)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Décembre 2011
    Messages : 265
    Points : 352
    Points
    352
    Par défaut
    donc j'ai resolu un "bug" mais la matrice n'est toujours pas alloué
    étant en 64 bits il faut redefinir les type de base mwIndex et mwSize en size_t ce que j'ai donc fait mais cela ne fonctionne toujours pas
    la matrice sparce alloué est parfaitement utilisable deriere elle n'est juste pas initialisé correctement...

    je joint un .zip avec de quoi tester le mex file
    le .zip comprend main.m qui charge les donnée et une fct interpo
    cette fonction utilise fichier mex ( le but etant de remplacé la commande au dessus de l'appel pa le fichier mex)
    Fichiers attachés Fichiers attachés

Discussions similaires

  1. Compréhension matrice sparse
    Par Frobenius dans le forum ALM
    Réponses: 0
    Dernier message: 10/04/2012, 23h44
  2. Multiplication matrices sparse - lent
    Par Mathieu999 dans le forum MATLAB
    Réponses: 8
    Dernier message: 25/01/2012, 23h45
  3. Remplissage matrice, toutes les solutions possibles
    Par corky813 dans le forum Algorithmes et structures de données
    Réponses: 18
    Dernier message: 16/06/2011, 12h25
  4. Matrices sparse avec Boost
    Par circe dans le forum Boost
    Réponses: 1
    Dernier message: 12/12/2009, 13h14
  5. Matrices sparse et fonction spconvert
    Par alain10025 dans le forum MATLAB
    Réponses: 8
    Dernier message: 17/11/2007, 00h19

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo