Discussion :

[epsilon68]

Bonjour,

j'aimerais votre aide pour comprendre pourquoi macosx n'est pas plus rapide en multi-thread qu'en single thread en rapport avec la discussion [MT] Benchmark de produits matriciels sur processeurs multi-cores

avec Windows XP avec VMWare c'est beaucoup plus rapide que sur macosx.
tout a l'air d'indiquer que le multithread est aussi rapide que le single thread !!!
quand j'ouvre le gestionnaire des taches, je vois que les deux cores sont utilisés...

Quelqu'un a-t-il une idée ?

Merci

[NiamorH]

Je n'ai pas eu ce probleme là lorsque je développais mon appli macosx multi core.

[epsilon68]

tu utilisais quoi ?
les pthreads ? les threads mac ?
t'as observé un speed-up par rapport au single thread ?

a+

[NiamorH]

J'utilisais l'API MP : http://developer.apple.com/documenta...rocessing.html

Et oui, mes calculs étant répartis sur 2 voire 4 processeurs, j'avais à peu près le même ratio de gain de perfs, si je me souviens bien.

[epsilon68]

c'est tres bizarre, peut-etre tu pourrais jeter un oeil sur l'appli de Charlemagne pour macosx ?

[NiamorH]

J'ai relu en diagonale les dernières pages, je ne sais pas trop ce que je pourrais faire de plus pour aider vu que je n'ai plus accès aux machines MacOsX.

Je pense qu'il faudrait trouver un autre dual core Mac pour tester...

Je veux bien essayer de retrouver des anciens bouts de code de test avec l'API Mac pour que tu fasse plus de tests.

[epsilon68]

pourquoi pas...

[PRomu@ld]

J'utilise les pthread pour un de mes projets, les performances sont bien là. Il doit y avoir un problème ailleurs .

[epsilon68]

tu voudrais tester le programme de Charlemagne pour macosx ?
j'ai un laptop macbook pro core duo et un autre tres recent macbook pro core 2 duo. les deux ont le meme symptome...

[PRomu@ld]

J'ai un MBP core 2 duo, je veux bien tester je récupère ça où ?

[Charlemagne]

Tu trouveras ci-joint les sources de ma librairie (GPL) ainsi que le programme de benchmark.
J'ai déconnecté le bench d'ATLAS pour que tu n'aies pas à l'installer.

Avec GCC 4.2 tu peux tester le multi-threading utilisant OpenMP (plutôt que mon implémentation) en ajoutant l'option "-DOPENMP".
Le multi-threading marche bien sous Windows avec GCC et ICL.
Mais même avec OpenMP les résultats ne sont pas au rendez-vous sur l'ordi d'Epsilon68, ce qui est très bizarre, car sinon j'aurais conclu à un bug de mon programme révélé uniquement sous MacOS.

