计算给定两点的法向量

VECTOR GetVector(float * v1,float * v2,float * v3) 
{
 //计算给定两向量的法向量

 VECTOR vec1;
 vec1.x = v2[0] - v1[0];
 vec1.y = v2[1] - v1[1];
 vec1.z = v2[2] - v1[2];

 VECTOR vec2;
 vec2.x = v3[0]- v1[0];
 vec2.y = v3[1] - v1[1];
 vec2.z = v3[2] - v1[2];

 VECTOR vec;
 vec.x = vec1.y * vec2.z - vec1.z * vec2.y;
 vec.y = vec1.z * vec2.x - vec1.x * vec2.z;
 vec.z = vec1.x * vec2.y - vec1.y * vec2.x;

 double d = sqrt (vec.x * vec.x + vec.y  * vec.y  + vec.z * vec.z);
 vec.x = float(vec.x/d);
 vec.y = float(vec.y/d);
 vec.z = float(vec.z/d);
 return vec;
}