Nektar++
Functions
Nektar::MatrixUnitTests Namespace Reference

Functions

 BOOST_AUTO_TEST_CASE (TestFullMatrixInversion)
 
 BOOST_AUTO_TEST_CASE (TestDiagonalMatrixInversion)
 
 BOOST_AUTO_TEST_CASE (TestNekMatrixConstruction)
 
 BOOST_AUTO_TEST_CASE (TestFullNekMatrixGetValue)
 
 BOOST_AUTO_TEST_CASE (TestDiagonalMatrixGetValue)
 
 BOOST_AUTO_TEST_CASE (TestFullNekMatrixSetValue)
 
 BOOST_AUTO_TEST_CASE (TestDiagonalNekMatrixSetValue)
 
 BOOST_AUTO_TEST_CASE (TestFullFullMatrixAddition)
 

Function Documentation

◆ BOOST_AUTO_TEST_CASE() [1/8]

Nektar::MatrixUnitTests::BOOST_AUTO_TEST_CASE ( TestDiagonalMatrixGetValue  )

Definition at line 241 of file testNekMatrix.cpp.

242{
243 UnitTests::RedirectCerrIfNeeded();
244 double data[] = {8.9, 3.4, 5.7};
245 Nektar::NekMatrix<double> m1(3, 3, data, eDIAGONAL);
246 std::shared_ptr<Nektar::Matrix<double>> m2(
247 new Nektar::NekMatrix<double>(3, 3, data, eDIAGONAL));
248 std::shared_ptr<Nektar::NekMatrix<double>> m3(
249 new Nektar::NekMatrix<double>(3, 3, data, eDIAGONAL));
250
251 BOOST_CHECK_EQUAL(m1(0, 0), 8.9);
252 BOOST_CHECK_EQUAL((*m2)(0, 0), 8.9);
253 BOOST_CHECK_EQUAL((*m3)(0, 0), 8.9);
254
255 BOOST_CHECK_EQUAL(m1(0, 1), 0.0);
256 BOOST_CHECK_EQUAL((*m2)(0, 1), 0.0);
257 BOOST_CHECK_EQUAL((*m3)(0, 1), 0.0);
258
259 BOOST_CHECK_EQUAL(m1(0, 2), 0.0);
260 BOOST_CHECK_EQUAL((*m2)(0, 2), 0.0);
261 BOOST_CHECK_EQUAL((*m3)(0, 2), 0.0);
262
263 BOOST_CHECK_EQUAL(m1(1, 0), 0.0);
264 BOOST_CHECK_EQUAL((*m2)(1, 0), 0.0);
265 BOOST_CHECK_EQUAL((*m3)(1, 0), 0.0);
266
267 BOOST_CHECK_EQUAL(m1(1, 1), 3.4);
268 BOOST_CHECK_EQUAL((*m2)(1, 1), 3.4);
269 BOOST_CHECK_EQUAL((*m3)(1, 1), 3.4);
270
271 BOOST_CHECK_EQUAL(m1(1, 2), 0.0);
272 BOOST_CHECK_EQUAL((*m2)(1, 2), 0.0);
273 BOOST_CHECK_EQUAL((*m3)(1, 2), 0.0);
274
275 BOOST_CHECK_EQUAL(m1(2, 0), 0.0);
276 BOOST_CHECK_EQUAL((*m2)(2, 0), 0.0);
277 BOOST_CHECK_EQUAL((*m3)(2, 0), 0.0);
278
279 BOOST_CHECK_EQUAL(m1(2, 1), 0.0);
280 BOOST_CHECK_EQUAL((*m2)(2, 1), 0.0);
281 BOOST_CHECK_EQUAL((*m3)(2, 1), 0.0);
282
283 BOOST_CHECK_EQUAL(m1(2, 2), 5.7);
284 BOOST_CHECK_EQUAL((*m2)(2, 2), 5.7);
285 BOOST_CHECK_EQUAL((*m3)(2, 2), 5.7);
286
287#ifdef NEKTAR_FULLDEBUG
288 BOOST_CHECK_THROW(m1(0, 4), ErrorUtil::NekError);
289 BOOST_CHECK_THROW((*m2)(0, 4), ErrorUtil::NekError);
290 BOOST_CHECK_THROW((*m3)(0, 4), ErrorUtil::NekError);
291
292 BOOST_CHECK_THROW(m1(4, 0), ErrorUtil::NekError);
293 BOOST_CHECK_THROW((*m2)(4, 0), ErrorUtil::NekError);
294 BOOST_CHECK_THROW((*m3)(4, 0), ErrorUtil::NekError);
295
296 BOOST_CHECK_THROW(m1(4, 4), ErrorUtil::NekError);
297 BOOST_CHECK_THROW((*m2)(4, 4), ErrorUtil::NekError);
298 BOOST_CHECK_THROW((*m3)(4, 4), ErrorUtil::NekError);
299#endif // NEKTAR_FULLDEBUG
300}
Nektar::UnitTests::RedirectCerrIfNeeded
void RedirectCerrIfNeeded()
Definition: util.cpp:41

