由于研究需要,经常要求向量积,自己写了个
double *cross_product(double *vector_a, double *vector_b)
{
        double result[3];
        result[0] = (*(vector_a+1))*(*(vector_b+2)) - (*(vector_a+2))*(*(vector_b+1));//X分量
        result[1] = (*(vector_a+3))*(*vector_b) - (*vector_a)*(*(vector_b+2));//Y分量
        result[2] = (*vector_a)*(*(vector_b+1)) - (*(vector_a+1))*(*vector_b);//Z分量
        return result;
}
然后在所需要的地方调用,返回的result值赋给一个数组c[3],如
c = cross_product(a, b);
其中a,b和c一样都是大小为3的数组

总感觉这样写不好,希望大家能给出更好的代码,谢谢!