OpenGL 求反射光线 反射向量

//myblbase.cpp //rorger, 2011 // last modified:2011-04-28 #include "myglbase.h" #include "math.h" Point2::Point2() { this->x=0; this->y=0; } Point2::Point2(double x, double y) { this->x=x; this->y=y; } void Point2::set(Point2& p) { this->x=p.x; this->y=p.y; } void Point2::set(double x, double y) { this->x=x; this->y=y; } Point2 Point2::operator /(double rhs) { Point2 result ; result.x=x/rhs; result.y=y/rhs; return result; } Point2 operator +(Point2& lhs,Vector2& rhs) { Point2 result ; result.x = lhs.x+rhs.x; result.y = lhs.y+rhs.y; return result; } Point2 operator +(Point2& lhs,Point2& rhs) { Point2 result ; result.x = lhs.x+rhs.x; result.y = lhs.y+rhs.y; return result; } Point2 operator - (Point2& lhs,Vector2& rhs) { return lhs + ((-1)*rhs); } Vector2 operator - (Point2& lhs,Point2& rhs) { Vector2 result ; result.x = lhs.x-rhs.x; result.y = lhs.y-rhs.y; return result; } Point2 operator* (Point2& lhs,double rhs) { Point2 result ; result.x=lhs.x * rhs ; result.y = lhs.y * rhs ; return result; } Point2 operator* (double lhs,Point2& rhs) { return rhs* lhs; } ostream& operator << (ostream& outputStream,Point2&rhs) { outputStream<<"("<<rhs.x<<","<<rhs.y<<")"; return outputStream ; } //implementation for class Vector2 Vector2::Vector2() { this->x=0; this->y=0; } Vector2::Vector2(double x, double y) { this->x=x; this->y=y; } void Vector2::set(Vector2&p) { this->x=p.x; this->y=p.y; } void Vector2::set(double x, double y) { this->x= x ; this->y=y ; } Vector2 Vector2::operator /(double rhs) { Vector2 result ; result.x = this->x / rhs ; result.y = this->y / rhs ; return result; } Vector2 operator +(Vector2&lhs,Vector2&rhs) { Vector2 result ; result.x = lhs.x + rhs.x ; result.y = lhs.y + rhs.y ; return result; } Vector2 operator - (Vector2&lhs,Vector2& rhs) { Vector2 result ; result.x = lhs.x - rhs.x ; result.y = lhs.y - rhs.y ; return result; } double operator *(Vector2&lhs,Vector2& rhs) { double result ; result = lhs.x * rhs.x + lhs.y * rhs.y ; return result; } Vector2 operator *(Vector2&lhs,double rhs) { Vector2 result ; result.x = lhs.x*rhs ; result.y = lhs.y*rhs ; return result; } Vector2 operator *(double lhs,Vector2& rhs) { return rhs*lhs; } double Vector2::GetLength() { return sqrt(this->x*this->x+this->y*this->y); } void Vector2::Normalized() { double length = GetLength(); this->x = x / length ; this->y = y / length ; } Vector2 Vector2::GetNormalVector() { Vector2 result ; result.x = -1*this->y; result.y = this->x ; return result ; } ostream& operator << (ostream& outputStream,Vector2&rhs) { outputStream<<"("<<rhs.x<<","<<rhs.y<<")"; return outputStream ; } //implementation for Vector3 Vector3::Vector3() { this->x= this->y=this->z=0; } void Vector3::set(Vector3&p) { this->x = p.x ; this->y =p.y ; this->z =p.z; } void Vector3::set(double x, double y, double z) { this->x=x; this->y=y; this->z=z; } double Vector3::GetLength() { return sqrt(this->x*this->x+this->y*this->y+this->z*this->z); } void Vector3::Normalized() { double length=GetLength(); this->x=x/length; this->y=y/length; this->z=z/length; } Vector3 Vector3::GetNormalVector() { throw "NotImplemented Implementation."; } Vector3 operator +(Vector3&lhs,Vector3&rhs) { Vector3 result ; result.x = lhs.x + rhs.x ; result.y = lhs.y + rhs.y ; result.z = lhs.z + rhs.z ; return result; } Vector3 operator - (Vector3&lhs,Vector3& rhs) { Vector3 result ; result.x = lhs.x - rhs.x ; result.y = lhs.y - rhs.y ; result.z = lhs.z - rhs.z ; return result; } double operator *(Vector3&lhs,Vector3& rhs) { double result ; result = lhs.x * rhs.x +lhs.y * rhs.y + lhs.z * rhs.z ; return result; } Vector3 operator *(Vector3&lhs,double rhs) { Vector3 result ; result.x = lhs.x * rhs ; result.y = lhs.y * rhs ; result.z = lhs.z * rhs ; return result; } Vector3 Vector3::operator /(double rhs) { Vector3 result ; result.x = this->x / rhs ; result.y = this->y / rhs ; result.z = this->z / rhs ; return result ; } ostream& operator << (ostream& outputStream,Vector3& rhs) { cout <<"("<< rhs.x<<","<<rhs.y<<","<<rhs.z<<")"; return outputStream; } //implementation for Point3 Point3::Point3() { this->x = this->y = this->z = 0 ; } void Point3::set(Point3& p) { this->x = p.x ; this->y = p.y ; this->z = p.z ; } void Point3::set(double x, double y, double z) { this->x = x ; this->y = y ; this->z = z ; } Point3 Point3::operator /(double rhs) { Point3 result ; result.x = this->x / rhs; result.y = this->y / rhs ; result.z = this->z / rhs ; return result ; } Point3 operator +(Point3& lhs,Vector3& rhs) { Point3 result ; result.x = lhs.x + rhs.x ; result.y = lhs.y + rhs.y ; result.z = lhs.z + rhs.z ; return result ; } Vector3 operator - (Point3& lhs,Point3& rhs) { Vector3 result ; result.x = lhs.x - rhs.x ; result.y = lhs.y - rhs.y ; result.z = lhs.z - rhs.z ; return result ; } Point3 operator* (Point3& lhs,double rhs) { Point3 result ; result.x = lhs.x * rhs ; result.y = lhs.y * rhs; result.z = lhs.z * rhs ; return result ; } ostream& operator << (ostream& outputStream,Point3& rhs) { cout <<"("<< rhs.x<<","<<rhs.y<<","<<rhs.z<<")"; return outputStream; }