Discussion :

[S4sha]

J'ai un petit souci dans un MEX ou je remplis une matrice creuse. Voici le code, en C

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
 
double *mxRES = NULL ; // sortie Matlab
mwIndex *ir = NULL ;
mwIndex *jc = NULL ;
mwIndex entry_counter = 0 ;  // position dans la matrice
 
[...]
 
// creation d'une matrice creuse lmax x kmax
plhs[0] = mxCreateSparse(lmax,kmax,Nelts, mxREAL) ;
 
mxRES = mxGetPr(plhs[0]) ;
ir = mxGetIr(plhs[0]) ;
jc = mxGetJc(plhs[0]) ;
 
for (k=0 ; k<kmax ; k++){
 
    jc[(mwIndex)(k)] = entry_counter ;   // position du 1er élément de la colonne k
 
    for (l = 0 ; l<lmax; l++){
 
        j = some_function( l,k ) ;
 
        // 'some_function' est une fonction qui renvoit un entier entre  0 et lmax-1, de telle sorte que j reste dans le bon intervalle
 
        mxRES[entry_counter] = f(l) ;   // valeur qui depend de l
        ir[entry_counter] = (mwIndex)(j) ;   //  position dans la matrice du j-eme element de la colonne k.
        entry_counter++ ;
 
    }
}
 
jc[ (mwIndex)(kmax) ] = entry_counter ;

Mon problème est que la fonction 'some_function' n'est pas croissante : pour un k donnée, il se peut que j = some_function( l,k ) soit plus petit que j' = some_function( l-1,k ). Du coup la matrice H (sur laquelle pointe prhs[0], la sortie) a des propriété étrange, comme [H'*H]' =/= H'*H

Je ne peux pas rendre 'some_function' croissante.

Est-il possible de remplir la matrice dans un ordre de ligne non croissant ?

Merci beaucoup

PS j'ai posté le même message sur Mathwork

EDIT - je crois que ma question est mal formulée. Reformulation : si, à k fixé (indice de colonne), les indices de lignes j ne sont pas croissants, est-il quand même possible de remplir une matrice creuse (avec une autre fonction peut-être...) ?

EDIT2 - correction dans le code :

Code c :

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

mxRES[entry_counter] = f(l) ;   // valeur qui depend de l

[Jerome Briot]

Pourquoi ne pas faire simplement ceci ?

• Stocker les "j" dans un tableau
• Trier le tableau des "j"
• Copier le tableau trié dans le tableau ir

Grosso modo :

Code c :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
for (k=0; k<kmax; k++){
 
    for (l=0; l<lmax; l++)
        j[l] = some_function(l, k);
 
    /* Tri de j ici */
 
    for (l=0; l<lmax; l++) {
        ir[entry_counter] = (mwIndex)(j[l]);
        entry_counter++;
    }
 
}

Ou alors j'ai mal compris ton problème.

[S4sha]

Salut Dut,

D'accord avec toi là dessus : c'est un problème de tri. Cependant il faut garder en mémoire la positions des j avant triage (le rang de chaque j). Sais-tu quelle fonction C de base (dans algorithm.h par exemple) fait ce travail ?

edit -
je vois en quoi mon problème n'était pas clair : les valeurs que je veux mettre dans la matrice dépendent de l. Je vais corriger ça

[Invité]

Bonjour,

Tu peux utiliser une structure contenant à la fois tes valeurs à trier et tes indices.
Tu tries alors avec qsort un tableau de cette structure en utilisant une fonction qui compare tes valeurs.

[Jerome Briot]

Tu tries alors avec qsort un tableau de cette structure en utilisant une fonction qui compare tes valeurs.

Effectivement, voir ici : Du bon usage de qsort()

[S4sha]

Merci je vais tenter qsort

j'ai implémenté mon propre algo de triage mais il a l'air de ramer... Je ne suis pas une brute en C !

[S4sha]

Bon en mettant les indices de lignes j avec les valeurs f(l) dans une structure, et en triant cette structure par rapport a j, ça a l'air de marcher.

j'ai créé une structure :

Code c :

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

1
2
3
4
5
struct cone_value_idx_type
{
    double value;
    mwIndex vol_index;
} ;

j'ai créé un vecteur du même type que cette structure, avec un nombre d’éléments égal au nombre d'indice j de ligne

Code :

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

std::vector<cone_value_idx_type> cone(lmax);

je remplis le vecteur de structure

Code c :

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

1
2
3
4
5
for (l = 0 ; l<lmax; l++){
    j = some_funtion(l,k)
    cone[n_cone].vol_index = j ;
    cone[n_cone].value = f(l) ;
}

(j'utilise un vecteur car la taille lmax varie tout le temps)

je trie

Code :

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

std::sort(cone.begin(), cone.end(), compare_function);

('compare_function' compare la partie 'value' de la structure)
et je remplis

Code c :

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

1
2
3
4
5
for (l = 0 ; l<lmax; l++){
    mxSYST[entry_counter] = cone[l].value ;
    ir[entry_counter] = cone[l].vol_index ;
    entry_counter++ ;
}

En gros, lorsque l'on rempli une matrice creuse avec mxCreateSparse il faut bien faire attention à remplir dans le même ordre que les indices.