#include <stdio.h>  /* pour les entrer/sortie... */
#include <stdlib.h> /* pour les allocations dynamique... */
//#include "nrutil.h"
#include <math.h>



/* Les prototypes de fonctions utilisé */
int Menu ( void );

float** Allocate_Mat_Memory ( int );
void Free_Mat_Memory ( float**, int );

void Affiche_Mat ( float**, int );
void Saisi_Mat ( float**, int );
void qrdcmp(float **A, int N, float *c, float *d, int *sing);

/* Cette fonction permet d'allouer de la memoire pour une matrice qulquonque */
float** Allocate_Mat_Memory(int N)
{
    int i;
    float **A = NULL;

    /* On fait l'allocation de memoire pour un vecteur de N cases,
      C'est un vecteur de pointeurs */
    A = (float**) malloc( sizeof(float*) * N );

    /* Si l'allocation à réussi, On fait une allocation de N cases
       pour chaque ligne */
    if( A != NULL )
        for( i=0; i<N; i++ )
            A[i] = (float*) malloc( sizeof(float) * N );

    return A;

}


/* cette fonction fait la désallocation d'une matrice (liberer la memoire)*/
void Free_Mat_Memory( float **A, int N )
{
    int i;
    /* librer chaque ligne: */
    for( i=0; i<N; i++ )
        free( A[i] );

    free( A );
}

/* Affichage des élement de la matrice */
void Affiche_Mat( float **A, int N )
{
    int i, j;

    for(i=0; i<N; i++)
    {
        for(j=0; j<N; j++)
		printf("%f\t", A[i][j]);
		printf("\n");
    }
	
}



/* Saisir les élément de la matrice */
void Saisi_Mat( float **A, int N )
{
    int i, j;
//char titre;
    FILE *fpt;
	
	

    fpt = fopen("notesma11.dat","r");
    if(fpt == NULL)
	{
		printf("\nError: File couldn't be opened (1).\n");
		exit(-1);
	}
//fscanf(fpt, "%c",&titre);
	for(i= 0; i < N; i++)
		for(j = 0; j < N; j++)
			fscanf(fpt, "%f\t",&A[i][j]);
	fclose(fpt);

//	return 0;

}


void qrdcmp(float **A, int N, float *c, float *d, int *sing)
{
	int i,j,k;
	float scale,sigma,sum,tau;
	
	*sing=0;
	for (k=0;k<N-1;k++) 
		{
		scale=0.0;
		for (i=k;i<N;i++) 
			{
			if (scale>=fabs(A[i][k])) 
				{
				scale=scale;
				} 
			else 
				{
				scale=fabs(A[i][k]);
				}
			}
		
		if (scale == 0.0) 
			{ 
			*sing=1;
			c[k]=d[k]=0.0;
			} 
		else 
			{ 
			for (i=k;i<N;i++) A[i][k] /= scale;
			for (sum=0.0,i=k;i<=N;i++) sum += (A[i][k])*(A[i][k]);
			if (A[k][k]>=0) 
			{
				sigma=sqrt(sum);
			}
			else 
			{
				sigma=-sqrt(sum);
			}
			A[k][k] += sigma;
			c[k]=sigma*A[k][k];
			d[k] = -scale*sigma;
			for (j=k+1;j<N;j++) 
				{
				for (sum=0.0,i=k;i<=N;i++) sum += A[i][k]*A[i][j];
				tau=sum/c[k];
				for (i=k;i<N;i++) A[i][j] -= tau*A[i][k];
				}
			}
		}
	d[N]=A[N-1][N-1];
	if (d[N-1] == 0.0) *sing=1;
}

void main()//void)
{ 
    FILE *fpt;
    float *c=0, *d=0;
    int *sing=0;
    float **matrice1;
    int N=1;
    char caractere;

    fpt = fopen("notesma11.dat","r");
	if(fpt == NULL)
	{
		printf("\nError: File couldn't be opened (1).\n");
		exit(-1);
	}
	do
        {
            caractere = fgetc(fpt); // On lit le caractère
          
			if(caractere == '\n')
			{
				N=N+1;
			}
			

        } while (caractere != EOF);
    fclose(fpt);

          
                printf("Nombre de lignes: %d  \n",N);     
                matrice1 = Allocate_Mat_Memory( N );
                         
                Saisi_Mat( matrice1, N );
                printf("La Matrice saisie:\n\n");
                Affiche_Mat( matrice1, N);
		  qrdcmp(matrice1, N, c, d, sing);
		 Affiche_Mat( matrice1, N);
        //       return 0;    
}