← Back to branch summary

~libmatrix-dev/libmatrix/trunk

« back to all changes in this revision

Viewing changes to matrix_inverse_test.cc

  • Committer: Jesse Barker
  • Date: 2011-06-16 21:35:35 UTC
  • Revision ID: jesse.barker@linaro.org-20110616213535-u48sdrp56pr8k10i
Add matrix_inverse_test and build it by default.

Show diffs side-by-side

added added

removed removed

Lines of Context:
matrix_inverse_test.cc
 
1
#include <iostream>
 
2
#include "mat.h"
 
3
 
 
4
using LibMatrix::mat4;
 
5
using LibMatrix::mat3;
 
6
using LibMatrix::mat2;
 
7
using std::cerr;
 
8
using std::cout;
 
9
using std::endl;
 
10
 
 
11
bool mat2OK()
 
12
{
 
13
    mat2 m;
 
14
    cout << "Starting with mat2 (should be identity): " << endl << endl;
 
15
    m.print();
 
16
 
 
17
    m[0][1] = -2.5;
 
18
    
 
19
    cout << endl << "Matrix should now have (0, 1) == -2.500000" << endl << endl;
 
20
    m.print();
 
21
    
 
22
    mat2 mi(m);
 
23
 
 
24
    cout << endl << "Copy of previous matrix (should have (0, 1) == -2.500000)" << endl << endl;
 
25
    mi.print();
 
26
 
 
27
    mi.inverse();
 
28
 
 
29
    cout << endl << "Inverse of copy: " << endl << endl;
 
30
    mi.print();
 
31
 
 
32
    mat2 i = m * mi;
 
33
 
 
34
    cout << endl << "Product of original and inverse (should be identity): " << endl << endl;
 
35
    i.print();
 
36
 
 
37
    mat2 ident;
 
38
    if (i != ident)
 
39
    {
 
40
        return false;
 
41
    }
 
42
 
 
43
    return true;
 
44
}
 
45
 
 
46
bool mat3OK()
 
47
{
 
48
    mat3 m;
 
49
    cout << "Starting with mat3 (should be identity): " << endl << endl;
 
50
    m.print();
 
51
 
 
52
    m[1][2] = -2.5;
 
53
    
 
54
    cout << endl << "Matrix should now have (1, 2) == -2.500000" << endl << endl;
 
55
    m.print();
 
56
    
 
57
    mat3 mi(m);
 
58
 
 
59
    cout << endl << "Copy of previous matrix (should have (1, 2) == -2.500000)" << endl << endl;
 
60
    mi.print();
 
61
 
 
62
    mi.inverse();
 
63
 
 
64
    cout << endl << "Inverse of copy: " << endl << endl;
 
65
    mi.print();
 
66
 
 
67
    mat3 i = m * mi;
 
68
 
 
69
    cout << endl << "Product of original and inverse (should be identity): " << endl << endl;
 
70
    i.print();
 
71
 
 
72
    mat3 ident;
 
73
    if (i != ident)
 
74
    {
 
75
        return false;
 
76
    }
 
77
 
 
78
    return true;
 
79
}
 
80
 
 
81
bool mat4OK()
 
82
{
 
83
    mat4 m;
 
84
    cout << "Starting with mat4 (should be identity): " << endl << endl;
 
85
    m.print();
 
86
 
 
87
    m[2][3] = -2.5;
 
88
    
 
89
    cout << endl << "Matrix should now have (2, 3) == -2.500000" << endl << endl;
 
90
    m.print();
 
91
    
 
92
    mat4 mi(m);
 
93
 
 
94
    cout << endl << "Copy of previous matrix (should have (2, 3) == -2.500000)" << endl << endl;
 
95
    mi.print();
 
96
 
 
97
    mi.inverse();
 
98
 
 
99
    cout << endl << "Inverse of copy: " << endl << endl;
 
100
    mi.print();
 
101
 
 
102
    mat4 i = m * mi;
 
103
 
 
104
    cout << endl <<  "Product of original and inverse (should be identity): " << endl << endl;
 
105
    i.print();
 
106
 
 
107
    mat4 ident;
 
108
    if (i != ident)
 
109
    {
 
110
        return false;
 
111
    }
 
112
 
 
113
    return true;
 
114
}
 
115
 
 
116
int
 
117
main(int argc, char** argv)
 
118
{
 
119
    if (!mat2OK())
 
120
    {
 
121
        cerr << "mat2::inverse() does not work!" << endl;
 
122
        return 1;
 
123
    }
 
124
    cout << "mat2::inverse() is okay!" << endl << endl;
 
125
 
 
126
    if (!mat3OK())
 
127
    {
 
128
        cerr << "mat3::inverse() does not work!" << endl;
 
129
        return 1;
 
130
    }
 
131
    cout << "mat3::inverse() is okay!" << endl << endl;
 
132
 
 
133
    if (!mat4OK())
 
134
    {
 
135
        cerr << "mat4::inverse() does not work!" << endl;
 
136
        return 1;
 
137
    }
 
138
    cout << "mat4::inverse() is okay!" << endl << endl;
 
139
 
 
140
    return 0;
 
141
}