Matrix3.cpp 2.0 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859
  1. //
  2. // Urho3D Engine
  3. // Copyright (c) 2008-2011 Lasse Öörni
  4. //
  5. // Permission is hereby granted, free of charge, to any person obtaining a copy
  6. // of this software and associated documentation files (the "Software"), to deal
  7. // in the Software without restriction, including without limitation the rights
  8. // to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
  9. // copies of the Software, and to permit persons to whom the Software is
  10. // furnished to do so, subject to the following conditions:
  11. //
  12. // The above copyright notice and this permission notice shall be included in
  13. // all copies or substantial portions of the Software.
  14. //
  15. // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
  16. // IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
  17. // FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
  18. // AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
  19. // LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
  20. // OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
  21. // THE SOFTWARE.
  22. //
  23. #include "Precompiled.h"
  24. #include "Matrix3.h"
  25. const Matrix3 Matrix3::sZero(
  26. 0.0f, 0.0f, 0.0f,
  27. 0.0f, 0.0f, 0.0f,
  28. 0.0f, 0.0f, 0.0f);
  29. const Matrix3 Matrix3::sIdentity(
  30. 1.0f, 0.0f, 0.0f,
  31. 0.0f, 1.0f, 0.0f,
  32. 0.0f, 0.0f, 1.0f);
  33. Matrix3 Matrix3::getInverse() const
  34. {
  35. float det = m00 * m11 * m22 +
  36. m10 * m21 * m02 +
  37. m20 * m01 * m12 -
  38. m20 * m11 * m02 -
  39. m10 * m01 * m22 -
  40. m00 * m21 * m12;
  41. float invDet = 1.0f / det;
  42. return Matrix3(
  43. (m11 * m22 - m21 * m12) * invDet,
  44. -(m01 * m22 - m21 * m02) * invDet,
  45. (m01 * m12 - m11 * m02) * invDet,
  46. -(m10 * m22 - m20 * m12) * invDet,
  47. (m00 * m22 - m20 * m02) * invDet,
  48. -(m00 * m12 - m10 * m02) * invDet,
  49. (m10 * m21 - m20 * m11) * invDet,
  50. -(m00 * m21 - m20 * m01) * invDet,
  51. (m00 * m11 - m10 * m01) * invDet
  52. );
  53. }