tests/Matrix44Test.cpp - platform/external/skia - Gitiles

 /*
  * Copyright 2011 Google Inc.
  *
  * Use of this source code is governed by a BSD-style license that can be
  * found in the LICENSE file.
  */

 #include "Test.h"
 #include "SkMatrix44.h"

 static bool nearly_equal_double(double a, double b) {
     const double tolerance = 1e-7;
     double diff = a - b;
     if (diff < 0)
         diff = -diff;
     return diff <= tolerance;
 }

 static bool nearly_equal_mscalar(SkMScalar a, SkMScalar b) {
     const SkMScalar tolerance = SK_MScalar1 / 200000;

     return SkTAbs<SkMScalar>(a - b) <= tolerance;
 }

 static bool nearly_equal_scalar(SkScalar a, SkScalar b) {
     // Note that we get more compounded error for multiple operations when
     // SK_SCALAR_IS_FIXED.
 #ifdef SK_SCALAR_IS_FLOAT
     const SkScalar tolerance = SK_Scalar1 / 200000;
 #else
     const SkScalar tolerance = SK_Scalar1 / 1024;
 #endif

     return SkScalarAbs(a - b) <= tolerance;
 }

 template <typename T> void assert16(skiatest::Reporter* reporter, const T data[],
                                     T m0,  T m1,  T m2,  T m3,
                                     T m4,  T m5,  T m6,  T m7,
                                     T m8,  T m9,  T m10, T m11,
                                     T m12, T m13, T m14, T m15) {
     REPORTER_ASSERT(reporter, data[0] == m0);
     REPORTER_ASSERT(reporter, data[1] == m1);
     REPORTER_ASSERT(reporter, data[2] == m2);
     REPORTER_ASSERT(reporter, data[3] == m3);

     REPORTER_ASSERT(reporter, data[4] == m4);
     REPORTER_ASSERT(reporter, data[5] == m5);
     REPORTER_ASSERT(reporter, data[6] == m6);
     REPORTER_ASSERT(reporter, data[7] == m7);

     REPORTER_ASSERT(reporter, data[8] == m8);
     REPORTER_ASSERT(reporter, data[9] == m9);
     REPORTER_ASSERT(reporter, data[10] == m10);
     REPORTER_ASSERT(reporter, data[11] == m11);

     REPORTER_ASSERT(reporter, data[12] == m12);
     REPORTER_ASSERT(reporter, data[13] == m13);
     REPORTER_ASSERT(reporter, data[14] == m14);
     REPORTER_ASSERT(reporter, data[15] == m15);
 }

 static bool nearly_equal(const SkMatrix44& a, const SkMatrix44& b) {
     for (int i = 0; i < 4; ++i) {
         for (int j = 0; j < 4; ++j) {
             if (!nearly_equal_mscalar(a.get(i, j), b.get(i, j))) {
                 SkDebugf("not equal %g %g\n", a.get(i, j), b.get(i, j));
                 return false;
             }
         }
     }
     return true;
 }

 static bool is_identity(const SkMatrix44& m) {
     SkMatrix44 identity(SkMatrix44::kIdentity_Constructor);
     return nearly_equal(m, identity);
 }

 ///////////////////////////////////////////////////////////////////////////////
 static bool bits_isonly(int value, int mask) {
     return 0 == (value & ~mask);
 }

 static void test_constructor(skiatest::Reporter* reporter) {
     // Allocate a matrix on the heap
     SkMatrix44* placeholderMatrix = new SkMatrix44(SkMatrix44::kUninitialized_Constructor);
     SkAutoTDelete<SkMatrix44> deleteMe(placeholderMatrix);

     for (int row = 0; row < 4; ++row) {
         for (int col = 0; col < 4; ++col) {
             placeholderMatrix->setDouble(row, col, row * col);
         }
     }

     // Use placement-new syntax to trigger the constructor on top of the heap
     // address we already initialized. This allows us to check that the
     // constructor did avoid initializing the matrix contents.
     SkMatrix44* testMatrix = new(placeholderMatrix) SkMatrix44(SkMatrix44::kUninitialized_Constructor);
     REPORTER_ASSERT(reporter, testMatrix == placeholderMatrix);
     REPORTER_ASSERT(reporter, !testMatrix->isIdentity());
     for (int row = 0; row < 4; ++row) {
         for (int col = 0; col < 4; ++col) {
             REPORTER_ASSERT(reporter, nearly_equal_double(row * col, testMatrix->getDouble(row, col)));
         }
     }

     // Verify that kIdentity_Constructor really does initialize to an identity matrix.
     testMatrix = 0;
     testMatrix = new(placeholderMatrix) SkMatrix44(SkMatrix44::kIdentity_Constructor);
     REPORTER_ASSERT(reporter, testMatrix == placeholderMatrix);
     REPORTER_ASSERT(reporter, testMatrix->isIdentity());
     REPORTER_ASSERT(reporter, *testMatrix == SkMatrix44::I());
 }

 static void test_translate(skiatest::Reporter* reporter) {
     SkMatrix44 mat(SkMatrix44::kUninitialized_Constructor);
     SkMatrix44 inverse(SkMatrix44::kUninitialized_Constructor);

     mat.setTranslate(0, 0, 0);
     REPORTER_ASSERT(reporter, bits_isonly(mat.getType(), SkMatrix44::kIdentity_Mask));
     mat.setTranslate(1, 2, 3);
     REPORTER_ASSERT(reporter, bits_isonly(mat.getType(), SkMatrix44::kTranslate_Mask));
     REPORTER_ASSERT(reporter, mat.invert(&inverse));
     REPORTER_ASSERT(reporter, bits_isonly(inverse.getType(), SkMatrix44::kTranslate_Mask));

     SkMatrix44 a(SkMatrix44::kUninitialized_Constructor);
     SkMatrix44 b(SkMatrix44::kUninitialized_Constructor);
     SkMatrix44 c(SkMatrix44::kUninitialized_Constructor);
     a.set3x3(1, 2, 3, 4, 5, 6, 7, 8, 9);
     b.setTranslate(10, 11, 12);

     c.setConcat(a, b);
     mat = a;
     mat.preTranslate(10, 11, 12);
     REPORTER_ASSERT(reporter, mat == c);

     c.setConcat(b, a);
     mat = a;
     mat.postTranslate(10, 11, 12);
     REPORTER_ASSERT(reporter, mat == c);
 }

 static void test_scale(skiatest::Reporter* reporter) {
     SkMatrix44 mat(SkMatrix44::kUninitialized_Constructor);
     SkMatrix44 inverse(SkMatrix44::kUninitialized_Constructor);

     mat.setScale(1, 1, 1);
     REPORTER_ASSERT(reporter, bits_isonly(mat.getType(), SkMatrix44::kIdentity_Mask));
     mat.setScale(1, 2, 3);
     REPORTER_ASSERT(reporter, bits_isonly(mat.getType(), SkMatrix44::kScale_Mask));
     REPORTER_ASSERT(reporter, mat.invert(&inverse));
     REPORTER_ASSERT(reporter, bits_isonly(inverse.getType(), SkMatrix44::kScale_Mask));

     SkMatrix44 a(SkMatrix44::kUninitialized_Constructor);
     SkMatrix44 b(SkMatrix44::kUninitialized_Constructor);
     SkMatrix44 c(SkMatrix44::kUninitialized_Constructor);
     a.set3x3(1, 2, 3, 4, 5, 6, 7, 8, 9);
     b.setScale(10, 11, 12);

     c.setConcat(a, b);
     mat = a;
     mat.preScale(10, 11, 12);
     REPORTER_ASSERT(reporter, mat == c);

     c.setConcat(b, a);
     mat = a;
     mat.postScale(10, 11, 12);
     REPORTER_ASSERT(reporter, mat == c);
 }

 static void make_i(SkMatrix44* mat) { mat->setIdentity(); }
 static void make_t(SkMatrix44* mat) { mat->setTranslate(1, 2, 3); }
 static void make_s(SkMatrix44* mat) { mat->setScale(1, 2, 3); }
 static void make_st(SkMatrix44* mat) {
     mat->setScale(1, 2, 3);
     mat->postTranslate(1, 2, 3);
 }
 static void make_a(SkMatrix44* mat) {
     mat->setRotateDegreesAbout(1, 2, 3, 45);
 }
 static void make_p(SkMatrix44* mat) {
     SkMScalar data[] = {
         1, 2, 3, 4, 5, 6, 7, 8,
         1, 2, 3, 4, 5, 6, 7, 8,
     };
     mat->setRowMajor(data);
 }

 typedef void (*Make44Proc)(SkMatrix44*);

 static const Make44Proc gMakeProcs[] = {
     make_i, make_t, make_s, make_st, make_a, make_p
 };

 static void test_map2(skiatest::Reporter* reporter, const SkMatrix44& mat) {
     SkMScalar src2[] = { 1, 2 };
     SkMScalar src4[] = { src2[0], src2[1], 0, 1 };
     SkMScalar dstA[4], dstB[4];

     for (int i = 0; i < 4; ++i) {
         dstA[i] = 123456789;
         dstB[i] = 987654321;
     }

     mat.map2(src2, 1, dstA);
     mat.mapMScalars(src4, dstB);

     for (int i = 0; i < 4; ++i) {
         REPORTER_ASSERT(reporter, dstA[i] == dstB[i]);
     }
 }

 static void test_map2(skiatest::Reporter* reporter) {
     SkMatrix44 mat(SkMatrix44::kUninitialized_Constructor);

     for (size_t i = 0; i < SK_ARRAY_COUNT(gMakeProcs); ++i) {
         gMakeProcs[i](&mat);
         test_map2(reporter, mat);
     }
 }

 static void test_gettype(skiatest::Reporter* reporter) {
     SkMatrix44 matrix(SkMatrix44::kIdentity_Constructor);

     REPORTER_ASSERT(reporter, matrix.isIdentity());
     REPORTER_ASSERT(reporter, SkMatrix44::kIdentity_Mask == matrix.getType());

     int expectedMask;

     matrix.set(1, 1, 0);
     expectedMask = SkMatrix44::kScale_Mask;
     REPORTER_ASSERT(reporter, matrix.getType() == expectedMask);

     matrix.set(0, 3, 1);    // translate-x
     expectedMask |= SkMatrix44::kTranslate_Mask;
     REPORTER_ASSERT(reporter, matrix.getType() == expectedMask);

     matrix.set(2, 0, 1);
     expectedMask |= SkMatrix44::kAffine_Mask;
     REPORTER_ASSERT(reporter, matrix.getType() == expectedMask);

     matrix.set(3, 2, 1);
     REPORTER_ASSERT(reporter, matrix.getType() & SkMatrix44::kPerspective_Mask);

     // ensure that negative zero is treated as zero
     SkMScalar dx = 0;
     SkMScalar dy = 0;
     SkMScalar dz = 0;
     matrix.setTranslate(-dx, -dy, -dz);
     REPORTER_ASSERT(reporter, matrix.isIdentity());
     matrix.preTranslate(-dx, -dy, -dz);
     REPORTER_ASSERT(reporter, matrix.isIdentity());
     matrix.postTranslate(-dx, -dy, -dz);
     REPORTER_ASSERT(reporter, matrix.isIdentity());
 }

 static void test_common_angles(skiatest::Reporter* reporter) {
     SkMatrix44 rot(SkMatrix44::kUninitialized_Constructor);
     // Test precision of rotation in common cases
     int common_angles[] = { 0, 90, -90, 180, -180, 270, -270, 360, -360 };
     for (int i = 0; i < 9; ++i) {
         rot.setRotateDegreesAbout(0, 0, -1, SkIntToScalar(common_angles[i]));

         SkMatrix rot3x3 = rot;
         REPORTER_ASSERT(reporter, rot3x3.rectStaysRect());
     }
 }

 static void test_concat(skiatest::Reporter* reporter) {
     int i;
     SkMatrix44 a(SkMatrix44::kUninitialized_Constructor);
     SkMatrix44 b(SkMatrix44::kUninitialized_Constructor);
     SkMatrix44 c(SkMatrix44::kUninitialized_Constructor);
     SkMatrix44 d(SkMatrix44::kUninitialized_Constructor);

     a.setTranslate(10, 10, 10);
     b.setScale(2, 2, 2);

     SkScalar src[8] = {
         0, 0, 0, 1,
         1, 1, 1, 1
     };
     SkScalar dst[8];

     c.setConcat(a, b);

     d = a;
     d.preConcat(b);
     REPORTER_ASSERT(reporter, d == c);

     c.mapScalars(src, dst); c.mapScalars(src + 4, dst + 4);
     for (i = 0; i < 3; ++i) {
         REPORTER_ASSERT(reporter, 10 == dst[i]);
         REPORTER_ASSERT(reporter, 12 == dst[i + 4]);
     }

     c.setConcat(b, a);

     d = a;
     d.postConcat(b);
     REPORTER_ASSERT(reporter, d == c);

     c.mapScalars(src, dst); c.mapScalars(src + 4, dst + 4);
     for (i = 0; i < 3; ++i) {
         REPORTER_ASSERT(reporter, 20 == dst[i]);
         REPORTER_ASSERT(reporter, 22 == dst[i + 4]);
     }
 }

 static void test_determinant(skiatest::Reporter* reporter) {
     SkMatrix44 a(SkMatrix44::kIdentity_Constructor);
     REPORTER_ASSERT(reporter, nearly_equal_double(1, a.determinant()));
     a.set(1, 1, 2);
     REPORTER_ASSERT(reporter, nearly_equal_double(2, a.determinant()));
     SkMatrix44 b(SkMatrix44::kUninitialized_Constructor);
     REPORTER_ASSERT(reporter, a.invert(&b));
     REPORTER_ASSERT(reporter, nearly_equal_double(0.5, b.determinant()));
     SkMatrix44 c = b = a;
     c.set(0, 1, 4);
     b.set(1, 0, 4);
     REPORTER_ASSERT(reporter,
                     nearly_equal_double(a.determinant(),
                                         b.determinant()));
     SkMatrix44 d = a;
     d.set(0, 0, 8);
     REPORTER_ASSERT(reporter, nearly_equal_double(16, d.determinant()));

     SkMatrix44 e = a;
     e.postConcat(d);
     REPORTER_ASSERT(reporter, nearly_equal_double(32, e.determinant()));
     e.set(0, 0, 0);
     REPORTER_ASSERT(reporter, nearly_equal_double(0, e.determinant()));
 }

 static void test_invert(skiatest::Reporter* reporter) {
     SkMatrix44 inverse(SkMatrix44::kUninitialized_Constructor);
     double inverseData[16];

     SkMatrix44 identity(SkMatrix44::kIdentity_Constructor);
     identity.invert(&inverse);
     inverse.asRowMajord(inverseData);
     assert16<double>(reporter, inverseData,
                      1, 0, 0, 0,
                      0, 1, 0, 0,
                      0, 0, 1, 0,
                      0, 0, 0, 1);

     SkMatrix44 translation(SkMatrix44::kUninitialized_Constructor);
     translation.setTranslate(2, 3, 4);
     translation.invert(&inverse);
     inverse.asRowMajord(inverseData);
     assert16<double>(reporter, inverseData,
                      1, 0, 0, -2,
                      0, 1, 0, -3,
                      0, 0, 1, -4,
                      0, 0, 0, 1);

     SkMatrix44 scale(SkMatrix44::kUninitialized_Constructor);
     scale.setScale(2, 4, 8);
     scale.invert(&inverse);
     inverse.asRowMajord(inverseData);
     assert16<double>(reporter, inverseData,
                      0.5, 0,    0,     0,
                      0,   0.25, 0,     0,
                      0,   0,    0.125, 0,
                      0,   0,    0,     1);

     SkMatrix44 scaleTranslation(SkMatrix44::kUninitialized_Constructor);
     scaleTranslation.setScale(10, 100, 1000);
     scaleTranslation.preTranslate(2, 3, 4);
     scaleTranslation.invert(&inverse);
     inverse.asRowMajord(inverseData);
     assert16<double>(reporter, inverseData,
                      0.1,  0,    0,   -2,
                      0,   0.01,  0,   -3,
                      0,    0,  0.001, -4,
                      0,    0,    0,   1);

     SkMatrix44 rotation(SkMatrix44::kUninitialized_Constructor);
     rotation.setRotateDegreesAbout(0, 0, 1, 90);
     rotation.invert(&inverse);
     SkMatrix44 expected(SkMatrix44::kUninitialized_Constructor);
     double expectedInverseRotation[16] =
             {0,  1, 0, 0,
              -1, 0, 0, 0,
              0,  0, 1, 0,
              0,  0, 0, 1};
     expected.setRowMajord(expectedInverseRotation);
     REPORTER_ASSERT(reporter, nearly_equal(expected, inverse));

     SkMatrix44 affine(SkMatrix44::kUninitialized_Constructor);
     affine.setRotateDegreesAbout(0, 0, 1, 90);
     affine.preScale(10, 20, 100);
     affine.preTranslate(2, 3, 4);
     affine.invert(&inverse);
     double expectedInverseAffine[16] =
             {0,    0.1,  0,   -2,
              -0.05, 0,   0,   -3,
              0,     0,  0.01, -4,
              0,     0,   0,   1};
     expected.setRowMajord(expectedInverseAffine);
     REPORTER_ASSERT(reporter, nearly_equal(expected, inverse));

     SkMatrix44 perspective(SkMatrix44::kIdentity_Constructor);
     perspective.setDouble(3, 2, 1.0);
     perspective.invert(&inverse);
     double expectedInversePerspective[16] =
             {1, 0,  0, 0,
              0, 1,  0, 0,
              0, 0,  1, 0,
              0, 0, -1, 1};
     expected.setRowMajord(expectedInversePerspective);
     REPORTER_ASSERT(reporter, nearly_equal(expected, inverse));

     SkMatrix44 affineAndPerspective(SkMatrix44::kIdentity_Constructor);
     affineAndPerspective.setDouble(3, 2, 1.0);
     affineAndPerspective.preScale(10, 20, 100);
     affineAndPerspective.preTranslate(2, 3, 4);
     affineAndPerspective.invert(&inverse);
     double expectedInverseAffineAndPerspective[16] =
             {0.1, 0,    2,   -2,
              0,  0.05,  3,   -3,
              0,   0,   4.01, -4,
              0,   0,   -1,    1};
     expected.setRowMajord(expectedInverseAffineAndPerspective);
     REPORTER_ASSERT(reporter, nearly_equal(expected, inverse));
 }

 static void test_transpose(skiatest::Reporter* reporter) {
     SkMatrix44 a(SkMatrix44::kUninitialized_Constructor);
     SkMatrix44 b(SkMatrix44::kUninitialized_Constructor);

     int i = 0;
     for (int row = 0; row < 4; ++row) {
         for (int col = 0; col < 4; ++col) {
             a.setDouble(row, col, i);
             b.setDouble(col, row, i++);
         }
     }

     a.transpose();
     REPORTER_ASSERT(reporter, nearly_equal(a, b));
 }

 static void test_get_set_double(skiatest::Reporter* reporter) {
     SkMatrix44 a(SkMatrix44::kUninitialized_Constructor);
     for (int row = 0; row < 4; ++row) {
         for (int col = 0; col < 4; ++col) {
             a.setDouble(row, col, 3.141592653589793);
             REPORTER_ASSERT(reporter,
                             nearly_equal_double(3.141592653589793,
                                                 a.getDouble(row, col)));
             a.setDouble(row, col, 0);
             REPORTER_ASSERT(reporter,
                             nearly_equal_double(0, a.getDouble(row, col)));
         }
     }
 }

 static void test_set_row_col_major(skiatest::Reporter* reporter) {
     SkMatrix44 a(SkMatrix44::kUninitialized_Constructor);
     SkMatrix44 b(SkMatrix44::kUninitialized_Constructor);

     for (int row = 0; row < 4; ++row) {
         for (int col = 0; col < 4; ++col) {
             a.setDouble(row, col, row * 4 + col);
         }
     }

     double bufferd[16];
     float bufferf[16];
     a.asColMajord(bufferd);
     b.setColMajord(bufferd);
     REPORTER_ASSERT(reporter, nearly_equal(a, b));
     b.setRowMajord(bufferd);
     b.transpose();
     REPORTER_ASSERT(reporter, nearly_equal(a, b));
     a.asColMajorf(bufferf);
     b.setColMajorf(bufferf);
     REPORTER_ASSERT(reporter, nearly_equal(a, b));
     b.setRowMajorf(bufferf);
     b.transpose();
     REPORTER_ASSERT(reporter, nearly_equal(a, b));
 }

 static void test_3x3_conversion(skiatest::Reporter* reporter) {
     SkMScalar values4x4[16] = { 1, 2, 3, 4,
                                 5, 6, 7, 8,
                                 9, 10, 11, 12,
                                 13, 14, 15, 16 };
     SkScalar values3x3[9] = { 1, 2, 4,
                               5, 6, 8,
                               13, 14, 16 };
     SkMScalar values4x4flattened[16] = { 1, 2, 0, 4,
                                          5, 6, 0, 8,
                                          0, 0, 1, 0,
                                          13, 14, 0, 16 };
     SkMatrix44 a44(SkMatrix44::kUninitialized_Constructor);
     a44.setRowMajor(values4x4);

     SkMatrix a33 = a44;
     SkMatrix expected33;
     for (int i = 0; i < 9; i++) expected33[i] = values3x3[i];
     REPORTER_ASSERT(reporter, expected33 == a33);

     SkMatrix44 a44flattened = a33;
     SkMatrix44 expected44flattened(SkMatrix44::kUninitialized_Constructor);
     expected44flattened.setRowMajor(values4x4flattened);
     REPORTER_ASSERT(reporter, nearly_equal(a44flattened, expected44flattened));

     // Test that a point with a Z value of 0 is transformed the same way.
     SkScalar vec4[4] = { 2, 4, 0, 8 };
     SkScalar vec3[3] = { 2, 4, 8 };

     SkScalar vec4transformed[4];
     SkScalar vec3transformed[3];
     SkScalar vec4transformed2[4];
     a44.mapScalars(vec4, vec4transformed);
     a33.mapHomogeneousPoints(vec3transformed, vec3, 1);
     a44flattened.mapScalars(vec4, vec4transformed2);
     REPORTER_ASSERT(reporter, nearly_equal_scalar(vec4transformed[0], vec3transformed[0]));
     REPORTER_ASSERT(reporter, nearly_equal_scalar(vec4transformed[1], vec3transformed[1]));
     REPORTER_ASSERT(reporter, nearly_equal_scalar(vec4transformed[3], vec3transformed[2]));
     REPORTER_ASSERT(reporter, nearly_equal_scalar(vec4transformed[0], vec4transformed2[0]));
     REPORTER_ASSERT(reporter, nearly_equal_scalar(vec4transformed[1], vec4transformed2[1]));
     REPORTER_ASSERT(reporter, !nearly_equal_scalar(vec4transformed[2], vec4transformed2[2]));
     REPORTER_ASSERT(reporter, nearly_equal_scalar(vec4transformed[3], vec4transformed2[3]));
 }

 static void TestMatrix44(skiatest::Reporter* reporter) {
     SkMatrix44 mat(SkMatrix44::kUninitialized_Constructor);
     SkMatrix44 inverse(SkMatrix44::kUninitialized_Constructor);
     SkMatrix44 iden1(SkMatrix44::kUninitialized_Constructor);
     SkMatrix44 iden2(SkMatrix44::kUninitialized_Constructor);
     SkMatrix44 rot(SkMatrix44::kUninitialized_Constructor);

     mat.setTranslate(1, 1, 1);
     mat.invert(&inverse);
     iden1.setConcat(mat, inverse);
     REPORTER_ASSERT(reporter, is_identity(iden1));

     mat.setScale(2, 2, 2);
     mat.invert(&inverse);
     iden1.setConcat(mat, inverse);
     REPORTER_ASSERT(reporter, is_identity(iden1));

     mat.setScale(SK_MScalar1/2, SK_MScalar1/2, SK_MScalar1/2);
     mat.invert(&inverse);
     iden1.setConcat(mat, inverse);
     REPORTER_ASSERT(reporter, is_identity(iden1));

     mat.setScale(3, 3, 3);
     rot.setRotateDegreesAbout(0, 0, -1, 90);
     mat.postConcat(rot);
     REPORTER_ASSERT(reporter, mat.invert(NULL));
     mat.invert(&inverse);
     iden1.setConcat(mat, inverse);
     REPORTER_ASSERT(reporter, is_identity(iden1));
     iden2.setConcat(inverse, mat);
     REPORTER_ASSERT(reporter, is_identity(iden2));

     // test tiny-valued matrix inverse
     mat.reset();
     mat.setScale(1.0e-12, 1.0e-12, 1.0e-12);
     rot.setRotateDegreesAbout(0, 0, -1, 90);
     mat.postConcat(rot);
     mat.postTranslate(1.0e-12, 1.0e-12, 1.0e-12);
     REPORTER_ASSERT(reporter, mat.invert(NULL));
     mat.invert(&inverse);
     iden1.setConcat(mat, inverse);
     REPORTER_ASSERT(reporter, is_identity(iden1));

     // test mixed-valued matrix inverse
     mat.reset();
     mat.setScale(1.0e-10, 3.0, 1.0e+10);
     rot.setRotateDegreesAbout(0, 0, -1, 90);
     mat.postConcat(rot);
     mat.postTranslate(1.0e+10, 3.0, 1.0e-10);
     REPORTER_ASSERT(reporter, mat.invert(NULL));
     mat.invert(&inverse);
     iden1.setConcat(mat, inverse);
     REPORTER_ASSERT(reporter, is_identity(iden1));

     // test degenerate matrix
     mat.reset();
     mat.set3x3(1.0, 1.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0);
     REPORTER_ASSERT(reporter, !mat.invert(NULL));

     // test rol/col Major getters
     {
         mat.setTranslate(2, 3, 4);
         float dataf[16];
         double datad[16];

         mat.asColMajorf(dataf);
         assert16<float>(reporter, dataf,
                  1, 0, 0, 0,
                  0, 1, 0, 0,
                  0, 0, 1, 0,
                  2, 3, 4, 1);
         mat.asColMajord(datad);
         assert16<double>(reporter, datad, 1, 0, 0, 0,
                         0, 1, 0, 0,
                         0, 0, 1, 0,
                         2, 3, 4, 1);
         mat.asRowMajorf(dataf);
         assert16<float>(reporter, dataf, 1, 0, 0, 2,
                         0, 1, 0, 3,
                         0, 0, 1, 4,
                         0, 0, 0, 1);
         mat.asRowMajord(datad);
         assert16<double>(reporter, datad, 1, 0, 0, 2,
                         0, 1, 0, 3,
                         0, 0, 1, 4,
                         0, 0, 0, 1);
     }

     test_concat(reporter);

     if (false) { // avoid bit rot, suppress warning (working on making this pass)
         test_common_angles(reporter);
     }

     test_constructor(reporter);
     test_gettype(reporter);
     test_determinant(reporter);
     test_invert(reporter);
     test_transpose(reporter);
     test_get_set_double(reporter);
     test_set_row_col_major(reporter);
     test_translate(reporter);
     test_scale(reporter);
     test_map2(reporter);
     test_3x3_conversion(reporter);
 }

 #include "TestClassDef.h"
 DEFINE_TESTCLASS("Matrix44", Matrix44TestClass, TestMatrix44)
	/*
	* Copyright 2011 Google Inc.
	*
	* Use of this source code is governed by a BSD-style license that can be
	* found in the LICENSE file.
	*/

	#include "Test.h"
	#include "SkMatrix44.h"

	static bool nearly_equal_double(double a, double b) {
	const double tolerance = 1e-7;
	double diff = a - b;
	if (diff < 0)
	diff = -diff;
	return diff <= tolerance;
	}

	static bool nearly_equal_mscalar(SkMScalar a, SkMScalar b) {
	const SkMScalar tolerance = SK_MScalar1 / 200000;

	return SkTAbs<SkMScalar>(a - b) <= tolerance;
	}

	static bool nearly_equal_scalar(SkScalar a, SkScalar b) {
	// Note that we get more compounded error for multiple operations when
	// SK_SCALAR_IS_FIXED.
	#ifdef SK_SCALAR_IS_FLOAT
	const SkScalar tolerance = SK_Scalar1 / 200000;
	#else
	const SkScalar tolerance = SK_Scalar1 / 1024;
	#endif

	return SkScalarAbs(a - b) <= tolerance;
	}

	template <typename T> void assert16(skiatest::Reporter* reporter, const T data[],
	T m0, T m1, T m2, T m3,
	T m4, T m5, T m6, T m7,
	T m8, T m9, T m10, T m11,
	T m12, T m13, T m14, T m15) {
	REPORTER_ASSERT(reporter, data[0] == m0);
	REPORTER_ASSERT(reporter, data[1] == m1);
	REPORTER_ASSERT(reporter, data[2] == m2);
	REPORTER_ASSERT(reporter, data[3] == m3);

	REPORTER_ASSERT(reporter, data[4] == m4);
	REPORTER_ASSERT(reporter, data[5] == m5);
	REPORTER_ASSERT(reporter, data[6] == m6);
	REPORTER_ASSERT(reporter, data[7] == m7);

	REPORTER_ASSERT(reporter, data[8] == m8);
	REPORTER_ASSERT(reporter, data[9] == m9);
	REPORTER_ASSERT(reporter, data[10] == m10);
	REPORTER_ASSERT(reporter, data[11] == m11);

	REPORTER_ASSERT(reporter, data[12] == m12);
	REPORTER_ASSERT(reporter, data[13] == m13);
	REPORTER_ASSERT(reporter, data[14] == m14);
	REPORTER_ASSERT(reporter, data[15] == m15);
	}

	static bool nearly_equal(const SkMatrix44& a, const SkMatrix44& b) {
	for (int i = 0; i < 4; ++i) {
	for (int j = 0; j < 4; ++j) {
	if (!nearly_equal_mscalar(a.get(i, j), b.get(i, j))) {
	SkDebugf("not equal %g %g\n", a.get(i, j), b.get(i, j));
	return false;
	}
	}
	}
	return true;
	}

	static bool is_identity(const SkMatrix44& m) {
	SkMatrix44 identity(SkMatrix44::kIdentity_Constructor);
	return nearly_equal(m, identity);
	}

	///////////////////////////////////////////////////////////////////////////////
	static bool bits_isonly(int value, int mask) {
	return 0 == (value & ~mask);
	}

	static void test_constructor(skiatest::Reporter* reporter) {
	// Allocate a matrix on the heap
	SkMatrix44* placeholderMatrix = new SkMatrix44(SkMatrix44::kUninitialized_Constructor);
	SkAutoTDelete<SkMatrix44> deleteMe(placeholderMatrix);

	for (int row = 0; row < 4; ++row) {
	for (int col = 0; col < 4; ++col) {
	placeholderMatrix->setDouble(row, col, row * col);
	}
	}

	// Use placement-new syntax to trigger the constructor on top of the heap
	// address we already initialized. This allows us to check that the
	// constructor did avoid initializing the matrix contents.
	SkMatrix44* testMatrix = new(placeholderMatrix) SkMatrix44(SkMatrix44::kUninitialized_Constructor);
	REPORTER_ASSERT(reporter, testMatrix == placeholderMatrix);
	REPORTER_ASSERT(reporter, !testMatrix->isIdentity());
	for (int row = 0; row < 4; ++row) {
	for (int col = 0; col < 4; ++col) {
	REPORTER_ASSERT(reporter, nearly_equal_double(row * col, testMatrix->getDouble(row, col)));
	}
	}

	// Verify that kIdentity_Constructor really does initialize to an identity matrix.
	testMatrix = 0;
	testMatrix = new(placeholderMatrix) SkMatrix44(SkMatrix44::kIdentity_Constructor);
	REPORTER_ASSERT(reporter, testMatrix == placeholderMatrix);
	REPORTER_ASSERT(reporter, testMatrix->isIdentity());
	REPORTER_ASSERT(reporter, *testMatrix == SkMatrix44::I());
	}

	static void test_translate(skiatest::Reporter* reporter) {
	SkMatrix44 mat(SkMatrix44::kUninitialized_Constructor);
	SkMatrix44 inverse(SkMatrix44::kUninitialized_Constructor);

	mat.setTranslate(0, 0, 0);
	REPORTER_ASSERT(reporter, bits_isonly(mat.getType(), SkMatrix44::kIdentity_Mask));
	mat.setTranslate(1, 2, 3);
	REPORTER_ASSERT(reporter, bits_isonly(mat.getType(), SkMatrix44::kTranslate_Mask));
	REPORTER_ASSERT(reporter, mat.invert(&inverse));
	REPORTER_ASSERT(reporter, bits_isonly(inverse.getType(), SkMatrix44::kTranslate_Mask));

	SkMatrix44 a(SkMatrix44::kUninitialized_Constructor);
	SkMatrix44 b(SkMatrix44::kUninitialized_Constructor);
	SkMatrix44 c(SkMatrix44::kUninitialized_Constructor);
	a.set3x3(1, 2, 3, 4, 5, 6, 7, 8, 9);
	b.setTranslate(10, 11, 12);

	c.setConcat(a, b);
	mat = a;
	mat.preTranslate(10, 11, 12);
	REPORTER_ASSERT(reporter, mat == c);

	c.setConcat(b, a);
	mat = a;
	mat.postTranslate(10, 11, 12);
	REPORTER_ASSERT(reporter, mat == c);
	}

	static void test_scale(skiatest::Reporter* reporter) {
	SkMatrix44 mat(SkMatrix44::kUninitialized_Constructor);
	SkMatrix44 inverse(SkMatrix44::kUninitialized_Constructor);

	mat.setScale(1, 1, 1);
	REPORTER_ASSERT(reporter, bits_isonly(mat.getType(), SkMatrix44::kIdentity_Mask));
	mat.setScale(1, 2, 3);
	REPORTER_ASSERT(reporter, bits_isonly(mat.getType(), SkMatrix44::kScale_Mask));
	REPORTER_ASSERT(reporter, mat.invert(&inverse));
	REPORTER_ASSERT(reporter, bits_isonly(inverse.getType(), SkMatrix44::kScale_Mask));

	SkMatrix44 a(SkMatrix44::kUninitialized_Constructor);
	SkMatrix44 b(SkMatrix44::kUninitialized_Constructor);
	SkMatrix44 c(SkMatrix44::kUninitialized_Constructor);
	a.set3x3(1, 2, 3, 4, 5, 6, 7, 8, 9);
	b.setScale(10, 11, 12);

	c.setConcat(a, b);
	mat = a;
	mat.preScale(10, 11, 12);
	REPORTER_ASSERT(reporter, mat == c);

	c.setConcat(b, a);
	mat = a;
	mat.postScale(10, 11, 12);
	REPORTER_ASSERT(reporter, mat == c);
	}

	static void make_i(SkMatrix44* mat) { mat->setIdentity(); }
	static void make_t(SkMatrix44* mat) { mat->setTranslate(1, 2, 3); }
	static void make_s(SkMatrix44* mat) { mat->setScale(1, 2, 3); }
	static void make_st(SkMatrix44* mat) {
	mat->setScale(1, 2, 3);
	mat->postTranslate(1, 2, 3);
	}
	static void make_a(SkMatrix44* mat) {
	mat->setRotateDegreesAbout(1, 2, 3, 45);
	}
	static void make_p(SkMatrix44* mat) {
	SkMScalar data[] = {
	1, 2, 3, 4, 5, 6, 7, 8,
	1, 2, 3, 4, 5, 6, 7, 8,
	};
	mat->setRowMajor(data);
	}

	typedef void (Make44Proc)(SkMatrix44);

	static const Make44Proc gMakeProcs[] = {
	make_i, make_t, make_s, make_st, make_a, make_p
	};

	static void test_map2(skiatest::Reporter* reporter, const SkMatrix44& mat) {
	SkMScalar src2[] = { 1, 2 };
	SkMScalar src4[] = { src2[0], src2[1], 0, 1 };
	SkMScalar dstA[4], dstB[4];

	for (int i = 0; i < 4; ++i) {
	dstA[i] = 123456789;
	dstB[i] = 987654321;
	}

	mat.map2(src2, 1, dstA);
	mat.mapMScalars(src4, dstB);

	for (int i = 0; i < 4; ++i) {
	REPORTER_ASSERT(reporter, dstA[i] == dstB[i]);
	}
	}

	static void test_map2(skiatest::Reporter* reporter) {
	SkMatrix44 mat(SkMatrix44::kUninitialized_Constructor);

	for (size_t i = 0; i < SK_ARRAY_COUNT(gMakeProcs); ++i) {
	gMakeProcs[i](&mat);
	test_map2(reporter, mat);
	}
	}

	static void test_gettype(skiatest::Reporter* reporter) {
	SkMatrix44 matrix(SkMatrix44::kIdentity_Constructor);

	REPORTER_ASSERT(reporter, matrix.isIdentity());
	REPORTER_ASSERT(reporter, SkMatrix44::kIdentity_Mask == matrix.getType());

	int expectedMask;

	matrix.set(1, 1, 0);
	expectedMask = SkMatrix44::kScale_Mask;
	REPORTER_ASSERT(reporter, matrix.getType() == expectedMask);

	matrix.set(0, 3, 1); // translate-x
	expectedMask \|= SkMatrix44::kTranslate_Mask;
	REPORTER_ASSERT(reporter, matrix.getType() == expectedMask);

	matrix.set(2, 0, 1);
	expectedMask \|= SkMatrix44::kAffine_Mask;
	REPORTER_ASSERT(reporter, matrix.getType() == expectedMask);

	matrix.set(3, 2, 1);
	REPORTER_ASSERT(reporter, matrix.getType() & SkMatrix44::kPerspective_Mask);

	// ensure that negative zero is treated as zero
	SkMScalar dx = 0;
	SkMScalar dy = 0;
	SkMScalar dz = 0;
	matrix.setTranslate(-dx, -dy, -dz);
	REPORTER_ASSERT(reporter, matrix.isIdentity());
	matrix.preTranslate(-dx, -dy, -dz);
	REPORTER_ASSERT(reporter, matrix.isIdentity());
	matrix.postTranslate(-dx, -dy, -dz);
	REPORTER_ASSERT(reporter, matrix.isIdentity());
	}

	static void test_common_angles(skiatest::Reporter* reporter) {
	SkMatrix44 rot(SkMatrix44::kUninitialized_Constructor);
	// Test precision of rotation in common cases
	int common_angles[] = { 0, 90, -90, 180, -180, 270, -270, 360, -360 };
	for (int i = 0; i < 9; ++i) {
	rot.setRotateDegreesAbout(0, 0, -1, SkIntToScalar(common_angles[i]));

	SkMatrix rot3x3 = rot;
	REPORTER_ASSERT(reporter, rot3x3.rectStaysRect());
	}
	}

	static void test_concat(skiatest::Reporter* reporter) {
	int i;
	SkMatrix44 a(SkMatrix44::kUninitialized_Constructor);
	SkMatrix44 b(SkMatrix44::kUninitialized_Constructor);
	SkMatrix44 c(SkMatrix44::kUninitialized_Constructor);
	SkMatrix44 d(SkMatrix44::kUninitialized_Constructor);

	a.setTranslate(10, 10, 10);
	b.setScale(2, 2, 2);

	SkScalar src[8] = {
	0, 0, 0, 1,
	1, 1, 1, 1
	};
	SkScalar dst[8];

	c.setConcat(a, b);

	d = a;
	d.preConcat(b);
	REPORTER_ASSERT(reporter, d == c);

	c.mapScalars(src, dst); c.mapScalars(src + 4, dst + 4);
	for (i = 0; i < 3; ++i) {
	REPORTER_ASSERT(reporter, 10 == dst[i]);
	REPORTER_ASSERT(reporter, 12 == dst[i + 4]);
	}

	c.setConcat(b, a);

	d = a;
	d.postConcat(b);
	REPORTER_ASSERT(reporter, d == c);

	c.mapScalars(src, dst); c.mapScalars(src + 4, dst + 4);
	for (i = 0; i < 3; ++i) {
	REPORTER_ASSERT(reporter, 20 == dst[i]);
	REPORTER_ASSERT(reporter, 22 == dst[i + 4]);
	}
	}

	static void test_determinant(skiatest::Reporter* reporter) {
	SkMatrix44 a(SkMatrix44::kIdentity_Constructor);
	REPORTER_ASSERT(reporter, nearly_equal_double(1, a.determinant()));
	a.set(1, 1, 2);
	REPORTER_ASSERT(reporter, nearly_equal_double(2, a.determinant()));
	SkMatrix44 b(SkMatrix44::kUninitialized_Constructor);
	REPORTER_ASSERT(reporter, a.invert(&b));
	REPORTER_ASSERT(reporter, nearly_equal_double(0.5, b.determinant()));
	SkMatrix44 c = b = a;
	c.set(0, 1, 4);
	b.set(1, 0, 4);
	REPORTER_ASSERT(reporter,
	nearly_equal_double(a.determinant(),
	b.determinant()));
	SkMatrix44 d = a;
	d.set(0, 0, 8);
	REPORTER_ASSERT(reporter, nearly_equal_double(16, d.determinant()));

	SkMatrix44 e = a;
	e.postConcat(d);
	REPORTER_ASSERT(reporter, nearly_equal_double(32, e.determinant()));
	e.set(0, 0, 0);
	REPORTER_ASSERT(reporter, nearly_equal_double(0, e.determinant()));
	}

	static void test_invert(skiatest::Reporter* reporter) {
	SkMatrix44 inverse(SkMatrix44::kUninitialized_Constructor);
	double inverseData[16];

	SkMatrix44 identity(SkMatrix44::kIdentity_Constructor);
	identity.invert(&inverse);
	inverse.asRowMajord(inverseData);
	assert16<double>(reporter, inverseData,
	1, 0, 0, 0,
	0, 1, 0, 0,
	0, 0, 1, 0,
	0, 0, 0, 1);

	SkMatrix44 translation(SkMatrix44::kUninitialized_Constructor);
	translation.setTranslate(2, 3, 4);
	translation.invert(&inverse);
	inverse.asRowMajord(inverseData);
	assert16<double>(reporter, inverseData,
	1, 0, 0, -2,
	0, 1, 0, -3,
	0, 0, 1, -4,
	0, 0, 0, 1);

	SkMatrix44 scale(SkMatrix44::kUninitialized_Constructor);
	scale.setScale(2, 4, 8);
	scale.invert(&inverse);
	inverse.asRowMajord(inverseData);
	assert16<double>(reporter, inverseData,
	0.5, 0, 0, 0,
	0, 0.25, 0, 0,
	0, 0, 0.125, 0,
	0, 0, 0, 1);

	SkMatrix44 scaleTranslation(SkMatrix44::kUninitialized_Constructor);
	scaleTranslation.setScale(10, 100, 1000);
	scaleTranslation.preTranslate(2, 3, 4);
	scaleTranslation.invert(&inverse);
	inverse.asRowMajord(inverseData);
	assert16<double>(reporter, inverseData,
	0.1, 0, 0, -2,
	0, 0.01, 0, -3,
	0, 0, 0.001, -4,
	0, 0, 0, 1);

	SkMatrix44 rotation(SkMatrix44::kUninitialized_Constructor);
	rotation.setRotateDegreesAbout(0, 0, 1, 90);
	rotation.invert(&inverse);
	SkMatrix44 expected(SkMatrix44::kUninitialized_Constructor);
	double expectedInverseRotation[16] =
	{0, 1, 0, 0,
	-1, 0, 0, 0,
	0, 0, 1, 0,
	0, 0, 0, 1};
	expected.setRowMajord(expectedInverseRotation);
	REPORTER_ASSERT(reporter, nearly_equal(expected, inverse));

	SkMatrix44 affine(SkMatrix44::kUninitialized_Constructor);
	affine.setRotateDegreesAbout(0, 0, 1, 90);
	affine.preScale(10, 20, 100);
	affine.preTranslate(2, 3, 4);
	affine.invert(&inverse);
	double expectedInverseAffine[16] =
	{0, 0.1, 0, -2,
	-0.05, 0, 0, -3,
	0, 0, 0.01, -4,
	0, 0, 0, 1};
	expected.setRowMajord(expectedInverseAffine);
	REPORTER_ASSERT(reporter, nearly_equal(expected, inverse));

	SkMatrix44 perspective(SkMatrix44::kIdentity_Constructor);
	perspective.setDouble(3, 2, 1.0);
	perspective.invert(&inverse);
	double expectedInversePerspective[16] =
	{1, 0, 0, 0,
	0, 1, 0, 0,
	0, 0, 1, 0,
	0, 0, -1, 1};
	expected.setRowMajord(expectedInversePerspective);
	REPORTER_ASSERT(reporter, nearly_equal(expected, inverse));

	SkMatrix44 affineAndPerspective(SkMatrix44::kIdentity_Constructor);
	affineAndPerspective.setDouble(3, 2, 1.0);
	affineAndPerspective.preScale(10, 20, 100);
	affineAndPerspective.preTranslate(2, 3, 4);
	affineAndPerspective.invert(&inverse);
	double expectedInverseAffineAndPerspective[16] =
	{0.1, 0, 2, -2,
	0, 0.05, 3, -3,
	0, 0, 4.01, -4,
	0, 0, -1, 1};
	expected.setRowMajord(expectedInverseAffineAndPerspective);
	REPORTER_ASSERT(reporter, nearly_equal(expected, inverse));
	}

	static void test_transpose(skiatest::Reporter* reporter) {
	SkMatrix44 a(SkMatrix44::kUninitialized_Constructor);
	SkMatrix44 b(SkMatrix44::kUninitialized_Constructor);

	int i = 0;
	for (int row = 0; row < 4; ++row) {
	for (int col = 0; col < 4; ++col) {
	a.setDouble(row, col, i);
	b.setDouble(col, row, i++);
	}
	}

	a.transpose();
	REPORTER_ASSERT(reporter, nearly_equal(a, b));
	}

	static void test_get_set_double(skiatest::Reporter* reporter) {
	SkMatrix44 a(SkMatrix44::kUninitialized_Constructor);
	for (int row = 0; row < 4; ++row) {
	for (int col = 0; col < 4; ++col) {
	a.setDouble(row, col, 3.141592653589793);
	REPORTER_ASSERT(reporter,
	nearly_equal_double(3.141592653589793,
	a.getDouble(row, col)));
	a.setDouble(row, col, 0);
	REPORTER_ASSERT(reporter,
	nearly_equal_double(0, a.getDouble(row, col)));
	}
	}
	}

	static void test_set_row_col_major(skiatest::Reporter* reporter) {
	SkMatrix44 a(SkMatrix44::kUninitialized_Constructor);
	SkMatrix44 b(SkMatrix44::kUninitialized_Constructor);

	for (int row = 0; row < 4; ++row) {
	for (int col = 0; col < 4; ++col) {
	a.setDouble(row, col, row * 4 + col);
	}
	}

	double bufferd[16];
	float bufferf[16];
	a.asColMajord(bufferd);
	b.setColMajord(bufferd);
	REPORTER_ASSERT(reporter, nearly_equal(a, b));
	b.setRowMajord(bufferd);
	b.transpose();
	REPORTER_ASSERT(reporter, nearly_equal(a, b));
	a.asColMajorf(bufferf);
	b.setColMajorf(bufferf);
	REPORTER_ASSERT(reporter, nearly_equal(a, b));
	b.setRowMajorf(bufferf);
	b.transpose();
	REPORTER_ASSERT(reporter, nearly_equal(a, b));
	}

	static void test_3x3_conversion(skiatest::Reporter* reporter) {
	SkMScalar values4x4[16] = { 1, 2, 3, 4,
	5, 6, 7, 8,
	9, 10, 11, 12,
	13, 14, 15, 16 };
	SkScalar values3x3[9] = { 1, 2, 4,
	5, 6, 8,
	13, 14, 16 };
	SkMScalar values4x4flattened[16] = { 1, 2, 0, 4,
	5, 6, 0, 8,
	0, 0, 1, 0,
	13, 14, 0, 16 };
	SkMatrix44 a44(SkMatrix44::kUninitialized_Constructor);
	a44.setRowMajor(values4x4);

	SkMatrix a33 = a44;
	SkMatrix expected33;
	for (int i = 0; i < 9; i++) expected33[i] = values3x3[i];
	REPORTER_ASSERT(reporter, expected33 == a33);

	SkMatrix44 a44flattened = a33;
	SkMatrix44 expected44flattened(SkMatrix44::kUninitialized_Constructor);
	expected44flattened.setRowMajor(values4x4flattened);
	REPORTER_ASSERT(reporter, nearly_equal(a44flattened, expected44flattened));

	// Test that a point with a Z value of 0 is transformed the same way.
	SkScalar vec4[4] = { 2, 4, 0, 8 };
	SkScalar vec3[3] = { 2, 4, 8 };

	SkScalar vec4transformed[4];
	SkScalar vec3transformed[3];
	SkScalar vec4transformed2[4];
	a44.mapScalars(vec4, vec4transformed);
	a33.mapHomogeneousPoints(vec3transformed, vec3, 1);
	a44flattened.mapScalars(vec4, vec4transformed2);
	REPORTER_ASSERT(reporter, nearly_equal_scalar(vec4transformed[0], vec3transformed[0]));
	REPORTER_ASSERT(reporter, nearly_equal_scalar(vec4transformed[1], vec3transformed[1]));
	REPORTER_ASSERT(reporter, nearly_equal_scalar(vec4transformed[3], vec3transformed[2]));
	REPORTER_ASSERT(reporter, nearly_equal_scalar(vec4transformed[0], vec4transformed2[0]));
	REPORTER_ASSERT(reporter, nearly_equal_scalar(vec4transformed[1], vec4transformed2[1]));
	REPORTER_ASSERT(reporter, !nearly_equal_scalar(vec4transformed[2], vec4transformed2[2]));
	REPORTER_ASSERT(reporter, nearly_equal_scalar(vec4transformed[3], vec4transformed2[3]));
	}

	static void TestMatrix44(skiatest::Reporter* reporter) {
	SkMatrix44 mat(SkMatrix44::kUninitialized_Constructor);
	SkMatrix44 inverse(SkMatrix44::kUninitialized_Constructor);
	SkMatrix44 iden1(SkMatrix44::kUninitialized_Constructor);
	SkMatrix44 iden2(SkMatrix44::kUninitialized_Constructor);
	SkMatrix44 rot(SkMatrix44::kUninitialized_Constructor);

	mat.setTranslate(1, 1, 1);
	mat.invert(&inverse);
	iden1.setConcat(mat, inverse);
	REPORTER_ASSERT(reporter, is_identity(iden1));

	mat.setScale(2, 2, 2);
	mat.invert(&inverse);
	iden1.setConcat(mat, inverse);
	REPORTER_ASSERT(reporter, is_identity(iden1));

	mat.setScale(SK_MScalar1/2, SK_MScalar1/2, SK_MScalar1/2);
	mat.invert(&inverse);
	iden1.setConcat(mat, inverse);
	REPORTER_ASSERT(reporter, is_identity(iden1));

	mat.setScale(3, 3, 3);
	rot.setRotateDegreesAbout(0, 0, -1, 90);
	mat.postConcat(rot);
	REPORTER_ASSERT(reporter, mat.invert(NULL));
	mat.invert(&inverse);
	iden1.setConcat(mat, inverse);
	REPORTER_ASSERT(reporter, is_identity(iden1));
	iden2.setConcat(inverse, mat);
	REPORTER_ASSERT(reporter, is_identity(iden2));

	// test tiny-valued matrix inverse
	mat.reset();
	mat.setScale(1.0e-12, 1.0e-12, 1.0e-12);
	rot.setRotateDegreesAbout(0, 0, -1, 90);
	mat.postConcat(rot);
	mat.postTranslate(1.0e-12, 1.0e-12, 1.0e-12);
	REPORTER_ASSERT(reporter, mat.invert(NULL));
	mat.invert(&inverse);
	iden1.setConcat(mat, inverse);
	REPORTER_ASSERT(reporter, is_identity(iden1));

	// test mixed-valued matrix inverse
	mat.reset();
	mat.setScale(1.0e-10, 3.0, 1.0e+10);
	rot.setRotateDegreesAbout(0, 0, -1, 90);
	mat.postConcat(rot);
	mat.postTranslate(1.0e+10, 3.0, 1.0e-10);
	REPORTER_ASSERT(reporter, mat.invert(NULL));
	mat.invert(&inverse);
	iden1.setConcat(mat, inverse);
	REPORTER_ASSERT(reporter, is_identity(iden1));

	// test degenerate matrix
	mat.reset();
	mat.set3x3(1.0, 1.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0);
	REPORTER_ASSERT(reporter, !mat.invert(NULL));

	// test rol/col Major getters
	{
	mat.setTranslate(2, 3, 4);
	float dataf[16];
	double datad[16];

	mat.asColMajorf(dataf);
	assert16<float>(reporter, dataf,
	1, 0, 0, 0,
	0, 1, 0, 0,
	0, 0, 1, 0,
	2, 3, 4, 1);
	mat.asColMajord(datad);
	assert16<double>(reporter, datad, 1, 0, 0, 0,
	0, 1, 0, 0,
	0, 0, 1, 0,
	2, 3, 4, 1);
	mat.asRowMajorf(dataf);
	assert16<float>(reporter, dataf, 1, 0, 0, 2,
	0, 1, 0, 3,
	0, 0, 1, 4,
	0, 0, 0, 1);
	mat.asRowMajord(datad);
	assert16<double>(reporter, datad, 1, 0, 0, 2,
	0, 1, 0, 3,
	0, 0, 1, 4,
	0, 0, 0, 1);
	}

	test_concat(reporter);

	if (false) { // avoid bit rot, suppress warning (working on making this pass)
	test_common_angles(reporter);
	}

	test_constructor(reporter);
	test_gettype(reporter);
	test_determinant(reporter);
	test_invert(reporter);
	test_transpose(reporter);
	test_get_set_double(reporter);
	test_set_row_col_major(reporter);
	test_translate(reporter);
	test_scale(reporter);
	test_map2(reporter);
	test_3x3_conversion(reporter);
	}

	#include "TestClassDef.h"
	DEFINE_TESTCLASS("Matrix44", Matrix44TestClass, TestMatrix44)