/* Kaksiulotteinen taulukko funktion argumenttina */

#include <stdio.h>
#define M 3
#define N 4

/* Kun kaksiulotteinen taulukko on funktion argumenttina, "vaakarivin"
 * pituus täytyy ilmoittaa, jotta ohjelma tietää, kuinka kaukana toisistaan
 * tietokoneen muistissa peräkkäsillä vaakariveillä olevat alkiot ovat. */
void print_matrix(double[][N]);

/* Matriin alkiot voidaan käydä läpi myös pointterin avulla (ks. toteutus). */
void print_matrix_using_pointers(double*);

main()
{
	double A[M][N];
	int m,n;

	for(m=0;m<M;m++)
		for(n=0;n<N;n++)
			A[m][n]=(m+n)*(m+n);

	print_matrix(A);

	printf("\n");

	/* Kääntäjä saattaa valittaa, että A on vääränlainen pointteri.
	 * Funktion pitäisi toimia kuitenkin oikein. Voitaisiin myös
	 * kirjoittaa print_matrix_using_pointers(&A[0][0]), jolloin
	 * valitusta ei tule. */
	print_matrix_using_pointers(A);
}

void print_matrix(double A[M][N])
{
	int i,j;

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

void print_matrix_using_pointers(double* p)
{
	int i,j;

	/* Miten päästää käsiksi alkioon [i][j], kun alkion [0][0] osoite
	 * tiedetään? Samalla pystyrivillä, peräkkäsillä vaakariveillä olevien
	 * alkioiden etäisyys on N*(yhden double-muuttujan viemä tila).
	 * Vastaavasti samalla vaakarivillä, peräkkäisillä pystyriveillä
	 * olevat alkiot ovat muistissa peräkkäin.
	 *
	 * Kun muistetaan pointtereiden laskusäännöt, nähdään että alkion
	 * [i][j] osoite on p+i*N+j. */

	for(i=0;i<M;i++)
	{
		for(j=0;j<N;j++)
		{
			printf("%6.2lf",*(p+i*N+j));
		}
		printf("\n");
	}

	/* Koska kaikki matriisin alkiot ovat muistissa peräkkäin, voitaisiin
	 * käyttää yhtä indeksiä, joka saa arvot 0...M*N-1:
	 *
	 * for(i=0;i<M*N;i++)
	 * {
	 * 		printf("%6.2lf",*(p+i));
	 * 		if((i-1)%N==0)
	 * 			printf("\n");  N alkion välein rivinsiirto 
	 * }
	 * */
}