Comment accélérer le remplacement de Mathematica des éléments matriciels

Peut être xxx pré>

maintenant effectuez le remplacement p> xxx pré>

modifier (1) strong> p> p> P>

Voici une vitesse rapide comparant les 4 méthodes ci-dessus, la boucle, puis Brett, moi et Verbeia. Peut-être que quelqu'un puisse les vérifier. J'ai utilisé les mêmes données pour tous. créé des données aléatoires une fois, puis utilisées pour chaque test. Même limite (appelée L) J'ai utilisé la taille de matrice de 2 000 par 2 000 par 2 000. p>

SE SE SE Vitesse Les numéros de synchronisation ci-dessous n'incluent pas l'allocation de données. P>

i Exécutez les tests une fois. P>

C'est ce que je vois: p>

pour 2 000 sur 2 000 matrices: h2>

1. Bill (boucle): 16 secondes li>
2. moi ( replacePart code>): 21 secondes li>
3. Brett ( SPARSARRAY CODE>): 7.27 secondes li>
4. verbeia ( cartthread code>): 32 secondes li> ol>

   pour 3 000 sur 3 000 matrices: h2>

   1. Bill (boucle): 37 secondes li>
   2. moi ( replactpart code>): 48 secondes li>
   3. Brett ( SPARSARRAY CODE>): 16 secondes LI>
   4. verbeia ( cartthread code>): 79 secondes li> ol>

      Donc, il semble que SPARSARRAY code> est le plus rapide. (Mais s'il vous plaît vérifier pour vous assurer que je n'ai pas cassé quelque chose) p>

      Code ci-dessous: P>

      Génération de données forte> p>

      %test, on same machine, 4GB ram, timing uses cpu timing using tic/toc
      %allocate data
      nRow = 3000;
      nCol = 3000;
      
      %generate a random matrix of real values
      %between 1 and 3
      xnow = 1 + (3-1).*rand(nRow,nRow);
      
      %allocate the other 3 matrices
      a=zeros(nRow,nCol);
      b=a;
      c=b;
      
      %set limit
      limit=2;
      
      %engine
      tstart=tic;
      
      for i=1:nRow
          for j=1:nCol
              if xnow(i,j) > limit
                  xnow(i,j) = 0;
                  a(i,j) = 0;
                  b(i,j) = 0;
                  c(i,j) = 0;
              end
          end
      end
      toc(tstart)
      

remplapparts est notoirement lent.

cartthread devrait faire ce que vous voulez - notez le troisième argument. xxx

et pour les quatre matrices xxx