向（矢）量叉乘法編程實現

在一個朋友的邀請下，我特意用c語言實現了三維向量的叉乘法，因為他在網上找不到關於向量乘法的相關代碼。既然寫出來了，我就想和大家一起分享，在寫的時候對數據結構考慮的不多，重在實現運算過程，有好的建議的朋友希望踴躍留言。
 
  ///////////////////////////////////main.c/////////////////////////////////////////
#include <stdio.h>
#include "vector_multiply.h"
int main(void)
{
 float a[3] = {1.0, 2.0, 3.0}, b[3] = {7.0, 8.0, 9.0}, c[3] = {0};
 
 vector_multiply( a,  b, c);
 printf("%.1f:%.1f:%.1f\n", c[0], c[1], c[2]);
 
 return(0);
}
 
///////////////////////////////////vector_multiply.c///////////////////////////////
void vector_multiply(float a[], float b[], float c[]) {
 int i, j, k, n = 0;
/*according the the formula we can get these marks*/
 int d[3][4] = { { 0, 1, -1, 0 }, { 0, -1, 1, 0 }, { 0, 1, -1, 0 } };
 
 for (k = 0; k < 3; k++) {
  for (i = 0; i < 3; i++) {
   for (j = 0; j < 3; j++)
    if (k != i && k != j)
     c[k] += a[i] * b[j] * d[k][n++];
  }
  n = 0;
 }
 return;
}
 
////////////////////////////////vector_multiply.h//////////////////////////////////
#ifndef VECTOR_MULTIPLY_H_
#define VECTOR_MULTIPLY_H_
void vector_multiply(float a[], float b[], float c[]);
#endif /*VECTOR_MULTIPLY_H_*/
///////////////////////////////////////////////////////////////////////////////////
gcc main.c vector_multiply.c -o  VectorOperations
 
 
最後再留下向量叉乘的公式：
向量：u=(u1,u2,u3) v=(v1,v2,v3)
叉積公式：u x v = | i    j     k|
                 |a1    b1   c1|
                 |a2    b2   c2|
               = { u2v3-v2u3 , u3v1-v3u1 , u1v2-u2v1 }