#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <math.h>
float *hilbert_t( float *M,int p, int q )
{
	int i,j;
	float a,b;
	M = malloc(sizeof (float)*q*p);
	for (i=0;i<p;i++)
	{
		for (j=0;j<q;j++)
		{
			a=i;
			b=j;
			M[i*q + j]=1/(float)(1+a+b);
		}
	}
return M;
}
float *aleatoire_t ( int p, int q)
{
	int i,j;
	float *M = malloc(sizeof (float)*q*p);
	for (i=0;i<p;i++)
	{
		for (j=0;j<q;j++)
		{
			M[i*q+j]= rand()%10  ;
		}
	}
return M;
}
float *stochastique_t ( int p,int q)
{
	int i,j,s=0;
	float *M = malloc(sizeof (float)*q*p);
	for (i=0;i<p;i++)
	{
		for (j=0;j<q;j++)
		{
			M[i*q+j]=rand();
			s=s+M[i*q+j];
		}
		for (j=0;j<q;j++)
		{
		
			M[i*q+j]= M[i*q+j]/s;
		}
	}
return M;
}
float *circulante_t(  int p, int q)
{
	int i,j;
	float *M = malloc(sizeof (float)*q*p);
	for (i=0;i<p;i++)
	{
		for (j=0;j<q;j++)
		{
			if (i==0)
				{
				M[i*p+j]=rand()%20;
				}
			else 
			{
				if (j==(q-1))
				{
					M[i*q+j]=M[(i-1)*q+0];
				}
				else
				{
					M[i*q+j]= M[(i-1)*q+(j+1)];
				}	
			}
		}
	}
return M;
}
float *dilation_t ( float *M, int p,int q,float s,float *N)
{
	int i,j;
	M = malloc(sizeof (float)*q*p);
	N = malloc(sizeof (float)*q*p);
	for (i=0;i<p;i++)
	{
		for (j=0;j<q;j++)
		{ 
			N[i*q+j]= s * M[i*q+j];
		}
	}
return N;
} 
float *addition_t ( float *A,int p,int q,float *B,int r,int s, float *C)
{
	int i,j;	
	A = malloc(sizeof (float)*q*p);
	B = malloc(sizeof (float)*r*s);
	C = malloc(sizeof (float)*q*p);
	if (( p==r) & (q==r))
	{
		for (i=0;i<p;i++)
		{
			for (j=0;j<q;j++)
			{ 
				C[r*q+s]= A[r*q+s]+ B[r*q+s];
			}
		}
	return C;
	}
	else
	{
	 return 0;
	}
}
float *multiplier_t ( float *A,int p,int q, float *B,int r,int s,float *C)
{
	int i,j,PLUS,w=0;
	A = malloc(sizeof (float)*q*p);
	B= malloc(sizeof (float)*r*s);
	C = malloc(sizeof (float)*q*s);
	if ( q==s )	
	{
		for ( i=0; i<p; i++)
		{
			for ( j=0;j<s;j++)
			{
			PLUS=0;
			for ( w=0;w<q;w++)
			{
			PLUS= A[i*q+w]*B[w*q+j];
			C[i*q+j]= C[i*q+j] + PLUS;
				}
			}
		
		}
	return C;
	}
	else
	{
	 return 0;
	}
}
void *affichage (float *A, int p,int q)
{
	int i,j;
	for (i=0;i<p;i++)
	{
		for (j=0;j<q;j++)
		{
			printf("%.2f\t",A[i*q+j]);
		}
	printf("\n");
	}
 }

int main()
{
float N[2][3]= { {1,1,1},{2,2,2}};
float A[2][3];
dilation_t(*N, 2, 3, 10,*A);

affichage ( *A, 2, 3);


return EXIT_SUCCESS;
}