References Nektar::eDIAGONAL, and Nektar::UnitTests::RedirectCerrIfNeeded().

◆ BOOST_AUTO_TEST_CASE() [2/8]

Nektar::MatrixUnitTests::BOOST_AUTO_TEST_CASE ( TestDiagonalMatrixInversion  )

Definition at line 100 of file testNekMatrix.cpp.

101{
102 double buf[] = {1.0, 2.0, 3.0, 4.0};
103 NekMatrix<double> m(4, 4, buf, eDIAGONAL);
104 m.Invert();
105
106 BOOST_CHECK_EQUAL(m(0, 0), 1.0 / 1.0);
107 BOOST_CHECK_EQUAL(m(1, 1), 1.0 / 2.0);
108 BOOST_CHECK_EQUAL(m(2, 2), 1.0 / 3.0);
109 BOOST_CHECK_EQUAL(m(3, 3), 1.0 / 4.0);
110}

References Nektar::eDIAGONAL.

◆ BOOST_AUTO_TEST_CASE() [3/8]

Nektar::MatrixUnitTests::BOOST_AUTO_TEST_CASE ( TestDiagonalNekMatrixSetValue  )

Definition at line 337 of file testNekMatrix.cpp.

338{
339 UnitTests::RedirectCerrIfNeeded();
340 double data[] = {8.9, 3.4, 5.7};
341 Nektar::NekMatrix<double> m1(3, 3, data, eDIAGONAL);
342 std::shared_ptr<Nektar::Matrix<double>> m2(
343 new Nektar::NekMatrix<double>(3, 3, data, eDIAGONAL));
344 std::shared_ptr<Nektar::NekMatrix<double>> m3(
345 new Nektar::NekMatrix<double>(3, 3, data, eDIAGONAL));
346
347 m1.SetValue(0, 0, 1.0);
348 m2->SetValue(1, 1, 2.0);
349 m3->SetValue(2, 2, 3.0);
350
351 BOOST_CHECK_EQUAL(m1.GetStorageType(), eDIAGONAL);
352 BOOST_CHECK_EQUAL(m2->GetStorageType(), eDIAGONAL);
353 BOOST_CHECK_EQUAL(m3->GetStorageType(), eDIAGONAL);
354
355 BOOST_CHECK_EQUAL(m1(0, 0), 1.0);
356 BOOST_CHECK_EQUAL(m1(1, 1), 3.4);
357 BOOST_CHECK_EQUAL(m1(2, 2), 5.7);
358
359 BOOST_CHECK_EQUAL((*m2)(0, 0), 8.9);
360 BOOST_CHECK_EQUAL((*m2)(1, 1), 2.0);
361 BOOST_CHECK_EQUAL((*m2)(2, 2), 5.7);
362
363 BOOST_CHECK_EQUAL((*m3)(0, 0), 8.9);
364 BOOST_CHECK_EQUAL((*m3)(1, 1), 3.4);
365 BOOST_CHECK_EQUAL((*m3)(2, 2), 3.0);
366
367#ifdef NEKTAR_FULLDEBUG
368 BOOST_CHECK_THROW(m1.SetValue(0, 3, 2.0), ErrorUtil::NekError);
369 BOOST_CHECK_THROW(m1.SetValue(3, 0, 2.0), ErrorUtil::NekError);
370 BOOST_CHECK_THROW(m1.SetValue(3, 3, 2.0), ErrorUtil::NekError);
371
372 BOOST_CHECK_THROW(m2->SetValue(0, 3, 2.0), ErrorUtil::NekError);
373 BOOST_CHECK_THROW(m2->SetValue(3, 0, 2.0), ErrorUtil::NekError);
374 BOOST_CHECK_THROW(m2->SetValue(3, 3, 2.0), ErrorUtil::NekError);
375
376 BOOST_CHECK_THROW(m3->SetValue(0, 3, 2.0), ErrorUtil::NekError);
377 BOOST_CHECK_THROW(m3->SetValue(3, 0, 2.0), ErrorUtil::NekError);
378 BOOST_CHECK_THROW(m3->SetValue(3, 3, 2.0), ErrorUtil::NekError);
379#endif // NEKTAR_FULLDEBUG
380}

