开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > Windows编程 > 查看源码
DUMB3D.CPP
资源名称：MSDN_VC98.zip　[点击查看]
上传用户：bangxh
上传日期：2007-01-31
资源大小：42235k
文件大小：8k
源码类别：

Windows编程

开发平台：

Visual C++

DUMB3D.CPP：源码内容
  1. /**************************************************************************
  2.   dumb3d.cpp - A simple linear algebra library for 3D.
  4.  **************************************************************************/
  5. /**************************************************************************
  7.     (C) Copyright 1995-1996 Microsoft Corp.  All rights reserved.
  9.     You have a royalty-free right to use, modify, reproduce and 
  10.     distribute the Sample Files (and/or any modified version) in 
  11.     any way you find useful, provided that you agree that 
  12.     Microsoft has no warranty obligations or liability for any 
  13.     Sample Application Files which are modified. 
  15.  **************************************************************************/
  17. #include<memory.h>
  18. #include<math.h>
  19. #include"dumb3d.h"
  21. /*--------------------------------------------------------------------------
  23. matrix multiplication.
  25. */
  27. matrix_4x4 operator*( matrix_4x4 const &Multiplicand,
  28.     matrix_4x4 const &Multiplier )
  29. {
  30.   matrix_4x4 ReturnMatrix;
  32.   for(int i = 0;i < 4;i++)
  33.   {
  34.     for(int j = 0;j < 4;j++)
  35.     {
  36.       real Value = 0;
  37.       
  38.       for(int k = 0;k < 4;k++)
  39.       {
  40.         Value += Multiplicand.GetElement(i,k) *
  41.             Multiplier.GetElement(k,j);
  42.       }
  44.       ReturnMatrix.SetElement(i,j,Value);
  45.     }
  46.   }
  48.   return ReturnMatrix;
  49. }
  51. vector_4 operator*( matrix_4x4 const &Multiplicand,
  52.     vector_4 const &Multiplier )
  53. {
  54.   vector_4 ReturnPoint;
  56.   for(int i = 0;i < 4;i++)
  57.   {
  58.     real Value = 0;
  59.     
  60.     for(int k = 0;k < 4;k++)
  61.     {
  62.       Value += Multiplicand.GetElement(i,k) *
  63.           Multiplier.GetElement(k);
  64.     }
  66.     ReturnPoint.SetElement(i,Value);
  67.   }
  69.   return ReturnPoint;
  70. }
  72. point_4 operator*( matrix_4x4 const &Multiplicand,
  73.     point_4 const &Multiplier )
  74. {
  75.   point_4 ReturnPoint;
  77.   for(int i = 0;i < 4;i++)
  78.   {
  79.     real Value = 0;
  80.     
  81.     for(int k = 0;k < 4;k++)
  82.     {
  83.       Value += Multiplicand.GetElement(i,k) *
  84.           Multiplier.GetElement(k);
  85.     }
  87.     ReturnPoint.SetElement(i,Value);
  88.   }
  90.   return ReturnPoint;
  91. }
  94. /*--------------------------------------------------------------------------
  96. constructor.
  98. */
  100. matrix_4x4::matrix_4x4( void )
  101. {
  102.   for(int Counter = 0;Counter < 16;Counter++)
  103.   {
  104.     aElements[0][Counter] = 0;
  105.   }
  107.   aElements[0][0] = aElements[1][1] = aElements[2][2] = aElements[3][3] = 1;
  108. }
  111. /*--------------------------------------------------------------------------
  113. Rotations.
  115. */
  117. matrix_4x4 &matrix_4x4::ConcatenateXRotation( real Degrees )
  118. {
  119.   real Temp01, Temp11, Temp21, Temp31;
  120.   real Temp02, Temp12, Temp22, Temp32;
  122.   real Radians = (Degrees/360) * M_PI * 2;
  124.   real Sin = sin(Radians), Cos = cos(Radians);
  126.   Temp01 = aElements[0][1] * Cos + aElements[0][2] * Sin;
  127.   Temp11 = aElements[1][1] * Cos + aElements[1][2] * Sin;
  128.   Temp21 = aElements[2][1] * Cos + aElements[2][2] * Sin;
  129.   Temp31 = aElements[3][1] * Cos + aElements[3][2] * Sin;
  131.   Temp02 = aElements[0][1] * -Sin + aElements[0][2] * Cos;
  132.   Temp12 = aElements[1][1] * -Sin + aElements[1][2] * Cos;
  133.   Temp22 = aElements[2][1] * -Sin + aElements[2][2] * Cos;
  134.   Temp32 = aElements[3][1] * -Sin + aElements[3][2] * Cos;
  136.   aElements[0][1] = Temp01;
  137.   aElements[1][1] = Temp11;
  138.   aElements[2][1] = Temp21;
  139.   aElements[3][1] = Temp31;
  140.   aElements[0][2] = Temp02;
  141.   aElements[1][2] = Temp12;
  142.   aElements[2][2] = Temp22;
  143.   aElements[3][2] = Temp32;
  145.   return *this;
  146. }
  148. matrix_4x4 &matrix_4x4::ConcatenateYRotation( real Degrees )
  149. {
  150.   real Temp00, Temp10, Temp20, Temp30;
  151.   real Temp02, Temp12, Temp22, Temp32;
  153.   real Radians = (Degrees/360) * M_PI * 2;
  155.   real Sin = sin(Radians), Cos = cos(Radians);
  157.   Temp00 = aElements[0][0] * Cos + aElements[0][2] * -Sin;
  158.   Temp10 = aElements[1][0] * Cos + aElements[1][2] * -Sin;
  159.   Temp20 = aElements[2][0] * Cos + aElements[2][2] * -Sin;
  160.   Temp30 = aElements[3][0] * Cos + aElements[3][2] * -Sin;
  162.   Temp02 = aElements[0][0] * Sin + aElements[0][2] * Cos;
  163.   Temp12 = aElements[1][0] * Sin + aElements[1][2] * Cos;
  164.   Temp22 = aElements[2][0] * Sin + aElements[2][2] * Cos;
  165.   Temp32 = aElements[3][0] * Sin + aElements[3][2] * Cos;
  167.   aElements[0][0] = Temp00;
  168.   aElements[1][0] = Temp10;
  169.   aElements[2][0] = Temp20;
  170.   aElements[3][0] = Temp30;
  171.   aElements[0][2] = Temp02;
  172.   aElements[1][2] = Temp12;
  173.   aElements[2][2] = Temp22;
  174.   aElements[3][2] = Temp32;
  176.   return *this;
  177. }
  179. matrix_4x4 &matrix_4x4::ConcatenateZRotation( real Degrees )
  180. {
  181.   real Temp00, Temp10, Temp20, Temp30;
  182.   real Temp01, Temp11, Temp21, Temp31;
  184.   real Radians = (Degrees/360) * M_PI * 2;
  186.   real Sin = sin(Radians), Cos = cos(Radians);
  188.   Temp00 = aElements[0][0] * Cos + aElements[0][1] * Sin;
  189.   Temp10 = aElements[1][0] * Cos + aElements[1][1] * Sin;
  190.   Temp20 = aElements[2][0] * Cos + aElements[2][1] * Sin;
  191.   Temp30 = aElements[3][0] * Cos + aElements[3][1] * Sin;
  193.   Temp01 = aElements[0][0] * -Sin + aElements[0][1] * Cos;
  194.   Temp11 = aElements[1][0] * -Sin + aElements[1][1] * Cos;
  195.   Temp21 = aElements[2][0] * -Sin + aElements[2][1] * Cos;
  196.   Temp31 = aElements[3][0] * -Sin + aElements[3][1] * Cos;
  198.   aElements[0][0] = Temp00;
  199.   aElements[1][0] = Temp10;
  200.   aElements[2][0] = Temp20;
  201.   aElements[3][0] = Temp30;
  202.   aElements[0][1] = Temp01;
  203.   aElements[1][1] = Temp11;
  204.   aElements[2][1] = Temp21;
  205.   aElements[3][1] = Temp31;
  207.   return *this;
  208. }
  210. /*--------------------------------------------------------------------------
  212. Translations.
  214. */
  216. matrix_4x4 &matrix_4x4::ConcatenateXTranslation( real Distance )
  217. {
  218.   aElements[0][3] = aElements[0][0] * Distance + aElements[0][3];
  219.   aElements[1][3] = aElements[1][0] * Distance + aElements[1][3];
  220.   aElements[2][3] = aElements[2][0] * Distance + aElements[2][3];
  221.   aElements[3][3] = aElements[3][0] * Distance + aElements[3][3];
  223.   return *this;
  224. }
  226. matrix_4x4 &matrix_4x4::ConcatenateYTranslation( real Distance )
  227. {
  228.   aElements[0][3] = aElements[0][1] * Distance + aElements[0][3];
  229.   aElements[1][3] = aElements[1][1] * Distance + aElements[1][3];
  230.   aElements[2][3] = aElements[2][1] * Distance + aElements[2][3];
  231.   aElements[3][3] = aElements[3][1] * Distance + aElements[3][3];
  233.   return *this;
  234. }
  236. matrix_4x4 &matrix_4x4::ConcatenateZTranslation( real Distance )
  237. {
  238.   aElements[0][3] = aElements[0][2] * Distance + aElements[0][3];
  239.   aElements[1][3] = aElements[1][2] * Distance + aElements[1][3];
  240.   aElements[2][3] = aElements[2][2] * Distance + aElements[2][3];
  241.   aElements[3][3] = aElements[3][2] * Distance + aElements[3][3];
  243.   return *this;
  244. }
  246. /*--------------------------------------------------------------------------
  248. vector normalize.
  250. */
  252. vector_4 &vector_4::Normalize( void )
  253. {
  254.   real Length = sqrt(GetX()*GetX() + GetY()*GetY() + GetZ()*GetZ());
  256.   SetX(GetX() / Length);
  257.   SetY(GetY() / Length);
  258.   SetZ(GetZ() / Length);
  260.   return *this;
  261. }
  262.   
  263. /*--------------------------------------------------------------------------
  265. view transform ctor.
  267. */
  269. view_transform::view_transform( point_4 const &Viewpoint,
  270.   vector_4 const &ViewDirection, vector_4 const &Up )
  271. {
  272.   // translate the viewpoint to the origin
  274.   this->ConcatenateXTranslation(-Viewpoint.GetX());
  275.   this->ConcatenateYTranslation(-Viewpoint.GetY());
  276.   this->ConcatenateZTranslation(-Viewpoint.GetZ());
  278.   // get view vectors set up
  280.   vector_4 Right = -CrossProduct(ViewDirection,Up);
  281.   vector_4 ReallyUp = CrossProduct(Right,ViewDirection);
  282.     
  283.   matrix_4x4 LookDownZ;
  285.   for(int Counter = 0;Counter < 3;Counter++)
  286.   {
  287.     LookDownZ.SetElement(0,Counter,Right.GetElement(Counter));
  288.   }
  290.   for(Counter = 0;Counter < 3;Counter++)
  291.   {
  292.     LookDownZ.SetElement(1,Counter,ReallyUp.GetElement(Counter));
  293.   }
  295.   for(Counter = 0;Counter < 3;Counter++)
  296.   {
  297.     LookDownZ.SetElement(2,Counter,ViewDirection.GetElement(Counter));
  298.   }
  300.   this->matrix_4x4::operator=(LookDownZ * *this);
  301. }