Une classe d'algèbre (calcul matriciel) performante ?
Sylphe67
Messages postés
4
Date d'inscription
jeudi 10 juillet 2008
Statut
Membre
Dernière intervention
11 juillet 2008
-
10 juil. 2008 à 11:32
Bahoumi88 Messages postés 1 Date d'inscription samedi 7 février 2009 Statut Membre Dernière intervention 9 mars 2010 - 9 mars 2010 à 11:07
Bahoumi88 Messages postés 1 Date d'inscription samedi 7 février 2009 Statut Membre Dernière intervention 9 mars 2010 - 9 mars 2010 à 11:07
A voir également:
- Une classe d'algèbre (calcul matriciel) performante ?
- Calcul qtc ✓ - Forum C++ & C++ .NET
- Module classe ✓ - Forum Visual Basic 6
- Calcul signe astrologique - Forum Java
- Calcul moyenne bac - Forum Visual Basic 6
- Calcul volume fond bombé elliptique - Forum Visual Basic
3 réponses
cs_Lucky92
Messages postés
180
Date d'inscription
mercredi 22 décembre 2004
Statut
Membre
Dernière intervention
16 août 2012
2
10 juil. 2008 à 19:41
10 juil. 2008 à 19:41
Bonjour,
la librairie uBLAS de boost répondra peut-être à ton besoin. En revanche, la documentation n'est pas géniale, et il faut un sacré background en algèbre et en c++, ou éventuellement une bonne connaissance de la librairie fortran originale, pour l'utiliser (à mon sens en tout cas). voici, toutefois - si cela t'intéresse - l'exemple de base qui permet d'inverser une matrice avec l'élimination de gauss-jordan :
/*_________________________________________________________________________________
*/
#include
#include
/*_________________________________________________________________________________
*/
using namespace boost::numeric;
/*_________________________________________________________________________________
*/
int main( int argc, char** argv )
{
ublas::matrix<double> A( 2 , 2 ) ; A( 0 , 0 ) 1 ; A( 0 , 1 ) 2 ; A( 1 , 0 ) 3 ; A( 1 , 1 ) 4 ;
std::cout << A << std::endl;
ublas::permutation_matrix< std::size_t > pm( A.size1() );
int res( 0 );
if ( res = ublas::lu_factorize( A , pm ) )
{
std::cout << "matrix rank = " << res - 1 << std::endl;
}
else
{
ublas::matrix<double> inv( ublas::identity_matrix<double>( A.size1() ) );
ublas::lu_substitute( A , pm , inv );
std::cout << inv << std::endl;
}
}
/*_________________________________________________________________________________
*/
la librairie uBLAS de boost répondra peut-être à ton besoin. En revanche, la documentation n'est pas géniale, et il faut un sacré background en algèbre et en c++, ou éventuellement une bonne connaissance de la librairie fortran originale, pour l'utiliser (à mon sens en tout cas). voici, toutefois - si cela t'intéresse - l'exemple de base qui permet d'inverser une matrice avec l'élimination de gauss-jordan :
/*_________________________________________________________________________________
*/
#include
#include
/*_________________________________________________________________________________
*/
using namespace boost::numeric;
/*_________________________________________________________________________________
*/
int main( int argc, char** argv )
{
ublas::matrix<double> A( 2 , 2 ) ; A( 0 , 0 ) 1 ; A( 0 , 1 ) 2 ; A( 1 , 0 ) 3 ; A( 1 , 1 ) 4 ;
std::cout << A << std::endl;
ublas::permutation_matrix< std::size_t > pm( A.size1() );
int res( 0 );
if ( res = ublas::lu_factorize( A , pm ) )
{
std::cout << "matrix rank = " << res - 1 << std::endl;
}
else
{
ublas::matrix<double> inv( ublas::identity_matrix<double>( A.size1() ) );
ublas::lu_substitute( A , pm , inv );
std::cout << inv << std::endl;
}
}
/*_________________________________________________________________________________
*/
Sylphe67
Messages postés
4
Date d'inscription
jeudi 10 juillet 2008
Statut
Membre
Dernière intervention
11 juillet 2008
11 juil. 2008 à 14:40
11 juil. 2008 à 14:40
Bonjour,
Merci pour ta réponse. A première vue de l'exemple que tu donnes, je dirais que c'est une solution qui ne me convient pas vraiment, car il faut visiblement choisir la méthode d'inversion. Alors que je recherche - dans l'idéal - une solution qui va surcharger une méthode "inversion" (ou opérateur ^(-1), etc) de sorte que l'utilisateur pourrait se contenter d'écrire INV A^(-1), la librairie se chargeant d'effectuer le travail suivant : détection automatique des propriétés de A, choix d'un algorithme d'inversion en fct de ces propriétés, puis execution. Idem pour un système linéaire A x b, l'utilisateur se contenterait d'écrire : x = b / A ou x = A^(-1) * b.
Merci pour ta réponse. A première vue de l'exemple que tu donnes, je dirais que c'est une solution qui ne me convient pas vraiment, car il faut visiblement choisir la méthode d'inversion. Alors que je recherche - dans l'idéal - une solution qui va surcharger une méthode "inversion" (ou opérateur ^(-1), etc) de sorte que l'utilisateur pourrait se contenter d'écrire INV A^(-1), la librairie se chargeant d'effectuer le travail suivant : détection automatique des propriétés de A, choix d'un algorithme d'inversion en fct de ces propriétés, puis execution. Idem pour un système linéaire A x b, l'utilisateur se contenterait d'écrire : x = b / A ou x = A^(-1) * b.
Bahoumi88
Messages postés
1
Date d'inscription
samedi 7 février 2009
Statut
Membre
Dernière intervention
9 mars 2010
9 mars 2010 à 11:07
9 mars 2010 à 11:07
Bienvenue ,
j'aime la programmation bouceaups
J'aime tous sur l'informatique et je veut apprendre la langue anglaise ( grace à vos aides) .
j'aime la programmation bouceaups
J'aime tous sur l'informatique et je veut apprendre la langue anglaise ( grace à vos aides) .