References Nektar::eDIAGONAL, and Nektar::UnitTests::RedirectCerrIfNeeded().

◆ BOOST_AUTO_TEST_CASE() [4/8]

Nektar::MatrixUnitTests::BOOST_AUTO_TEST_CASE ( TestFullFullMatrixAddition  )

Definition at line 382 of file testNekMatrix.cpp.

383{
384 double lhs_data[] = {1.0, 4.0, 2.0, 5.0, 3.0, 6.0};
385 double rhs_data[] = {10.0, 13.0, 11.0, 14.0, 12.0, 15.0};
386
387 Nektar::NekMatrix<double> lhs(2, 3, lhs_data);
388 Nektar::NekMatrix<double> rhs(2, 3, rhs_data);
389 Nektar::NekMatrix<double> result = lhs + rhs;
390
391 BOOST_CHECK_EQUAL(result.GetRows(), 2);
392 BOOST_CHECK_EQUAL(result.GetColumns(), 3);
393 BOOST_CHECK_EQUAL(result.GetStorageType(), eFULL);
394 BOOST_CHECK_EQUAL(result(0, 0), 11.0);
395 BOOST_CHECK_EQUAL(result(0, 1), 13.0);
396 BOOST_CHECK_EQUAL(result(0, 2), 15.0);
397 BOOST_CHECK_EQUAL(result(1, 0), 17.0);
398 BOOST_CHECK_EQUAL(result(1, 1), 19.0);
399 BOOST_CHECK_EQUAL(result(1, 2), 21.0);
400}

References Nektar::eFULL.

◆ BOOST_AUTO_TEST_CASE() [5/8]

Nektar::MatrixUnitTests::BOOST_AUTO_TEST_CASE ( TestFullMatrixInversion  )

Definition at line 45 of file testNekMatrix.cpp.

