gistfile1.txt

#include "Matrix.h"

float const Matrix::identityMatrix[] = {
    1,0,0,0,
    0,1,0,0,
    0,0,1,0,
    0,0,0,1
};

Matrix Matrix::MultiplyMatrices(const Matrix* m1, const Matrix* m2) {
    Matrix result;
    
    result.m = Matrix::Identity();

    unsigned int row;
    unsigned int column;
    unsigned int rowOffset;

    for(row = 0, rowOffset = row*4; row < 4; ++row, rowOffset = row*4)
        for(column = 0; column < 4; ++column)
            result.m[rowOffset + column] = (m1->m[rowOffset + 0] * m2->m[column + 0]) +
                                           (m1->m[rowOffset + 1] * m2->m[column + 4]) +
                                           (m1->m[rowOffset + 2] * m2->m[column + 8]) +
                                           (m1->m[rowOffset + 3] * m2->m[column + 12]);
    return result;
}