DGtal 1.3.0
Loading...
Searching...
No Matches
testSimpleMatrix.cpp
Go to the documentation of this file.
1
31#include <iostream>
32#include "DGtal/base/Common.h"
33#include "DGtal/math/linalg/SimpleMatrix.h"
34#include "DGtal/math/linalg/CStaticMatrix.h"
35#include "DGtal/math/linalg/CDenseMatrix.h"
36#include "DGtal/math/linalg/CStaticVector.h"
37#include "DGtal/math/linalg/CDenseVector.h"
38#include "DGtal/math/linalg/CLinearAlgebra.h"
39#include "DGtal/helpers/StdDefs.h"
41
42using namespace std;
43using namespace DGtal;
44
46// Functions for testing class SimpleMatrix.
48
52bool testSimpleMatrix()
53{
54 unsigned int nbok = 0;
55 unsigned int nb = 0;
56
57 trace.beginBlock ( "Testing create ..." );
58
59 typedef SimpleMatrix<double,3,4> M34d;
60
61 M34d m34d;
62 trace.info() << m34d<<std::endl;
63
64 m34d.setComponent(1,2, 0.5);
65 trace.info() << m34d<<std::endl;
66
67 nbok += (m34d(1,2) == 0.5) ? 1 : 0;
68 nb++;
69 trace.info() << "(" << nbok << "/" << nb << ") "
70 << "true == true" << std::endl;
71
72 M34d matrix;
73 bool res=true;
74
75 matrix.constant(12.3);
76 trace.info() << matrix;
77 for(DGtal::Dimension i = 0; i< 3; ++i)
78 for(DGtal::Dimension j = 0; j< 4; ++j)
79 res = res && (matrix(i,j) == 12.3);
80 nbok += res ? 1 : 0;
81 nb++;
82 trace.info() << "(" << nbok << "/" << nb << ") "
83 << "all equals to 12.3" << std::endl;
84
85
86 trace.endBlock();
87
88 return nbok == nb;
89}
90
91bool testArithm()
92{
93 unsigned int nbok = 0;
94 unsigned int nb = 0;
95
96
97 typedef SimpleMatrix<double,3,4> M34d;
98 typedef SimpleMatrix<double,4,3> M43d;
99 typedef SimpleMatrix<double,3,3> M33d;
100
101 M34d m34d, two,four;
102 M34d m34dbis, resadd, ressub;
103
104 two.constant(2);
105 four.constant(4);
106
107 for(DGtal::Dimension i = 0; i< 3; ++i)
108 for(DGtal::Dimension j = 0; j< 4; ++j)
109 {
110 m34d.setComponent(i,j,i*j);
111 m34dbis.setComponent(i,j,i+j);
112 resadd.setComponent(i,j,i*j+i+j);
113 ressub.setComponent(i,j,i*j-(double)i-(double)j);
114 }
115
116
117 trace.info() << m34d <<std::endl;
118 trace.info() << m34dbis<<std::endl;
119
120 trace.beginBlock ( "Testing add ..." );
121 nbok += ((m34d + m34dbis) == resadd) ? 1 : 0;
122 nb++;
123 trace.info() << "(" << nbok << "/" << nb << ") "
124 << "ok" << std::endl;
125 nbok += ((m34dbis + m34d) == resadd) ? 1 : 0;
126 nb++;
127 trace.info() << "(" << nbok << "/" << nb << ") "
128 << "ok commutative" << std::endl;
129
130 M34d other;
131 other += m34d;
132 nbok += (other == m34d) ? 1 : 0;
133 nb++;
134 trace.info() << "(" << nbok << "/" << nb << ") "
135 << "ok +=" << std::endl;
136
137 trace.endBlock();
138
139 trace.beginBlock ( "Testing substraction ..." );
140 nbok += ((m34d - m34dbis) == ressub) ? 1 : 0;
141 nb++;
142 trace.info()<<ressub<<std::endl;
143 trace.info()<<m34d - m34dbis<<std::endl;
144
145 trace.info() << "(" << nbok << "/" << nb << ") "
146 << "ok simple" << std::endl;
147 trace.endBlock();
148
149 trace.beginBlock ( "Testing scalar product/divide ..." );
150 nbok += ( (two*2.0) == four) ? 1 : 0;
151 nb++;
152 trace.info()<<ressub<<std::endl;
153 trace.info() << "(" << nbok << "/" << nb << ") "
154 << " [2]*2 == [4]" << std::endl;
155
156 nbok += ( two == four/2.0) ? 1 : 0;
157 nb++;
158 trace.info()<<ressub<<std::endl;
159 trace.info() << "(" << nbok << "/" << nb << ") "
160 << " [2]= [4]/2" << std::endl;
161 trace.endBlock();
162
163
164 trace.beginBlock ( "Testing transpose ..." );
165 M43d transp = m34d.transpose();
166 nbok += (transp.transpose() == m34d) ? 1 : 0;
167 nb++;
168 trace.info() << "(" << nbok << "/" << nb << ") "
169 << "ok idem potent" << std::endl;
170 trace.endBlock();
171
172 trace.beginBlock ( "Testing product ..." );
173
174 M43d one;
175 M33d eight33;
176
177 one.constant(1);
178 eight33.constant(8);
179 trace.info() << two * one<<std::endl;
180 nbok += (two * one == eight33) ? 1 : 0;
181 nb++;
182 trace.info() << "(" << nbok << "/" << nb << ") "
183 << " [2]*[1] = [8]" << std::endl;
184 trace.endBlock();
185
186
187
188 return nbok == nb;
189
190}
191
192bool testColRow()
193{
194 unsigned int nbok = 0;
195 unsigned int nb = 0;
196
197 DGtal::SimpleMatrix<double,3,4> mat;
198 for(DGtal::Dimension i = 0; i< 3; ++i)
199 for(DGtal::Dimension j = 0; j< 4; ++j)
200 mat.setComponent(i,j,i+j);
201
202 trace.beginBlock("Get Row");
203 trace.info() << mat <<std::endl;
204 DGtal::SimpleMatrix<double,3,4>::RowVector row;
205 row = mat.row(1);
206 trace.info() << row << std::endl;
207 nbok += (row[1] == 2 ) ? 1 : 0;
208 nb++;
209 trace.info() << "(" << nbok << "/" << nb << ") "
210 << " row value" << std::endl;
211 trace.endBlock();
212
213 trace.beginBlock("Get Col");
214 DGtal::SimpleMatrix<double,3,4>::ColumnVector col;
215 col = mat.column(1);
216 trace.info() << row << std::endl;
217 nbok += (col[1] == 2 ) ? 1 : 0;
218 nb++;
219 trace.info() << "(" << nbok << "/" << nb << ") "
220 << " col value" << std::endl;
221 trace.endBlock();
222
223
224
225 trace.beginBlock("Prod Matrix x Row^t");
226 //Row vector is a dim 4 vector
227 DGtal::SimpleMatrix<double,3,4>::RowVector r(1,2,3,4);
228 DGtal::SimpleMatrix<double,3,4>::ColumnVector c;
229 c = mat*r;
230 DGtal::SimpleMatrix<double,3,4>::ColumnVector expected(20,30,40);
231
232 trace.info() << c << std::endl;
233 nbok += (c == expected) ? 1 : 0;
234 nb++;
235 trace.info() << "(" << nbok << "/" << nb << ") "
236 << " mat*row^t" << std::endl;
237 trace.endBlock();
238
239 return nbok == nb;
240}
241
242bool testDetCofactor()
243{
244 unsigned int nbok = 0;
245 unsigned int nb = 0;
246
247 typedef DGtal::SimpleMatrix<double,2,2> MAT2;
248 MAT2 mat2;
249 mat2.setComponent(0,0,1);
250 mat2.setComponent(1,1,2);
251
252 trace.beginBlock("det2x2 tests...");
253 trace.info() << mat2<<std::endl;
254 trace.info() << mat2.determinant() << std::endl;
255 nbok += (mat2.determinant() == 2) ? 1 : 0;
256 nb++;
257 trace.info() << "(" << nbok << "/" << nb << ") "
258 << " 2" << std::endl;
259 trace.endBlock();
260
261 typedef DGtal::SimpleMatrix<double,3,3> MAT;
262 MAT mat;
263 mat.setComponent(0,0,1);
264 mat.setComponent(1,1,2);
265 mat.setComponent(2,2,4);
266
267 trace.beginBlock("det3x3 tests...");
268 trace.info() << mat<<std::endl;
269 nbok += (mat.determinant() == 8) ? 1 : 0;
270 nb++;
271 trace.info() << "(" << nbok << "/" << nb << ") "
272 << " 8" << std::endl;
273 trace.endBlock();
274
275
276 typedef DGtal::SimpleMatrix<double,4,4> MAT44;
277 MAT44 mat44;
278 mat44.setComponent(0,0,1);
279 mat44.setComponent(1,1,2);
280 mat44.setComponent(2,2,4);
281 mat44.setComponent(3,3,4);
282
283 trace.beginBlock("det4x4 tests...");
284 trace.info() << mat44 <<std::endl;
285 trace.info() << mat44.determinant() << std::endl;
286 nbok += (mat44.determinant() == 32) ? 1 : 0;
287 nb++;
288 trace.info() << "(" << nbok << "/" << nb << ") "
289 << " 32" << std::endl;
290 trace.endBlock();
291
292
293 return nbok == nb;
294}
295
296bool testM1Matrix()
297{
298 trace.beginBlock("Mx1 matrix test");
299 SimpleMatrix<double, 3,1> mat;
300 trace.info() << mat<<std::endl;
301 trace.endBlock();
302 return true;
303}
304
305bool testInverse()
306{
307 unsigned int nbok = 0;
308 unsigned int nb = 0;
309
310 trace.beginBlock("Inverse tests 2X2...");
311
312 typedef DGtal::SimpleMatrix<double,2,2> MAT2;
313 MAT2 mat2;
314 mat2.setComponent(0,0,1);
315 mat2.setComponent(1,1,2);
316
317 MAT2 Id2;
318 Id2.identity();
319
320 trace.info() << mat2<<std::endl;
321 trace.info() << mat2.inverse() << std::endl;
322 nbok += (( mat2 * mat2.inverse() )== Id2 ) ? 1 : 0;
323 nb++;
324 trace.info() << "(" << nbok << "/" << nb << ") "
325 << " M*M^-1=Id" << std::endl;
326
327 trace.endBlock();
328
329 trace.beginBlock("Inverse tests 6x6 random...");
330
331 typedef DGtal::SimpleMatrix<double,6,6> MAT6;
332 MAT6 mat;
333
334 for(unsigned int i=0; i< 6; i++)
335 for(unsigned int j=0; j< 6; j++)
336 mat.setComponent(i,j, rand() % 10);
337
338 MAT6 Id6;
339 Id6.identity();
340
341 trace.info() << "M= "<<mat<<std::endl;
342 trace.info() << "M^-1=" <<mat.inverse() << std::endl;
343 trace.info() << "det(M)= "<<mat.determinant() <<std::endl;
344 trace.info() << "M*M^-1= "<<mat.inverse()*mat << std::endl;
345
346
347
348 trace.endBlock();
349
350
351 return nbok == nb;
352}
353
354bool testConstructor()
355{
356 unsigned int nbok = 0;
357 unsigned int nb = 0;
358
359 trace.beginBlock( "Initilizer-list constructor test" );
360 SimpleMatrix<double, 3, 3> mat = {1, 2, 3, 4, 5, 6, 7, 8, 9};
361 trace.info() << mat << std::endl;
362
363 trace.info() << "Testing values: ";
364 trace.info() << mat( 0, 0 );
365 nbok += ( mat( 0, 0 ) == 1 ) ? 1 : 0;
366 nb++;
367 trace.info() << "(" << nbok << "/" << nb << ") ";
368
369 trace.info() << mat( 0, 1 );
370 nbok += ( mat( 0, 1 ) == 2 ) ? 1 : 0;
371 nb++;
372 trace.info() << "(" << nbok << "/" << nb << ") ";
373
374 trace.info() << mat( 2, 2 );
375 nbok += ( mat( 2, 2 ) == 9 ) ? 1 : 0;
376 nb++;
377 trace.info() << "(" << nbok << "/" << nb << ") ";
378
379 trace.info() << std::endl;
380 trace.endBlock();
381 return nbok == nb;
382}
383
384bool testConcepts()
385{
386 typedef DGtal::SimpleMatrix<double,3,3> Matrix;
387 typedef Matrix::ColumnVector Vector;
388
389 BOOST_CONCEPT_ASSERT(( concepts::CStaticVector<Vector> ));
390 BOOST_CONCEPT_ASSERT(( concepts::CDenseVector<Vector> ));
391 BOOST_CONCEPT_ASSERT(( concepts::CStaticMatrix<Matrix> ));
392 BOOST_CONCEPT_ASSERT(( concepts::CDenseMatrix<Matrix> ));
393 BOOST_CONCEPT_ASSERT(( concepts::CLinearAlgebra<Vector, Matrix> ));
394 BOOST_CONCEPT_ASSERT(( concepts::CStaticVector<DGtal::Z3i::Space::Vector> ));
395 BOOST_CONCEPT_ASSERT(( concepts::CDenseVector<DGtal::Z3i::Space::Vector> ));
396
397 return true;
398}
399
401// Standard services - public :
402
403int main( int argc, char** argv )
404{
405 trace.beginBlock ( "Testing class SimpleMatrix" );
406 trace.info() << "Args:";
407 for ( int i = 0; i < argc; ++i )
408 trace.info() << " " << argv[ i ];
409 trace.info() << endl;
410
411 bool res = testSimpleMatrix() && testArithm() && testColRow() &&
412 testDetCofactor() && testM1Matrix() && testInverse() &&
413 testConcepts() && testConstructor();
414 trace.emphase() << ( res ? "Passed." : "Error." ) << endl;
415 trace.endBlock();
416 return res ? 0 : 1;
417}
418// //
DGtal::PointVector
Aim: Implements basic operations that will be used in Point and Vector classes.
Definition: PointVector.h:593
DGtal::SimpleMatrix
Aim: implements basic MxN Matrix services (M,N>=1).
Definition: SimpleMatrix.h:76
DGtal::SimpleMatrix::column
ColumnVector column(const DGtal::Dimension j) const
DGtal::SimpleMatrix::row
RowVector row(const DGtal::Dimension i) const
DGtal::SimpleMatrix::setComponent
void setComponent(const DGtal::Dimension i, const DGtal::Dimension j, const Component &aValue)
DGtal::SimpleMatrix::constant
void constant(const Component &aScalar)
DGtal::Trace::beginBlock
void beginBlock(const std::string &keyword="")
DGtal::Trace::emphase
std::ostream & emphase()
DGtal::Trace::info
std::ostream & info()
DGtal::Trace::endBlock
double endBlock()
DGtal
DGtal is the top-level namespace which contains all DGtal functions and types.
Definition: ClosedIntegerHalfPlane.h:49
DGtal::Dimension
DGtal::uint32_t Dimension
Definition: Common.h:137
DGtal::trace
Trace trace
Definition: Common.h:154
std
STL namespace.
DGtal::concepts::CDenseMatrix
Aim: Represent any dynamic or static sized matrix having dense representation.
Definition: CDenseMatrix.h:89
DGtal::concepts::CDenseVector
Aim: Represent any dynamic or static sized matrix having dense representation.
Definition: CDenseVector.h:89
DGtal::concepts::CLinearAlgebra
Aim: Check right multiplication between matrix and vector and internal matrix multiplication....
Definition: CLinearAlgebra.h:94
DGtal::concepts::CStaticMatrix
Aim: Represent any static sized matrix having sparse or dense representation.
Definition: CStaticMatrix.h:90
DGtal::concepts::CStaticVector
Aim: Represent any static sized column vector having sparse or dense representation.
Definition: CStaticVector.h:91
main
int main()
Definition: testBits.cpp:56
Vector
FreemanChain< int >::Vector Vector
Definition: testCombinDSS.cpp:60
testInverse
bool testInverse()
Definition: testSimpleMatrix.cpp:305
testM1Matrix
bool testM1Matrix()
Definition: testSimpleMatrix.cpp:296
testDetCofactor
bool testDetCofactor()
Definition: testSimpleMatrix.cpp:242
testColRow
bool testColRow()
Definition: testSimpleMatrix.cpp:192
testConstructor
bool testConstructor()
Definition: testSimpleMatrix.cpp:354
testConcepts
bool testConcepts()
Definition: testSimpleMatrix.cpp:384
testArithm
bool testArithm()
Definition: testSimpleMatrix.cpp:91
testSimpleMatrix
bool testSimpleMatrix()
Definition: testSimpleMatrix.cpp:52
Generated on Sun Nov 27 2022 15:15:57 for DGtal by doxygen 1.9.5