46{
47 {
48 double buf[] = {1.0, 3.0, 2.0, 4.0};
49 NekMatrix<double> m(2, 2, buf);
50 m.Invert();
51
52 BOOST_CHECK_CLOSE(*m(0, 0), -2.0, .00001);
53 BOOST_CHECK_CLOSE(*m(0, 1), 1.0, .00001);
54 BOOST_CHECK_CLOSE(*m(1, 0), 3.0 / 2.0, .00001);
55 BOOST_CHECK_CLOSE(*m(1, 1), -1.0 / 2.0, .00001);
56 }
57
58 {
59 double buf[] = {1.7, 5.6, 9.1, -3.4, -2.0, 4.5, 7.8, 8.2, -2.5,
60 5.0, 9.0, 3.6, 7.3, .98, 1.0, -12.6, 2.9, 6.4,
61 .02, -4.0, -1, 1.7, 5.6, 2.98, 5.0};
62 NekMatrix<double> m(5, 5, buf);
63 m.Invert();
64
65 BOOST_CHECK_CLOSE(*m(0, 0), 0.0005010995978, .000001);
66 BOOST_CHECK_CLOSE(*m(0, 1), -0.4704403712, .000001);
67 BOOST_CHECK_CLOSE(*m(0, 2), 0.2719063614, .000001);
68 BOOST_CHECK_CLOSE(*m(0, 3), -0.3941557805, .000001);
69 BOOST_CHECK_CLOSE(*m(0, 4), 0.09043166650, .000001);
70
71 BOOST_CHECK_CLOSE(*m(1, 0), -0.01087166322, .000001);
72 BOOST_CHECK_CLOSE(*m(1, 1), 0.3242048735, .000001);
73 BOOST_CHECK_CLOSE(*m(1, 2), -0.1605116333, .000001);
74 BOOST_CHECK_CLOSE(*m(1, 3), 0.09133673974, .000001);
75 BOOST_CHECK_CLOSE(*m(1, 4), 0.01293234281, .000001);
76
77 BOOST_CHECK_CLOSE(*m(2, 0), 0.06465598241, .000001);
78 BOOST_CHECK_CLOSE(*m(2, 1), 0.2593895742, .000001);
79 BOOST_CHECK_CLOSE(*m(2, 2), -0.09057233795, .000001);
80 BOOST_CHECK_CLOSE(*m(2, 3), 0.3106562459, .000001);
81 BOOST_CHECK_CLOSE(*m(2, 4), -0.1589713628, .000001);
82
83 BOOST_CHECK_CLOSE(*m(3, 0), -0.03464982687, .000001);
84 BOOST_CHECK_CLOSE(*m(3, 1), 0.2655177914, .000001);
85 BOOST_CHECK_CLOSE(*m(3, 2), -0.1016866584, .000001);
86 BOOST_CHECK_CLOSE(*m(3, 3), 0.2125363445, .000001);
87 BOOST_CHECK_CLOSE(*m(3, 4), -0.1099886183, .000001);
88
89 BOOST_CHECK_CLOSE(*m(4, 0), -0.02957895492, .000001);
90 BOOST_CHECK_CLOSE(*m(4, 1), -0.3518447037, .000001);
91 BOOST_CHECK_CLOSE(*m(4, 2), 0.2060393188, .000001);
92 BOOST_CHECK_CLOSE(*m(4, 3), -0.1411012255, .000001);
93 BOOST_CHECK_CLOSE(*m(4, 4), 0.1670437017, .000001);
94 }
95
96 {
97 }
98}

◆ BOOST_AUTO_TEST_CASE() [6/8]

Nektar::MatrixUnitTests::BOOST_AUTO_TEST_CASE ( TestFullNekMatrixGetValue  )

Definition at line 161 of file testNekMatrix.cpp.