Voici une ligne de commande qui devrait marcher (changer éventuellement les chemins d'accès).

Code :

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

1
2
 
/usr/local/bin/g++ main.cpp -O2 -I../genial -msse3 -Winline -DNDEBUG -DBLAS_THREADING -DMAX_THREADS=2  -framework Carbon -I/System/Library/Frameworks/CoreServices.framework/Frameworks/CarbonCore.framework/Headers

Si y'a des problèmes, je suis là...

[PRomu@ld]

Moi je ne suis pas contrariant mais :

In file included from ./genial/array/vectoriterator.h:21,
from ./genial/array/vector.h:49,
from ./genial/array/matrix.h:21,
from ./genial/blas/blas.h:22,
from ./genial/blas/copy.h:24,
from ./genial/blas/gemv.h:24,
from main.cpp:7:
./genial/array/arrayiterator.h:128: error: declaration of 'typedef struct shift_array_iterator_traits<typename array_iterator<A>::array_type::iterator> shift_array_iterator<A>::shift_array_iterator_traits'
./genial/array/arrayiterator.h:113: error: changes meaning of 'shift_array_iterator_traits' from 'struct shift_array_iterator_traits<typename array_iterator<A>::array_type::iterator>'
In file included from ./genial/array/matrix.h:21,
from ./genial/blas/blas.h:22,
from ./genial/blas/copy.h:24,
from ./genial/blas/gemv.h:24,
from main.cpp:7:
./genial/array/vector.h:94: error: declaration of 'typedef class std::reverse_iterator<typename generator_traits<G>::generator_type::iterator_rebind<Vector<G> >::other> Vector<G>::reverse_iterator'
/usr/local/ada-4.3/lib/gcc/i686-apple-darwin8/4.3.0/../../../../include/c++/4.3.0/bits/stl_iterator.h:100: error: changes meaning of 'reverse_iterator' from 'class std::reverse_iterator<typename generator_traits<G>::generator_type::iterator_rebind<Vector<G> >::other>'
./genial/array/vector.h: In function 'Vector<tiny_vector_generator<4, std::complex<float> > > shuffle(const Vector<tiny_vector_generator<4, std::complex<float> > >&, const Vector<tiny_vector_generator<4, std::complex<float> > >&)':
./genial/array/vector.h:1137: error: cannot convert 'const m256cf' to 'float __vector__' for argument '1' to 'float __vector__ _mm_shuffle_ps(float __vector__, float __vector__, int)'
In file included from ./genial/blas/blas.h:22,
from ./genial/blas/copy.h:24,
from ./genial/blas/gemv.h:24,
from main.cpp:7:
./genial/array/matrix.h: At global scope:
./genial/array/matrix.h:65: error: declaration of 'typedef class std::reverse_iterator<typename generator_traits<G>::generator_type::iterator_rebind<Matrix<G> >::other> Matrix<G>::reverse_iterator'
/usr/local/ada-4.3/lib/gcc/i686-apple-darwin8/4.3.0/../../../../include/c++/4.3.0/bits/stl_iterator.h:100: error: changes meaning of 'reverse_iterator' from 'class std::reverse_iterator<typename generator_traits<G>::generator_type::iterator_rebind<Matrix<G> >::other>'
In file included from ./genial/blas/gemm.h:27,
from main.cpp:8:
./genial/threads.h:381: error: declaration of 'typedef class gmt::scoped_lock<gmt::mutex> gmt::mutex::scoped_lock'
./genial/threads.h:130: error: changes meaning of 'scoped_lock' from 'class gmt::scoped_lock<gmt::mutex>'
./genial/threads.h:425: error: declaration of 'typedef class gmt::scoped_lock<gmt::condition_mutex> gmt::condition_mutex::scoped_lock'
./genial/threads.h:130: error: changes meaning of 'scoped_lock' from 'class gmt::scoped_lock<gmt::condition_mutex>'
./genial/threads.h:447: error: declaration of 'typedef class gmt::scoped_lock<gmt::semaphore> gmt::semaphore::scoped_lock'
./genial/threads.h:130: error: changes meaning of 'scoped_lock' from 'class gmt::scoped_lock<gmt::semaphore>'
In file included from main.cpp:10:
./genial/blas/trsm.h: In member function 'void simd_trsm010_function<N, Alpha, Unit, F>::operator()(const typename Matrix<G2>::value_type&, const Matrix<G>&, Matrix<G2>&) const':
./genial/blas/trsm.h:143: error: 'base_trsm010' was not declared in this scope
./genial/blas/trsm.h:143: error: missing template arguments before '>' token
./genial/blas/trsm.h: In member function 'void simd_trsm011_function<N, Alpha, Unit, F>::operator()(const typename Matrix<G2>::value_type&, const Matrix<G>&, Matrix<G2>&) const':
./genial/blas/trsm.h:144: error: 'base_trsm011' was not declared in this scope
./genial/blas/trsm.h:144: error: missing template arguments before '>' token
./genial/blas/trsm.h: In member function 'void simd_trsm100_function<N, Alpha, Unit, F>::operator()(const typename Matrix<G2>::value_type&, const Matrix<G>&, Matrix<G2>&) const':
./genial/blas/trsm.h:145: error: 'base_trsm100' was not declared in this scope
./genial/blas/trsm.h:145: error: missing template arguments before '>' token
./genial/blas/trsm.h: In member function 'void simd_trsm101_function<N, Alpha, Unit, F>::operator()(const typename Matrix<G2>::value_type&, const Matrix<G>&, Matrix<G2>&) const':
./genial/blas/trsm.h:146: error: 'base_trsm101' was not declared in this scope
./genial/blas/trsm.h:146: error: missing template arguments before '>' token
./genial/blas/trsm.h: In member function 'void simd_trsm110_function<N, Alpha, Unit, F>::operator()(const typename Matrix<G2>::value_type&, const Matrix<G>&, Matrix<G2>&) const':
./genial/blas/trsm.h:147: error: 'base_trsm110' was not declared in this scope
./genial/blas/trsm.h:147: error: missing template arguments before '>' token
./genial/blas/trsm.h: In member function 'void simd_trsm111_function<N, Alpha, Unit, F>::operator()(const typename Matrix<G2>::value_type&, const Matrix<G>&, Matrix<G2>&) const':
./genial/blas/trsm.h:148: error: 'base_trsm111' was not declared in this scope
./genial/blas/trsm.h:148: error: missing template arguments before '>' token
./genial/blas/trsm.h: In member function 'void simd_trsm010_function<0, Alpha, Unit, F>::operator()(const typename Matrix<G2>::value_type&, const Matrix<G>&, Matrix<G2>&) const':
./genial/blas/trsm.h:152: error: 'base_trsm010' was not declared in this scope
./genial/blas/trsm.h:152: error: missing template arguments before '>' token
./genial/blas/trsm.h: In member function 'void simd_trsm011_function<0, Alpha, Unit, F>::operator()(const typename Matrix<G2>::value_type&, const Matrix<G>&, Matrix<G2>&) const':
./genial/blas/trsm.h:153: error: 'base_trsm011' was not declared in this scope
./genial/blas/trsm.h:153: error: missing template arguments before '>' token
./genial/blas/trsm.h: In member function 'void simd_trsm100_function<0, Alpha, Unit, F>::operator()(const typename Matrix<G2>::value_type&, const Matrix<G>&, Matrix<G2>&) const':
./genial/blas/trsm.h:154: error: 'base_trsm100' was not declared in this scope
./genial/blas/trsm.h:154: error: missing template arguments before '>' token
./genial/blas/trsm.h: In member function 'void simd_trsm101_function<0, Alpha, Unit, F>::operator()(const typename Matrix<G2>::value_type&, const Matrix<G>&, Matrix<G2>&) const':
./genial/blas/trsm.h:155: error: 'base_trsm101' was not declared in this scope
./genial/blas/trsm.h:155: error: missing template arguments before '>' token
./genial/blas/trsm.h: In member function 'void simd_trsm110_function<0, Alpha, Unit, F>::operator()(const typename Matrix<G2>::value_type&, const Matrix<G>&, Matrix<G2>&) const':
./genial/blas/trsm.h:156: error: 'base_trsm110' was not declared in this scope
./genial/blas/trsm.h:156: error: missing template arguments before '>' token
./genial/blas/trsm.h: In member function 'void simd_trsm111_function<0, Alpha, Unit, F>::operator()(const typename Matrix<G2>::value_type&, const Matrix<G>&, Matrix<G2>&) const':
./genial/blas/trsm.h:157: error: 'base_trsm111' was not declared in this scope
./genial/blas/trsm.h:157: error: missing template arguments before '>' token
./genial/blas/trsm.h: In member function 'void simd_trsm010_function<1, Alpha, Unit, F>::operator()(const typename Matrix<G2>::value_type&, const Matrix<G>&, Matrix<G2>&) const':
./genial/blas/trsm.h:161: error: 'base_trsm010' was not declared in this scope
./genial/blas/trsm.h:161: error: missing template arguments before '>' token
./genial/blas/trsm.h: In member function 'void simd_trsm011_function<1, Alpha, Unit, F>::operator()(const typename Matrix<G2>::value_type&, const Matrix<G>&, Matrix<G2>&) const':
./genial/blas/trsm.h:162: error: 'base_trsm011' was not declared in this scope
./genial/blas/trsm.h:162: error: missing template arguments before '>' token
./genial/blas/trsm.h: In member function 'void simd_trsm100_function<1, Alpha, Unit, F>::operator()(const typename Matrix<G2>::value_type&, const Matrix<G>&, Matrix<G2>&) const':
./genial/blas/trsm.h:163: error: 'base_trsm100' was not declared in this scope
./genial/blas/trsm.h:163: error: missing template arguments before '>' token
./genial/blas/trsm.h: In member function 'void simd_trsm101_function<1, Alpha, Unit, F>::operator()(const typename Matrix<G2>::value_type&, const Matrix<G>&, Matrix<G2>&) const':
./genial/blas/trsm.h:164: error: 'base_trsm101' was not declared in this scope
./genial/blas/trsm.h:164: error: missing template arguments before '>' token
./genial/blas/trsm.h: In member function 'void simd_trsm110_function<1, Alpha, Unit, F>::operator()(const typename Matrix<G2>::value_type&, const Matrix<G>&, Matrix<G2>&) const':
./genial/blas/trsm.h:165: error: 'base_trsm110' was not declared in this scope
./genial/blas/trsm.h:165: error: missing template arguments before '>' token
./genial/blas/trsm.h: In member function 'void simd_trsm111_function<1, Alpha, Unit, F>::operator()(const typename Matrix<G2>::value_type&, const Matrix<G>&, Matrix<G2>&) const':
./genial/blas/trsm.h:166: error: 'base_trsm111' was not declared in this scope
./genial/blas/trsm.h:166: error: missing template arguments before '>' token
main.cpp: In member function 'void atlas_gemv_function<float>::operator()(CBLAS_TRANSPOSE, const Matrix<G>&, const Vector<G>&, Vector<G2>&)':
main.cpp:84: error: 'AtlasRowMajor' was not declared in this scope
main.cpp: In member function 'void atlas_gemv_function<double>::operator()(CBLAS_TRANSPOSE, const Matrix<G>&, const Vector<G>&, Vector<G2>&)':
main.cpp:85: error: 'AtlasRowMajor' was not declared in this scope
main.cpp: In member function 'void atlas_gemv_function<std::complex<float> >::operator()(CBLAS_TRANSPOSE, const Matrix<G>&, const Vector<G>&, Vector<G2>&)':
main.cpp:86: error: 'AtlasRowMajor' was not declared in this scope
main.cpp:86: error: there are no arguments to 'cblas_cgemv' that depend on a template parameter, so a declaration of 'cblas_cgemv' must be available
main.cpp:86: error: (if you use '-fpermissive', G++ will accept your code, but allowing the use of an undeclared name is deprecated)
main.cpp: In member function 'void atlas_gemv_function<std::complex<double> >::operator()(CBLAS_TRANSPOSE, const Matrix<G>&, const Vector<G>&, Vector<G2>&)':
main.cpp:87: error: 'AtlasRowMajor' was not declared in this scope
main.cpp:87: error: there are no arguments to 'cblas_zgemv' that depend on a template parameter, so a declaration of 'cblas_zgemv' must be available
main.cpp: In function 'void atlas_gemv(const Matrix<G>&, const Vector<G2>&, Vector<G2>&)':
main.cpp:92: error: 'AtlasNoTrans' was not declared in this scope
main.cpp:93: error: 'AtlasTrans' was not declared in this scope
main.cpp: In member function 'void atlas_gemm_function<float>::operator()(CBLAS_TRANSPOSE, CBLAS_TRANSPOSE, const Matrix<G>&, const Matrix<G>&, Matrix<G>&)':
main.cpp:141: error: 'AtlasRowMajor' was not declared in this scope
main.cpp: In member function 'void atlas_gemm_function<double>::operator()(CBLAS_TRANSPOSE, CBLAS_TRANSPOSE, const Matrix<G>&, const Matrix<G>&, Matrix<G>&)':
main.cpp:142: error: 'AtlasRowMajor' was not declared in this scope
main.cpp: In member function 'void atlas_gemm_function<std::complex<float> >::operator()(CBLAS_TRANSPOSE, CBLAS_TRANSPOSE, const Matrix<G>&, const Matrix<G>&, Matrix<G>&)':
main.cpp:143: error: 'AtlasRowMajor' was not declared in this scope
main.cpp:143: error: there are no arguments to 'cblas_cgemm' that depend on a template parameter, so a declaration of 'cblas_cgemm' must be available
main.cpp: In member function 'void atlas_gemm_function<std::complex<double> >::operator()(CBLAS_TRANSPOSE, CBLAS_TRANSPOSE, const Matrix<G>&, const Matrix<G>&, Matrix<G>&)':
main.cpp:144: error: 'AtlasRowMajor' was not declared in this scope
main.cpp:144: error: there are no arguments to 'cblas_zgemm' that depend on a template parameter, so a declaration of 'cblas_zgemm' must be available
main.cpp: In member function 'void atlas_trsm_function<std::complex<float> >::operator()(CBLAS_SIDE, CBLAS_UPLO, CBLAS_TRANSPOSE, CBLAS_DIAG, const Matrix<G>&, Matrix<G>&)':
main.cpp:271: error: there are no arguments to 'cblas_ctrsm' that depend on a template parameter, so a declaration of 'cblas_ctrsm' must be available
main.cpp: In member function 'void atlas_trsm_function<std::complex<double> >::operator()(CBLAS_SIDE, CBLAS_UPLO, CBLAS_TRANSPOSE, CBLAS_DIAG, const Matrix<G>&, Matrix<G>&)':
main.cpp:272: error: there are no arguments to 'cblas_ztrsm' that depend on a template parameter, so a declaration of 'cblas_ztrsm' must be available

EDIT: j'ai gcc 4.3

[epsilon68]

il faudrait deja enlevé ATLAS je pense pour un test simple...
Charlemagne: peux-tu enlever les appels a ATLAS dans le main.cpp (moi j'enleve toujours ce qu'il ne faut pas )

EDIT: j'ai rien dit c'est deja fait

bizarre les erreurs

[Charlemagne]

C'est quoi le compilo que tu utilises? Normalement ça passe sans problème avec GCC, VC2005, ICC et ICL (rien que ça...)
J'ai vérifié vite fait pour la 1ère erreur de compilation, et y'a rien d'anormal dans mon code. Je vais voir ce que je peux faire quand-même, mais je comprends pas le problème de ton compilo...

[PRomu@ld]

gcc 4.3

[Charlemagne]

A première vue la série des 1ères erreurs est due à un conflit entre des définitions locale de types et des définition globales. Ca me paraît pas normal de la part d'un compilo. Je vais voir si je peux renommer quelques définitions.

[Charlemagne]

gcc 4.3

Il est déjà sorti? C'st pas une vesion test?
Bon je vais essayer de faire les changements qui s'imposent, faudra que je les fasse tôt ou tard...

[Charlemagne]

Essaye avec ceci.
Toutes les erreurs de compilation ne seront peut-être pas parties, mais on devrait y voir beaucoup plus clair.

[epsilon68]

GCC 4.3 rentre juste en stage 3.
Honnetement je ne pense pas qu'il faille modifier le code pour l'instant pour GCC 4.3, il y a encore beaucoup trop de regressions.

[epsilon68]

j'ai aussi ouvert la discussion du coté mac pour elargir le nombre potentiel de testeurs...

[MACOSX] APPLICATION MULTITHREAD