determinant static method

double determinant(Float32List matrix, int offset)

Compute the determinant of the 4x4 matrix starting at offset.

Implementation

static double determinant(Float32List matrix, int offset) {
  final m0 = matrix[0 + offset];
  final m1 = matrix[1 + offset];
  final m2 = matrix[2 + offset];
  final m3 = matrix[3 + offset];
  final m4 = matrix[4 + offset];
  final m5 = matrix[5 + offset];
  final m6 = matrix[6 + offset];
  final m7 = matrix[7 + offset];

  final det2_01_01 = m0 * m5 - m1 * m4;
  final det2_01_02 = m0 * m6 - m2 * m4;
  final det2_01_03 = m0 * m7 - m3 * m4;
  final det2_01_12 = m1 * m6 - m2 * m5;
  final det2_01_13 = m1 * m7 - m3 * m5;
  final det2_01_23 = m2 * m7 - m3 * m6;

  final m8 = matrix[8 + offset];
  final m9 = matrix[9 + offset];
  final m10 = matrix[10 + offset];
  final m11 = matrix[11 + offset];

  final det3_201_012 = m8 * det2_01_12 - m9 * det2_01_02 + m10 * det2_01_01;
  final det3_201_013 = m8 * det2_01_13 - m9 * det2_01_03 + m11 * det2_01_01;
  final det3_201_023 = m8 * det2_01_23 - m10 * det2_01_03 + m11 * det2_01_02;
  final det3_201_123 = m9 * det2_01_23 - m10 * det2_01_13 + m11 * det2_01_12;

  final m12 = matrix[12 + offset];
  final m13 = matrix[13 + offset];
  final m14 = matrix[14 + offset];
  final m15 = matrix[15 + offset];

  return -det3_201_123 * m12 +
      det3_201_023 * m13 -
      det3_201_013 * m14 +
      det3_201_012 * m15;
}