/* atk1 h7t3 */

#include <stdio.h>

#define N 3

/* argumentit ovat reaalilukutaulukkoja (tai, yhtäpitävästi,
osoittimia reaalilukutaulukkojen ensimmäisiin alkioihin) */
double dot_product(double*,double*);
/* kolmanteen argumenttiin lasketaan itse ristitulovektori */
void cross_product(double*,double*,double*);

main()
{
	double x[N]={1.0,-2.3,3.0};
	double y[N]={-0.2,3.2,1.1};
	double w[N];

	printf("x = (%.2lf,%.2lf,%.2lf)\n",x[0],x[1],x[2]);
	printf("y = (%.2lf,%.2lf,%.2lf)\n\n",y[0],y[1],y[2]);

	printf("x . y = %.2lf\n",dot_product(x,y));

	cross_product(x,y,w);

	printf("x x y = (%.2lf,%.2lf,%.2lf)\n",w[0],w[1],w[2]);
}

double dot_product(double* x,double* y)
{
	int i;
	double sum;
	for(i=0,sum=0.0;i<N;i++)
		sum += (*(x+i))*(*(y+i));
	return sum;
}

void cross_product(double* x,double* y,double* w)
{
	/* kahden vektorin A ja B ristitulo määritellään kaavalla
	A x B = i(A_yB_z-A_zB_y)+j(A_zB_x-A_xB_z)+k(A_xB_y-A_yB_x) 
	nyt x-komponenttia vastaa A[0], y:tä A[1] ja z:aa A[2] */

	/* muista, että jos double* p = taulu, niin
	p[i] = taulu[i] */
	w[0]=x[1]*y[2]-x[2]*y[1];
	w[1]=x[2]*y[0]-x[0]*y[2];
	w[2]=x[0]*y[1]-x[1]*y[0];
}