#include <stdio.h>
#define N 100

/* Vektorien summa, vakio kertaa vektori, ja vektorien ristitulo
 * antavat tulokseksi vektorin. Laskun suorittavalle funktiolle
 * annetaan argumentiksi taulukko (aina nimeltään z), johon tulos lasketaan.
 * Viimeinen argumentti on aina vektorin pituus. */
void vec_sum(double x[],double y[],double z[],int n);
void scal_mult(double x[],double a,double z[],int n);
double dot_prod(double x[],double y[],int n);
/* Ristitulon laskeva funktio palauttaa arvon 1, jos lasku onnistui
 * (ts. jos n oli 3), ja arvon 0 muuten. Tietenkin voitaisiin myös
 * pääohjelmassa tarkistaa, onko n==3. */
int cross_prod(double x[],double y[],double z[],int n);

void prn_vec(double x[],int n); /* Tulostaa n-komponenttisen vektorin. */

int main()
{
	/* Taulukkoon za lasketaan x+y, zb:hen ax ja zd:hen ristitulo. */
	double x[N], y[N], za[N], zb[N], zd[N];
	double a;
	int n,i;

	/* Kysytään käyttäjältä vektorien x ja y pituus ja niiden
	 * komponentit, sekä vakion a arvo. */

	printf("Enter dimension of vectors: "); scanf("%i",&n);

	printf("Enter components of x:\n");
	for(i=0;i<n;i++)
	{
		printf("x[%2i] = ",i); scanf("%lf",&x[i]);
	}

	printf("Enter components of y:\n");
	for(i=0;i<n;i++)
	{
		printf("y[%2i] = ",i); scanf("%lf",&y[i]);
	}

	printf("Enter scalar a: "); scanf("%lf",&a);

	/* Lasketaan x+y ja ax. */

	vec_sum(x,y,za,n);
	scal_mult(x,a,zb,n);

	/* Tulostetaan tuloksen prn_vec-funktiolla. */

	printf("%10s","x = "); prn_vec(x,n);
	printf("%10s","y = "); prn_vec(y,n);
	printf("%10s","x+y = "); prn_vec(za,n);
	printf("%10s","ax = "); prn_vec(zb,n);
	printf("%10s%6.2lf\n","x.y = ",dot_prod(x,y,n));

	/* Lasketaan ristitulo. Jos cross_prod palauttaa arvon 1,
	 * ristitulon lasku onnistui, joten tulos voidaan tulostaa. */

	if(cross_prod(x,y,zd,n))
	{
		printf("%10s","x x y = ");
		prn_vec(zd,n);
	}
	else
	{
		printf("Cross product of x and y not defined!\n");
	}
}

void vec_sum(double x[],double y[],double z[],int n)
{
	int i;
	for(i=0;i<n;i++)
	{
		z[i]=x[i]+y[i];
	}
}

void scal_mult(double x[],double a,double z[],int n)
{
	int i;
	for(i=0;i<n;i++)
	{
		z[i]=a*x[i];
	}
}

double dot_prod(double x[],double y[],int n)
{
	double result;
	int i;

	for(i=0,result=0.0;i<n;i++)
	{
		result += x[i]*y[i];
	}
	return result;
}

int cross_prod(double x[],double y[],double z[],int n)
{
	/* Tehdään laskutoimitukset vain jos n==3. */
	if(n==3)
	{
		z[0]=x[1]*y[2]-x[2]*y[1];
		z[1]=x[2]*y[0]-x[0]*y[2];
		z[2]=x[0]*y[1]-x[1]*y[0];
		return 1;
	}
	else
	{
		return 0;
	}
}

void prn_vec(double x[],int n)
{
	/* Tulostetaan vektori muodossa (x[0], x[1], x[2]). */
	int i;
	printf("(");
	for(i=0;i<n;i++)
	{
		printf("%6.2lf",x[i]);
		if(i<n-1)
			printf(",");
	}
	printf(")\n");
}