计算射线与平面的交点

#include "stdafx.h"#include<iostream>using namespace std;struct Point3{float x;float y;float z;};struct Vector{float x;float y;float z;};struct Ray{       //一点，和一个方向向量（两点求差）确定一条射线，Point3 p0;Vector u;};struct Plane{     //一点，和一个法向量确定一个平面Point3 p1;Vector n;};float operator*(Vector v1,Vector v2){float x = v1.x*v2.x+v1.y*v2.y+v1.z*v2.z;return x;}float operator*(Vector v1, Point3 v2){float x = v1.x*v2.x + v1.y*v2.y + v1.z*v2.z;return x;}Point3 operator*(float t, Vector v){Point3 p;p = { t*v.x, t*v.y, t*v.z };return p;}Vector operator-(Vector v1, Vector v2){Vector v3 = { v2.x-v1.x, v2.y-v1.y, v2.z-v1.z };return v3;}Point3 operator+(Point3 v1, Point3 v2){Point3 v3 = { v1.x + v2.x, v1.y + v2.y, v1.z + v2.z };return v3;}Point3 operator-(Point3 v1, Point3 v2){Point3 v3 = { v2.x - v1.x, v2.y - v1.y, v2.z - v1.z };return v3;}Point3 RayPlaneIntersection(Ray ray,Plane plane){Point3 p;float t;t = (plane.n*plane.p1 - plane.n*ray.p0) / (plane.n*ray.u);p = ray.p0 + t*ray.u;return p;}int _tmain(int argc, _TCHAR* argv[]){Ray r1;Plane p1;r1 = { { 0, 0, 0 }, { 0, 1, 0 } };p1 = { { 0, 3, 0 }, { 0, -1, 0 } };Point3 pt = RayPlaneIntersection(r1, p1);cout << pt.x << " " << pt.y << " " << pt.z << endl;return 0;}