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

  1. #21
    Rédacteur/Modérateur

    Citation Envoyé par Dut Voir le message
    je rappelle que les structures itératives peuvent aussi être intégrées dans un fichier MEX...

    Quelqu'un veut essayer ?
    Bon ben je me lance...

    Le fichier Dutmex.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
    #include "mex.h"
     
    void mexFunction( int nlhs, mxArray *plhs[],
                      int nrhs, const mxArray *prhs[])
    {
    int i, j, m ,n;
    double *M, *X, *P;
     
    m = mxGetM(prhs[0]);
    n = mxGetN(prhs[0]);
     
    M = mxGetPr(prhs[0]);
    X = mxGetPr(prhs[1]);
     
    plhs[0] = mxCreateDoubleMatrix(m, n, mxREAL);
     
    P = mxGetPr(plhs[0]);
     
    j=0;
    for (i = 0; i < m*n; i++) 
        {
        j = (i / m);     
        P[i]=M[i]-X[j];
        }
    }


    A compiler comme ceci :

    Code :Sélectionner tout -Visualiser dans une fenêtre à part
    mex Dutmex.c


    Et le fichier m qui va avec :

    Code :Sélectionner tout -Visualiser dans une fenêtre à part
    1
    2
    3
    function C = Dut(A,B);
     
    C = Dutmex(A,B);
    Mes compétences :
    • conception mécanique 3D (Autodesk Fusion 360)
    • développement informatique (MATLAB, Python, C, VBA)
    • impression 3D (Ultimaker)
    • programmation de microcontrôleur (Microchip PIC et Arduino)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  2. #22
    Rédacteur/Modérateur

    Merci de faire tourner les nouveaux code pour pouvoir mettre les résultats à jour...

    Mes compétences :
    • conception mécanique 3D (Autodesk Fusion 360)
    • développement informatique (MATLAB, Python, C, VBA)
    • impression 3D (Ultimaker)
    • programmation de microcontrôleur (Microchip PIC et Arduino)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  3. #23
    Membre éprouvé
    Voici mes résultats:
    Code :Sélectionner tout -Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
       0.000078200000000   0.000134600000000   0.000062400000000   0.000084200000000   0.000327200000000
       0.000036400000000   0.000131000000000   0.000026000000000   0.000062200000000   0.000052400000000
       0.000050000000000   0.000142400000000   0.000039000000000   0.000082400000000   0.000072600000000
       0.000222000000000   0.000157200000000   0.000055000000000   0.000192800000000   0.000162600000000
       0.002152600000000   0.001704200000000   0.000903200000000   0.001920200000000   0.001592400000000
       0.023429800000000   0.027307800000000   0.028054800000000   0.019134200000000   0.018003600000000
       0.247296600000000   0.259066800000000   0.275122400000000   0.186231000000000   0.174070600000000


    Matlab 2008a Linux 32 bits 3.2GHz. Ram: 2Go

    [Edit=Dut] Premier post mis à jour
    ZiK un lecteur audio et son blog.

  4. #24
    Membre expert
    Donc après avoir modifié microdefi2.m pour rajouter l'utilisateur Dut et apres avoir mexé ...j'arrive premier sur la même Linux 64-bits 8 cores (multithreading on) en 8b

    [Edit=Dut] Oui mais tout le monde n'a pas un 8 cœurs
    Résultats mis à jour dans le premier post
    MATLAB 8.2 (R2013b)

    le blog des MATLAB geeks

###raw>template_hook.ano_emploi###