fant3
17-05-2004, 21:12
Devo calcolare il determinante di una matrice quadrata, usando l'allocazione dinamica per la matrice e funzioni void.
Questo è quello che ho scritto:
#include <stdio.h>
#include <stdlib.h>
void scambia_righe(float **matrice, int n, int r1, int r2, int start)
{
int i;
float tmp;
for(i=start; i<n; ++i)
{
tmp = matrice[r1][i];
matrice[r1][i] = matrice[r2][i];
matrice[r2][i] = tmp;
}
}
void det(float **matrice, int n, float *risultato, int start)
{
int r;
int sgn = 1;
float det_min;
*risultato = 0;
if(start == (n - 1))
return matrice[n-1][n-1];
for(r=start; r<n; ++r)
{
if(r != start)
{
sgn = -1; /*scambiando due righe il determinante del minore cambia di segno*/
scambia_righe(matrice, n, start, r, start);
}
det(matrice, n, &det_min, start+1);
*risultato = sgn * matrice[start][start] * det_min;
if(r != start)
scambia_righe(matrice, n, start, r, start);
}
}
int main()
{
int n, i, j;
float **mat, risultato;
printf("inserisci l'ordine della matrice quadrata\n");
scanf("%d", &n);
mat = (float **)calloc(n, sizeof(float *));
for(i=0;i<n;i++)
{
mat[i] = (float *)calloc(n, sizeof(float));
}
printf("inserisci gli elementi della matrice \n");
for(i=0;i<n;i++)
{
for(j=0;j<n;j++)
{
scanf("%f", &mat[i][j]);
}
}
det(mat, n, &risultato, 0);
printf("Determinante = %f", risultato);
free(mat);
}
Come faccio a sostituire quel return nella funzione det?
Questo è quello che ho scritto:
#include <stdio.h>
#include <stdlib.h>
void scambia_righe(float **matrice, int n, int r1, int r2, int start)
{
int i;
float tmp;
for(i=start; i<n; ++i)
{
tmp = matrice[r1][i];
matrice[r1][i] = matrice[r2][i];
matrice[r2][i] = tmp;
}
}
void det(float **matrice, int n, float *risultato, int start)
{
int r;
int sgn = 1;
float det_min;
*risultato = 0;
if(start == (n - 1))
return matrice[n-1][n-1];
for(r=start; r<n; ++r)
{
if(r != start)
{
sgn = -1; /*scambiando due righe il determinante del minore cambia di segno*/
scambia_righe(matrice, n, start, r, start);
}
det(matrice, n, &det_min, start+1);
*risultato = sgn * matrice[start][start] * det_min;
if(r != start)
scambia_righe(matrice, n, start, r, start);
}
}
int main()
{
int n, i, j;
float **mat, risultato;
printf("inserisci l'ordine della matrice quadrata\n");
scanf("%d", &n);
mat = (float **)calloc(n, sizeof(float *));
for(i=0;i<n;i++)
{
mat[i] = (float *)calloc(n, sizeof(float));
}
printf("inserisci gli elementi della matrice \n");
for(i=0;i<n;i++)
{
for(j=0;j<n;j++)
{
scanf("%f", &mat[i][j]);
}
}
det(mat, n, &risultato, 0);
printf("Determinante = %f", risultato);
free(mat);
}
Come faccio a sostituire quel return nella funzione det?