preCICE v3.1.2
Loading...
Searching...
No Matches
AccelerationIntraCommTest.cpp
Go to the documentation of this file.
1#include <Eigen/Core>
2#include <algorithm>
3#include <cmath>
4#include <map>
5#include <memory>
6#include <utility>
7#include <vector>
8#include "acceleration/Acceleration.hpp"
9#include "acceleration/BaseQNAcceleration.hpp"
10#include "acceleration/IQNILSAcceleration.hpp"
11#include "acceleration/IQNIMVJAcceleration.hpp"
12#include "acceleration/SharedPointer.hpp"
13#include "acceleration/impl/ConstantPreconditioner.hpp"
14#include "acceleration/impl/QRFactorization.hpp"
15#include "acceleration/impl/ResidualSumPreconditioner.hpp"
16#include "acceleration/impl/SharedPointer.hpp"
17#include "acceleration/test/helper.hpp"
18#include "cplscheme/Constants.hpp"
19#include "cplscheme/CouplingData.hpp"
20#include "cplscheme/SharedPointer.hpp"
21#include "mesh/Mesh.hpp"
22#include "mesh/SharedPointer.hpp"
23#include "testing/TestContext.hpp"
24#include "testing/Testing.hpp"
25#include "utils/EigenHelperFunctions.hpp"
26
27using namespace precice;
28using namespace precice::cplscheme;
29using namespace precice::acceleration;
30using namespace precice::acceleration::impl;
31
32using precice::testing::makeCouplingData;
33
34#ifndef PRECICE_NO_MPI
35
36BOOST_AUTO_TEST_SUITE(AccelerationTests)
37
38using DataMap = std::map<int, PtrCouplingData>;
39
40BOOST_AUTO_TEST_SUITE(AccelerationIntraCommTests)
41
43BOOST_AUTO_TEST_CASE(testVIQNILSppWithoutSubsteps)
44{
45 PRECICE_TEST(""_on(4_ranks).setupIntraComm());
46 double initialRelaxation = 0.01;
47 int maxIterationsUsed = 50;
48 int timeWindowsReused = 6;
49 int filter = BaseQNAcceleration::QR1FILTER;
50 double singularityLimit = 1e-10;
51 bool enforceInitialRelaxation = false;
52 std::vector<int> dataIDs;
53 dataIDs.push_back(0);
54 dataIDs.push_back(1);
55 std::vector<double> factors;
56 factors.resize(2, 1.0);
57 PtrPreconditioner prec(new ConstantPreconditioner(factors));
58 std::vector<int> vertexOffsets{4, 8, 8, 10};
59
60 mesh::PtrMesh dummyMesh(new mesh::Mesh("DummyMesh", 3, testing::nextMeshID()));
61 dummyMesh->setVertexOffsets(vertexOffsets);
62
63 IQNILSAcceleration pp(initialRelaxation, enforceInitialRelaxation, maxIterationsUsed,
64 timeWindowsReused, filter, singularityLimit, dataIDs, prec);
65
66 Eigen::VectorXd dcol1;
67 Eigen::VectorXd fcol1;
68
69 DataMap data;
70 mesh::PtrData displacements(new mesh::Data("dvalues", -1, 1));
71 mesh::PtrData forces(new mesh::Data("fvalues", -1, 1));
72
73 bool exchangeSubsteps = false; // @todo add testVIQNILSppWithSubsteps, where exchangeSubsteps = true as soon as acceleration scheme supports subcycling.
74
75 if (context.isPrimary()) { // Primary
80 // init displacements
81 Eigen::VectorXd insert(4);
82 insert << 1.0, 1.0, 1.0, 1.0;
83 utils::append(displacements->values(), insert);
84 // init forces
85 insert << 0.2, 0.2, 0.2, 0.2;
86 utils::append(forces->values(), insert);
87
88 PtrCouplingData dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
89 dpcd->setSampleAtTime(0, dpcd->sample());
90 PtrCouplingData fpcd = makeCouplingData(forces, dummyMesh, exchangeSubsteps);
91 fpcd->setSampleAtTime(0, fpcd->sample());
92
93 data.insert(std::pair<int, PtrCouplingData>(0, dpcd));
94 data.insert(std::pair<int, PtrCouplingData>(1, fpcd));
95
96 dpcd->storeIteration();
97 fpcd->storeIteration();
98
99 pp.initialize(data);
100
101 insert << 1.0, 2.0, 3.0, 4.0;
102 displacements->values() = insert;
103 insert << 0.1, 0.1, 0.1, 0.1;
104 forces->values() = insert;
105
106 } else if (context.isRank(1)) { // SecondaryRank1
107
112 // init displacements
113 Eigen::VectorXd insert(4);
114 insert << 1.0, 1.0, 1.0, 1.0;
115 utils::append(displacements->values(), insert);
116 // init forces
117 insert << 0.2, 0.2, 0.2, 0.2;
118 utils::append(forces->values(), insert);
119
120 PtrCouplingData dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
121 dpcd->setSampleAtTime(0, dpcd->sample());
122 PtrCouplingData fpcd = makeCouplingData(forces, dummyMesh, exchangeSubsteps);
123 fpcd->setSampleAtTime(0, fpcd->sample());
124
125 data.insert(std::pair<int, PtrCouplingData>(0, dpcd));
126 data.insert(std::pair<int, PtrCouplingData>(1, fpcd));
127
128 dpcd->storeIteration();
129 fpcd->storeIteration();
130
131 pp.initialize(data);
132
133 insert << 5.0, 6.0, 7.0, 8.0;
134 displacements->values() = insert;
135 insert << 0.1, 0.1, 0.1, 0.1;
136 forces->values() = insert;
137
138 } else if (context.isRank(2)) { // Secondary rank 2
139
144 // init displacements
145 PtrCouplingData dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
146 dpcd->setSampleAtTime(0, dpcd->sample());
147
148 // init forces
149 PtrCouplingData fpcd = makeCouplingData(forces, dummyMesh, exchangeSubsteps);
150 fpcd->setSampleAtTime(0, fpcd->sample());
151
152 data.insert(std::pair<int, PtrCouplingData>(0, dpcd));
153 data.insert(std::pair<int, PtrCouplingData>(1, fpcd));
154
155 dpcd->storeIteration();
156 fpcd->storeIteration();
157
158 pp.initialize(data);
159
160 } else if (context.isRank(3)) { // Secondary rank 3
161
166 // init displacements
167 Eigen::VectorXd insert(2);
168 insert << 1.0, 1.0;
169 utils::append(displacements->values(), insert);
170 // init forces
171 insert << 0.2, 0.2;
172 utils::append(forces->values(), insert);
173
174 PtrCouplingData dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
175 dpcd->setSampleAtTime(0, dpcd->sample());
176 PtrCouplingData fpcd = makeCouplingData(forces, dummyMesh, exchangeSubsteps);
177 fpcd->setSampleAtTime(0, fpcd->sample());
178
179 data.insert(std::pair<int, PtrCouplingData>(0, dpcd));
180 data.insert(std::pair<int, PtrCouplingData>(1, fpcd));
181
182 dpcd->storeIteration();
183 fpcd->storeIteration();
184
185 pp.initialize(data);
186
187 insert << 1.0, 2.0;
188 displacements->values() = insert;
189 insert << 0.1, 0.1;
190 forces->values() = insert;
191 }
192
193 pp.performAcceleration(data);
194
195 Eigen::VectorXd newdvalues;
196 if (context.isPrimary()) { // Primary
197
198 BOOST_TEST(testing::equals(data.at(0)->values()(0), 1.00), data.at(0)->values()(0));
199 BOOST_TEST(testing::equals(data.at(0)->values()(1), 1.01), data.at(0)->values()(1));
200 BOOST_TEST(testing::equals(data.at(0)->values()(2), 1.02), data.at(0)->values()(2));
201 BOOST_TEST(testing::equals(data.at(0)->values()(3), 1.03), data.at(0)->values()(3));
202 BOOST_TEST(testing::equals(data.at(1)->values()(0), 0.199), data.at(1)->values()(0));
203 BOOST_TEST(testing::equals(data.at(1)->values()(1), 0.199), data.at(1)->values()(1));
204 BOOST_TEST(testing::equals(data.at(1)->values()(2), 0.199), data.at(1)->values()(2));
205 BOOST_TEST(testing::equals(data.at(1)->values()(3), 0.199), data.at(1)->values()(3));
206
207 utils::append(newdvalues, 10.0);
208 utils::append(newdvalues, 10.0);
209 utils::append(newdvalues, 10.0);
210 utils::append(newdvalues, 10.0);
211
212 } else if (context.isRank(1)) { // SecondaryRank1
213
214 BOOST_TEST(testing::equals(data.at(0)->values()(0), 1.04), data.at(0)->values()(0));
215 BOOST_TEST(testing::equals(data.at(0)->values()(1), 1.05), data.at(0)->values()(1));
216 BOOST_TEST(testing::equals(data.at(0)->values()(2), 1.06), data.at(0)->values()(2));
217 BOOST_TEST(testing::equals(data.at(0)->values()(3), 1.07), data.at(0)->values()(3));
218 BOOST_TEST(testing::equals(data.at(1)->values()(0), 0.199), data.at(1)->values()(0));
219 BOOST_TEST(testing::equals(data.at(1)->values()(1), 0.199), data.at(1)->values()(1));
220 BOOST_TEST(testing::equals(data.at(1)->values()(2), 0.199), data.at(1)->values()(2));
221 BOOST_TEST(testing::equals(data.at(1)->values()(3), 0.199), data.at(1)->values()(3));
222
223 utils::append(newdvalues, 10.0);
224 utils::append(newdvalues, 10.0);
225 utils::append(newdvalues, 10.0);
226 utils::append(newdvalues, 10.0);
227
228 } else if (context.isRank(2)) { // Secondary rank 2
229 // empty proc
230 } else if (context.isRank(3)) { // Secondary rank 3
231
232 BOOST_TEST(testing::equals(data.at(0)->values()(0), 1.00), data.at(0)->values()(0));
233 BOOST_TEST(testing::equals(data.at(0)->values()(1), 1.01), data.at(0)->values()(1));
234 BOOST_TEST(testing::equals(data.at(1)->values()(0), 0.199), data.at(1)->values()(0));
235 BOOST_TEST(testing::equals(data.at(1)->values()(1), 0.199), data.at(1)->values()(1));
236
237 utils::append(newdvalues, 10.0);
238 utils::append(newdvalues, 10.0);
239 }
240
241 data.begin()->second->values() = newdvalues;
242
243 pp.performAcceleration(data);
244
245 if (context.isPrimary()) { // Primary
246 BOOST_TEST(testing::equals(data.at(0)->values()(0), -1.51483105223442748866e+00), data.at(0)->values()(0));
247 BOOST_TEST(testing::equals(data.at(0)->values()(1), -2.35405379763935940218e-01), data.at(0)->values()(1));
248 BOOST_TEST(testing::equals(data.at(0)->values()(2), 1.04402029270655560822e+00), data.at(0)->values()(2));
249 BOOST_TEST(testing::equals(data.at(0)->values()(3), 2.32344596517704804484e+00), data.at(0)->values()(3));
250 BOOST_TEST(testing::equals(data.at(1)->values()(0), 7.23368584254212854123e-02), data.at(1)->values()(0));
251 BOOST_TEST(testing::equals(data.at(1)->values()(1), 7.23368584254212854123e-02), data.at(1)->values()(1));
252 BOOST_TEST(testing::equals(data.at(1)->values()(2), 7.23368584254212854123e-02), data.at(1)->values()(2));
253 BOOST_TEST(testing::equals(data.at(1)->values()(3), 7.23368584254212854123e-02), data.at(1)->values()(3));
254 } else if (context.isRank(1)) { // SecondaryRank1
255 BOOST_TEST(testing::equals(data.at(0)->values()(0), 3.60287163764754048145e+00), data.at(0)->values()(0));
256 BOOST_TEST(testing::equals(data.at(0)->values()(1), 4.88229731011803202989e+00), data.at(0)->values()(1));
257 BOOST_TEST(testing::equals(data.at(0)->values()(2), 6.16172298258852357833e+00), data.at(0)->values()(2));
258 BOOST_TEST(testing::equals(data.at(0)->values()(3), 7.44114865505901601495e+00), data.at(0)->values()(3));
259 BOOST_TEST(testing::equals(data.at(1)->values()(0), 7.23368584254212854123e-02), data.at(1)->values()(0));
260 BOOST_TEST(testing::equals(data.at(1)->values()(1), 7.23368584254212854123e-02), data.at(1)->values()(1));
261 BOOST_TEST(testing::equals(data.at(1)->values()(2), 7.23368584254212854123e-02), data.at(1)->values()(2));
262 BOOST_TEST(testing::equals(data.at(1)->values()(3), 7.23368584254212854123e-02), data.at(1)->values()(3));
263 } else if (context.isRank(2)) { // Secondary rank 2
264 // empty proc
265 } else if (context.isRank(3)) { // Secondary rank 3
266 BOOST_TEST(testing::equals(data.at(0)->values()(0), -1.51483105223442748866e+00), data.at(0)->values()(0));
267 BOOST_TEST(testing::equals(data.at(0)->values()(1), -2.35405379763935940218e-01), data.at(0)->values()(1));
268 BOOST_TEST(testing::equals(data.at(1)->values()(0), 7.23368584254212854123e-02), data.at(1)->values()(0));
269 BOOST_TEST(testing::equals(data.at(1)->values()(1), 7.23368584254212854123e-02), data.at(1)->values()(1));
270 }
271}
272
274BOOST_AUTO_TEST_CASE(testVIQNIMVJppWithoutSubsteps)
275{
276 PRECICE_TEST(""_on(4_ranks).setupIntraComm());
277 double initialRelaxation = 0.01;
278 int maxIterationsUsed = 50;
279 int timeWindowsReused = 6;
280 int filter = BaseQNAcceleration::QR1FILTER;
281 int restartType = IQNIMVJAcceleration::NO_RESTART;
282 int chunkSize = 0;
283 int reusedTimeWindowsAtRestart = 0;
284 double singularityLimit = 1e-10;
285 double svdTruncationEps = 0.0;
286 bool enforceInitialRelaxation = false;
287 bool alwaysBuildJacobian = false;
288
289 std::vector<int> dataIDs;
290 dataIDs.push_back(0);
291 dataIDs.push_back(1);
292 std::vector<double> factors;
293 factors.resize(2, 1.0);
294 PtrPreconditioner prec(new ConstantPreconditioner(factors));
295 std::vector<int> vertexOffsets{4, 8, 8, 10};
296
297 mesh::PtrMesh dummyMesh(new mesh::Mesh("DummyMesh", 3, testing::nextMeshID()));
298 dummyMesh->setVertexOffsets(vertexOffsets);
299
300 IQNIMVJAcceleration pp(initialRelaxation, enforceInitialRelaxation, maxIterationsUsed,
301 timeWindowsReused, filter, singularityLimit, dataIDs, prec, alwaysBuildJacobian,
302 restartType, chunkSize, reusedTimeWindowsAtRestart, svdTruncationEps);
303
304 Eigen::VectorXd dcol1;
305 Eigen::VectorXd fcol1;
306
307 mesh::PtrData displacements(new mesh::Data("dvalues", -1, 1));
308 mesh::PtrData forces(new mesh::Data("fvalues", -1, 1));
309
310 DataMap data;
311
312 bool exchangeSubsteps = false; // @todo add testVIQNIMVJppWithSubsteps, where exchangeSubsteps = true as soon as acceleration scheme supports subcycling.
313
314 if (context.isPrimary()) { // Primary
315
320 // init displacements
321 Eigen::VectorXd insert(4);
322 insert << 1.0, 1.0, 1.0, 1.0;
323 utils::append(displacements->values(), insert);
324 // init forces
325 insert << 0.2, 0.2, 0.2, 0.2;
326 utils::append(forces->values(), insert);
327
328 PtrCouplingData dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
329 dpcd->setSampleAtTime(0, dpcd->sample());
330 PtrCouplingData fpcd = makeCouplingData(forces, dummyMesh, exchangeSubsteps);
331 fpcd->setSampleAtTime(0, fpcd->sample());
332
333 dpcd->storeIteration();
334 fpcd->storeIteration();
335
336 data.insert(std::pair<int, PtrCouplingData>(0, dpcd));
337 data.insert(std::pair<int, PtrCouplingData>(1, fpcd));
338
339 pp.initialize(data);
340
341 // update displacement
342 insert << 1.0, 2.0, 3.0, 4.0;
343 dpcd->values() = insert;
344
345 // update forces
346 insert << 0.1, 0.1, 0.1, 0.1;
347 fpcd->values() = insert;
348
349 // check for correct initial data
350 BOOST_TEST(testing::equals(data.at(0)->previousIteration()(0), 1.0), data.at(0)->previousIteration()(0));
351 BOOST_TEST(testing::equals(data.at(0)->previousIteration()(1), 1.0), data.at(0)->previousIteration()(1));
352 BOOST_TEST(testing::equals(data.at(0)->previousIteration()(2), 1.0), data.at(0)->previousIteration()(2));
353 BOOST_TEST(testing::equals(data.at(0)->previousIteration()(3), 1.0), data.at(0)->previousIteration()(3));
354 BOOST_TEST(testing::equals(data.at(0)->values()(0), 1.0), data.at(0)->values()(0));
355 BOOST_TEST(testing::equals(data.at(0)->values()(1), 2.0), data.at(0)->values()(1));
356 BOOST_TEST(testing::equals(data.at(0)->values()(2), 3.0), data.at(0)->values()(2));
357 BOOST_TEST(testing::equals(data.at(0)->values()(3), 4.0), data.at(0)->values()(3));
358 BOOST_TEST(testing::equals(data.at(1)->previousIteration()(0), 0.2), data.at(1)->previousIteration()(0));
359 BOOST_TEST(testing::equals(data.at(1)->previousIteration()(1), 0.2), data.at(1)->previousIteration()(1));
360 BOOST_TEST(testing::equals(data.at(1)->previousIteration()(2), 0.2), data.at(1)->previousIteration()(2));
361 BOOST_TEST(testing::equals(data.at(1)->previousIteration()(3), 0.2), data.at(1)->previousIteration()(3));
362 BOOST_TEST(testing::equals(data.at(1)->values()(0), 0.1), data.at(1)->values()(0));
363 BOOST_TEST(testing::equals(data.at(1)->values()(1), 0.1), data.at(1)->values()(1));
364 BOOST_TEST(testing::equals(data.at(1)->values()(2), 0.1), data.at(1)->values()(2));
365 BOOST_TEST(testing::equals(data.at(1)->values()(3), 0.1), data.at(1)->values()(3));
366 } else if (context.isRank(1)) { // SecondaryRank1
367
372 // init displacements
373 Eigen::VectorXd insert(4);
374 insert << 1.0, 1.0, 1.0, 1.0;
375 utils::append(displacements->values(), insert);
376 // init forces
377 insert << 0.2, 0.2, 0.2, 0.2;
378 utils::append(forces->values(), insert);
379
380 PtrCouplingData dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
381 dpcd->setSampleAtTime(0, dpcd->sample());
382 PtrCouplingData fpcd = makeCouplingData(forces, dummyMesh, exchangeSubsteps);
383 fpcd->setSampleAtTime(0, fpcd->sample());
384
385 dpcd->storeIteration();
386 fpcd->storeIteration();
387
388 data.insert(std::pair<int, PtrCouplingData>(0, dpcd));
389 data.insert(std::pair<int, PtrCouplingData>(1, fpcd));
390
391 pp.initialize(data);
392
393 // update displacements
394 insert << 5.0, 6.0, 7.0, 8.0;
395 dpcd->values() = insert;
396 // update forces
397 insert << 0.1, 0.1, 0.1, 0.1;
398 fpcd->values() = insert;
399
400 // check for correct initial data
401 BOOST_TEST(testing::equals(data.at(0)->previousIteration()(0), 1.0), data.at(0)->previousIteration()(0));
402 BOOST_TEST(testing::equals(data.at(0)->previousIteration()(1), 1.0), data.at(0)->previousIteration()(1));
403 BOOST_TEST(testing::equals(data.at(0)->previousIteration()(2), 1.0), data.at(0)->previousIteration()(2));
404 BOOST_TEST(testing::equals(data.at(0)->previousIteration()(3), 1.0), data.at(0)->previousIteration()(3));
405 BOOST_TEST(testing::equals(data.at(0)->values()(0), 5.0), data.at(0)->values()(0));
406 BOOST_TEST(testing::equals(data.at(0)->values()(1), 6.0), data.at(0)->values()(1));
407 BOOST_TEST(testing::equals(data.at(0)->values()(2), 7.0), data.at(0)->values()(2));
408 BOOST_TEST(testing::equals(data.at(0)->values()(3), 8.0), data.at(0)->values()(3));
409 BOOST_TEST(testing::equals(data.at(1)->previousIteration()(0), 0.2), data.at(1)->previousIteration()(0));
410 BOOST_TEST(testing::equals(data.at(1)->previousIteration()(1), 0.2), data.at(1)->previousIteration()(1));
411 BOOST_TEST(testing::equals(data.at(1)->previousIteration()(2), 0.2), data.at(1)->previousIteration()(2));
412 BOOST_TEST(testing::equals(data.at(1)->previousIteration()(3), 0.2), data.at(1)->previousIteration()(3));
413 BOOST_TEST(testing::equals(data.at(1)->values()(0), 0.1), data.at(1)->values()(0));
414 BOOST_TEST(testing::equals(data.at(1)->values()(1), 0.1), data.at(1)->values()(1));
415 BOOST_TEST(testing::equals(data.at(1)->values()(2), 0.1), data.at(1)->values()(2));
416 BOOST_TEST(testing::equals(data.at(1)->values()(3), 0.1), data.at(1)->values()(3));
417 } else if (context.isRank(2)) { // Secondary rank 2
418
423 // init displacements
424 PtrCouplingData dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
425 dpcd->setSampleAtTime(0, dpcd->sample());
426
427 // init forces
428 PtrCouplingData fpcd = makeCouplingData(forces, dummyMesh, exchangeSubsteps);
429 fpcd->setSampleAtTime(0, fpcd->sample());
430
431 dpcd->storeIteration();
432 fpcd->storeIteration();
433
434 data.insert(std::pair<int, PtrCouplingData>(0, dpcd));
435 data.insert(std::pair<int, PtrCouplingData>(1, fpcd));
436
437 pp.initialize(data);
438
439 // check for correct initial data
440 // empty proc
441 } else if (context.isRank(3)) { // Secondary rank 3
442
447 // init displacements
448 Eigen::VectorXd insert(2);
449 insert << 1.0, 1.0;
450 utils::append(displacements->values(), insert);
451 // init forces
452 insert << 0.2, 0.2;
453 utils::append(forces->values(), insert);
454
455 PtrCouplingData dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
456 dpcd->setSampleAtTime(0, dpcd->sample());
457 PtrCouplingData fpcd = makeCouplingData(forces, dummyMesh, exchangeSubsteps);
458 fpcd->setSampleAtTime(0, fpcd->sample());
459
460 dpcd->storeIteration();
461 fpcd->storeIteration();
462
463 data.insert(std::pair<int, PtrCouplingData>(0, dpcd));
464 data.insert(std::pair<int, PtrCouplingData>(1, fpcd));
465
466 pp.initialize(data);
467
468 // update displacements
469 insert << 1.0, 2.0;
470 displacements->values() = insert;
471
472 // update forces
473 insert << 0.1, 0.1;
474 forces->values() = insert;
475
476 // check for correct initial data
477 BOOST_TEST(testing::equals(data.at(0)->previousIteration()(0), 1.0), data.at(0)->previousIteration()(0));
478 BOOST_TEST(testing::equals(data.at(0)->previousIteration()(1), 1.0), data.at(0)->previousIteration()(1));
479 BOOST_TEST(testing::equals(data.at(0)->values()(0), 1.0), data.at(0)->values()(0));
480 BOOST_TEST(testing::equals(data.at(0)->values()(1), 2.0), data.at(0)->values()(1));
481 BOOST_TEST(testing::equals(data.at(1)->previousIteration()(0), 0.2), data.at(1)->previousIteration()(0));
482 BOOST_TEST(testing::equals(data.at(1)->previousIteration()(1), 0.2), data.at(1)->previousIteration()(1));
483 BOOST_TEST(testing::equals(data.at(1)->values()(0), 0.1), data.at(1)->values()(0));
484 BOOST_TEST(testing::equals(data.at(1)->values()(1), 0.1), data.at(1)->values()(1));
485 }
486
487 pp.performAcceleration(data);
488
489 Eigen::VectorXd newdvalues;
490 if (context.isPrimary()) { // Primary
491 BOOST_TEST(testing::equals(data.at(0)->values()(0), 1.00000000000000000000e+00), data.at(0)->values()(0));
492 BOOST_TEST(testing::equals(data.at(0)->values()(1), 1.01000000000000000888e+00), data.at(0)->values()(1));
493 BOOST_TEST(testing::equals(data.at(0)->values()(2), 1.02000000000000001776e+00), data.at(0)->values()(2));
494 BOOST_TEST(testing::equals(data.at(0)->values()(3), 1.03000000000000002665e+00), data.at(0)->values()(3));
495 BOOST_TEST(testing::equals(data.at(1)->values()(0), 1.99000000000000010214e-01), data.at(1)->values()(0));
496 BOOST_TEST(testing::equals(data.at(1)->values()(1), 1.99000000000000010214e-01), data.at(1)->values()(1));
497 BOOST_TEST(testing::equals(data.at(1)->values()(2), 1.99000000000000010214e-01), data.at(1)->values()(2));
498 BOOST_TEST(testing::equals(data.at(1)->values()(3), 1.99000000000000010214e-01), data.at(1)->values()(3));
499 utils::append(newdvalues, 10.0);
500 utils::append(newdvalues, 10.0);
501 utils::append(newdvalues, 10.0);
502 utils::append(newdvalues, 10.0);
503 } else if (context.isRank(1)) { // SecondaryRank1
504 BOOST_TEST(testing::equals(data.at(0)->values()(0), 1.04000000000000003553e+00), data.at(0)->values()(0));
505 BOOST_TEST(testing::equals(data.at(0)->values()(1), 1.05000000000000004441e+00), data.at(0)->values()(1));
506 BOOST_TEST(testing::equals(data.at(0)->values()(2), 1.06000000000000005329e+00), data.at(0)->values()(2));
507 BOOST_TEST(testing::equals(data.at(0)->values()(3), 1.07000000000000006217e+00), data.at(0)->values()(3));
508 BOOST_TEST(testing::equals(data.at(1)->values()(0), 1.99000000000000010214e-01), data.at(1)->values()(0));
509 BOOST_TEST(testing::equals(data.at(1)->values()(1), 1.99000000000000010214e-01), data.at(1)->values()(1));
510 BOOST_TEST(testing::equals(data.at(1)->values()(2), 1.99000000000000010214e-01), data.at(1)->values()(2));
511 BOOST_TEST(testing::equals(data.at(1)->values()(3), 1.99000000000000010214e-01), data.at(1)->values()(3));
512 utils::append(newdvalues, 10.0);
513 utils::append(newdvalues, 10.0);
514 utils::append(newdvalues, 10.0);
515 utils::append(newdvalues, 10.0);
516 } else if (context.isRank(2)) { // Secondary rank 2
517 // empty proc
518 } else if (context.isRank(3)) { // Secondary rank 3
519 BOOST_TEST(testing::equals(data.at(0)->values()(0), 1.00000000000000000000e+00), data.at(0)->values()(0));
520 BOOST_TEST(testing::equals(data.at(0)->values()(1), 1.01000000000000000888e+00), data.at(0)->values()(1));
521 BOOST_TEST(testing::equals(data.at(1)->values()(0), 1.99000000000000010214e-01), data.at(1)->values()(0));
522 BOOST_TEST(testing::equals(data.at(1)->values()(1), 1.99000000000000010214e-01), data.at(1)->values()(1));
523 utils::append(newdvalues, 10.0);
524 utils::append(newdvalues, 10.0);
525 }
526
527 data.begin()->second->values() = newdvalues;
528 pp.performAcceleration(data);
529
530 if (context.isPrimary()) { // Primary
531 BOOST_TEST(testing::equals(data.at(0)->values()(0), -1.51483105223442748866e+00), data.at(0)->values()(0));
532 BOOST_TEST(testing::equals(data.at(0)->values()(1), -2.35405379763935940218e-01), data.at(0)->values()(1));
533 BOOST_TEST(testing::equals(data.at(0)->values()(2), 1.04402029270655738458e+00), data.at(0)->values()(2));
534 BOOST_TEST(testing::equals(data.at(0)->values()(3), 2.32344596517704893301e+00), data.at(0)->values()(3));
535 BOOST_TEST(testing::equals(data.at(1)->values()(0), 7.23368584254213131679e-02), data.at(1)->values()(0));
536 BOOST_TEST(testing::equals(data.at(1)->values()(1), 7.23368584254213131679e-02), data.at(1)->values()(1));
537 BOOST_TEST(testing::equals(data.at(1)->values()(2), 7.23368584254213131679e-02), data.at(1)->values()(2));
538 BOOST_TEST(testing::equals(data.at(1)->values()(3), 7.23368584254213131679e-02), data.at(1)->values()(3));
539 } else if (context.isRank(1)) { // SecondaryRank1
540 BOOST_TEST(testing::equals(data.at(0)->values()(0), 3.60287163764754048145e+00), data.at(0)->values()(0));
541 BOOST_TEST(testing::equals(data.at(0)->values()(1), 4.88229731011803202989e+00), data.at(0)->values()(1));
542 BOOST_TEST(testing::equals(data.at(0)->values()(2), 6.16172298258852446651e+00), data.at(0)->values()(2));
543 BOOST_TEST(testing::equals(data.at(0)->values()(3), 7.44114865505901601495e+00), data.at(0)->values()(3));
544 BOOST_TEST(testing::equals(data.at(1)->values()(0), 7.23368584254213131679e-02), data.at(1)->values()(0));
545 BOOST_TEST(testing::equals(data.at(1)->values()(1), 7.23368584254213131679e-02), data.at(1)->values()(1));
546 BOOST_TEST(testing::equals(data.at(1)->values()(2), 7.23368584254213131679e-02), data.at(1)->values()(2));
547 BOOST_TEST(testing::equals(data.at(1)->values()(3), 7.23368584254213131679e-02), data.at(1)->values()(3));
548 } else if (context.isRank(2)) { // Secondary rank 2
549 // empty proc
550 } else if (context.isRank(3)) { // Secondary rank 3
551 BOOST_TEST(testing::equals(data.at(0)->values()(0), -1.51483105223442748866e+00), data.at(0)->values()(0));
552 BOOST_TEST(testing::equals(data.at(0)->values()(1), -2.35405379763935940218e-01), data.at(0)->values()(1));
553 BOOST_TEST(testing::equals(data.at(1)->values()(0), 7.23368584254213131679e-02), data.at(1)->values()(0));
554 BOOST_TEST(testing::equals(data.at(1)->values()(1), 7.23368584254213131679e-02), data.at(1)->values()(1));
555 }
556}
557
559BOOST_AUTO_TEST_CASE(testIMVJ_effUpdate_ppWithoutSubsteps)
560{
561 PRECICE_TEST(""_on(4_ranks).setupIntraComm());
562 // config:
563 double initialRelaxation = 0.1;
564 int maxIterationsUsed = 30;
565 int timeWindowsReused = 0;
566 int filter = BaseQNAcceleration::QR2FILTER;
567 int restartType = IQNIMVJAcceleration::NO_RESTART;
568 int chunkSize = 0;
569 int reusedTimeWindowsAtRestart = 0;
570 double singularityLimit = 1e-2;
571 double svdTruncationEps = 0.0;
572 bool enforceInitialRelaxation = false;
573 bool alwaysBuildJacobian = false;
574
575 std::vector<int> dataIDs;
576 dataIDs.push_back(4);
577 dataIDs.push_back(5);
578 PtrPreconditioner _preconditioner = PtrPreconditioner(new ResidualSumPreconditioner(-1));
579 std::vector<int> vertexOffsets{0, 11, 22, 22};
580
581 mesh::PtrMesh dummyMesh(new mesh::Mesh("dummyMesh", 2, testing::nextMeshID()));
582 dummyMesh->setVertexOffsets(vertexOffsets);
583
584 IQNIMVJAcceleration pp(initialRelaxation, enforceInitialRelaxation, maxIterationsUsed,
585 timeWindowsReused, filter, singularityLimit, dataIDs, _preconditioner, alwaysBuildJacobian,
586 restartType, chunkSize, reusedTimeWindowsAtRestart, svdTruncationEps);
587
588 mesh::PtrData displacements(new mesh::Data("dvalues", -1, 2));
589 mesh::PtrData forces(new mesh::Data("fvalues", -1, 2));
590
591 Eigen::VectorXd dref;
592 Eigen::VectorXd fref;
593 double drefNorm = 0., frefNorm = 0.;
594
595 DataMap data;
596 PtrCouplingData dpcd;
597 PtrCouplingData fpcd;
598
599 bool exchangeSubsteps = false; // @todo add testIMVJ_effUpdate_ppWithSubsteps, where exchangeSubsteps = true as soon as acceleration scheme supports subcycling.
600
601 if (context.isPrimary()) { // Primary
606 displacements->values().resize(0);
607 forces->values().resize(0);
608
609 // init displacements
610 dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
611 dpcd->setSampleAtTime(0, dpcd->sample());
612
613 // init forces
614 fpcd = makeCouplingData(forces, dummyMesh, exchangeSubsteps);
615 fpcd->setSampleAtTime(0, fpcd->sample());
616
617 dpcd->storeIteration();
618 fpcd->storeIteration();
619
620 data.insert(std::pair<int, PtrCouplingData>(4, dpcd));
621 data.insert(std::pair<int, PtrCouplingData>(5, fpcd));
622
623 pp.initialize(data);
624 } else if (context.isRank(1)) { // SecondaryRank1
629 // init displacements
630 displacements->values().resize(22);
631 forces->values().resize(22);
632 displacements->values() << 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0;
633 forces->values() << 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0;
634
635 dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
636 dpcd->setSampleAtTime(0, dpcd->sample());
637 fpcd = makeCouplingData(forces, dummyMesh, exchangeSubsteps);
638 fpcd->setSampleAtTime(0, fpcd->sample());
639
640 dpcd->storeIteration();
641 fpcd->storeIteration();
642
643 data.insert(std::pair<int, PtrCouplingData>(4, dpcd));
644 data.insert(std::pair<int, PtrCouplingData>(5, fpcd));
645
646 pp.initialize(data);
647
648 forces->values() << -0.01765744149520443, -0.000534499502588083, 0.05397520666020422, 0.0005546984205735067, 0.05213823386543703, 0.0007618478879228568, -0.01944857239806249, -0.0009206665792022876, -0.02459872346309381, -0.001296931976456198, 0.04688718434761113, 0.001346643628716769, -0.01063536095060684, -0.01905148710330257, 0.02514593936525903, -0.01643393169986981, -0.02189723835016068, -0.000912218689367709, 0.04985117008772211, 0.0009615805506705544, 0.05534647415570375, 0.0004068469082890895;
649 } else if (context.isRank(2)) { // Secondary rank 2
654 // init displacements
655 displacements->values().resize(22);
656 forces->values().resize(22);
657 displacements->values() << 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0;
658 forces->values() << 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0;
659
660 dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
661 dpcd->setSampleAtTime(0, dpcd->sample());
662 fpcd = makeCouplingData(forces, dummyMesh, exchangeSubsteps);
663 fpcd->setSampleAtTime(0, fpcd->sample());
664
665 dpcd->storeIteration();
666 fpcd->storeIteration();
667
668 data.insert(std::pair<int, PtrCouplingData>(4, dpcd));
669 data.insert(std::pair<int, PtrCouplingData>(5, fpcd));
670
671 pp.initialize(data);
672
673 forces->values() << -0.01465589151503364, -0.0002670111835650672, 0.05711438689366102, 0.0002383730129847531, -0.01536098575916998, -0.000285287812066552, 0.05638274807579218, 0.000283961973555227, -0.006856432131857973, -0.006815594391460808, 0.02901925611525407, -0.02907380915674757, 0.05800715138289463, 9.667376010126116e-05, -0.01376443700165205, -9.547563271960956e-05, 0.05768190311116184, 0.0001311583226994801, -0.01408147387131287, -0.0001216961377915992, -0.0163823504288376, -0.0003874626690545313;
674 } else if (context.isRank(3)) { // Secondary rank 3
679 displacements->values().resize(0);
680 forces->values().resize(0);
681
682 dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
683 dpcd->setSampleAtTime(0, dpcd->sample());
684 fpcd = makeCouplingData(forces, dummyMesh, exchangeSubsteps);
685 fpcd->setSampleAtTime(0, fpcd->sample());
686
687 dpcd->storeIteration();
688 fpcd->storeIteration();
689
690 data.insert(std::pair<int, PtrCouplingData>(4, dpcd));
691 data.insert(std::pair<int, PtrCouplingData>(5, fpcd));
692
693 pp.initialize(data);
694 }
695
696 pp.performAcceleration(data);
697
698 // necessary because acceleration does not directly work on storage, but on CouplingData::values. See and https://github.com/precice/precice/issues/1645 current implementation in BaseCouplingScheme::doImplicitStep()
699 for (auto &pair : data) {
700 pair.second->setSampleAtTime(1, pair.second->sample());
701 }
702
703 // underrelaxation, first iteration
704
705 if (context.isPrimary()) { // Primary
706
707 } else if (context.isRank(1)) { // SecondaryRank1
708
709 dref = Eigen::VectorXd::Zero(22);
710 fref = Eigen::VectorXd::Zero(22);
711
712 dref << 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0;
713 fref << -0.001765744149520443, -5.344995025880831e-05, 0.005397520666020422, 5.546984205735068e-05, 0.005213823386543704, 7.618478879228569e-05, -0.001944857239806249, -9.206665792022876e-05, -0.002459872346309381, -0.0001296931976456198, 0.004688718434761114, 0.000134664362871677, -0.001063536095060684, -0.001905148710330257, 0.002514593936525903, -0.001643393169986981, -0.002189723835016068, -9.12218689367709e-05, 0.004985117008772211, 9.615805506705544e-05, 0.005534647415570375, 4.068469082890896e-05;
714 frefNorm = 0.01286041129960619;
715 drefNorm = 0.0;
716
717 // validate values
718 for (int i = 0; i < data.at(4)->values().size(); i++)
719 BOOST_TEST(testing::equals(data.at(4)->values()(i), dref(i)));
720
721 for (int i = 0; i < data.at(5)->values().size(); i++)
722 BOOST_TEST(testing::equals(data.at(5)->values()(i), fref(i)));
723
724 // validate norm
725 BOOST_TEST(testing::equals(data.at(4)->values().norm(), drefNorm), data.at(4)->values().norm());
726 BOOST_TEST(testing::equals(data.at(5)->values().norm(), frefNorm), data.at(5)->values().norm());
727
728 // update cplData
729 fpcd->storeIteration();
730 displacements->values() << 1.790053057185293e-06, -2.44566429072041e-08, 1.889281703254964e-06, -1.972492834475447e-07, 1.681634609242917e-06, -2.373356532433882e-07, 1.585003447958184e-06, -5.301850772916681e-08, 1.274187257620066e-06, -2.137488936999111e-07, 1.362955262700412e-06, -2.762153471191986e-07, 1.249747540920782e-06, -3.196338173465977e-07, 1.333501893726392e-06, -3.161541101487353e-07, 1.394538527892028e-06, -1.166536323805688e-07, 1.488382850875808e-06, -2.605379508545059e-07, 2.056077021837937e-06, -1.341692715765341e-07;
731 forces->values() << -0.01765744187144705, -0.000534499502451157, 0.05397520721069472, 0.0005546984181272257, 0.05213823442800309, 0.000761847881060882, -0.01944857277019029, -0.0009206665773591249, -0.02459872381892192, -0.001296931982922439, 0.04688718490326162, 0.001346643636856003, -0.01063536111298416, -0.01905148734069537, 0.02514593966043068, -0.01643393192020026, -0.02189723869781963, -0.0009122186870252733, 0.04985117065739809, 0.0009615805515004192, 0.05534647470527156, 0.0004068469091761907;
732 } else if (context.isRank(2)) { // Secondary rank 2
733
734 dref = Eigen::VectorXd::Zero(22);
735 fref = Eigen::VectorXd::Zero(22);
736
737 dref << 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0;
738 fref << -0.001465589151503364, -2.670111835650672e-05, 0.005711438689366103, 2.383730129847531e-05, -0.001536098575916998, -2.85287812066552e-05, 0.005638274807579218, 2.83961973555227e-05, -0.0006856432131857974, -0.0006815594391460808, 0.002901925611525407, -0.002907380915674757, 0.005800715138289463, 9.667376010126117e-06, -0.001376443700165206, -9.547563271960956e-06, 0.005768190311116184, 1.311583226994801e-05, -0.001408147387131287, -1.216961377915992e-05, -0.00163823504288376, -3.874626690545313e-05;
739 frefNorm = 0.01265754337961098;
740 drefNorm = 0.0;
741
742 // validate values
743 for (int i = 0; i < data.at(4)->values().size(); i++)
744 BOOST_TEST(testing::equals(data.at(4)->values()(i), dref(i)));
745
746 for (int i = 0; i < data.at(5)->values().size(); i++)
747 BOOST_TEST(testing::equals(data.at(5)->values()(i), fref(i)));
748
749 // validate norm
750 BOOST_TEST(testing::equals(data.at(4)->values().norm(), drefNorm), data.at(4)->values().norm());
751 BOOST_TEST(testing::equals(data.at(5)->values().norm(), frefNorm), data.at(5)->values().norm());
752
753 // update cplData
754 fpcd->storeIteration();
755 displacements->values() << 1.848184969639987e-06, -1.983566187932991e-07, 1.952383060128974e-06, 1.050101286643166e-07, 2.020975712018586e-06, -9.297459906882382e-08, 2.123910878481957e-06, -3.349554682884977e-08, 0, 0, 0, 0, 7.715047421278781e-07, 2.958323850532032e-07, 6.5137785527863e-07, -3.40165313149562e-07, 1.498023570500414e-06, 2.492038233690158e-07, 1.395223018993416e-06, -3.150663149441921e-07, 1.954718171910318e-06, -3.415637300374603e-08;
756 forces->values() << -0.0146558918972568, -0.000267011181975166, 0.05711438744699839, 0.0002383730136872111, -0.0153609861368436, -0.0002852878106683293, 0.05638274862725741, 0.0002839619744993407, -0.00685643232676097, -0.006815594586569211, 0.02901925639144463, -0.02907380943293575, 0.05800715193585099, 9.667375963025685e-05, -0.01376443739049903, -9.547563172575954e-05, 0.05768190366530584, 0.0001311583223016465, -0.01408147425699792, -0.0001216961368213471, -0.01638235080508845, -0.0003874626694560972;
757 } else if (context.isRank(3)) { // Secondary rank 3
758 // Dummy Secondary rank to be able to reuse the 4 proc Intra-participant communication Fixture
759 }
760
761 // QN- Update, 2. iteration
762 pp.performAcceleration(data);
763
764 // necessary because acceleration does not directly work on storage, but on CouplingData::values. See and https://github.com/precice/precice/issues/1645 current implementation in BaseCouplingScheme::doImplicitStep()
765 for (auto &pair : data) {
766 pair.second->setSampleAtTime(1, pair.second->sample());
767 }
768
769 if (context.isPrimary()) { // Primary
770
771 } else if (context.isRank(1)) { // SecondaryRank1
772
773 dref = Eigen::VectorXd::Zero(22);
774 fref = Eigen::VectorXd::Zero(22);
775
776 dref << 2.182991240484406e-07, -2.982517027848927e-09, 2.30400176824822e-07, -2.405478743936717e-08, 2.050773638768581e-07, -2.89433684663868e-08, 1.932930774951963e-07, -6.465670807451687e-09, 1.553886691210811e-07, -2.606693476135542e-08, 1.662140341429634e-07, -3.368479391313121e-08, 1.524082162646531e-07, -3.897972859683727e-08, 1.62622160359393e-07, -3.85553741174048e-08, 1.700656513328811e-07, -1.422605082208521e-08, 1.815100794307223e-07, -3.17729165761968e-08, 2.507410666078866e-07, -1.636210409619326e-08;
777 fref << -0.01765744154108767, -0.0005344995025713848, 0.0539752067273372, 0.0005546984202751798, 0.05213823393404264, 0.0007618478870860308, -0.01944857244344392, -0.0009206665789775117, -0.02459872350648748, -0.001296931977244764, 0.04688718441537337, 0.001346643629709359, -0.01063536097040895, -0.01905148713225291, 0.02514593940125557, -0.01643393172673937, -0.0218972383925581, -0.0009122186890820461, 0.04985117015719479, 0.0009615805507717574, 0.05534647422272421, 0.0004068469083972726;
778 drefNorm = 6.435143632392166e-07;
779 frefNorm = 0.1286041131776192;
780
781 // validate values
782 for (int i = 0; i < data.at(4)->values().size(); i++)
783 BOOST_TEST(testing::equals(data.at(4)->values()(i), dref(i)));
784
785 for (int i = 0; i < data.at(5)->values().size(); i++)
786 BOOST_TEST(testing::equals(data.at(5)->values()(i), fref(i)));
787
788 // validate norm
789 BOOST_TEST(testing::equals(data.at(4)->values().norm(), drefNorm), data.at(4)->values().norm());
790 BOOST_TEST(testing::equals(data.at(5)->values().norm(), frefNorm), data.at(5)->values().norm());
791
792 // update cplData
793 dpcd->storeIteration();
794 fpcd->storeIteration();
795 displacements->values() << 1.790034504773721e-05, -2.446591076368466e-07, 1.889267115021718e-05, -1.972643201602028e-06, 1.681613350812527e-05, -2.373460013995369e-06, 1.584978895355817e-05, -5.302446869164338e-07, 1.274157692078479e-05, -2.137546278211264e-06, 1.362926508984742e-05, -2.762211725309514e-06, 1.249719424608544e-05, -3.19640295598053e-06, 1.333474052315949e-05, -3.16159193819195e-06, 1.394510078525391e-05, -1.166587691625877e-06, 1.488356439901566e-05, -2.605456452904905e-06, 2.056070000286195e-05, -1.341920935569228e-06;
796 forces->values() << -0.01845221261910751, -0.000473660279052688, 0.05115217408647257, 0.0004997712466226923, 0.04953299847638116, 0.0006834967349236343, -0.02003343430934222, -0.000812620519068711, -0.02461452082338934, -0.00115320030035888, 0.04486850176987132, 0.00121794049154945, -0.01080003301957858, -0.01839753254507924, 0.02398381340216066, -0.01606849579606463, -0.02220297183992202, -0.0008082032560853196, 0.04750952330271016, 0.0008672013640382038, 0.05236940113731539, 0.0003710183715860552;
797 } else if (context.isRank(2)) { // Secondary rank 2
798
799 dref = Eigen::VectorXd::Zero(22);
800 fref = Eigen::VectorXd::Zero(22);
801
802 dref << 2.253883807144274e-07, -2.418982831708627e-08, 2.380954632167933e-07, 1.280611153486133e-08, 2.464604196428374e-07, -1.133836421999323e-08, 2.590134870407767e-07, -4.084822236363769e-09, 0, 0, 0, 0, 9.408593154806163e-08, 3.607711529166728e-08, 7.943631316463149e-08, -4.148356921273439e-08, 1.826857767882953e-07, 3.039070609254422e-08, 1.701491258462493e-07, -3.842271618341636e-08, 2.38380232908047e-07, -4.165410783473636e-09;
803 fref << -0.01465589156164622, -0.0002670111833711768, 0.05711438696114118, 0.0002383730130704187, -0.01536098580522773, -0.0002852878118960371, 0.05638274814304403, 0.0002839619736703628, -0.006856432155626628, -0.006815594415254513, 0.02901925614893584, -0.02907380919042905, 0.05800715145032832, 9.667376004382162e-05, -0.01376443704907241, -9.547563259840837e-05, 0.05768190317874038, 0.0001311583226509638, -0.01408147391834763, -0.0001216961376732758, -0.01638235047472185, -0.0003874626691035027;
804 drefNorm = 6.131610103923933e-07;
805 frefNorm = 0.1265754339635572;
806
807 // validate values
808 for (int i = 0; i < data.at(4)->values().size(); i++)
809 BOOST_TEST(testing::equals(data.at(4)->values()(i), dref(i)));
810
811 for (int i = 0; i < data.at(5)->values().size(); i++)
812 BOOST_TEST(testing::equals(data.at(5)->values()(i), fref(i)));
813
814 // validate norm
815 BOOST_TEST(testing::equals(data.at(4)->values().norm(), drefNorm), data.at(4)->values().norm());
816 BOOST_TEST(testing::equals(data.at(5)->values().norm(), frefNorm), data.at(5)->values().norm());
817
818 // update cplData
819 dpcd->storeIteration();
820 fpcd->storeIteration();
821 displacements->values() << 1.848182952307335e-05, -1.983938722952872e-06, 1.952389995095743e-05, 1.049689886611777e-06, 2.020972044646931e-05, -9.30012125294331e-07, 2.123911759834233e-05, -3.352823479948144e-07, 0, 0, 0, 0, 7.715124780435689e-06, 2.958056858428718e-06, 6.513639301665504e-06, -3.401886529062288e-06, 1.498034283416962e-05, 2.491634858078641e-06, 1.39521486945152e-05, -3.151050708450101e-06, 1.954707223943552e-05, -3.417246252999375e-07;
822 forces->values() << -0.01568208277628194, -0.0002595395446636614, 0.0540328986967421, 0.0002362571305830931, -0.01637736854863682, -0.0002699645831085989, 0.05331751790879287, 0.0002707054191427001, -0.007277539612331946, -0.007235194100552225, 0.02757151633202504, -0.02762772092892902, 0.05505877464319012, 0.0001052840945529276, -0.01465499974491537, -0.0001017767294585529, 0.05464614037258596, 0.0001424559420056945, -0.01506072500921042, -0.0001315030046882618, -0.0173164149989076, -0.0003474184175392483;
823 } else if (context.isRank(3)) { // Secondary rank 3
824 // Dummy Secondary rank to be able to reuse the 4 proc Intra-participant communication Fixture
825 }
826
827 // QN- Update, 3. iteration
828 pp.performAcceleration(data);
829
830 // necessary because acceleration does not directly work on storage, but on CouplingData::values. See and https://github.com/precice/precice/issues/1645 current implementation in BaseCouplingScheme::doImplicitStep()
831 for (auto &pair : data) {
832 pair.second->setSampleAtTime(1, pair.second->sample());
833 }
834
835 if (context.isPrimary()) { // Primary
836
837 } else if (context.isRank(1)) { // SecondaryRank1
838
839 dref = Eigen::VectorXd::Zero(22);
840 fref = Eigen::VectorXd::Zero(22);
841
842 dref << 1.717512702194643e-05, -2.347247264125579e-07, 1.812723818498285e-05, -1.892683301884884e-06, 1.613485037277307e-05, -2.277271314691775e-06, 1.520766651015439e-05, -5.087470614753933e-07, 1.222540039026145e-05, -2.050926769187082e-06, 1.30771205492404e-05, -2.650282676556395e-06, 1.199091589196657e-05, -3.066880426987482e-06, 1.279452699120483e-05, -3.033483087437419e-06, 1.338015894263769e-05, -1.11930952734218e-06, 1.428059443561221e-05, -2.499874235762116e-06, 1.972766650768916e-05, -1.287497600071812e-06;
843 fref << -0.01823457911070198, -0.0004903200524407003, 0.05192521422826112, 0.0005148121016499705, 0.05024639863781041, 0.0007049518329474636, -0.01987328081455485, -0.0008422070677716292, -0.02461019581850909, -0.001192558759870676, 0.04542128452828625, 0.001253183648371612, -0.01075494078801318, -0.01857660729355974, 0.02430204227390975, -0.01616856463954442, -0.02211925279965755, -0.0008366860858679643, 0.04815074428053346, 0.0008930454810143837, 0.05318462260012811, 0.0003808294014608455;
844 drefNorm = 5.062970166651817e-05;
845 frefNorm = 0.1247902554601672;
846
847 // validate values
848 for (int i = 0; i < data.at(4)->values().size(); i++)
849 BOOST_TEST(testing::equals(data.at(4)->values()(i), dref(i)));
850
851 for (int i = 0; i < data.at(5)->values().size(); i++)
852 BOOST_TEST(testing::equals(data.at(5)->values()(i), fref(i)));
853
854 // validate norm
855 BOOST_TEST(testing::equals(data.at(4)->values().norm(), drefNorm), data.at(4)->values().norm());
856 BOOST_TEST(testing::equals(data.at(5)->values().norm(), frefNorm), data.at(5)->values().norm());
857
858 // update cplData
859 dpcd->storeIteration();
860 fpcd->storeIteration();
861 displacements->values() << 1.659080663925766e-05, -2.839283676791931e-07, 1.756292801739508e-05, -1.881726812992964e-06, 1.564798101471437e-05, -2.265931706775091e-06, 1.470124517392331e-05, -5.705378156142171e-07, 1.186047603634431e-05, -2.115667271562722e-06, 1.273027556448604e-05, -2.674541973319838e-06, 1.165645170777486e-05, -3.135385366949176e-06, 1.247728214631633e-05, -3.082564251671268e-06, 1.295443089215965e-05, -1.185450561958201e-06, 1.387356342346108e-05, -2.500933334689963e-06, 1.911143938064833e-05, -1.289577439500651e-06;
862 forces->values() << -0.08018882094302014, 0.004252226533647444, -0.1681454949889014, -0.00376664315520894, -0.1528473285523723, -0.005402482563102148, -0.06546415562528343, 0.007580544006107058, -0.02584070375141325, 0.01001310706105286, -0.1119519940946071, -0.008779400784032571, -0.02359260557541131, 0.03241153400001811, -0.06629773012631451, 0.01232457880930212, -0.04595115217779744, 0.007271575191182346, -0.1343966169086029, -0.006463593204557456, -0.1788931815052193, -0.00241194798896828;
863 } else if (context.isRank(2)) { // Secondary rank 2
864
865 dref = Eigen::VectorXd::Zero(22);
866 fref = Eigen::VectorXd::Zero(22);
867
868 dref << 1.773301313815872e-05, -1.903469331632798e-06, 1.873284151703964e-05, 1.007255996407046e-06, 1.939089962749955e-05, -8.922690925673692e-07, 2.037857848281338e-05, -3.216215773080572e-07, 0, 0, 0, 0, 7.402516024117914e-06, 2.838268713863215e-06, 6.249761203424372e-06, -3.263999161516882e-06, 1.437336512054573e-05, 2.390776372597577e-06, 1.338687393172981e-05, -3.023290389839287e-06, 1.875511666670207e-05, -3.278415622262174e-07;
869 fref << -0.01540107886229656, -0.0002615855206049247, 0.05487671246695029, 0.0002368365301820234, -0.0160990505051311, -0.0002741605822288765, 0.05415687969310452, 0.0002743355004921105, -0.007162227035058467, -0.007120294400987128, 0.02796795558583163, -0.02802370793272083, 0.05586613813214418, 0.0001029263016352541, -0.0144111353804976, -0.0001000512803033142, 0.05547743301530601, 0.0001393622825766435, -0.01479257484776803, -0.0001288175608001121, -0.01706063837893593, -0.0003583838471874983;
870 drefNorm = 4.824205753272403e-05;
871 frefNorm = 0.1225851627302816;
872
873 // validate values
874 for (int i = 0; i < data.at(4)->values().size(); i++)
875 BOOST_TEST(testing::equals(data.at(4)->values()(i), dref(i)));
876
877 for (int i = 0; i < data.at(5)->values().size(); i++)
878 BOOST_TEST(testing::equals(data.at(5)->values()(i), fref(i)));
879
880 // validate norm
881 BOOST_TEST(testing::equals(data.at(4)->values().norm(), drefNorm), data.at(4)->values().norm());
882 BOOST_TEST(testing::equals(data.at(5)->values().norm(), frefNorm), data.at(5)->values().norm());
883
884 // update cplData
885 dpcd->storeIteration();
886 fpcd->storeIteration();
887 displacements->values() << 1.716650969972045e-05, -1.856138836171773e-06, 1.818701485070425e-05, 9.439657883607802e-07, 1.874709534954619e-05, -8.85448704675396e-07, 1.975527973304359e-05, -3.501096287428596e-07, 0, 0, 0, 0, 7.228951427433641e-06, 2.745909101918556e-06, 6.052367643912141e-06, -3.179587143921995e-06, 1.398276918926419e-05, 2.29762824040882e-06, 1.297587398676e-05, -2.941551341709183e-06, 1.811863361465251e-05, -3.546317448342288e-07;
888 forces->values() << -0.09539527385890252, 0.0003208855941258066, -0.1853399184726223, 7.203155656644242e-05, -0.09532865072058605, 0.0009202649288056726, -0.1847925968312873, -0.0007589246108979722, -0.03998875591551594, -0.03982927597079221, -0.08489044889406808, 0.08470593806523596, -0.1739740974580442, 0.0007742373134568178, -0.08383286811708256, -0.0005911288917162662, -0.1811747642897668, 0.001020161732709184, -0.09112767929864005, -0.0008931566039992005, -0.08987332323372975, 0.002763113283891189;
889 } else if (context.isRank(3)) { // Secondary rank 3
890 // Dummy Secondary rank to be able to reuse the 4 proc Intra-participant communication Fixture
891 }
892
893 // QN- Update, 4. iteration
894 pp.performAcceleration(data);
895
896 // necessary because acceleration does not directly work on storage, but on CouplingData::values. See and https://github.com/precice/precice/issues/1645 current implementation in BaseCouplingScheme::doImplicitStep()
897 for (auto &pair : data) {
898 pair.second->setSampleAtTime(1, pair.second->sample());
899 }
900
901 if (context.isPrimary()) { // Primary
902
903 } else if (context.isRank(1)) { // SecondaryRank1
904
905 dref = Eigen::VectorXd::Zero(22);
906 fref = Eigen::VectorXd::Zero(22);
907
908 dref << 1.633368119919659e-05, -2.237360618564265e-07, 1.723960679736894e-05, -1.800451405634617e-06, 1.534489463502018e-05, -2.166297690952519e-06, 1.446268777631733e-05, -4.845152271528126e-07, 1.162685558347312e-05, -1.951619069146136e-06, 1.243725578988821e-05, -2.521440338918289e-06, 1.140405450734714e-05, -2.91813193235939e-06, 1.216867994540041e-05, -2.886191281742917e-06, 1.272486031358505e-05, -1.065380276784719e-06, 1.35816114080243e-05, -2.3781612091675e-06, 1.876166803225929e-05, -1.224865187618123e-06;
909 fref << -0.01890549129941078, -0.0004389621840029381, 0.04954204963698766, 0.0004684469509178089, 0.04804708246519945, 0.0006388129115723739, -0.02036699235959898, -0.0007509963669595499, -0.02462352250264116, -0.00107121393884661, 0.04371708893453349, 0.001144538793507394, -0.01089395922122854, -0.01802447191842139, 0.02332094230979109, -0.01586002180275038, -0.02237733336181922, -0.0007488803918274927, 0.04617393036901896, 0.00081337820788483, 0.05067142953803216, 0.0003505856246817933;
910 drefNorm = 4.815113092042934e-05;
911 frefNorm = 0.1204042970796453;
912
913 // validate values
914 for (int i = 0; i < data.at(4)->values().size(); i++)
915 BOOST_TEST(testing::equals(data.at(4)->values()(i), dref(i)));
916
917 for (int i = 0; i < data.at(5)->values().size(); i++)
918 BOOST_TEST(testing::equals(data.at(5)->values()(i), fref(i)));
919
920 // validate norm
921 BOOST_TEST(testing::equals(data.at(4)->values().norm(), drefNorm));
922 BOOST_TEST(testing::equals(data.at(5)->values().norm(), frefNorm));
923
924 // update cplData
925 dpcd->storeIteration();
926 fpcd->storeIteration();
927 displacements->values() << 1.506845042291629e-05, -3.295713481574521e-07, 1.601708402767785e-05, -1.776032790440438e-06, 1.428998106709373e-05, -2.140925300298825e-06, 1.336604025915203e-05, -6.173668108734595e-07, 1.083614857997936e-05, -2.09020895349816e-06, 1.168515223592441e-05, -2.572621503366579e-06, 1.067901778881212e-05, -3.064421860106172e-06, 1.14804152344671e-05, -2.990689693425248e-06, 1.180274523671064e-05, -1.207361572630013e-06, 1.26994053163659e-05, -2.379420351559266e-06, 1.742665917249236e-05, -1.228726437307901e-06;
928 forces->values() << -0.07712271523301416, 0.004018032584755363, -0.1572746631951394, -0.003554957358717858, -0.1428154301136939, -0.005100471254411636, -0.06320679903752099, 0.00716455247252809, -0.02577587645775314, 0.009459639059138256, -0.1041793821597049, -0.008283348947413409, -0.02295652849878388, 0.02989440569403293, -0.06182269085435656, 0.01091909076145766, -0.04476922496853244, 0.00687112729333933, -0.1253800696382085, -0.006099818450316078, -0.1674291774005812, -0.002273881806091619;
929 } else if (context.isRank(2)) { // Secondary rank 2
930
931 dref = Eigen::VectorXd::Zero(22);
932 fref = Eigen::VectorXd::Zero(22);
933
934 dref << 1.686458723791834e-05, -1.810446627874213e-06, 1.781592254220799e-05, 9.575750096228968e-07, 1.844107125032693e-05, -8.488149913661512e-07, 1.938083795505493e-05, -3.062726369792662e-07, 0, 0, 0, 0, 7.040556208586022e-06, 2.699202297106055e-06, 5.943687851414432e-06, -3.104374082089934e-06, 1.367008353382051e-05, 2.273466384432095e-06, 1.27314174291305e-05, -2.875442657958392e-06, 1.783629704415144e-05, -3.12140240279707e-07;
935 fref << -0.01626734826572977, -0.0002552779342950638, 0.05227538136126555, 0.0002350515292044683, -0.01695703999572442, -0.0002612258237118909, 0.05156927126324851, 0.0002631458821373055, -0.007517710152292339, -0.007474504607240594, 0.02674580053052369, -0.02680294718917284, 0.05337717466458169, 0.0001101957844225173, -0.01516291356738642, -0.0001053694865272746, 0.05291470170150131, 0.0001489003323108161, -0.01561921932892748, -0.0001370950096012615, -0.01784913805041409, -0.000324580522802835;
936 drefNorm = 4.587992970636867e-05;
937 frefNorm = 0.1180354804938519;
938
939 // validate values
940 for (int i = 0; i < data.at(4)->values().size(); i++)
941 BOOST_TEST(testing::equals(data.at(4)->values()(i), dref(i)));
942
943 for (int i = 0; i < data.at(5)->values().size(); i++)
944 BOOST_TEST(testing::equals(data.at(5)->values()(i), fref(i)));
945
946 // validate norm
947 BOOST_TEST(testing::equals(data.at(4)->values().norm(), drefNorm));
948 BOOST_TEST(testing::equals(data.at(5)->values().norm(), frefNorm));
949
950 // update cplData
951 dpcd->storeIteration();
952 fpcd->storeIteration();
953 displacements->values() << 1.563743909676446e-05, -1.707572586404205e-06, 1.663287551913161e-05, 8.210579991784308e-07, 1.704678071734513e-05, -8.336427145015805e-07, 1.803030552728031e-05, -3.673472962716038e-07, 0, 0, 0, 0, 6.663771864888832e-06, 2.499283366425937e-06, 5.516134032932667e-06, -2.921164340377279e-06, 1.282308279788757e-05, 2.07209067754735e-06, 1.184094543159743e-05, -2.698009821996337e-06, 1.645805878055576e-05, -3.696322259193852e-07;
954 forces->values() << -0.09143825207871489, 0.0002922798859936043, -0.1734744585823354, 8.018501629471556e-05, -0.09140909287296797, 0.0008614262538692869, -0.1729893406976169, -0.0007078492982043917, -0.03836482525057584, -0.03821118279703851, -0.07931609050916776, 0.07913795276507131, -0.1626218046186428, 0.0007411076261039719, -0.08039872451576649, -0.0005667402343291361, -0.1694857588798654, 0.0009766586358261331, -0.0873517838382746, -0.0008552683303008771, -0.08627064033821233, 0.002609015553424872;
955 } else if (context.isRank(3)) { // Secondary rank 3
956 // Dummy Secondary rank to be able to reuse the 4 proc Intra-participant communication Fixture
957 }
958
959 // QN- Update, 5. iteration
960 pp.performAcceleration(data);
961
962 // necessary because acceleration does not directly work on storage, but on CouplingData::values. See and https://github.com/precice/precice/issues/1645 current implementation in BaseCouplingScheme::doImplicitStep()
963 for (auto &pair : data) {
964 pair.second->setSampleAtTime(1, pair.second->sample());
965 }
966
967 if (context.isPrimary()) { // Primary
968
969 } else if (context.isRank(1)) { // SecondaryRank1
970
971 dref = Eigen::VectorXd::Zero(22);
972 fref = Eigen::VectorXd::Zero(22);
973
974 dref << 1.275776729912441e-06, -7.411719120649928e-07, 2.009844043916140e-06, -8.166362562036991e-07,
975 1.98429549697312266297e-06, -1.006128466370864e-06, 1.26860076890459333335e-06, -1.038973060940335e-06,
976 1.553926578960109e-06, -1.85501566670221e-06, 2.21293251735091021427e-06, -1.645357854619679e-06,
977 1.820890655907848e-06, -2.41565468511731e-06, 2.44472211085250353034e-06, -2.153167028992053e-06,
978 1.367377699262639e-06, -1.402360838056698e-06, 2.05877363493188315457e-06, -1.275601072931189e-06,
979 2.16056223288899190422e-06, -6.758668464219906e-07;
980
981 fref << -0.02494062387644205, 2.83457783442084623737e-05, 2.78926996928726099456e-02, 4.927873553950413e-05,
982 0.02806450388830189, 4.133693821366594e-05, -0.02479702682061509, 7.828236584438153e-05,
983 -0.02470268458079963, 3.17420733502629e-05, 0.0282235973672153, 0.0001624809599050109,
984 -0.01213681562464092, -0.01299265397817983, 0.01440494173854578, -0.01303642911469127,
985 -0.02467306876921166, 4.955128897708022e-05, 0.02820821342518329, 9.359794339403125e-05,
986 0.02784148826297314, 7.700134509804057e-05;
987
988 drefNorm = 7.910283453653413e-06;
989 frefNorm = 0.08415800797163485;
990
991 // validate values
992 for (int i = 0; i < data.at(4)->values().size(); i++)
993 BOOST_TEST(testing::equals(data.at(4)->values()(i), dref(i)));
994
995 for (int i = 0; i < data.at(5)->values().size(); i++)
996 BOOST_TEST(testing::equals(data.at(5)->values()(i), fref(i), 1e-8));
997
998 // validate norm
999 BOOST_TEST(testing::equals(data.at(4)->values().norm(), drefNorm));
1000 BOOST_TEST(testing::equals(data.at(5)->values().norm(), frefNorm));
1001
1002 } else if (context.isRank(2)) { // Secondary rank 2
1003
1004 dref = Eigen::VectorXd::Zero(22);
1005 fref = Eigen::VectorXd::Zero(22);
1006
1007 dref << 1.782695896956317e-06, -3.609665456581163e-07, 2.549165865006206e-06, -2.9152246796606e-07,
1008 1.641256623136569e-06, -3.63485032734166e-07, 2.399541139942214e-06, -5.221960935703056e-07,
1009 0, 0, 0, 0,
1010 1.540562829986484e-06, 2.646900958343007e-07, 6.572995607065404e-07, -5.789378440843989e-07,
1011 2.312916074266845e-06, 2.909615579313369e-08, 1.55690286980020939237e-06, -4.907962897313283e-07,
1012 1.41306136610511245499e-06, -5.042593080795539e-07;
1013
1014 fref << -0.02407925914468757, -0.0001962517462601011, 2.86388026412357221684e-02, 0.0002189694547629064,
1015 -0.02469057752137756, -0.0001420431374215673, 0.02806138636325711, 0.0001618659175134536,
1016 -0.01072215615486178, -0.01066770812537758, 0.01563815124261267, -0.01570783185165062,
1017 0.03075510899502418, 0.0001765574924817709, -0.02194129299700746, -0.0001523229943385313,
1018 0.02962481369149697, 2.35608194343154197722e-04, -0.02307508702308544, -0.0002109167490697018,
1019 -0.02495036462961023, -1.65224214831653608282e-05;
1020
1021 drefNorm = 5.676684399367158e-06;
1022 frefNorm = 0.08353026170200345;
1023
1024 // validate values
1025 for (int i = 0; i < data.at(4)->values().size(); i++)
1026 BOOST_TEST(testing::equals(data.at(4)->values()(i), dref(i)));
1027
1028 for (int i = 0; i < data.at(5)->values().size(); i++)
1029 BOOST_TEST(testing::equals(data.at(5)->values()(i), fref(i)));
1030
1031 // validate norm
1032 BOOST_TEST(testing::equals(data.at(4)->values().norm(), drefNorm));
1033 BOOST_TEST(testing::equals(data.at(5)->values().norm(), frefNorm));
1034 } else if (context.isRank(3)) { // Secondary rank 3
1035 // Dummy Secondary rank to be able to reuse the 4 proc Intra-participant communication Fixture
1036 }
1037}
1038
1040BOOST_AUTO_TEST_CASE(testColumnsLoggingWithoutSubsteps)
1041{
1042 PRECICE_TEST(""_on(4_ranks).setupIntraComm());
1043 double initialRelaxation = 0.1;
1044 int maxIterationsUsed = 3;
1045 int timeWindowsReused = 1;
1046 int filter = BaseQNAcceleration::QR1FILTER;
1047 double singularityLimit = 0.1;
1048 bool enforceInitialRelaxation = false;
1049 std::vector<int> dataIDs;
1050 dataIDs.push_back(0);
1051 std::vector<double> factors;
1052 factors.resize(1.0, 1.0);
1053 PtrPreconditioner prec(new ConstantPreconditioner(factors));
1054 std::vector<int> vertexOffsets{2, 3, 3, 4};
1055
1056 mesh::PtrMesh dummyMesh(new mesh::Mesh("DummyMesh", 3, testing::nextMeshID()));
1057 dummyMesh->setVertexOffsets(vertexOffsets);
1058
1059 IQNILSAcceleration acc(initialRelaxation, enforceInitialRelaxation, maxIterationsUsed,
1060 timeWindowsReused, filter, singularityLimit, dataIDs, prec);
1061
1062 mesh::PtrData displacements(new mesh::Data("dvalues", -1, 1));
1063
1064 DataMap data;
1065
1066 Eigen::VectorXd insert;
1067 // init displacements
1068 if (context.isPrimary()) { // 2 vertices
1069 insert.resize(2);
1070 insert << 0.5, 0.5;
1071 } else if (context.isRank(1)) { // 1 vertex
1072 insert.resize(1);
1073 insert << 0.5;
1074 } else if (context.isRank(2)) { // no vertices
1075 } else { // 1 vertex
1076 insert.resize(1);
1077 insert << 0.5;
1078 }
1079
1080 utils::append(displacements->values(), insert);
1081
1082 bool exchangeSubsteps = false; // @todo add testColumnsLoggingWithSubsteps, where exchangeSubsteps = true as soon as acceleration scheme supports subcycling.
1083
1084 PtrCouplingData dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
1085 data.insert(std::pair<int, PtrCouplingData>(0, dpcd));
1086 dpcd->setSampleAtTime(0, dpcd->sample());
1087 dpcd->storeIteration();
1088
1089 acc.initialize(data);
1090
1091 // update displacements
1092 if (context.isPrimary()) { // 2 vertices
1093 insert << 1.0, 1.0;
1094 } else if (context.isRank(1)) { // 1 vertex
1095 insert << 1.0;
1096 } else if (context.isRank(2)) { // no vertices
1097 } else { // 1 vertex
1098 insert << 1.0;
1099 }
1100
1101 displacements->values() = insert;
1102
1103 acc.performAcceleration(data);
1104
1105 Eigen::VectorXd newdvalues1;
1106 if (context.isPrimary()) {
1107 utils::append(newdvalues1, 1.1);
1108 utils::append(newdvalues1, 1.0);
1109 } else if (context.isRank(1)) {
1110 utils::append(newdvalues1, 1.0);
1111 } else if (context.isRank(2)) {
1112 } else if (context.isRank(3)) {
1113 utils::append(newdvalues1, 1.0);
1114 }
1115 data.begin()->second->values() = newdvalues1;
1116
1117 acc.performAcceleration(data);
1118
1119 Eigen::VectorXd newdvalues2;
1120 if (context.isPrimary()) {
1121 utils::append(newdvalues2, 1.0);
1122 utils::append(newdvalues2, 2.0);
1123 } else if (context.isRank(1)) {
1124 utils::append(newdvalues2, 1.0);
1125 } else if (context.isRank(2)) {
1126 } else if (context.isRank(3)) {
1127 utils::append(newdvalues2, 1.0);
1128 }
1129 data.begin()->second->values() = newdvalues2;
1130
1131 acc.iterationsConverged(data);
1132
1133 BOOST_TEST(acc.getLSSystemCols() == 2);
1134 BOOST_TEST(acc.getDeletedColumns() == 0);
1135 BOOST_TEST(acc.getDroppedColumns() == 0);
1136
1137 Eigen::VectorXd newdvalues3;
1138 if (context.isPrimary()) {
1139 utils::append(newdvalues3, 1.1);
1140 utils::append(newdvalues3, 2.0);
1141 } else if (context.isRank(1)) {
1142 utils::append(newdvalues3, 1.0);
1143 } else if (context.isRank(2)) {
1144 } else if (context.isRank(3)) {
1145 utils::append(newdvalues3, 1.0);
1146 }
1147 data.begin()->second->values() = newdvalues3;
1148
1149 acc.performAcceleration(data);
1150
1151 Eigen::VectorXd newdvalues4;
1152 if (context.isPrimary()) {
1153 utils::append(newdvalues4, 1.0);
1154 utils::append(newdvalues4, 1.0);
1155 } else if (context.isRank(1)) {
1156 utils::append(newdvalues4, 1.0);
1157 } else if (context.isRank(2)) {
1158 } else if (context.isRank(3)) {
1159 utils::append(newdvalues4, 5.0);
1160 }
1161 data.begin()->second->values() = newdvalues4;
1162
1163 acc.iterationsConverged(data);
1164
1165 BOOST_TEST(acc.getLSSystemCols() == 1);
1166 BOOST_TEST(acc.getDeletedColumns() == 1);
1167 BOOST_TEST(acc.getDroppedColumns() == 1);
1168}
1169
1170BOOST_AUTO_TEST_SUITE_END()
1171BOOST_AUTO_TEST_SUITE_END()
1172
1173#endif // PRECICE_NO_MPI
BOOST_AUTO_TEST_CASE
BOOST_AUTO_TEST_CASE(testVIQNILSppWithoutSubsteps)
Test that runs on 4 processors.
Definition AccelerationIntraCommTest.cpp:43
BOOST_AUTO_TEST_SUITE
BOOST_AUTO_TEST_SUITE(PreProcess)
BOOST_AUTO_TEST_SUITE_END
BOOST_AUTO_TEST_SUITE_END()
PRECICE_TEST
#define PRECICE_TEST(...)
Definition Testing.hpp:27
precice::acceleration::BaseQNAcceleration::getLSSystemCols
virtual int getLSSystemCols() const
: computes number of cols in least squares system, i.e, number of cols in _matrixV,...
Definition BaseQNAcceleration.cpp:600
precice::acceleration::BaseQNAcceleration::iterationsConverged
virtual void iterationsConverged(const DataMap &cplData)
Marks a iteration sequence as converged.
Definition BaseQNAcceleration.cpp:462
precice::acceleration::BaseQNAcceleration::getDroppedColumns
virtual int getDroppedColumns() const
how many QN columns were dropped (went out of scope) in this time window
Definition BaseQNAcceleration.cpp:595
precice::acceleration::BaseQNAcceleration::performAcceleration
virtual void performAcceleration(DataMap &cplData)
Performs one acceleration step.
Definition BaseQNAcceleration.cpp:275
precice::acceleration::BaseQNAcceleration::getDeletedColumns
virtual int getDeletedColumns() const
how many QN columns were deleted in this time window
Definition BaseQNAcceleration.cpp:590
precice::acceleration::IQNILSAcceleration
Interface quasi-Newton with interface least-squares approximation.
Definition IQNILSAcceleration.hpp:26
precice::acceleration::IQNILSAcceleration::initialize
virtual void initialize(const DataMap &cplData)
Initializes the acceleration.
Definition IQNILSAcceleration.cpp:41
precice::acceleration::IQNIMVJAcceleration
Multi vector quasi-Newton update scheme.
Definition IQNIMVJAcceleration.hpp:33
precice::acceleration::IQNIMVJAcceleration::initialize
virtual void initialize(const DataMap &cplData)
Initializes the acceleration.
Definition IQNIMVJAcceleration.cpp:72
precice::acceleration::impl::ConstantPreconditioner
Preconditioner that uses the constant user-defined factors to scale the quasi-Newton system.
Definition ConstantPreconditioner.hpp:15
precice::acceleration::impl::ResidualSumPreconditioner
Preconditioner that uses the residuals of all iterations of the current time window summed up to scal...
Definition ResidualSumPreconditioner.hpp:18
precice::cplscheme::CouplingData::values
Eigen::VectorXd & values()
Returns a reference to the data values.
Definition CouplingData.cpp:42
precice::cplscheme::CouplingData::setSampleAtTime
void setSampleAtTime(double time, time::Sample sample)
Add sample at given time to _timeStepsStorage.
Definition CouplingData.cpp:84
precice::cplscheme::CouplingData::storeIteration
void storeIteration()
store _data->values() in read-only variable _previousIteration for convergence checks etc.
Definition CouplingData.cpp:102
precice::cplscheme::CouplingData::sample
time::Sample & sample()
Returns a reference to the gradient data Sample.
Definition CouplingData.cpp:163
precice::mesh::Data
Describes a set of data values belonging to the vertices of a mesh.
Definition Data.hpp:29
precice::mesh::Mesh
Container and creator for meshes.
Definition Mesh.hpp:39
precice::acceleration::impl::PtrPreconditioner
std::shared_ptr< Preconditioner > PtrPreconditioner
Definition SharedPointer.hpp:12
precice::acceleration
contains implementations of acceleration schemes.
Definition Acceleration.cpp:5
precice::cplscheme
contains implementations of coupling schemes for coupled simulations.
Definition BaseCouplingScheme.cpp:31
precice::testing::equals
boost::test_tools::predicate_result equals(const std::vector< float > &VectorA, const std::vector< float > &VectorB, float tolerance)
equals to be used in tests. Compares two std::vectors using a given tolerance. Prints both operands o...
Definition Testing.cpp:65
precice::testing::makeCouplingData
cplscheme::PtrCouplingData makeCouplingData(mesh::PtrData data, mesh::PtrMesh mesh, bool exchangeSubsteps)
Definition helper.hpp:10
precice::utils::append
void append(Eigen::VectorXd &v, double value)
Definition EigenHelperFunctions.cpp:48
precice
Main namespace of the precice library.
Definition Acceleration.cpp:5
std
STL namespace.
std::vector::push_back
T push_back(T... args)
std::vector::resize
T resize(T... args)