preCICE v3.2.0
Loading...
Searching...
No Matches
AccelerationIntraCommTest.cpp
Go to the documentation of this file.
1#include <Eigen/Core>
2#include <algorithm>
3#include <boost/test/data/test_case.hpp>
4#include <cmath>
5#include <map>
6#include <memory>
7#include <utility>
8#include <vector>
9#include "acceleration/Acceleration.hpp"
10#include "acceleration/BaseQNAcceleration.hpp"
11#include "acceleration/IQNILSAcceleration.hpp"
12#include "acceleration/IQNIMVJAcceleration.hpp"
13#include "acceleration/SharedPointer.hpp"
14#include "acceleration/impl/ConstantPreconditioner.hpp"
15#include "acceleration/impl/QRFactorization.hpp"
16#include "acceleration/impl/ResidualSumPreconditioner.hpp"
17#include "acceleration/impl/SharedPointer.hpp"
18#include "acceleration/test/helper.hpp"
19#include "cplscheme/Constants.hpp"
20#include "cplscheme/CouplingData.hpp"
21#include "cplscheme/SharedPointer.hpp"
22#include "mesh/Mesh.hpp"
23#include "mesh/SharedPointer.hpp"
24#include "testing/Meshes.hpp"
25#include "testing/TestContext.hpp"
26#include "testing/Testing.hpp"
27#include "utils/EigenHelperFunctions.hpp"
28
29using namespace precice;
30using namespace precice::cplscheme;
31using namespace precice::acceleration;
32using namespace precice::acceleration::impl;
33
34using precice::testing::makeCouplingData;
35
36#ifndef PRECICE_NO_MPI
37
38BOOST_AUTO_TEST_SUITE(AccelerationTests)
39
40using DataMap = std::map<int, PtrCouplingData>;
41
42BOOST_AUTO_TEST_SUITE(AccelerationIntraCommTests)
43
45PRECICE_TEST_SETUP(""_on(4_ranks).setupIntraComm())
46BOOST_DATA_TEST_CASE(testVIQNILSppWithoutSubsteps, boost::unit_test::data::make({true, false}), exchangeSubsteps)
47{
48 PRECICE_TEST()
49 double initialRelaxation = 0.01;
50 int maxIterationsUsed = 50;
51 int timeWindowsReused = 6;
52 int filter = BaseQNAcceleration::QR1FILTER;
53 double singularityLimit = 1e-10;
54 bool enforceInitialRelaxation = false;
55 std::vector<int> dataIDs;
56 dataIDs.push_back(0);
57 dataIDs.push_back(1);
58 std::vector<double> factors;
59 factors.resize(2, 1.0);
60 PtrPreconditioner prec(new ConstantPreconditioner(factors));
61 std::vector<int> vertexOffsets{4, 8, 8, 10};
62
63 mesh::PtrMesh dummyMesh(new mesh::Mesh("DummyMesh", 3, testing::nextMeshID()));
64 dummyMesh->setVertexOffsets(vertexOffsets);
65
66 IQNILSAcceleration pp(initialRelaxation, enforceInitialRelaxation, maxIterationsUsed,
67 timeWindowsReused, filter, singularityLimit, dataIDs, prec, !exchangeSubsteps);
68
69 Eigen::VectorXd dcol1;
70 Eigen::VectorXd fcol1;
71
72 DataMap data;
73 mesh::PtrData displacements(new mesh::Data("dvalues", -1, 1));
74 mesh::PtrData forces(new mesh::Data("fvalues", -1, 1));
75 PtrCouplingData dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
76 PtrCouplingData fpcd = makeCouplingData(forces, dummyMesh, exchangeSubsteps);
77
78 const double windowStart = 0;
79 const double windowEnd = 1;
80
81 data.insert(std::pair<int, PtrCouplingData>(0, dpcd));
82 data.insert(std::pair<int, PtrCouplingData>(1, fpcd));
83
84 if (context.isPrimary()) { // Primary
88 dummyMesh->createVertex(Eigen::Vector3d{0, 0, 0});
89 dummyMesh->createVertex(Eigen::Vector3d{1, 0, 0});
90 dummyMesh->createVertex(Eigen::Vector3d{2, 0, 0});
91 dummyMesh->createVertex(Eigen::Vector3d{3, 0, 0});
92
93 // init displacements & forces
94 dpcd->emplaceSampleAtTime(windowStart, {1.0, 1.0, 1.0, 1.0});
95 dpcd->emplaceSampleAtTime(windowEnd, {1.0, 1.0, 1.0, 1.0});
96
97 fpcd->emplaceSampleAtTime(windowStart, {0.2, 0.2, 0.2, 0.2});
98 fpcd->emplaceSampleAtTime(windowEnd, {0.2, 0.2, 0.2, 0.2});
99
100 dpcd->storeIteration();
101 fpcd->storeIteration();
102
103 pp.initialize(data);
104
105 dpcd->emplaceSampleAtTime(windowEnd, {1.0, 2.0, 3.0, 4.0});
106 fpcd->emplaceSampleAtTime(windowEnd, {0.1, 0.1, 0.1, 0.1});
107
108 } else if (context.isRank(1)) { // SecondaryRank1
109
113 dummyMesh->createVertex(Eigen::Vector3d{0, 0, 0});
114 dummyMesh->createVertex(Eigen::Vector3d{1, 0, 0});
115 dummyMesh->createVertex(Eigen::Vector3d{2, 0, 0});
116 dummyMesh->createVertex(Eigen::Vector3d{3, 0, 0});
117
118 // init displacements & forces
119 dpcd->emplaceSampleAtTime(windowStart, {1.0, 1.0, 1.0, 1.0});
120 dpcd->emplaceSampleAtTime(windowEnd, {1.0, 1.0, 1.0, 1.0});
121
122 fpcd->emplaceSampleAtTime(windowStart, {0.2, 0.2, 0.2, 0.2});
123 fpcd->emplaceSampleAtTime(windowEnd, {0.2, 0.2, 0.2, 0.2});
124
125 dpcd->storeIteration();
126 fpcd->storeIteration();
127
128 pp.initialize(data);
129
130 dpcd->emplaceSampleAtTime(windowEnd, {5.0, 6.0, 7.0, 8.0});
131 fpcd->emplaceSampleAtTime(windowEnd, {0.1, 0.1, 0.1, 0.1});
132
133 } else if (context.isRank(2)) { // Secondary rank 2
134
138
139 // init displacements
140 dpcd->emplaceSampleAtTime(windowStart);
141 dpcd->emplaceSampleAtTime(windowEnd);
142
143 // init forces
144 fpcd->emplaceSampleAtTime(windowStart);
145 fpcd->emplaceSampleAtTime(windowEnd);
146
147 data.insert(std::pair<int, PtrCouplingData>(0, dpcd));
148 data.insert(std::pair<int, PtrCouplingData>(1, fpcd));
149
150 dpcd->storeIteration();
151 fpcd->storeIteration();
152
153 pp.initialize(data);
154
155 dpcd->emplaceSampleAtTime(windowEnd);
156 fpcd->emplaceSampleAtTime(windowEnd);
157
158 } else if (context.isRank(3)) { // Secondary rank 3
159
163 dummyMesh->createVertex(Eigen::Vector3d{0, 0, 0});
164 dummyMesh->createVertex(Eigen::Vector3d{1, 0, 0});
165
166 // init displacements & forces
167 dpcd->emplaceSampleAtTime(windowStart, {1.0, 1.0});
168 dpcd->emplaceSampleAtTime(windowEnd, {1.0, 1.0});
169
170 fpcd->emplaceSampleAtTime(windowStart, {0.2, 0.2});
171 fpcd->emplaceSampleAtTime(windowEnd, {0.2, 0.2});
172
173 data.insert(std::pair<int, PtrCouplingData>(0, dpcd));
174 data.insert(std::pair<int, PtrCouplingData>(1, fpcd));
175
176 dpcd->storeIteration();
177 fpcd->storeIteration();
178
179 pp.initialize(data);
180
181 dpcd->emplaceSampleAtTime(windowEnd, {1.0, 2.0});
182 fpcd->emplaceSampleAtTime(windowEnd, {0.1, 0.1});
183 }
184
185 pp.performAcceleration(data, windowStart, windowEnd);
186
187 Eigen::VectorXd newdvalues;
188 if (context.isPrimary()) { // Primary
189
190 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(0), 1.00), dpcd->timeStepsStorage().sample(windowEnd)(0));
191 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(1), 1.01), dpcd->timeStepsStorage().sample(windowEnd)(1));
192 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(2), 1.02), dpcd->timeStepsStorage().sample(windowEnd)(2));
193 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(3), 1.03), dpcd->timeStepsStorage().sample(windowEnd)(3));
194 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(0), 0.199), fpcd->timeStepsStorage().sample(windowEnd)(0));
195 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(1), 0.199), fpcd->timeStepsStorage().sample(windowEnd)(1));
196 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(2), 0.199), fpcd->timeStepsStorage().sample(windowEnd)(2));
197 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(3), 0.199), fpcd->timeStepsStorage().sample(windowEnd)(3));
198
199 utils::append(newdvalues, 10.0);
200 utils::append(newdvalues, 10.0);
201 utils::append(newdvalues, 10.0);
202 utils::append(newdvalues, 10.0);
203
204 } else if (context.isRank(1)) { // SecondaryRank1
205
206 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(0), 1.04), dpcd->timeStepsStorage().sample(windowEnd)(0));
207 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(1), 1.05), dpcd->timeStepsStorage().sample(windowEnd)(1));
208 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(2), 1.06), dpcd->timeStepsStorage().sample(windowEnd)(2));
209 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(3), 1.07), dpcd->timeStepsStorage().sample(windowEnd)(3));
210 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(0), 0.199), fpcd->timeStepsStorage().sample(windowEnd)(0));
211 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(1), 0.199), fpcd->timeStepsStorage().sample(windowEnd)(1));
212 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(2), 0.199), fpcd->timeStepsStorage().sample(windowEnd)(2));
213 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(3), 0.199), fpcd->timeStepsStorage().sample(windowEnd)(3));
214
215 utils::append(newdvalues, 10.0);
216 utils::append(newdvalues, 10.0);
217 utils::append(newdvalues, 10.0);
218 utils::append(newdvalues, 10.0);
219
220 } else if (context.isRank(2)) { // Secondary rank 2
221 // empty proc
222 } else if (context.isRank(3)) { // Secondary rank 3
223
224 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(0), 1.00), dpcd->timeStepsStorage().sample(windowEnd)(0));
225 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(1), 1.01), dpcd->timeStepsStorage().sample(windowEnd)(1));
226 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(0), 0.199), fpcd->timeStepsStorage().sample(windowEnd)(0));
227 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(1), 0.199), fpcd->timeStepsStorage().sample(windowEnd)(1));
228
229 utils::append(newdvalues, 10.0);
230 utils::append(newdvalues, 10.0);
231 }
232
233 dpcd->setSampleAtTime(windowEnd, time::Sample(dpcd->getDimensions(), newdvalues));
234 fpcd->setSampleAtTime(windowEnd, fpcd->sample());
235
236 pp.performAcceleration(data, windowStart, windowEnd);
237
238 if (context.isPrimary()) { // Primary
239 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(0), -1.51483105223442748866e+00), dpcd->timeStepsStorage().sample(windowEnd)(0));
240 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(1), -2.35405379763935940218e-01), dpcd->timeStepsStorage().sample(windowEnd)(1));
241 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(2), 1.04402029270655560822e+00), dpcd->timeStepsStorage().sample(windowEnd)(2));
242 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(3), 2.32344596517704804484e+00), dpcd->timeStepsStorage().sample(windowEnd)(3));
243 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(0), 7.23368584254212854123e-02), fpcd->timeStepsStorage().sample(windowEnd)(0));
244 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(1), 7.23368584254212854123e-02), fpcd->timeStepsStorage().sample(windowEnd)(1));
245 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(2), 7.23368584254212854123e-02), fpcd->timeStepsStorage().sample(windowEnd)(2));
246 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(3), 7.23368584254212854123e-02), fpcd->timeStepsStorage().sample(windowEnd)(3));
247 } else if (context.isRank(1)) { // SecondaryRank1
248 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(0), 3.60287163764754048145e+00), dpcd->timeStepsStorage().sample(windowEnd)(0));
249 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(1), 4.88229731011803202989e+00), dpcd->timeStepsStorage().sample(windowEnd)(1));
250 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(2), 6.16172298258852357833e+00), dpcd->timeStepsStorage().sample(windowEnd)(2));
251 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(3), 7.44114865505901601495e+00), dpcd->timeStepsStorage().sample(windowEnd)(3));
252 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(0), 7.23368584254212854123e-02), fpcd->timeStepsStorage().sample(windowEnd)(0));
253 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(1), 7.23368584254212854123e-02), fpcd->timeStepsStorage().sample(windowEnd)(1));
254 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(2), 7.23368584254212854123e-02), fpcd->timeStepsStorage().sample(windowEnd)(2));
255 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(3), 7.23368584254212854123e-02), fpcd->timeStepsStorage().sample(windowEnd)(3));
256 } else if (context.isRank(2)) { // Secondary rank 2
257 // empty proc
258 } else if (context.isRank(3)) { // Secondary rank 3
259 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(0), -1.51483105223442748866e+00), dpcd->timeStepsStorage().sample(windowEnd)(0));
260 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(1), -2.35405379763935940218e-01), dpcd->timeStepsStorage().sample(windowEnd)(1));
261 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(0), 7.23368584254212854123e-02), fpcd->timeStepsStorage().sample(windowEnd)(0));
262 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(1), 7.23368584254212854123e-02), fpcd->timeStepsStorage().sample(windowEnd)(1));
263 }
264}
265
267PRECICE_TEST_SETUP(""_on(4_ranks).setupIntraComm())
268BOOST_AUTO_TEST_CASE(testVIQNIMVJppWithoutSubsteps)
269{
270 PRECICE_TEST();
271 double initialRelaxation = 0.01;
272 int maxIterationsUsed = 50;
273 int timeWindowsReused = 6;
274 int filter = BaseQNAcceleration::QR1FILTER;
275 int restartType = IQNIMVJAcceleration::NO_RESTART;
276 int chunkSize = 0;
277 int reusedTimeWindowsAtRestart = 0;
278 double singularityLimit = 1e-10;
279 double svdTruncationEps = 0.0;
280 bool enforceInitialRelaxation = false;
281 bool alwaysBuildJacobian = false;
282 bool exchangeSubsteps = false; // @todo add testVIQNIMVJppWithSubsteps, where exchangeSubsteps = true as soon as acceleration scheme supports subcycling.
283
284 std::vector<int> dataIDs;
285 dataIDs.push_back(0);
286 dataIDs.push_back(1);
287 std::vector<double> factors;
288 factors.resize(2, 1.0);
289 PtrPreconditioner prec(new ConstantPreconditioner(factors));
290 std::vector<int> vertexOffsets{4, 8, 8, 10};
291
292 mesh::PtrMesh dummyMesh(new mesh::Mesh("DummyMesh", 3, testing::nextMeshID()));
293 dummyMesh->setVertexOffsets(vertexOffsets);
294
295 IQNIMVJAcceleration pp(initialRelaxation, enforceInitialRelaxation, maxIterationsUsed,
296 timeWindowsReused, filter, singularityLimit, dataIDs, prec, alwaysBuildJacobian,
297 restartType, chunkSize, reusedTimeWindowsAtRestart, svdTruncationEps, !exchangeSubsteps);
298
299 Eigen::VectorXd dcol1;
300 Eigen::VectorXd fcol1;
301
302 mesh::PtrData displacements(new mesh::Data("dvalues", -1, 1));
303 mesh::PtrData forces(new mesh::Data("fvalues", -1, 1));
304 PtrCouplingData dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
305 PtrCouplingData fpcd = makeCouplingData(forces, dummyMesh, exchangeSubsteps);
306
307 DataMap data;
308
309 const double windowStart = 0;
310 const double windowEnd = 1;
311
312 if (context.isPrimary()) { // Primary
313
317 dummyMesh->createVertex(Eigen::Vector3d{0, 0, 0});
318 dummyMesh->createVertex(Eigen::Vector3d{1, 0, 0});
319 dummyMesh->createVertex(Eigen::Vector3d{2, 0, 0});
320 dummyMesh->createVertex(Eigen::Vector3d{3, 0, 0});
321
322 // init displacements & forces
323 // Need to store 2 values in the waveform iteration
324 dpcd->emplaceSampleAtTime(windowStart, {1.0, 1.0, 1.0, 1.0});
325 dpcd->emplaceSampleAtTime(windowEnd, {1.0, 1.0, 1.0, 1.0});
326 fpcd->emplaceSampleAtTime(windowStart, {0.2, 0.2, 0.2, 0.2});
327 fpcd->emplaceSampleAtTime(windowEnd, {0.2, 0.2, 0.2, 0.2});
328
329 dpcd->storeIteration();
330 fpcd->storeIteration();
331
332 data.insert(std::pair<int, PtrCouplingData>(0, dpcd));
333 data.insert(std::pair<int, PtrCouplingData>(1, fpcd));
334
335 pp.initialize(data);
336
337 // update displacement & forces
338 dpcd->emplaceSampleAtTime(windowEnd, {1.0, 2.0, 3.0, 4.0});
339 fpcd->emplaceSampleAtTime(windowEnd, {0.1, 0.1, 0.1, 0.1});
340
341 // check for correct initial data
342
343 BOOST_TEST(testing::equals(dpcd->getPreviousValuesAtTime(windowEnd)(0), 1.0), dpcd->getPreviousValuesAtTime(windowEnd)(0));
344 BOOST_TEST(testing::equals(dpcd->getPreviousValuesAtTime(windowEnd)(1), 1.0), dpcd->getPreviousValuesAtTime(windowEnd)(1));
345 BOOST_TEST(testing::equals(dpcd->getPreviousValuesAtTime(windowEnd)(2), 1.0), dpcd->getPreviousValuesAtTime(windowEnd)(2));
346 BOOST_TEST(testing::equals(dpcd->getPreviousValuesAtTime(windowEnd)(3), 1.0), dpcd->getPreviousValuesAtTime(windowEnd)(3));
347 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(0), 1.0), dpcd->timeStepsStorage().sample(windowEnd)(0));
348 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(1), 2.0), dpcd->timeStepsStorage().sample(windowEnd)(1));
349 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(2), 3.0), dpcd->timeStepsStorage().sample(windowEnd)(2));
350 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(3), 4.0), dpcd->timeStepsStorage().sample(windowEnd)(3));
351 BOOST_TEST(testing::equals(fpcd->getPreviousValuesAtTime(windowEnd)(0), 0.2), fpcd->getPreviousValuesAtTime(windowEnd)(0));
352 BOOST_TEST(testing::equals(fpcd->getPreviousValuesAtTime(windowEnd)(1), 0.2), fpcd->getPreviousValuesAtTime(windowEnd)(1));
353 BOOST_TEST(testing::equals(fpcd->getPreviousValuesAtTime(windowEnd)(2), 0.2), fpcd->getPreviousValuesAtTime(windowEnd)(2));
354 BOOST_TEST(testing::equals(fpcd->getPreviousValuesAtTime(windowEnd)(3), 0.2), fpcd->getPreviousValuesAtTime(windowEnd)(3));
355 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(0), 0.1), fpcd->timeStepsStorage().sample(windowEnd)(0));
356 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(1), 0.1), fpcd->timeStepsStorage().sample(windowEnd)(1));
357 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(2), 0.1), fpcd->timeStepsStorage().sample(windowEnd)(2));
358 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(3), 0.1), fpcd->timeStepsStorage().sample(windowEnd)(3));
359 } else if (context.isRank(1)) { // SecondaryRank1
360
364 dummyMesh->createVertex(Eigen::Vector3d{0, 0, 0});
365 dummyMesh->createVertex(Eigen::Vector3d{1, 0, 0});
366 dummyMesh->createVertex(Eigen::Vector3d{2, 0, 0});
367 dummyMesh->createVertex(Eigen::Vector3d{3, 0, 0});
368
369 // init displacements & forces
370 dpcd->emplaceSampleAtTime(windowStart, {1.0, 1.0, 1.0, 1.0});
371 dpcd->emplaceSampleAtTime(windowEnd, {1.0, 1.0, 1.0, 1.0});
372 fpcd->emplaceSampleAtTime(windowStart, {0.2, 0.2, 0.2, 0.2});
373 fpcd->emplaceSampleAtTime(windowEnd, {0.2, 0.2, 0.2, 0.2});
374
375 dpcd->storeIteration();
376 fpcd->storeIteration();
377
378 data.insert(std::pair<int, PtrCouplingData>(0, dpcd));
379 data.insert(std::pair<int, PtrCouplingData>(1, fpcd));
380
381 pp.initialize(data);
382
383 // update displacements & forces
384 dpcd->emplaceSampleAtTime(windowEnd, {5.0, 6.0, 7.0, 8.0});
385 fpcd->emplaceSampleAtTime(windowEnd, {0.1, 0.1, 0.1, 0.1});
386
387 // check for correct initial data
388 BOOST_TEST(testing::equals(dpcd->getPreviousValuesAtTime(windowEnd)(0), 1.0), dpcd->getPreviousValuesAtTime(windowEnd)(0));
389 BOOST_TEST(testing::equals(dpcd->getPreviousValuesAtTime(windowEnd)(1), 1.0), dpcd->getPreviousValuesAtTime(windowEnd)(1));
390 BOOST_TEST(testing::equals(dpcd->getPreviousValuesAtTime(windowEnd)(2), 1.0), dpcd->getPreviousValuesAtTime(windowEnd)(2));
391 BOOST_TEST(testing::equals(dpcd->getPreviousValuesAtTime(windowEnd)(3), 1.0), dpcd->getPreviousValuesAtTime(windowEnd)(3));
392 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(0), 5.0), dpcd->timeStepsStorage().sample(windowEnd)(0));
393 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(1), 6.0), dpcd->timeStepsStorage().sample(windowEnd)(1));
394 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(2), 7.0), dpcd->timeStepsStorage().sample(windowEnd)(2));
395 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(3), 8.0), dpcd->timeStepsStorage().sample(windowEnd)(3));
396 BOOST_TEST(testing::equals(fpcd->getPreviousValuesAtTime(windowEnd)(0), 0.2), fpcd->getPreviousValuesAtTime(windowEnd)(0));
397 BOOST_TEST(testing::equals(fpcd->getPreviousValuesAtTime(windowEnd)(1), 0.2), fpcd->getPreviousValuesAtTime(windowEnd)(1));
398 BOOST_TEST(testing::equals(fpcd->getPreviousValuesAtTime(windowEnd)(2), 0.2), fpcd->getPreviousValuesAtTime(windowEnd)(2));
399 BOOST_TEST(testing::equals(fpcd->getPreviousValuesAtTime(windowEnd)(3), 0.2), fpcd->getPreviousValuesAtTime(windowEnd)(3));
400 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(0), 0.1), fpcd->timeStepsStorage().sample(windowEnd)(0));
401 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(1), 0.1), fpcd->timeStepsStorage().sample(windowEnd)(1));
402 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(2), 0.1), fpcd->timeStepsStorage().sample(windowEnd)(2));
403 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(3), 0.1), fpcd->timeStepsStorage().sample(windowEnd)(3));
404 } else if (context.isRank(2)) { // Secondary rank 2
405
409
410 // init displacements
411 dpcd->emplaceSampleAtTime(windowStart, {});
412 dpcd->emplaceSampleAtTime(windowEnd, {});
413
414 // init forces
415 fpcd->emplaceSampleAtTime(windowStart, {});
416 fpcd->emplaceSampleAtTime(windowEnd, {});
417
418 dpcd->storeIteration();
419 fpcd->storeIteration();
420
421 data.insert(std::pair<int, PtrCouplingData>(0, dpcd));
422 data.insert(std::pair<int, PtrCouplingData>(1, fpcd));
423
424 pp.initialize(data);
425
426 // check for correct initial data
427 // empty proc
428 } else if (context.isRank(3)) { // Secondary rank 3
429
433 dummyMesh->createVertex(Eigen::Vector3d{0, 0, 0});
434 dummyMesh->createVertex(Eigen::Vector3d{1, 0, 0});
435
436 // init displacements & forces
437 dpcd->emplaceSampleAtTime(windowStart, {1.0, 1.0});
438 dpcd->emplaceSampleAtTime(windowEnd, {1.0, 1.0});
439 fpcd->emplaceSampleAtTime(windowStart, {0.2, 0.2});
440 fpcd->emplaceSampleAtTime(windowEnd, {0.2, 0.2});
441
442 dpcd->storeIteration();
443 fpcd->storeIteration();
444
445 data.insert(std::pair<int, PtrCouplingData>(0, dpcd));
446 data.insert(std::pair<int, PtrCouplingData>(1, fpcd));
447
448 pp.initialize(data);
449
450 // update displacements & forces
451 dpcd->emplaceSampleAtTime(windowEnd, {1.0, 2.0});
452 fpcd->emplaceSampleAtTime(windowEnd, {0.1, 0.1});
453
454 // check for correct initial data
455 BOOST_TEST(testing::equals(dpcd->getPreviousValuesAtTime(windowEnd)(0), 1.0), dpcd->getPreviousValuesAtTime(windowEnd)(0));
456 BOOST_TEST(testing::equals(dpcd->getPreviousValuesAtTime(windowEnd)(1), 1.0), dpcd->getPreviousValuesAtTime(windowEnd)(1));
457 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(0), 1.0), dpcd->timeStepsStorage().sample(windowEnd)(0));
458 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(1), 2.0), dpcd->timeStepsStorage().sample(windowEnd)(1));
459 BOOST_TEST(testing::equals(fpcd->getPreviousValuesAtTime(windowEnd)(0), 0.2), fpcd->getPreviousValuesAtTime(windowEnd)(0));
460 BOOST_TEST(testing::equals(fpcd->getPreviousValuesAtTime(windowEnd)(1), 0.2), fpcd->getPreviousValuesAtTime(windowEnd)(1));
461 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(0), 0.1), fpcd->timeStepsStorage().sample(windowEnd)(0));
462 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(1), 0.1), fpcd->timeStepsStorage().sample(windowEnd)(1));
463 }
464
465 pp.performAcceleration(data, windowStart, windowEnd);
466
467 Eigen::VectorXd newdvalues;
468 if (context.isPrimary()) { // Primary
469 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(0), 1.00000000000000000000e+00), dpcd->timeStepsStorage().sample(windowEnd)(0));
470 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(1), 1.01000000000000000888e+00), dpcd->timeStepsStorage().sample(windowEnd)(1));
471 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(2), 1.02000000000000001776e+00), dpcd->timeStepsStorage().sample(windowEnd)(2));
472 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(3), 1.03000000000000002665e+00), dpcd->timeStepsStorage().sample(windowEnd)(3));
473 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(0), 1.99000000000000010214e-01), fpcd->timeStepsStorage().sample(windowEnd)(0));
474 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(1), 1.99000000000000010214e-01), fpcd->timeStepsStorage().sample(windowEnd)(1));
475 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(2), 1.99000000000000010214e-01), fpcd->timeStepsStorage().sample(windowEnd)(2));
476 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(3), 1.99000000000000010214e-01), fpcd->timeStepsStorage().sample(windowEnd)(3));
477 utils::append(newdvalues, 10.0);
478 utils::append(newdvalues, 10.0);
479 utils::append(newdvalues, 10.0);
480 utils::append(newdvalues, 10.0);
481 } else if (context.isRank(1)) { // SecondaryRank1
482 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(0), 1.04000000000000003553e+00), dpcd->timeStepsStorage().sample(windowEnd)(0));
483 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(1), 1.05000000000000004441e+00), dpcd->timeStepsStorage().sample(windowEnd)(1));
484 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(2), 1.06000000000000005329e+00), dpcd->timeStepsStorage().sample(windowEnd)(2));
485 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(3), 1.07000000000000006217e+00), dpcd->timeStepsStorage().sample(windowEnd)(3));
486 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(0), 1.99000000000000010214e-01), fpcd->timeStepsStorage().sample(windowEnd)(0));
487 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(1), 1.99000000000000010214e-01), fpcd->timeStepsStorage().sample(windowEnd)(1));
488 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(2), 1.99000000000000010214e-01), fpcd->timeStepsStorage().sample(windowEnd)(2));
489 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(3), 1.99000000000000010214e-01), fpcd->timeStepsStorage().sample(windowEnd)(3));
490 utils::append(newdvalues, 10.0);
491 utils::append(newdvalues, 10.0);
492 utils::append(newdvalues, 10.0);
493 utils::append(newdvalues, 10.0);
494 } else if (context.isRank(2)) { // Secondary rank 2
495 // empty proc
496 } else if (context.isRank(3)) { // Secondary rank 3
497 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(0), 1.00000000000000000000e+00), dpcd->timeStepsStorage().sample(windowEnd)(0));
498 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(1), 1.01000000000000000888e+00), dpcd->timeStepsStorage().sample(windowEnd)(1));
499 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(0), 1.99000000000000010214e-01), fpcd->timeStepsStorage().sample(windowEnd)(0));
500 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(1), 1.99000000000000010214e-01), fpcd->timeStepsStorage().sample(windowEnd)(1));
501 utils::append(newdvalues, 10.0);
502 utils::append(newdvalues, 10.0);
503 }
504
505 dpcd->setSampleAtTime(windowEnd, time::Sample(dpcd->getDimensions(), newdvalues));
506 fpcd->setSampleAtTime(windowEnd, fpcd->sample());
507
508 pp.performAcceleration(data, windowStart, windowEnd);
509
510 if (context.isPrimary()) { // Primary
511 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(0), -1.51483105223442748866e+00), dpcd->timeStepsStorage().sample(windowEnd)(0));
512 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(1), -2.35405379763935940218e-01), dpcd->timeStepsStorage().sample(windowEnd)(1));
513 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(2), 1.04402029270655738458e+00), dpcd->timeStepsStorage().sample(windowEnd)(2));
514 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(3), 2.32344596517704893301e+00), dpcd->timeStepsStorage().sample(windowEnd)(3));
515 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(0), 7.23368584254213131679e-02), fpcd->timeStepsStorage().sample(windowEnd)(0));
516 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(1), 7.23368584254213131679e-02), fpcd->timeStepsStorage().sample(windowEnd)(1));
517 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(2), 7.23368584254213131679e-02), fpcd->timeStepsStorage().sample(windowEnd)(2));
518 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(3), 7.23368584254213131679e-02), fpcd->timeStepsStorage().sample(windowEnd)(3));
519 } else if (context.isRank(1)) { // SecondaryRank1
520 BOOST_TEST(testing::equals(dpcd->values()(0), 3.60287163764754048145e+00), dpcd->values()(0));
521 BOOST_TEST(testing::equals(dpcd->values()(1), 4.88229731011803202989e+00), dpcd->values()(1));
522 BOOST_TEST(testing::equals(dpcd->values()(2), 6.16172298258852446651e+00), dpcd->values()(2));
523 BOOST_TEST(testing::equals(dpcd->values()(3), 7.44114865505901601495e+00), dpcd->values()(3));
524 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(0), 7.23368584254213131679e-02), fpcd->timeStepsStorage().sample(windowEnd)(0));
525 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(1), 7.23368584254213131679e-02), fpcd->timeStepsStorage().sample(windowEnd)(1));
526 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(2), 7.23368584254213131679e-02), fpcd->timeStepsStorage().sample(windowEnd)(2));
527 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(3), 7.23368584254213131679e-02), fpcd->timeStepsStorage().sample(windowEnd)(3));
528 } else if (context.isRank(2)) { // Secondary rank 2
529 // empty proc
530 } else if (context.isRank(3)) { // Secondary rank 3
531 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(0), -1.51483105223442748866e+00), dpcd->timeStepsStorage().sample(windowEnd)(0));
532 BOOST_TEST(testing::equals(dpcd->timeStepsStorage().sample(windowEnd)(1), -2.35405379763935940218e-01), dpcd->timeStepsStorage().sample(windowEnd)(1));
533 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(0), 7.23368584254213131679e-02), fpcd->timeStepsStorage().sample(windowEnd)(0));
534 BOOST_TEST(testing::equals(fpcd->timeStepsStorage().sample(windowEnd)(1), 7.23368584254213131679e-02), fpcd->timeStepsStorage().sample(windowEnd)(1));
535 }
536}
537
539PRECICE_TEST_SETUP(""_on(4_ranks).setupIntraComm())
540BOOST_AUTO_TEST_CASE(testIMVJ_effUpdate_ppWithoutSubsteps)
541{
542 PRECICE_TEST();
543 // config:
544 double initialRelaxation = 0.1;
545 int maxIterationsUsed = 30;
546 int timeWindowsReused = 0;
547 int filter = BaseQNAcceleration::QR2FILTER;
548 int restartType = IQNIMVJAcceleration::NO_RESTART;
549 int chunkSize = 0;
550 int reusedTimeWindowsAtRestart = 0;
551 double singularityLimit = 1e-2;
552 double svdTruncationEps = 0.0;
553 bool enforceInitialRelaxation = false;
554 bool alwaysBuildJacobian = false;
555 bool exchangeSubsteps = false; // @todo add testIMVJ_effUpdate_ppWithSubsteps, where exchangeSubsteps = true as soon as acceleration scheme supports subcycling.
556
557 std::vector<int> dataIDs;
558 dataIDs.push_back(4);
559 dataIDs.push_back(5);
560 PtrPreconditioner _preconditioner = PtrPreconditioner(new ResidualSumPreconditioner(-1, true));
561 std::vector<int> vertexOffsets{0, 11, 22, 22};
562
563 mesh::PtrMesh dummyMesh(new mesh::Mesh("dummyMesh", 2, testing::nextMeshID()));
564 dummyMesh->setVertexOffsets(vertexOffsets);
565
566 IQNIMVJAcceleration pp(initialRelaxation, enforceInitialRelaxation, maxIterationsUsed,
567 timeWindowsReused, filter, singularityLimit, dataIDs, _preconditioner, alwaysBuildJacobian,
568 restartType, chunkSize, reusedTimeWindowsAtRestart, svdTruncationEps, !exchangeSubsteps);
569
570 mesh::PtrData displacements(new mesh::Data("dvalues", -1, 2));
571 mesh::PtrData forces(new mesh::Data("fvalues", -1, 2));
572
573 Eigen::VectorXd dref;
574 Eigen::VectorXd fref;
575 double drefNorm = 0., frefNorm = 0.;
576
577 DataMap data;
578 PtrCouplingData dpcd;
579 PtrCouplingData fpcd;
580
581 const double windowStart = 0;
582 const double windowEnd = 1;
583
584 if (context.isPrimary()) { // Primary
588
589 displacements->values().resize(0);
590 forces->values().resize(0);
591
592 // init displacements & forces
593 dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
594 dpcd->emplaceSampleAtTime(windowStart, {});
595 dpcd->emplaceSampleAtTime(windowEnd, {});
596
597 // init forces
598 fpcd = makeCouplingData(forces, dummyMesh, exchangeSubsteps);
599 fpcd->emplaceSampleAtTime(windowStart, {});
600 fpcd->emplaceSampleAtTime(windowEnd, {});
601
602 dpcd->storeIteration();
603 fpcd->storeIteration();
604
605 data.insert(std::pair<int, PtrCouplingData>(4, dpcd));
606 data.insert(std::pair<int, PtrCouplingData>(5, fpcd));
607
608 pp.initialize(data);
609 } else if (context.isRank(1)) { // SecondaryRank1
613 testing::addDummyVertices(11, *dummyMesh);
614
615 // init displacements & forces
616 dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
617 dpcd->emplaceSampleAtTime(windowStart, {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0});
618 dpcd->emplaceSampleAtTime(windowEnd, {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0});
619
620 fpcd = makeCouplingData(forces, dummyMesh, exchangeSubsteps);
621 fpcd->emplaceSampleAtTime(windowStart, {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0});
622 fpcd->emplaceSampleAtTime(windowEnd, {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0});
623
624 dpcd->storeIteration();
625 fpcd->storeIteration();
626
627 data.insert(std::pair<int, PtrCouplingData>(4, dpcd));
628 data.insert(std::pair<int, PtrCouplingData>(5, fpcd));
629
630 pp.initialize(data);
631
632 fpcd->emplaceSampleAtTime(windowEnd, {-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});
633 dpcd->setSampleAtTime(windowEnd, dpcd->sample());
634
635 } else if (context.isRank(2)) { // Secondary rank 2
639 testing::addDummyVertices(11, *dummyMesh);
640
641 // init displacements & forces
642 dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
643 dpcd->emplaceSampleAtTime(windowStart, {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0});
644 dpcd->emplaceSampleAtTime(windowEnd, {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0});
645
646 fpcd = makeCouplingData(forces, dummyMesh, exchangeSubsteps);
647 fpcd->emplaceSampleAtTime(windowStart, {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0});
648 fpcd->emplaceSampleAtTime(windowEnd, {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0});
649
650 dpcd->storeIteration();
651 fpcd->storeIteration();
652
653 data.insert(std::pair<int, PtrCouplingData>(4, dpcd));
654 data.insert(std::pair<int, PtrCouplingData>(5, fpcd));
655
656 pp.initialize(data);
657
658 fpcd->emplaceSampleAtTime(windowEnd, {-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});
659 dpcd->setSampleAtTime(windowEnd, dpcd->sample());
660
661 } else if (context.isRank(3)) { // Secondary rank 3
665
666 displacements->values().resize(0);
667 forces->values().resize(0);
668
669 dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
670 dpcd->setSampleAtTime(windowStart, dpcd->sample());
671 dpcd->setSampleAtTime(windowEnd, dpcd->sample());
672
673 fpcd = makeCouplingData(forces, dummyMesh, exchangeSubsteps);
674
675 fpcd->setSampleAtTime(windowStart, fpcd->sample());
676 fpcd->setSampleAtTime(windowEnd, fpcd->sample());
677
678 dpcd->storeIteration();
679 fpcd->storeIteration();
680
681 data.insert(std::pair<int, PtrCouplingData>(4, dpcd));
682 data.insert(std::pair<int, PtrCouplingData>(5, fpcd));
683
684 pp.initialize(data);
685 }
686
687 pp.performAcceleration(data, windowStart, windowEnd);
688
689 // 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()
690 for (auto &pair : data) {
691 pair.second->setSampleAtTime(windowEnd, pair.second->sample());
692 }
693
694 // underrelaxation, first iteration
695
696 if (context.isPrimary()) { // Primary
697
698 } else if (context.isRank(1)) { // SecondaryRank1
699
700 dref = Eigen::VectorXd::Zero(22);
701 fref = Eigen::VectorXd::Zero(22);
702
703 dref << 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0;
704 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;
705 frefNorm = 0.01286041129960619;
706 drefNorm = 0.0;
707
708 // validate values
709 for (int i = 0; i < data.at(4)->values().size(); i++)
710 BOOST_TEST(testing::equals(data.at(4)->values()(i), dref(i)));
711
712 for (int i = 0; i < data.at(5)->values().size(); i++)
713 BOOST_TEST(testing::equals(data.at(5)->values()(i), fref(i)));
714
715 // validate norm
716 BOOST_TEST(testing::equals(data.at(4)->values().norm(), drefNorm), data.at(4)->values().norm());
717 BOOST_TEST(testing::equals(data.at(5)->values().norm(), frefNorm), data.at(5)->values().norm());
718
719 // update cplData
720 fpcd->storeIteration();
721
722 fpcd->emplaceSampleAtTime(windowEnd, {-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});
723 dpcd->emplaceSampleAtTime(windowEnd, {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});
724
725 } else if (context.isRank(2)) { // Secondary rank 2
726
727 dref = Eigen::VectorXd::Zero(22);
728 fref = Eigen::VectorXd::Zero(22);
729
730 dref << 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0;
731 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;
732 frefNorm = 0.01265754337961098;
733 drefNorm = 0.0;
734
735 // validate values
736 for (int i = 0; i < data.at(4)->values().size(); i++)
737 BOOST_TEST(testing::equals(data.at(4)->values()(i), dref(i)));
738
739 for (int i = 0; i < data.at(5)->values().size(); i++)
740 BOOST_TEST(testing::equals(data.at(5)->values()(i), fref(i)));
741
742 // validate norm
743 BOOST_TEST(testing::equals(data.at(4)->values().norm(), drefNorm), data.at(4)->values().norm());
744 BOOST_TEST(testing::equals(data.at(5)->values().norm(), frefNorm), data.at(5)->values().norm());
745
746 // update cplData
747 fpcd->storeIteration();
748
749 fpcd->emplaceSampleAtTime(windowEnd, {-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});
750 dpcd->emplaceSampleAtTime(windowEnd, {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});
751
752 } else if (context.isRank(3)) { // Secondary rank 3
753 // Dummy Secondary rank to be able to reuse the 4 proc Intra-participant communication Fixture
754 }
755
756 // QN- Update, 2. iteration
757 pp.performAcceleration(data, windowStart, windowEnd);
758
759 // 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()
760 for (auto &pair : data) {
761 pair.second->setSampleAtTime(windowEnd, pair.second->sample());
762 }
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(windowEnd, 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
796 fpcd->emplaceSampleAtTime(windowEnd, {-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 dpcd->emplaceSampleAtTime(windowEnd, {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});
798
799 } else if (context.isRank(2)) { // Secondary rank 2
800
801 dref = Eigen::VectorXd::Zero(22);
802 fref = Eigen::VectorXd::Zero(22);
803
804 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;
805 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;
806 drefNorm = 6.131610103923933e-07;
807 frefNorm = 0.1265754339635572;
808
809 // validate values
810 for (int i = 0; i < data.at(4)->values().size(); i++)
811 BOOST_TEST(testing::equals(data.at(4)->values()(i), dref(i)));
812
813 for (int i = 0; i < data.at(5)->values().size(); i++)
814 BOOST_TEST(testing::equals(data.at(5)->values()(i), fref(i)));
815
816 // validate norm
817 BOOST_TEST(testing::equals(data.at(4)->values().norm(), drefNorm), data.at(4)->values().norm());
818 BOOST_TEST(testing::equals(data.at(5)->values().norm(), frefNorm), data.at(5)->values().norm());
819
820 // update cplData
821 dpcd->storeIteration();
822 fpcd->storeIteration();
823
824 fpcd->emplaceSampleAtTime(windowEnd, {-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});
825 dpcd->emplaceSampleAtTime(windowEnd, {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});
826
827 } else if (context.isRank(3)) { // Secondary rank 3
828 // Dummy Secondary rank to be able to reuse the 4 proc Intra-participant communication Fixture
829 }
830
831 // QN- Update, 3. iteration
832 pp.performAcceleration(data, windowStart, windowEnd);
833
834 // 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()
835 for (auto &pair : data) {
836 pair.second->setSampleAtTime(windowEnd, pair.second->sample());
837 }
838
839 if (context.isPrimary()) { // Primary
840
841 } else if (context.isRank(1)) { // SecondaryRank1
842
843 dref = Eigen::VectorXd::Zero(22);
844 fref = Eigen::VectorXd::Zero(22);
845
846 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;
847 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;
848 drefNorm = 5.062970166651817e-05;
849 frefNorm = 0.1247902554601672;
850
851 // validate values
852 for (int i = 0; i < data.at(4)->values().size(); i++)
853 BOOST_TEST(testing::equals(data.at(4)->values()(i), dref(i)));
854
855 for (int i = 0; i < data.at(5)->values().size(); i++)
856 BOOST_TEST(testing::equals(data.at(5)->values()(i), fref(i)));
857
858 // validate norm
859 BOOST_TEST(testing::equals(data.at(4)->values().norm(), drefNorm), data.at(4)->values().norm());
860 BOOST_TEST(testing::equals(data.at(5)->values().norm(), frefNorm), data.at(5)->values().norm());
861
862 // update cplData
863 dpcd->storeIteration();
864 fpcd->storeIteration();
865
866 fpcd->emplaceSampleAtTime(windowEnd, {-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});
867 dpcd->emplaceSampleAtTime(windowEnd, {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});
868
869 } else if (context.isRank(2)) { // Secondary rank 2
870
871 dref = Eigen::VectorXd::Zero(22);
872 fref = Eigen::VectorXd::Zero(22);
873
874 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;
875 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;
876 drefNorm = 4.824205753272403e-05;
877 frefNorm = 0.1225851627302816;
878
879 // validate values
880 for (int i = 0; i < data.at(4)->values().size(); i++)
881 BOOST_TEST(testing::equals(data.at(4)->values()(i), dref(i)));
882
883 for (int i = 0; i < data.at(5)->values().size(); i++)
884 BOOST_TEST(testing::equals(data.at(5)->values()(i), fref(i)));
885
886 // validate norm
887 BOOST_TEST(testing::equals(data.at(4)->values().norm(), drefNorm), data.at(4)->values().norm());
888 BOOST_TEST(testing::equals(data.at(5)->values().norm(), frefNorm), data.at(5)->values().norm());
889
890 // update cplData
891 dpcd->storeIteration();
892 fpcd->storeIteration();
893
894 fpcd->emplaceSampleAtTime(windowEnd, {-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});
895 dpcd->emplaceSampleAtTime(windowEnd, {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});
896
897 } else if (context.isRank(3)) { // Secondary rank 3
898 // Dummy Secondary rank to be able to reuse the 4 proc Intra-participant communication Fixture
899 }
900
901 // QN- Update, 4. iteration
902 pp.performAcceleration(data, windowStart, windowEnd);
903
904 // 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()
905 for (auto &pair : data) {
906 pair.second->setSampleAtTime(windowEnd, pair.second->sample());
907 }
908
909 if (context.isPrimary()) { // Primary
910
911 } else if (context.isRank(1)) { // SecondaryRank1
912
913 dref = Eigen::VectorXd::Zero(22);
914 fref = Eigen::VectorXd::Zero(22);
915
916 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;
917 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;
918 drefNorm = 4.815113092042934e-05;
919 frefNorm = 0.1204042970796453;
920
921 // validate values
922 for (int i = 0; i < data.at(4)->values().size(); i++)
923 BOOST_TEST(testing::equals(data.at(4)->values()(i), dref(i)));
924
925 for (int i = 0; i < data.at(5)->values().size(); i++)
926 BOOST_TEST(testing::equals(data.at(5)->values()(i), fref(i)));
927
928 // validate norm
929 BOOST_TEST(testing::equals(data.at(4)->values().norm(), drefNorm));
930 BOOST_TEST(testing::equals(data.at(5)->values().norm(), frefNorm));
931
932 // update cplData
933 dpcd->storeIteration();
934 fpcd->storeIteration();
935
936 fpcd->emplaceSampleAtTime(windowEnd, {-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});
937 dpcd->emplaceSampleAtTime(windowEnd, {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});
938
939 } else if (context.isRank(2)) { // Secondary rank 2
940
941 dref = Eigen::VectorXd::Zero(22);
942 fref = Eigen::VectorXd::Zero(22);
943
944 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;
945 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;
946 drefNorm = 4.587992970636867e-05;
947 frefNorm = 0.1180354804938519;
948
949 // validate values
950 for (int i = 0; i < data.at(4)->values().size(); i++)
951 BOOST_TEST(testing::equals(data.at(4)->values()(i), dref(i)));
952
953 for (int i = 0; i < data.at(5)->values().size(); i++)
954 BOOST_TEST(testing::equals(data.at(5)->values()(i), fref(i)));
955
956 // validate norm
957 BOOST_TEST(testing::equals(data.at(4)->values().norm(), drefNorm));
958 BOOST_TEST(testing::equals(data.at(5)->values().norm(), frefNorm));
959
960 // update cplData
961 dpcd->storeIteration();
962 fpcd->storeIteration();
963
964 fpcd->emplaceSampleAtTime(windowEnd, {-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});
965 dpcd->emplaceSampleAtTime(windowEnd, {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});
966
967 } else if (context.isRank(3)) { // Secondary rank 3
968 // Dummy Secondary rank to be able to reuse the 4 proc Intra-participant communication Fixture
969 }
970
971 // QN- Update, 5. iteratidataon
972 pp.performAcceleration(data, windowStart, windowEnd);
973
974 // 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()
975 for (auto &pair : data) {
976 pair.second->setSampleAtTime(windowEnd, pair.second->sample());
977 }
978
979 if (context.isPrimary()) { // Primary
980
981 } else if (context.isRank(1)) { // SecondaryRank1
982
983 dref = Eigen::VectorXd::Zero(22);
984 fref = Eigen::VectorXd::Zero(22);
985
986 dref << 1.275776729912441e-06, -7.411719120649928e-07, 2.009844043916140e-06, -8.166362562036991e-07,
987 1.98429549697312266297e-06, -1.006128466370864e-06, 1.26860076890459333335e-06, -1.038973060940335e-06,
988 1.553926578960109e-06, -1.85501566670221e-06, 2.21293251735091021427e-06, -1.645357854619679e-06,
989 1.820890655907848e-06, -2.41565468511731e-06, 2.44472211085250353034e-06, -2.153167028992053e-06,
990 1.367377699262639e-06, -1.402360838056698e-06, 2.05877363493188315457e-06, -1.275601072931189e-06,
991 2.16056223288899190422e-06, -6.758668464219906e-07;
992
993 fref << -0.02494062387644205, 2.83457783442084623737e-05, 2.78926996928726099456e-02, 4.927873553950413e-05,
994 0.02806450388830189, 4.133693821366594e-05, -0.02479702682061509, 7.828236584438153e-05,
995 -0.02470268458079963, 3.17420733502629e-05, 0.0282235973672153, 0.0001624809599050109,
996 -0.01213681562464092, -0.01299265397817983, 0.01440494173854578, -0.01303642911469127,
997 -0.02467306876921166, 4.955128897708022e-05, 0.02820821342518329, 9.359794339403125e-05,
998 0.02784148826297314, 7.700134509804057e-05;
999
1000 drefNorm = 7.910283453653413e-06;
1001 frefNorm = 0.08415800797163485;
1002
1003 // validate values
1004 for (int i = 0; i < data.at(4)->values().size(); i++)
1005 BOOST_TEST(testing::equals(data.at(4)->values()(i), dref(i)));
1006
1007 for (int i = 0; i < data.at(5)->values().size(); i++)
1008 BOOST_TEST(testing::equals(data.at(5)->values()(i), fref(i), 1e-8));
1009
1010 // validate norm
1011 BOOST_TEST(testing::equals(data.at(4)->values().norm(), drefNorm));
1012 BOOST_TEST(testing::equals(data.at(5)->values().norm(), frefNorm));
1013
1014 } else if (context.isRank(2)) { // Secondary rank 2
1015
1016 dref = Eigen::VectorXd::Zero(22);
1017 fref = Eigen::VectorXd::Zero(22);
1018
1019 dref << 1.782695896956317e-06, -3.609665456581163e-07, 2.549165865006206e-06, -2.9152246796606e-07,
1020 1.641256623136569e-06, -3.63485032734166e-07, 2.399541139942214e-06, -5.221960935703056e-07,
1021 0, 0, 0, 0,
1022 1.540562829986484e-06, 2.646900958343007e-07, 6.572995607065404e-07, -5.789378440843989e-07,
1023 2.312916074266845e-06, 2.909615579313369e-08, 1.55690286980020939237e-06, -4.907962897313283e-07,
1024 1.41306136610511245499e-06, -5.042593080795539e-07;
1025
1026 fref << -0.02407925914468757, -0.0001962517462601011, 2.86388026412357221684e-02, 0.0002189694547629064,
1027 -0.02469057752137756, -0.0001420431374215673, 0.02806138636325711, 0.0001618659175134536,
1028 -0.01072215615486178, -0.01066770812537758, 0.01563815124261267, -0.01570783185165062,
1029 0.03075510899502418, 0.0001765574924817709, -0.02194129299700746, -0.0001523229943385313,
1030 0.02962481369149697, 2.35608194343154197722e-04, -0.02307508702308544, -0.0002109167490697018,
1031 -0.02495036462961023, -1.65224214831653608282e-05;
1032
1033 drefNorm = 5.676684399367158e-06;
1034 frefNorm = 0.08353026170200345;
1035
1036 // validate values
1037 for (int i = 0; i < data.at(4)->values().size(); i++)
1038 BOOST_TEST(testing::equals(data.at(4)->values()(i), dref(i)));
1039
1040 for (int i = 0; i < data.at(5)->values().size(); i++)
1041 BOOST_TEST(testing::equals(data.at(5)->values()(i), fref(i)));
1042
1043 // validate norm
1044 BOOST_TEST(testing::equals(data.at(4)->values().norm(), drefNorm));
1045 BOOST_TEST(testing::equals(data.at(5)->values().norm(), frefNorm));
1046 } else if (context.isRank(3)) { // Secondary rank 3
1047 // Dummy Secondary rank to be able to reuse the 4 proc Intra-participant communication Fixture
1048 }
1049}
1050
1052PRECICE_TEST_SETUP(""_on(4_ranks).setupIntraComm())
1053BOOST_AUTO_TEST_CASE(testColumnsLoggingWithoutSubsteps)
1054{
1055 PRECICE_TEST();
1056 double initialRelaxation = 0.1;
1057 int maxIterationsUsed = 3;
1058 int timeWindowsReused = 1;
1059 int filter = BaseQNAcceleration::QR1FILTER;
1060 double singularityLimit = 0.1;
1061 bool enforceInitialRelaxation = false;
1062 bool exchangeSubsteps = false; // @todo add testColumnsLoggingWithSubsteps, where exchangeSubsteps = true as soon as acceleration scheme supports subcycling.
1063 std::vector<int> dataIDs;
1064 dataIDs.push_back(0);
1065 std::vector<double> factors;
1066 factors.resize(1.0, 1.0);
1067 PtrPreconditioner prec(new ConstantPreconditioner(factors));
1068 std::vector<int> vertexOffsets{2, 3, 3, 4};
1069
1070 mesh::PtrMesh dummyMesh(new mesh::Mesh("DummyMesh", 3, testing::nextMeshID()));
1071 dummyMesh->setVertexOffsets(vertexOffsets);
1072 if (context.isPrimary()) { // 2 vertices
1073 dummyMesh->createVertex(Eigen::Vector3d{0, 0, 0});
1074 dummyMesh->createVertex(Eigen::Vector3d{1, 0, 0});
1075 } else if (context.isRank(1)) { // 1 vertex
1076 dummyMesh->createVertex(Eigen::Vector3d{0, 0, 0});
1077 } else if (context.isRank(2)) { // no vertices
1078 } else { // 1 vertex
1079 dummyMesh->createVertex(Eigen::Vector3d{0, 0, 0});
1080 }
1081
1082 IQNILSAcceleration acc(initialRelaxation, enforceInitialRelaxation, maxIterationsUsed,
1083 timeWindowsReused, filter, singularityLimit, dataIDs, prec, !exchangeSubsteps);
1084
1085 mesh::PtrData displacements(new mesh::Data("dvalues", -1, 1));
1086
1087 DataMap data;
1088
1089 const double windowStart = 0;
1090 const double windowEnd = 1;
1091
1092 Eigen::VectorXd insert;
1093 // init displacements
1094 if (context.isPrimary()) { // 2 vertices
1095 insert.resize(2);
1096 insert << 0.5, 0.5;
1097 } else if (context.isRank(1)) { // 1 vertex
1098 insert.resize(1);
1099 insert << 0.5;
1100 } else if (context.isRank(2)) { // no vertices
1101 } else { // 1 vertex
1102 insert.resize(1);
1103 insert << 0.5;
1104 }
1105
1106 utils::append(displacements->values(), insert);
1107 PtrCouplingData dpcd = makeCouplingData(displacements, dummyMesh, exchangeSubsteps);
1108 data.insert(std::pair<int, PtrCouplingData>(0, dpcd));
1109
1110 // Need at least two samples to get a valid waveform
1111 dpcd->setSampleAtTime(windowStart, dpcd->sample());
1112 dpcd->setSampleAtTime(windowEnd, dpcd->sample());
1113 dpcd->storeIteration();
1114
1115 acc.initialize(data);
1116
1117 // update displacements
1118 if (context.isPrimary()) { // 2 vertices
1119 insert << 1.0, 1.0;
1120 } else if (context.isRank(1)) { // 1 vertex
1121 insert << 1.0;
1122 } else if (context.isRank(2)) { // no vertices
1123 } else { // 1 vertex
1124 insert << 1.0;
1125 }
1126
1127 displacements->values() = insert;
1128
1129 // Need to update the samples
1130 dpcd->setSampleAtTime(windowStart, dpcd->sample());
1131 dpcd->setSampleAtTime(windowEnd, dpcd->sample());
1132
1133 acc.performAcceleration(data, windowStart, windowEnd);
1134
1135 Eigen::VectorXd newdvalues1;
1136 if (context.isPrimary()) {
1137 utils::append(newdvalues1, 1.1);
1138 utils::append(newdvalues1, 1.0);
1139 } else if (context.isRank(1)) {
1140 utils::append(newdvalues1, 1.0);
1141 } else if (context.isRank(2)) {
1142 } else if (context.isRank(3)) {
1143 utils::append(newdvalues1, 1.0);
1144 }
1145
1146 dpcd->setSampleAtTime(windowEnd, time::Sample(dpcd->getDimensions(), newdvalues1));
1147
1148 acc.performAcceleration(data, windowStart, windowEnd);
1149
1150 Eigen::VectorXd newdvalues2;
1151 if (context.isPrimary()) {
1152 utils::append(newdvalues2, 1.0);
1153 utils::append(newdvalues2, 2.0);
1154 } else if (context.isRank(1)) {
1155 utils::append(newdvalues2, 1.0);
1156 } else if (context.isRank(2)) {
1157 } else if (context.isRank(3)) {
1158 utils::append(newdvalues2, 1.0);
1159 }
1160
1161 dpcd->setSampleAtTime(windowEnd, time::Sample(dpcd->getDimensions(), newdvalues2));
1162
1163 acc.iterationsConverged(data, windowStart);
1164
1165 BOOST_TEST(acc.getLSSystemCols() == 2);
1166 BOOST_TEST(acc.getDeletedColumns() == 0);
1167 BOOST_TEST(acc.getDroppedColumns() == 0);
1168
1169 Eigen::VectorXd newdvalues3;
1170 if (context.isPrimary()) {
1171 utils::append(newdvalues3, 1.1);
1172 utils::append(newdvalues3, 2.0);
1173 } else if (context.isRank(1)) {
1174 utils::append(newdvalues3, 1.0);
1175 } else if (context.isRank(2)) {
1176 } else if (context.isRank(3)) {
1177 utils::append(newdvalues3, 1.0);
1178 }
1179
1180 dpcd->setSampleAtTime(windowEnd, time::Sample(dpcd->getDimensions(), newdvalues3));
1181
1182 acc.performAcceleration(data, windowStart, windowEnd);
1183
1184 Eigen::VectorXd newdvalues4;
1185 if (context.isPrimary()) {
1186 utils::append(newdvalues4, 1.0);
1187 utils::append(newdvalues4, 1.0);
1188 } else if (context.isRank(1)) {
1189 utils::append(newdvalues4, 1.0);
1190 } else if (context.isRank(2)) {
1191 } else if (context.isRank(3)) {
1192 utils::append(newdvalues4, 5.0);
1193 }
1194
1195 dpcd->setSampleAtTime(windowEnd, time::Sample(dpcd->getDimensions(), newdvalues4));
1196
1197 acc.iterationsConverged(data, windowStart);
1198
1199 BOOST_TEST(acc.getLSSystemCols() == 1);
1200 BOOST_TEST(acc.getDeletedColumns() == 1);
1201 BOOST_TEST(acc.getDroppedColumns() == 1);
1202}
1203
1204BOOST_AUTO_TEST_SUITE_END()
1205BOOST_AUTO_TEST_SUITE_END()
1206
1207#endif // PRECICE_NO_MPI
makeCouplingData
cplscheme::PtrCouplingData makeCouplingData(mesh::PtrData data, mesh::PtrMesh mesh, bool exchangeSubsteps)
Definition helper.hpp:10
DataMap
std::map< int, PtrCouplingData > DataMap
Definition AccelerationIntraCommTest.cpp:40
BOOST_AUTO_TEST_CASE
BOOST_AUTO_TEST_CASE(testIQNIMVJPPWithSubsteps)
Definition AccelerationSerialTest.cpp:122
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:39
PRECICE_TEST_SETUP
#define PRECICE_TEST_SETUP(...)
Creates and attaches a TestSetup to a Boost test case.
Definition Testing.hpp:29
precice::acceleration::BaseQNAcceleration
Base Class for quasi-Newton acceleration schemes.
Definition BaseQNAcceleration.hpp:64
precice::acceleration::IQNILSAcceleration
Interface quasi-Newton with interface least-squares approximation.
Definition IQNILSAcceleration.hpp:25
precice::acceleration::IQNIMVJAcceleration
Multi vector quasi-Newton update scheme.
Definition IQNIMVJAcceleration.hpp:32
precice::acceleration::impl::ConstantPreconditioner
Preconditioner that uses the constant user-defined factors to scale the quasi-Newton system.
Definition ConstantPreconditioner.hpp:13
precice::acceleration::impl::ResidualSumPreconditioner
Preconditioner that uses the residuals of all iterations of the current time window summed up to scal...
Definition ResidualSumPreconditioner.hpp:16
precice::cplscheme::CouplingData::values
const Eigen::VectorXd & values() const
Returns a const reference to the data values.
Definition CouplingData.cpp:47
precice::cplscheme::CouplingData::timeStepsStorage
time::Storage & timeStepsStorage()
Returns a reference to the time step storage of the data.
Definition CouplingData.cpp:77
precice::cplscheme::CouplingData::emplaceSampleAtTime
void emplaceSampleAtTime(double time)
Creates an empty sample at given time.
Definition CouplingData.cpp:124
precice::cplscheme::CouplingData::setSampleAtTime
void setSampleAtTime(double time, time::Sample sample)
Add sample at given time to _timeStepsStorage.
Definition CouplingData.cpp:99
precice::cplscheme::CouplingData::sample
const time::Sample & sample() const
Returns a const reference to the data Sample.
Definition CouplingData.cpp:71
precice::cplscheme::CouplingData::storeIteration
void storeIteration()
store _data->values() in read-only variable _previousIteration for convergence checks etc.
Definition CouplingData.cpp:165
precice::cplscheme::CouplingData::getPreviousValuesAtTime
time::SampleResult getPreviousValuesAtTime(double relativeDt)
returns previous data interpolated to the relativeDt time
Definition CouplingData.cpp:89
precice::mesh::Data
Describes a set of data values belonging to the vertices of a mesh.
Definition Data.hpp:27
precice::mesh::Mesh
Container and creator for meshes.
Definition Mesh.hpp:38
precice::time::Storage::sample
SampleResult sample(double time) const
Need to use interpolation for the case with changing time grids.
Definition Storage.cpp:198
vector
STL class.
std::data
T data(T... args)
precice::acceleration::impl::PtrPreconditioner
std::shared_ptr< Preconditioner > PtrPreconditioner
Definition SharedPointer.hpp:10
precice::acceleration
contains implementations of acceleration schemes.
Definition Acceleration.cpp:5
precice::cplscheme
contains implementations of coupling schemes for coupled simulations.
Definition BaseCouplingScheme.cpp:33
precice::cplscheme::PtrCouplingData
std::shared_ptr< CouplingData > PtrCouplingData
Definition SharedPointer.hpp:14
precice::mesh::PtrData
std::shared_ptr< Data > PtrData
Definition SharedPointer.hpp:13
precice::mesh::PtrMesh
std::shared_ptr< Mesh > PtrMesh
Definition SharedPointer.hpp:15
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:93
precice::testing::makeCouplingData
cplscheme::PtrCouplingData makeCouplingData(mesh::PtrData data, mesh::PtrMesh mesh, bool exchangeSubsteps)
Definition helper.hpp:10
precice::testing::addDummyVertices
void addDummyVertices(size_t nVertices, mesh::Mesh &mesh)
Definition Meshes.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)
std::size
T size(T... args)