зеркало из https://github.com/mozilla/gecko-dev.git
420 строки
14 KiB
C++
420 строки
14 KiB
C++
/* -*- Mode: C++; tab-width: 20; indent-tabs-mode: nil; c-basic-offset: 4 -*-
|
|
* ***** BEGIN LICENSE BLOCK *****
|
|
* Version: MPL 1.1/GPL 2.0/LGPL 2.1
|
|
*
|
|
* The contents of this file are subject to the Mozilla Public License Version
|
|
* 1.1 (the "License"); you may not use this file except in compliance with
|
|
* the License. You may obtain a copy of the License at
|
|
* http://www.mozilla.org/MPL/
|
|
*
|
|
* Software distributed under the License is distributed on an "AS IS" basis,
|
|
* WITHOUT WARRANTY OF ANY KIND, either express or implied. See the License
|
|
* for the specific language governing rights and limitations under the
|
|
* License.
|
|
*
|
|
* The Original Code is Oracle Corporation code.
|
|
*
|
|
* The Initial Developer of the Original Code is Oracle Corporation.
|
|
* Portions created by the Initial Developer are Copyright (C) 2005
|
|
* the Initial Developer. All Rights Reserved.
|
|
*
|
|
* Contributor(s):
|
|
* Bas Schouten <bschouten@mozilla.com>
|
|
* Matt Woodrow <mwoodrow@mozilla.com>
|
|
*
|
|
* Alternatively, the contents of this file may be used under the terms of
|
|
* either the GNU General Public License Version 2 or later (the "GPL"), or
|
|
* the GNU Lesser General Public License Version 2.1 or later (the "LGPL"),
|
|
* in which case the provisions of the GPL or the LGPL are applicable instead
|
|
* of those above. If you wish to allow use of your version of this file only
|
|
* under the terms of either the GPL or the LGPL, and not to allow others to
|
|
* use your version of this file under the terms of the MPL, indicate your
|
|
* decision by deleting the provisions above and replace them with the notice
|
|
* and other provisions required by the GPL or the LGPL. If you do not delete
|
|
* the provisions above, a recipient may use your version of this file under
|
|
* the terms of any one of the MPL, the GPL or the LGPL.
|
|
*
|
|
* ***** END LICENSE BLOCK ***** */
|
|
|
|
#include "gfxMatrix.h"
|
|
#include "gfx3DMatrix.h"
|
|
#include <math.h>
|
|
#include <algorithm>
|
|
using namespace std;
|
|
|
|
gfx3DMatrix::gfx3DMatrix(void)
|
|
{
|
|
_11 = _22 = _33 = _44 = 1.0f;
|
|
_12 = _13 = _14 = 0.0f;
|
|
_21 = _23 = _24 = 0.0f;
|
|
_31 = _32 = _34 = 0.0f;
|
|
_41 = _42 = _43 = 0.0f;
|
|
}
|
|
|
|
gfx3DMatrix
|
|
gfx3DMatrix::operator*(const gfx3DMatrix &aMatrix) const
|
|
{
|
|
gfx3DMatrix matrix;
|
|
|
|
matrix._11 = _11 * aMatrix._11 + _12 * aMatrix._21 + _13 * aMatrix._31 + _14 * aMatrix._41;
|
|
matrix._21 = _21 * aMatrix._11 + _22 * aMatrix._21 + _23 * aMatrix._31 + _24 * aMatrix._41;
|
|
matrix._31 = _31 * aMatrix._11 + _32 * aMatrix._21 + _33 * aMatrix._31 + _34 * aMatrix._41;
|
|
matrix._41 = _41 * aMatrix._11 + _42 * aMatrix._21 + _43 * aMatrix._31 + _44 * aMatrix._41;
|
|
matrix._12 = _11 * aMatrix._12 + _12 * aMatrix._22 + _13 * aMatrix._32 + _14 * aMatrix._42;
|
|
matrix._22 = _21 * aMatrix._12 + _22 * aMatrix._22 + _23 * aMatrix._32 + _24 * aMatrix._42;
|
|
matrix._32 = _31 * aMatrix._12 + _32 * aMatrix._22 + _33 * aMatrix._32 + _34 * aMatrix._42;
|
|
matrix._42 = _41 * aMatrix._12 + _42 * aMatrix._22 + _43 * aMatrix._32 + _44 * aMatrix._42;
|
|
matrix._13 = _11 * aMatrix._13 + _12 * aMatrix._23 + _13 * aMatrix._33 + _14 * aMatrix._43;
|
|
matrix._23 = _21 * aMatrix._13 + _22 * aMatrix._23 + _23 * aMatrix._33 + _24 * aMatrix._43;
|
|
matrix._33 = _31 * aMatrix._13 + _32 * aMatrix._23 + _33 * aMatrix._33 + _34 * aMatrix._43;
|
|
matrix._43 = _41 * aMatrix._13 + _42 * aMatrix._23 + _43 * aMatrix._33 + _44 * aMatrix._43;
|
|
matrix._14 = _11 * aMatrix._14 + _12 * aMatrix._24 + _13 * aMatrix._34 + _14 * aMatrix._44;
|
|
matrix._24 = _21 * aMatrix._14 + _22 * aMatrix._24 + _23 * aMatrix._34 + _24 * aMatrix._44;
|
|
matrix._34 = _31 * aMatrix._14 + _32 * aMatrix._24 + _33 * aMatrix._34 + _34 * aMatrix._44;
|
|
matrix._44 = _41 * aMatrix._14 + _42 * aMatrix._24 + _43 * aMatrix._34 + _44 * aMatrix._44;
|
|
|
|
return matrix;
|
|
}
|
|
|
|
gfx3DMatrix&
|
|
gfx3DMatrix::operator*=(const gfx3DMatrix &aMatrix)
|
|
{
|
|
return *this = *this * aMatrix;
|
|
}
|
|
|
|
bool
|
|
gfx3DMatrix::operator==(const gfx3DMatrix& o) const
|
|
{
|
|
// XXX would be nice to memcmp here, but that breaks IEEE 754 semantics
|
|
return _11 == o._11 && _12 == o._12 && _13 == o._13 && _14 == o._14 &&
|
|
_21 == o._21 && _22 == o._22 && _23 == o._23 && _24 == o._24 &&
|
|
_31 == o._31 && _32 == o._32 && _33 == o._33 && _34 == o._34 &&
|
|
_41 == o._41 && _42 == o._42 && _43 == o._43 && _44 == o._44;
|
|
}
|
|
|
|
gfx3DMatrix&
|
|
gfx3DMatrix::operator/=(const gfxFloat scalar)
|
|
{
|
|
_11 /= scalar;
|
|
_12 /= scalar;
|
|
_13 /= scalar;
|
|
_14 /= scalar;
|
|
_21 /= scalar;
|
|
_22 /= scalar;
|
|
_23 /= scalar;
|
|
_24 /= scalar;
|
|
_31 /= scalar;
|
|
_32 /= scalar;
|
|
_33 /= scalar;
|
|
_34 /= scalar;
|
|
_41 /= scalar;
|
|
_42 /= scalar;
|
|
_43 /= scalar;
|
|
_44 /= scalar;
|
|
return *this;
|
|
}
|
|
|
|
gfx3DMatrix
|
|
gfx3DMatrix::From2D(const gfxMatrix &aMatrix)
|
|
{
|
|
gfx3DMatrix matrix;
|
|
matrix._11 = (float)aMatrix.xx;
|
|
matrix._12 = (float)aMatrix.yx;
|
|
matrix._21 = (float)aMatrix.xy;
|
|
matrix._22 = (float)aMatrix.yy;
|
|
matrix._41 = (float)aMatrix.x0;
|
|
matrix._42 = (float)aMatrix.y0;
|
|
return matrix;
|
|
}
|
|
|
|
PRBool
|
|
gfx3DMatrix::IsIdentity() const
|
|
{
|
|
return _11 == 1.0f && _12 == 0.0f && _13 == 0.0f && _14 == 0.0f &&
|
|
_21 == 0.0f && _22 == 1.0f && _23 == 0.0f && _24 == 0.0f &&
|
|
_31 == 0.0f && _32 == 0.0f && _33 == 1.0f && _34 == 0.0f &&
|
|
_41 == 0.0f && _42 == 0.0f && _43 == 0.0f && _44 == 1.0f;
|
|
}
|
|
|
|
gfx3DMatrix
|
|
gfx3DMatrix::Translation(float aX, float aY, float aZ)
|
|
{
|
|
gfx3DMatrix matrix;
|
|
|
|
matrix._41 = aX;
|
|
matrix._42 = aY;
|
|
matrix._43 = aZ;
|
|
return matrix;
|
|
}
|
|
|
|
gfx3DMatrix
|
|
gfx3DMatrix::Translation(const gfxPoint3D& aPoint)
|
|
{
|
|
gfx3DMatrix matrix;
|
|
|
|
matrix._41 = aPoint.x;
|
|
matrix._42 = aPoint.y;
|
|
matrix._43 = aPoint.z;
|
|
return matrix;
|
|
}
|
|
|
|
gfx3DMatrix
|
|
gfx3DMatrix::Scale(float aFactor)
|
|
{
|
|
gfx3DMatrix matrix;
|
|
|
|
matrix._11 = matrix._22 = matrix._33 = aFactor;
|
|
return matrix;
|
|
}
|
|
|
|
gfx3DMatrix
|
|
gfx3DMatrix::Scale(float aX, float aY, float aZ)
|
|
{
|
|
gfx3DMatrix matrix;
|
|
|
|
matrix._11 = aX;
|
|
matrix._22 = aY;
|
|
matrix._33 = aZ;
|
|
|
|
return matrix;
|
|
}
|
|
|
|
gfxFloat
|
|
gfx3DMatrix::Determinant() const
|
|
{
|
|
return _14 * _23 * _32 * _41
|
|
- _13 * _24 * _32 * _41
|
|
- _14 * _22 * _33 * _41
|
|
+ _12 * _24 * _33 * _41
|
|
+ _13 * _22 * _34 * _41
|
|
- _12 * _23 * _34 * _41
|
|
- _14 * _23 * _31 * _42
|
|
+ _13 * _24 * _31 * _42
|
|
+ _14 * _21 * _33 * _42
|
|
- _11 * _24 * _33 * _42
|
|
- _13 * _21 * _34 * _42
|
|
+ _11 * _23 * _34 * _42
|
|
+ _14 * _22 * _31 * _43
|
|
- _12 * _24 * _31 * _43
|
|
- _14 * _21 * _32 * _43
|
|
+ _11 * _24 * _32 * _43
|
|
+ _12 * _21 * _34 * _43
|
|
- _11 * _22 * _34 * _43
|
|
- _13 * _22 * _31 * _44
|
|
+ _12 * _23 * _31 * _44
|
|
+ _13 * _21 * _32 * _44
|
|
- _11 * _23 * _32 * _44
|
|
- _12 * _21 * _33 * _44
|
|
+ _11 * _22 * _33 * _44;
|
|
}
|
|
|
|
PRBool
|
|
gfx3DMatrix::IsSingular() const
|
|
{
|
|
return Determinant() == 0.0;
|
|
}
|
|
|
|
gfx3DMatrix&
|
|
gfx3DMatrix::Invert()
|
|
{
|
|
gfxFloat det = Determinant();
|
|
if (det == 0.0) {
|
|
return *this;
|
|
}
|
|
|
|
gfx3DMatrix temp = *this;
|
|
|
|
_11 = temp._23*temp._34*temp._42 - temp._24*temp._33*temp._42
|
|
+ temp._24*temp._32*temp._43 - temp._22*temp._34*temp._43
|
|
- temp._23*temp._32*temp._44 + temp._22*temp._33*temp._44;
|
|
_12 = temp._14*temp._33*temp._42 - temp._13*temp._34*temp._42
|
|
- temp._14*temp._32*temp._43 + temp._12*temp._34*temp._43
|
|
+ temp._13*temp._32*temp._44 - temp._12*temp._33*temp._44;
|
|
_13 = temp._13*temp._24*temp._42 - temp._14*temp._23*temp._42
|
|
+ temp._14*temp._22*temp._43 - temp._12*temp._24*temp._43
|
|
- temp._13*temp._22*temp._44 + temp._12*temp._23*temp._44;
|
|
_14 = temp._14*temp._23*temp._32 - temp._13*temp._24*temp._32
|
|
- temp._14*temp._22*temp._33 + temp._12*temp._24*temp._33
|
|
+ temp._13*temp._22*temp._34 - temp._12*temp._23*temp._34;
|
|
_21 = temp._24*temp._33*temp._41 - temp._23*temp._34*temp._41
|
|
- temp._24*temp._31*temp._43 + temp._21*temp._34*temp._43
|
|
+ temp._23*temp._31*temp._44 - temp._21*temp._33*temp._44;
|
|
_22 = temp._13*temp._34*temp._41 - temp._14*temp._33*temp._41
|
|
+ temp._14*temp._31*temp._43 - temp._11*temp._34*temp._43
|
|
- temp._13*temp._31*temp._44 + temp._11*temp._33*temp._44;
|
|
_23 = temp._14*temp._23*temp._41 - temp._13*temp._24*temp._41
|
|
- temp._14*temp._21*temp._43 + temp._11*temp._24*temp._43
|
|
+ temp._13*temp._21*temp._44 - temp._11*temp._23*temp._44;
|
|
_24 = temp._13*temp._24*temp._31 - temp._14*temp._23*temp._31
|
|
+ temp._14*temp._21*temp._33 - temp._11*temp._24*temp._33
|
|
- temp._13*temp._21*temp._34 + temp._11*temp._23*temp._34;
|
|
_31 = temp._22*temp._34*temp._41 - temp._24*temp._32*temp._41
|
|
+ temp._24*temp._31*temp._42 - temp._21*temp._34*temp._42
|
|
- temp._22*temp._31*temp._44 + temp._21*temp._32*temp._44;
|
|
_32 = temp._14*temp._32*temp._41 - temp._12*temp._34*temp._41
|
|
- temp._14*temp._31*temp._42 + temp._11*temp._34*temp._42
|
|
+ temp._12*temp._31*temp._44 - temp._11*temp._32*temp._44;
|
|
_33 = temp._12*temp._24*temp._41 - temp._14*temp._22*temp._41
|
|
+ temp._14*temp._21*temp._42 - temp._11*temp._24*temp._42
|
|
- temp._12*temp._21*temp._44 + temp._11*temp._22*temp._44;
|
|
_34 = temp._14*temp._22*temp._31 - temp._12*temp._24*temp._31
|
|
- temp._14*temp._21*temp._32 + temp._11*temp._24*temp._32
|
|
+ temp._12*temp._21*temp._34 - temp._11*temp._22*temp._34;
|
|
_41 = temp._23*temp._32*temp._41 - temp._22*temp._33*temp._41
|
|
- temp._23*temp._31*temp._42 + temp._21*temp._33*temp._42
|
|
+ temp._22*temp._31*temp._43 - temp._21*temp._32*temp._43;
|
|
_42 = temp._12*temp._33*temp._41 - temp._13*temp._32*temp._41
|
|
+ temp._13*temp._31*temp._42 - temp._11*temp._33*temp._42
|
|
- temp._12*temp._31*temp._43 + temp._11*temp._32*temp._43;
|
|
_43 = temp._13*temp._22*temp._41 - temp._12*temp._23*temp._41
|
|
- temp._13*temp._21*temp._42 + temp._11*temp._23*temp._42
|
|
+ temp._12*temp._21*temp._43 - temp._11*temp._22*temp._43;
|
|
_44 = temp._12*temp._23*temp._31 - temp._13*temp._22*temp._31
|
|
+ temp._13*temp._21*temp._32 - temp._11*temp._23*temp._32
|
|
- temp._12*temp._21*temp._33 + temp._11*temp._22*temp._33;
|
|
|
|
*this /= det;
|
|
return *this;
|
|
}
|
|
|
|
gfxPoint
|
|
gfx3DMatrix::Transform(const gfxPoint& point) const
|
|
{
|
|
gfxPoint3D vec3d(point.x, point.y, 0);
|
|
vec3d = Transform3D(vec3d);
|
|
return gfxPoint(vec3d.x, vec3d.y);
|
|
}
|
|
|
|
gfxPoint3D
|
|
gfx3DMatrix::Transform3D(const gfxPoint3D& point) const
|
|
{
|
|
gfxFloat x = point.x * _11 + point.y * _21 + point.z * _31 + _41;
|
|
gfxFloat y = point.x * _12 + point.y * _22 + point.z * _32 + _42;
|
|
gfxFloat z = point.x * _13 + point.y * _23 + point.z * _33 + _43;
|
|
gfxFloat w = point.x * _14 + point.y * _24 + point.z * _34 + _44;
|
|
|
|
x /= w;
|
|
y /= w;
|
|
z /= w;
|
|
|
|
return gfxPoint3D(x, y, z);
|
|
}
|
|
|
|
gfxRect
|
|
gfx3DMatrix::TransformBounds(const gfxRect& rect) const
|
|
{
|
|
gfxPoint points[4];
|
|
|
|
points[0] = Transform(rect.TopLeft());
|
|
points[1] = Transform(gfxPoint(rect.X() + rect.Width(), rect.Y()));
|
|
points[2] = Transform(gfxPoint(rect.X(), rect.Y() + rect.Height()));
|
|
points[3] = Transform(gfxPoint(rect.X() + rect.Width(),
|
|
rect.Y() + rect.Height()));
|
|
|
|
gfxFloat min_x, max_x;
|
|
gfxFloat min_y, max_y;
|
|
|
|
min_x = max_x = points[0].x;
|
|
min_y = max_y = points[0].y;
|
|
|
|
for (int i=1; i<4; i++) {
|
|
min_x = min(points[i].x, min_x);
|
|
max_x = max(points[i].x, max_x);
|
|
min_y = min(points[i].y, min_y);
|
|
max_y = max(points[i].y, max_y);
|
|
}
|
|
|
|
return gfxRect(min_x, min_y, max_x - min_x, max_y - min_y);
|
|
}
|
|
|
|
PRBool
|
|
gfx3DMatrix::Is2D(gfxMatrix* aMatrix) const
|
|
{
|
|
if (_13 != 0.0f || _14 != 0.0f ||
|
|
_23 != 0.0f || _24 != 0.0f ||
|
|
_31 != 0.0f || _32 != 0.0f || _33 != 1.0f || _34 != 0.0f ||
|
|
_43 != 0.0f || _44 != 1.0f) {
|
|
return PR_FALSE;
|
|
}
|
|
if (aMatrix) {
|
|
aMatrix->xx = _11;
|
|
aMatrix->yx = _12;
|
|
aMatrix->xy = _21;
|
|
aMatrix->yy = _22;
|
|
aMatrix->x0 = _41;
|
|
aMatrix->y0 = _42;
|
|
}
|
|
return PR_TRUE;
|
|
}
|
|
|
|
gfxPoint gfx3DMatrix::ProjectPoint(const gfxPoint& aPoint) const
|
|
{
|
|
// Define a ray of the form P + Ut where t is a real number
|
|
// w is assumed to always be 1 when transforming 3d points with our
|
|
// 4x4 matrix.
|
|
// p is our click point, q is another point on the same ray.
|
|
//
|
|
// Note: since the transformation is a general projective transformation and is not
|
|
// necessarily affine, we can't just take a unit vector u, back-transform it, and use
|
|
// it as unit vector on the back-transformed ray. Instead, we really must take two points
|
|
// on the ray and back-transform them.
|
|
gfxPoint3D p(aPoint.x, aPoint.y, 0);
|
|
gfxPoint3D q(aPoint.x, aPoint.y, 1);
|
|
|
|
// Back transform the vectors (using w = 1) and normalize
|
|
// back into 3d vectors by dividing by the w component.
|
|
gfxPoint3D pback = Transform3D(p);
|
|
gfxPoint3D qback = Transform3D(q);
|
|
gfxPoint3D uback = qback - pback;
|
|
|
|
// Find the point where the back transformed line intersects z=0
|
|
// and find t.
|
|
|
|
float t = -pback.z / uback.z;
|
|
|
|
gfxPoint result(pback.x + t*uback.x, pback.y + t*uback.y);
|
|
|
|
return result;
|
|
}
|
|
|
|
gfxRect gfx3DMatrix::ProjectRectBounds(const gfxRect& aRect) const
|
|
{
|
|
gfxPoint points[4];
|
|
|
|
points[0] = ProjectPoint(aRect.TopLeft());
|
|
points[1] = ProjectPoint(gfxPoint(aRect.X() + aRect.Width(), aRect.Y()));
|
|
points[2] = ProjectPoint(gfxPoint(aRect.X(), aRect.Y() + aRect.Height()));
|
|
points[3] = ProjectPoint(gfxPoint(aRect.X() + aRect.Width(),
|
|
aRect.Y() + aRect.Height()));
|
|
|
|
gfxFloat min_x, max_x;
|
|
gfxFloat min_y, max_y;
|
|
|
|
min_x = max_x = points[0].x;
|
|
min_y = max_y = points[0].y;
|
|
|
|
for (int i=1; i<4; i++) {
|
|
min_x = min(points[i].x, min_x);
|
|
max_x = max(points[i].x, max_x);
|
|
min_y = min(points[i].y, min_y);
|
|
max_y = max(points[i].y, max_y);
|
|
}
|
|
|
|
return gfxRect(min_x, min_y, max_x - min_x, max_y - min_y);
|
|
}
|
|
|
|
gfxPoint3D gfx3DMatrix::GetNormalVector() const
|
|
{
|
|
// Define a plane in transformed space as the transformations
|
|
// of 3 points on the z=0 screen plane.
|
|
gfxPoint3D a = Transform3D(gfxPoint3D(0, 0, 0));
|
|
gfxPoint3D b = Transform3D(gfxPoint3D(0, 1, 0));
|
|
gfxPoint3D c = Transform3D(gfxPoint3D(1, 0, 0));
|
|
|
|
// Convert to two vectors on the surface of the plane.
|
|
gfxPoint3D ab = b - a;
|
|
gfxPoint3D ac = c - a;
|
|
|
|
return ac.CrossProduct(ab);
|
|
}
|