Dokumentacja pliku detMac.c

#include "MyLib.h"

Idź do kodu źródłowego tego pliku.

Funkcje
void	WyswietlMacierzInd (int n, int m, double A[][matrix_M], int P[])
void	WyswietlWiersz (int m, double w[])
double	ObliczWyznacznik (int *row, int P[], int n, int m, double A[][matrix_M])
int	main ()
Zmienne
const char	det_naF [] = "A.dat"

---

Dokumentacja funkcji

int main

(

)

Definicja w linii 106 pliku detMac.c.

Odwołuje się do det_naF, matrix_M, matrix_N, ObliczWyznacznik(), OdczytMacierzy(), WyswietlMacierz() i WyswietlMacierzInd().

{
   double det,A[matrix_N][matrix_M];
   int n,m,P[matrix_N],row;

   OdczytMacierzy(det_naF, &n, &m, A);
   printf("wczytana macierz: \"%s\"\n", det_naF);
   WyswietlMacierz(n, m, A);
   det = ObliczWyznacznik(&row, P, n, m, A);
   if( row!=n )
   {
      printf("Macierz \"%s\" jest osobliwa\n", det_naF);
      exit(-1);
   }
   puts("-----------------------------");
   printf("wyznacznik: %16le\n", det);
   puts("-----------------------------");
   puts("Wynikowa macierz gornotrojkatna:");
   WyswietlMacierzInd(n, m, A, P);
   puts("\nBez uwzglednienia przestawien wierszy:");
   WyswietlMacierz(n, m, A);
   puts("\nPermutacja wierszy:");
   for(row=0; row<n; row=row+1)
   {
      printf("%2i ",row+1);
   }
   puts("");
   for(row=0; row<n; row=row+1)
   {
      printf("%2i ",P[row]+1);
   }
   puts("");
   return 0;
}

double ObliczWyznacznik	(	int *	row,
		int	P[],
		int	n,
		int	m,
		double	A[][matrix_M]
	)

obliczenie wyznacznika metoda eleminacji Gaussa (Gaussian elimination),

Parametry:

	P	opisuje zamiane wierszy,
	row	okresla wiersz do ktorego metoda "doszla".

Definicja w linii 52 pliku detMac.c.

Odwołuje się do matrix_N.

Odwołania w main().

{
   double d[matrix_M],det;
   int i,j;

   for(i=0; i<matrix_N; i=i+1)
   {
      P[i] = i;
   }
   det = 1.0;
   for(*row=0; *row<n; *row=*row+1)
   {
      if( fabs(A[P[*row]][*row])<=1.0e-100 )
      {
         for(i=*row+1; i<n; i=i+1 )
         {
            if( 1.0e-100<fabs( A[P[i]][*row] ) )
            {
               j       = P[*row];
               P[*row] = P[i];
               P[i]    = j;
               det     = -det;
               break;
            }
         }
         if( i==n )
         {
            for(i=0; i<*row-1; i=i+1)
            {
               det = A[P[i]][i]*det;
            }
            return det;
         }
      }
      for(i=*row+1; i<m; i=i+1)
      {
         d[i] = A[P[*row]][i]/A[P[*row]][*row];
      }
      for(j=*row+1; j<n; j=j+1)
      {
         for(i=*row+1; i<m; i=i+1)
         {
            A[P[j]][i] = A[P[j]][i]-d[i]*A[P[j]][*row];
         }
         A[P[j]][*row] = 0.0;
      }
   }
   for(i=0; i<n; i=i+1)
   {
      det = A[P[i]][i]*det;
   }
   return det;
}

void WyswietlMacierzInd	(	int	n,
		int	m,
		double	A[][matrix_M],
		int	P[]
	)

wyswietlenie macierzy z uwzglednieniem permutacji wierszy

Parametry:

P

permutacja wierszy

Definicja w linii 22 pliku detMac.c.

Odwołania w main().

{
    int i,j;

    for(i=0; i<n; i=i+1)
    {
        for(j=0; j<m; j=j+1)
        {
            printf("%6.2lf ",A[P[i]][j]);
        }
        printf("\n");
    }
}

void WyswietlWiersz	(	int	m,
		double	w[]
	)

Definicja w linii 36 pliku detMac.c.

{
    int j;

    for(j=0; j<m; j=j+1)
    {
        printf("%6.2lf ",w[j]);
    }
    printf("\n");
}

---

Dokumentacja zmiennych

const char det_naF[] = "A.dat"

nazwa pliku z macierza

Definicja w linii 16 pliku detMac.c.

Odwołania w main().