[MATLAB] alternativa a maxNumCompThreads()

[DanteCpp]

Si si, il mio algoritmo è in c

Pero non regge il confronto.

Comunque sia malgrado loro ci lavorino da tempo, penso che le performance di algoritmi del genere sono strettamente dipendenti dall'architettura. Quindi con il cambio generazionale anche loro hanno dovuto "reiniziare".

Oltretutto ci sono tre livelli: BLAS-LAPACK-APP

Credo proprio che matlab utilizzi le blas della intel, e chi meglio di loro conosce le architetture.

[dimaios]

Se vuoi ottenere performance degne di nota devi impiegare librerie in C++ mirate per il calcolo.
La Armadillo è una di queste.
Il segreto primario consiste nell'utilizzo dei template metaprogramming quindi è anche una questione di stile di programmazione e capacità nel software engineering.

[DanteCpp]

Various matrix decompositions are provided through integration with LAPACK, or one of its high performance drop-in replacements (eg. multi-threaded Intel MKL, or AMD ACML, or OpenBLAS)

Cosi non sarebbe più divertente, se avessi la necessità di fattorizzare matrici; userei matlab, mathematica o armadillo. Ma essendo un gioco...

Io penso, che molto si possa ottenere ottimizzando dal punto di vista dell'architettura di memoria. Per esempio sul mio sandy bridge a tre livelli di cache.

Per fattorizzare una matrice di ordine 10k, se implemento la matrice in un unico vettore di memoria è necessario circa un decimo del tempo di un'implementazione come vettore di vettori.

Il classico

Codice: Seleziona tutto

float *

contro

Codice: Seleziona tutto

float **

per intenderci.

Questo perche i dati sono molto più localizzati nel primo caso, quindi diminuiscono i cache miss.

[dimaios]

Se vuoi giocare puoi fare qualsiasi cosa ma ti assicuro che il "C" puro non è più il linguaggio principe per questo tipo di applicazioni ormai da molto tempo.
Se guardi il codice sorgente di alcune librerie matematiche moderne prendi paura, le ottimizzazioni sono incredibili e fanno spesso riferimento ad articoli specialistici recenti.

Comunque puoi fare qualche esperimento e confrontarlo con le librerie in rete.
Facci sapere dove riesci ad arrivare.
Buon divertimento.

[DanteCpp]

Facci sapere dove riesci ad arrivare.
Buon divertimento.

Senz'altro, grazie!