162{
163 UnitTests::RedirectCerrIfNeeded();
164 double data[] = {1.0, 89.0, 0.0, 45.12, 2.0, -12.3,
165 892.2532, 76.12, 3.0, -56.7, 211.22, 45.23};
166
167 Nektar::NekMatrix<double> m(4, 3, data);
168
169 BOOST_CHECK_EQUAL(m(0, 0), 1.0);
170 BOOST_CHECK_EQUAL(m(0, 1), 2.0);
171 BOOST_CHECK_EQUAL(m(0, 2), 3.0);
172
173 BOOST_CHECK_EQUAL(m(1, 0), 89.0);
174 BOOST_CHECK_EQUAL(m(1, 1), -12.3);
175 BOOST_CHECK_EQUAL(m(1, 2), -56.7);
176
177 BOOST_CHECK_EQUAL(m(2, 0), 0.0);
178 BOOST_CHECK_EQUAL(m(2, 1), 892.2532);
179 BOOST_CHECK_EQUAL(m(2, 2), 211.22);
180
181 BOOST_CHECK_EQUAL(m(3, 0), 45.12);
182 BOOST_CHECK_EQUAL(m(3, 1), 76.12);
183 BOOST_CHECK_EQUAL(m(3, 2), 45.23);
184
185#ifdef NEKTAR_FULLDEBUG
186 BOOST_CHECK_THROW(m(0, 3), ErrorUtil::NekError);
187 BOOST_CHECK_THROW(m(4, 1), ErrorUtil::NekError);
188 BOOST_CHECK_THROW(m(4, 3), ErrorUtil::NekError);
189#endif // NEKTAR_FULLDEBUG
190
191 std::shared_ptr<Nektar::Matrix<double>> m1(
192 new Nektar::NekMatrix<double>(4, 3, data));
193 std::shared_ptr<Nektar::NekMatrix<double>> m2(
194 new Nektar::NekMatrix<double>(4, 3, data));
195
196 BOOST_CHECK_EQUAL((*m1)(0, 0), 1.0);
197 BOOST_CHECK_EQUAL((*m1)(0, 1), 2.0);
198 BOOST_CHECK_EQUAL((*m1)(0, 2), 3.0);
199
200 BOOST_CHECK_EQUAL((*m1)(1, 0), 89.0);
201 BOOST_CHECK_EQUAL((*m1)(1, 1), -12.3);
202 BOOST_CHECK_EQUAL((*m1)(1, 2), -56.7);
203
204 BOOST_CHECK_EQUAL((*m1)(2, 0), 0.0);
205 BOOST_CHECK_EQUAL((*m1)(2, 1), 892.2532);
206 BOOST_CHECK_EQUAL((*m1)(2, 2), 211.22);
207
208 BOOST_CHECK_EQUAL((*m1)(3, 0), 45.12);
209 BOOST_CHECK_EQUAL((*m1)(3, 1), 76.12);
210 BOOST_CHECK_EQUAL((*m1)(3, 2), 45.23);
211
212#ifdef NEKTAR_FULLDEBUG
213 BOOST_CHECK_THROW((*m1)(0, 3), ErrorUtil::NekError);
214 BOOST_CHECK_THROW((*m1)(4, 1), ErrorUtil::NekError);
215 BOOST_CHECK_THROW((*m1)(4, 3), ErrorUtil::NekError);
216#endif // NEKTAR_FULLDEBUG
217
218 BOOST_CHECK_EQUAL((*m2)(0, 0), 1.0);
219 BOOST_CHECK_EQUAL((*m2)(0, 1), 2.0);
220 BOOST_CHECK_EQUAL((*m2)(0, 2), 3.0);
221
222 BOOST_CHECK_EQUAL((*m2)(1, 0), 89.0);
223 BOOST_CHECK_EQUAL((*m2)(1, 1), -12.3);
224 BOOST_CHECK_EQUAL((*m2)(1, 2), -56.7);
225
226 BOOST_CHECK_EQUAL((*m2)(2, 0), 0.0);
227 BOOST_CHECK_EQUAL((*m2)(2, 1), 892.2532);
228 BOOST_CHECK_EQUAL((*m2)(2, 2), 211.22);
229
230 BOOST_CHECK_EQUAL((*m2)(3, 0), 45.12);
231 BOOST_CHECK_EQUAL((*m2)(3, 1), 76.12);
232 BOOST_CHECK_EQUAL((*m2)(3, 2), 45.23);
233
234#ifdef NEKTAR_FULLDEBUG
235 BOOST_CHECK_THROW((*m2)(0, 3), ErrorUtil::NekError);
236 BOOST_CHECK_THROW((*m2)(4, 1), ErrorUtil::NekError);
237 BOOST_CHECK_THROW((*m2)(4, 3), ErrorUtil::NekError);
238#endif // NEKTAR_FULLDEBUG
239}

References Nektar::UnitTests::RedirectCerrIfNeeded().

◆ BOOST_AUTO_TEST_CASE() [7/8]

Nektar::MatrixUnitTests::BOOST_AUTO_TEST_CASE ( TestFullNekMatrixSetValue  )

Definition at line 302 of file testNekMatrix.cpp.

