//méthode du produit matriciel multiple
public Matrice FaireProduitMatriciel(out int operations, params Matrice[] matrices)
{
operations = 0;
Matrice temp = this;
for (int i = 0; i < matrices.Length; i++)
{
operations += temp.GetLength(0) * temp.GetLength(1) * matrices[i].GetLength(1);
temp = temp.FaireProduitMatriciel(matrices[i]);
}
return(temp);
}
//méthode du produit matriciel en commencant par la fin
public Matrice FaireProduitMatricielParFin(out int operations, params Matrice[] matrices)
{
operations = 0;
Matrice matriceA = this;
Matrice temp = this;
for (int i = matrices.Length - 1; i == 1; i--)
{
operations += matrices[i - 1].GetLength(0) * matrices[i - 1].GetLength(1) * matrices[i].GetLength(1);
temp = matrices[i - 1].FaireProduitMatriciel(matrices[i]);
}
operations += matriceA.GetLength(0) * matriceA.GetLength(1) * temp.GetLength(1);
temp = matriceA.FaireProduitMatriciel(temp);
return(temp);
}
static void Main(string[] args)
{
Matrice matrice1 = new Matrice(new double[3, 3] {
{ 3, 4, 1 }, { 4, 1, 2 }, { 4, 0, 2 }
});
Matrice matrice2 = new Matrice(new double[2, 3] {
{ 1, 2, 1 }, { 1, 2, 1 }
});
Matrice matrice3 = new Matrice(new double[3, 2] {
{ 3, 2 }, { 1, 2 }, { 1, 1 }
});
Matrice matrice4 = new Matrice(new double[2, 2] {
{ 3, 2 }, { 4, 5 }
});
Matrice matrice5 = new Matrice(new double[4, 4] {
{ 4, 2, 8, 3 }, { 5, 1, 7, 5 }, { 8, 0, 8, 5 }, { 3, 2, 3, 8 }
});
Matrice matrice6 = new Matrice(new double[3, 3] {
{ 1, 2, -1 }, { -2, 1, 1 }, { 0, 3, -3 }
});
Matrice matrice7 = new Matrice(new double[4, 4] {
{ 0, 2, 8, 3 }, { 0, 0, 7, 5 }, { 0, 0, 0, 5 }, { 0, 0, 0, 0 }
});
Matrice matriceA = new Matrice(new double[3, 3] {
{ 1, 3, 4 }, { 3, 5, -4 }, { 4, 7, -2 }
});
Matrice matriceB = new Matrice(new double[3, 1] {
{ 50 }, { 2 }, { 31 }
});
Systeme systeme1 = new Systeme(matriceA, matriceB);
Matrice matriceC = new Matrice(new double[3, 3] {
{ 1, 1, -2 }, { 4, -1, 2 }, { 2, -6, 3 }
});
Matrice matriceD = new Matrice(new double[3, 1] {
{ 0 }, { 0 }, { 0 }
});
Systeme systeme2 = new Systeme(matriceC, matriceD);
Matrice matriceE = new Matrice(new double[3, 3] {
{ 4, -1, 0 }, { -1, 4, -1 }, { 0, -1, 4 }
});
Matrice matriceF = new Matrice(new double[3, 1] {
{ 100 }, { 100 }, { 100 }
});
Systeme systeme3 = new Systeme(matriceE, matriceF);
////Additionner
//Console.WriteLine("--------------Addition----------------");
//(matrice1.Additionner(matrice6)).AfficheMatrice();
//Console.WriteLine();
////Produit Scalaire
//Console.WriteLine("----------Produit Scalaire----------------");
//matrice1.FaireProduitScalaire(6).AfficheMatrice();
//Console.WriteLine();
////Produit de matrices
Console.WriteLine("--------------Produit----------------");
int a = 0;
Console.WriteLine((matrice2.FaireProduitMatriciel(out a, matrice3, matrice4)).AfficheMatrice());
Console.WriteLine("Nombre de produits effectues : " + a);
Console.WriteLine();
////Triangulaire ?
//Console.WriteLine("--------------Triangulaire?----------------");
//bool m = matrice7.EstTriangulaire(1, true);
//matrice7.AfficheMatrice();
//Console.WriteLine("La matrice 7 est elle superieure strict ? " +m);
//Console.WriteLine();
//Console.WriteLine();
////Trace
//Console.WriteLine("--------------Trace----------------");
//double n = matrice5.Trace;
//matrice5.AfficheMatrice();
//Console.WriteLine("La matrice 5 a une trace equivalent à " + n);
//Console.WriteLine();
//Determinant
//Console.WriteLine("--------------Determinant----------------");
//double c = matrice5.Determinant;
//Console.WriteLine("Le determinant de la matrice 5 est de "+c);
//Console.WriteLine();
////Transposee
//Console.WriteLine("------------Transposee----------------");
//Console.WriteLine("matrice 6 avant");
//matrice6.AfficheMatrice();
//Console.WriteLine("matrice 6 apres");
//matrice6.Transposee.AfficheMatrice();
//Console.WriteLine();
////Comatrice
//Console.WriteLine("------------Comatrice----------------");
//Console.WriteLine("matrice 4 avant operation");
//matrice4.AfficheMatrice();
//Console.WriteLine("Apres");
//matrice4.Comatrice.AfficheMatrice();
//Console.WriteLine();
////Matrice inverse
//Console.WriteLine("------------Matrice Inverse----------------");
//Console.WriteLine("matrice 5 avant operation");
//matrice5.AfficheMatrice();
//Console.WriteLine("Apres operation");
//matrice5.MatriceInverse.AfficheMatrice();
//Console.WriteLine();
////Est reguliere ?
//Console.WriteLine("--------------Matrice Reguliere?----------------");
//bool o = matrice3.EstReguliere;
//double p = matrice3.Determinant;
//Console.WriteLine("Le determinant de la matrice 3 est de " + p);
//Console.WriteLine("La matrice 3 est elle reguliere ? " + o);
//Console.WriteLine();
//Console.WriteLine("--------------Méthode de CRAMER----------------");
//systeme1.TrouverXParCramer().AfficheMatrice();
//Console.WriteLine();
//systeme2.TrouverXParCramer().AfficheMatrice();
//Console.WriteLine("--------------Méthode d'inversion matricielle----------------");
//systeme1.TrouverXParInversionMatricielle().AfficheMatrice();
//Console.WriteLine("--------------Méthode de Jacobi----------------");
//systeme3.TrouverXParJacobi(0.02).AfficheMatrice();
Console.WriteLine();
Console.WriteLine("Appuyer sur Enter pour fermer");
Console.ReadLine();
//MODIFIER POUR QUE L'UTILISATEUR PUISSE ENTRER SES PROPRES MATRICES
}