303{
304 double data[] = {1.0, 3.0, 5.0, 2.0, 4.0, 6.0};
305 Nektar::NekMatrix<double> m1(3, 2, data);
306 std::shared_ptr<Nektar::Matrix<double>> m2(
307 new Nektar::NekMatrix<double>(3, 2, data));
308 std::shared_ptr<Nektar::NekMatrix<double>> m3(
309 new Nektar::NekMatrix<double>(3, 2, data));
310
311 m1.SetValue(0, 0, -1.0);
312 m2->SetValue(1, 1, -2.0);
313 m3->SetValue(2, 1, -3.0);
314
315 BOOST_CHECK_EQUAL(m1(0, 0), -1.0);
316 BOOST_CHECK_EQUAL(m1(0, 1), 2.0);
317 BOOST_CHECK_EQUAL(m1(1, 0), 3.0);
318 BOOST_CHECK_EQUAL(m1(1, 1), 4.0);
319 BOOST_CHECK_EQUAL(m1(2, 0), 5.0);
320 BOOST_CHECK_EQUAL(m1(2, 1), 6.0);
321
322 BOOST_CHECK_EQUAL((*m2)(0, 0), 1.0);
323 BOOST_CHECK_EQUAL((*m2)(0, 1), 2.0);
324 BOOST_CHECK_EQUAL((*m2)(1, 0), 3.0);
325 BOOST_CHECK_EQUAL((*m2)(1, 1), -2.0);
326 BOOST_CHECK_EQUAL((*m2)(2, 0), 5.0);
327 BOOST_CHECK_EQUAL((*m2)(2, 1), 6.0);
328
329 BOOST_CHECK_EQUAL((*m3)(0, 0), 1.0);
330 BOOST_CHECK_EQUAL((*m3)(0, 1), 2.0);
331 BOOST_CHECK_EQUAL((*m3)(1, 0), 3.0);
332 BOOST_CHECK_EQUAL((*m3)(1, 1), 4.0);
333 BOOST_CHECK_EQUAL((*m3)(2, 0), 5.0);
334 BOOST_CHECK_EQUAL((*m3)(2, 1), -3.0);
335}

◆ BOOST_AUTO_TEST_CASE() [8/8]

Nektar::MatrixUnitTests::BOOST_AUTO_TEST_CASE ( TestNekMatrixConstruction  )

Definition at line 112 of file testNekMatrix.cpp.

113{
114 {
115 std::shared_ptr<Nektar::Matrix<double>> a(
116 new Nektar::NekMatrix<double>(3, 4));
117 std::shared_ptr<Nektar::NekMatrix<double>> b(
118 new Nektar::NekMatrix<double>(5, 6));
119
120 BOOST_CHECK_EQUAL(a->GetRows(), 3u);
121 BOOST_CHECK_EQUAL(a->GetColumns(), 4u);
122 BOOST_CHECK_EQUAL(a->GetSize()[0], 3u);
123 BOOST_CHECK_EQUAL(a->GetSize()[1], 4u);
124
125 BOOST_CHECK_EQUAL(b->GetRows(), 5u);
126 BOOST_CHECK_EQUAL(b->GetColumns(), 6u);
127 BOOST_CHECK_EQUAL(b->GetSize()[0], 5u);
128 BOOST_CHECK_EQUAL(b->GetSize()[1], 6u);
129
130 BOOST_CHECK_EQUAL(a->GetStorageType(), eFULL);
131 BOOST_CHECK_EQUAL(b->GetStorageType(), eFULL);
132
133 BOOST_CHECK_EQUAL(a->GetStorageSize(), 3 * 4);
134 BOOST_CHECK_EQUAL(b->GetStorageSize(), 5 * 6);
135 }
136
137 {
138 std::shared_ptr<Nektar::Matrix<double>> a(
139 new Nektar::NekMatrix<double>(3, 3, eDIAGONAL));
140 std::shared_ptr<Nektar::NekMatrix<double>> b(
141 new Nektar::NekMatrix<double>(5, 5, eDIAGONAL));
142
143 BOOST_CHECK_EQUAL(a->GetRows(), 3);
144 BOOST_CHECK_EQUAL(a->GetColumns(), 3);
145 BOOST_CHECK_EQUAL(a->GetSize()[0], 3);
146 BOOST_CHECK_EQUAL(a->GetSize()[1], 3);
147
148 BOOST_CHECK_EQUAL(b->GetRows(), 5);
149 BOOST_CHECK_EQUAL(b->GetColumns(), 5);
150 BOOST_CHECK_EQUAL(b->GetSize()[0], 5);
151 BOOST_CHECK_EQUAL(b->GetSize()[1], 5);
152
153 BOOST_CHECK_EQUAL(a->GetStorageType(), eDIAGONAL);
154 BOOST_CHECK_EQUAL(b->GetStorageType(), eDIAGONAL);
155
156 BOOST_CHECK_EQUAL(a->GetStorageSize(), 3);
157 BOOST_CHECK_EQUAL(b->GetStorageSize(), 5);
158 }
159}

References Nektar::eDIAGONAL, and Nektar::eFULL.

Generated on Wed Jan 31 2024 09:36:04 for Nektar++ by doxygen 1.9.4