@@ -68,7 +68,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
68
68
Eigen::Vector3d sphereCenter = Eigen::Vector3d::Zero ();
69
69
SphericalDiffuse<> gf (eps1, eps2, width, sphereRadius, sphereCenter, maxL);
70
70
THEN (" the value of the Green's function inside the droplet is" ) {
71
- double value = 0.0204749331147813275 ;
71
+ double value = 0.0204749331147992819 ;
72
72
double gf_value = gf.kernelS (source1, probe1);
73
73
INFO (" ref_value = " << std::setprecision (
74
74
std::numeric_limits<long double >::digits10)
@@ -79,7 +79,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
79
79
REQUIRE (value == Approx (gf_value));
80
80
}
81
81
AND_THEN (" the value of the Green's function outside the droplet is" ) {
82
- double value = 0.488060362514614432 ;
82
+ double value = 0.487473605103865837 ;
83
83
double gf_value = gf.kernelS (source2, probe2);
84
84
INFO (" ref_value = " << std::setprecision (
85
85
std::numeric_limits<long double >::digits10)
@@ -91,7 +91,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
91
91
}
92
92
THEN (" the value of the Green's function directional derivative wrt the probe "
93
93
" point inside the droplet is" ) {
94
- double derProbe = -0.0124999765635167015 ;
94
+ double derProbe = -0.0124999765637075211 ;
95
95
double gf_derProbe = gf.derivativeProbe (probeNormal1, source1, probe1);
96
96
INFO (" ref_derProbe = "
97
97
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -103,7 +103,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
103
103
}
104
104
AND_THEN (" the value of the Green's function directional derivative wrt the "
105
105
" probe point outside the droplet is" ) {
106
- double derProbe = -0.290662079075743041 ;
106
+ double derProbe = -0.290662955233778053 ;
107
107
double gf_derProbe = gf.derivativeProbe (probeNormal2, source2, probe2);
108
108
INFO (" ref_derProbe = "
109
109
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -115,7 +115,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
115
115
}
116
116
THEN (" the value of the Green's function directional derivative wrt the source "
117
117
" point inside the droplet is" ) {
118
- double derSource = 0.00937504723103055326 ;
118
+ double derSource = 0.00937504723101320603 ;
119
119
double gf_derSource = gf.derivativeSource (sourceNormal1, source1, probe1);
120
120
INFO (" ref_derSource = "
121
121
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -127,7 +127,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
127
127
}
128
128
AND_THEN (" the value of the Green's function directional derivative wrt the "
129
129
" source point outside the droplet is" ) {
130
- double derSource = -0.938830970128590181 ;
130
+ double derSource = -0.901954472843635724 ;
131
131
double gf_derSource = gf.derivativeSource (sourceNormal2, source2, probe2);
132
132
INFO (" ref_derSource = "
133
133
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -144,7 +144,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
144
144
(Eigen::Vector3d () << 25.0 , 0.0 , 0.0 ).finished ();
145
145
SphericalDiffuse<> gf (eps1, eps2, width, sphereRadius, sphereCenter, maxL);
146
146
THEN (" the value of the Green's function inside the droplet is" ) {
147
- double value = 0.0173845118865096904 ;
147
+ double value = 0.0173845118865097424 ;
148
148
double gf_value = gf.kernelS (source1, probe1);
149
149
INFO (" ref_value = " << std::setprecision (
150
150
std::numeric_limits<long double >::digits10)
@@ -155,7 +155,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
155
155
REQUIRE (value == Approx (gf_value));
156
156
}
157
157
AND_THEN (" the value of the Green's function outside the droplet is" ) {
158
- double value = 0.501844654162973303 ;
158
+ double value = 0.501385603574595162 ;
159
159
double gf_value = gf.kernelS (source2, probe2);
160
160
INFO (" ref_value = " << std::setprecision (
161
161
std::numeric_limits<long double >::digits10)
@@ -167,7 +167,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
167
167
}
168
168
THEN (" the value of the Green's function directional derivative wrt the probe "
169
169
" point inside the droplet is" ) {
170
- double derProbe = -0.012483438527558649 ;
170
+ double derProbe = -0.0124834385275933435 ;
171
171
double gf_derProbe = gf.derivativeProbe (probeNormal1, source1, probe1);
172
172
INFO (" ref_derProbe = "
173
173
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -179,7 +179,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
179
179
}
180
180
AND_THEN (" the value of the Green's function directional derivative wrt the "
181
181
" probe point outside the droplet is" ) {
182
- double derProbe = -0.291257922496734878 ;
182
+ double derProbe = -0.291262383059764929 ;
183
183
double gf_derProbe = gf.derivativeProbe (probeNormal2, source2, probe2);
184
184
INFO (" ref_derProbe = "
185
185
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -191,7 +191,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
191
191
}
192
192
THEN (" the value of the Green's function directional derivative wrt the source "
193
193
" point inside the droplet is" ) {
194
- double derSource = 0.0124835646475064677 ;
194
+ double derSource = 0.012483564647454426 ;
195
195
double gf_derSource = gf.derivativeSource (sourceNormal1, source1, probe1);
196
196
INFO (" ref_derSource = "
197
197
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -203,7 +203,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
203
203
}
204
204
AND_THEN (" the value of the Green's function directional derivative wrt the "
205
205
" source point outside the droplet is" ) {
206
- double derSource = 5.01504540067754245 ;
206
+ double derSource = 2.720604738008503 ;
207
207
double gf_derSource = gf.derivativeSource (sourceNormal2, source2, probe2);
208
208
INFO (" ref_derSource = "
209
209
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -221,7 +221,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
221
221
SphericalDiffuse<OneLayerTanh> gf (
222
222
eps1, eps2, width, sphereRadius, sphereCenter, maxL);
223
223
THEN (" the value of the Green's function inside the droplet is" ) {
224
- double value = 0.0204265162808782985 ;
224
+ double value = 0.0204265162808963327 ;
225
225
double gf_value = gf.kernelS (source1, probe1);
226
226
INFO (" ref_value = " << std::setprecision (
227
227
std::numeric_limits<long double >::digits10)
@@ -232,7 +232,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
232
232
REQUIRE (value == Approx (gf_value));
233
233
}
234
234
AND_THEN (" the value of the Green's function outside the droplet is" ) {
235
- double value = 0.428395676083591359 ;
235
+ double value = 0.42831756921770181 ;
236
236
double gf_value = gf.kernelS (source2, probe2);
237
237
INFO (" ref_value = " << std::setprecision (
238
238
std::numeric_limits<long double >::digits10)
@@ -244,7 +244,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
244
244
}
245
245
THEN (" the value of the Green's function directional derivative wrt the probe "
246
246
" point inside the droplet is" ) {
247
- double derProbe = -0.0124999769449823939 ;
247
+ double derProbe = -0.0124999769448609632 ;
248
248
double gf_derProbe = gf.derivativeProbe (probeNormal1, source1, probe1);
249
249
INFO (" ref_derProbe = "
250
250
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -256,7 +256,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
256
256
}
257
257
AND_THEN (" the value of the Green's function directional derivative wrt the "
258
258
" probe point outside the droplet is" ) {
259
- double derProbe = -0.288385047571004804 ;
259
+ double derProbe = -0.288392889131738883 ;
260
260
double gf_derProbe = gf.derivativeProbe (probeNormal2, source2, probe2);
261
261
INFO (" ref_derProbe = "
262
262
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -280,7 +280,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
280
280
}
281
281
AND_THEN (" the value of the Green's function directional derivative wrt the "
282
282
" source point outside the droplet is" ) {
283
- double derSource = 3.22957321418654297 ;
283
+ double derSource = - 1.51576208806569745 ;
284
284
double gf_derSource = gf.derivativeSource (sourceNormal2, source2, probe2);
285
285
INFO (" ref_derSource = "
286
286
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -298,7 +298,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
298
298
SphericalDiffuse<OneLayerTanh> gf (
299
299
eps1, eps2, width, sphereRadius, sphereCenter, maxL);
300
300
THEN (" the value of the Green's function inside the droplet is" ) {
301
- double value = 0.0173358023608628994 ;
301
+ double value = 0.017335802360862948 ;
302
302
double gf_value = gf.kernelS (source1, probe1);
303
303
INFO (" ref_value = " << std::setprecision (
304
304
std::numeric_limits<long double >::digits10)
@@ -309,7 +309,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
309
309
REQUIRE (value == Approx (gf_value));
310
310
}
311
311
AND_THEN (" the value of the Green's function outside the droplet is" ) {
312
- double value = 0.471193749836353093 ;
312
+ double value = 0.470725113611886292 ;
313
313
double gf_value = gf.kernelS (source2, probe2);
314
314
INFO (" ref_value = " << std::setprecision (
315
315
std::numeric_limits<long double >::digits10)
@@ -321,7 +321,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
321
321
}
322
322
THEN (" the value of the Green's function directional derivative wrt the probe "
323
323
" point inside the droplet is" ) {
324
- double derProbe = -0.0124834504454732209 ;
324
+ double derProbe = -0.0124834504460630269 ;
325
325
double gf_derProbe = gf.derivativeProbe (probeNormal1, source1, probe1);
326
326
INFO (" ref_derProbe = "
327
327
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -333,7 +333,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
333
333
}
334
334
AND_THEN (" the value of the Green's function directional derivative wrt the "
335
335
" probe point outside the droplet is" ) {
336
- double derProbe = -0.290867341330436346 ;
336
+ double derProbe = -0.290870450536662162 ;
337
337
double gf_derProbe = gf.derivativeProbe (probeNormal2, source2, probe2);
338
338
INFO (" ref_derProbe = "
339
339
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -345,7 +345,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
345
345
}
346
346
THEN (" the value of the Green's function directional derivative wrt the source "
347
347
" point inside the droplet is" ) {
348
- double derSource = 0.0124835522706361057 ;
348
+ double derSource = 0.012483552270584064 ;
349
349
double gf_derSource = gf.derivativeSource (sourceNormal1, source1, probe1);
350
350
INFO (" ref_derSource = "
351
351
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -357,7 +357,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
357
357
}
358
358
AND_THEN (" the value of the Green's function directional derivative wrt the "
359
359
" source point outside the droplet is" ) {
360
- double derSource = 4.9830666203676266 ;
360
+ double derSource = 3.75908425086779463 ;
361
361
double gf_derSource = gf.derivativeSource (sourceNormal2, source2, probe2);
362
362
INFO (" ref_derSource = "
363
363
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -376,7 +376,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
376
376
SphericalDiffuse<OneLayerErf> gf (
377
377
eps1, eps2, width, sphereRadius, sphereCenter, maxL);
378
378
THEN (" the value of the Green's function inside the droplet is" ) {
379
- double value = 0.0204465601345453427 ;
379
+ double value = 0.0204465601345278428 ;
380
380
double gf_value = gf.kernelS (source1, probe1);
381
381
INFO (" ref_value = " << std::setprecision (
382
382
std::numeric_limits<long double >::digits10)
@@ -387,7 +387,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
387
387
REQUIRE (value == Approx (gf_value));
388
388
}
389
389
AND_THEN (" the value of the Green's function outside the droplet is" ) {
390
- double value = 0.546459075558889285 ;
390
+ double value = 0.546440841068159378 ;
391
391
double gf_value = gf.kernelS (source2, probe2);
392
392
INFO (" ref_value = " << std::setprecision (
393
393
std::numeric_limits<long double >::digits10)
@@ -399,7 +399,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
399
399
}
400
400
THEN (" the value of the Green's function directional derivative wrt the probe "
401
401
" point inside the droplet is" ) {
402
- double derProbe = -0.0124999769317464537 ;
402
+ double derProbe = -0.0124999769314168563 ;
403
403
double gf_derProbe = gf.derivativeProbe (probeNormal1, source1, probe1);
404
404
INFO (" ref_derProbe = "
405
405
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -411,7 +411,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
411
411
}
412
412
AND_THEN (" the value of the Green's function directional derivative wrt the "
413
413
" probe point outside the droplet is" ) {
414
- double derProbe = -0.292047321249766512 ;
414
+ double derProbe = -0.292045737975143993 ;
415
415
double gf_derProbe = gf.derivativeProbe (probeNormal2, source2, probe2);
416
416
INFO (" ref_derProbe = "
417
417
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -423,7 +423,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
423
423
}
424
424
THEN (" the value of the Green's function directional derivative wrt the source "
425
425
" point inside the droplet is" ) {
426
- double derSource = 0.00937504649521289646 ;
426
+ double derSource = 0.00937504649523024369 ;
427
427
double gf_derSource = gf.derivativeSource (sourceNormal1, source1, probe1);
428
428
INFO (" ref_derSource = "
429
429
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -435,7 +435,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
435
435
}
436
436
AND_THEN (" the value of the Green's function directional derivative wrt the "
437
437
" source point outside the droplet is" ) {
438
- double derSource = -2.1857233870498094 ;
438
+ double derSource = -1.60173403890240262 ;
439
439
double gf_derSource = gf.derivativeSource (sourceNormal2, source2, probe2);
440
440
INFO (" ref_derSource = "
441
441
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -453,7 +453,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
453
453
SphericalDiffuse<OneLayerErf> gf (
454
454
eps1, eps2, width, sphereRadius, sphereCenter, maxL);
455
455
THEN (" the value of the Green's function inside the droplet is" ) {
456
- double value = 0.0173558529449408631 ;
456
+ double value = 0.0173558563009676217 ;
457
457
double gf_value = gf.kernelS (source1, probe1);
458
458
INFO (" ref_value = " << std::setprecision (
459
459
std::numeric_limits<long double >::digits10)
@@ -464,7 +464,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
464
464
REQUIRE (value == Approx (gf_value));
465
465
}
466
466
AND_THEN (" the value of the Green's function outside the droplet is" ) {
467
- double value = 0.372262125530776311 ;
467
+ double value = 0.37238664692877399 ;
468
468
double gf_value = gf.kernelS (source2, probe2);
469
469
INFO (" ref_value = " << std::setprecision (
470
470
std::numeric_limits<long double >::digits10)
@@ -476,7 +476,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
476
476
}
477
477
THEN (" the value of the Green's function directional derivative wrt the probe "
478
478
" point inside the droplet is" ) {
479
- double derProbe = -0.0124834509424368023 ;
479
+ double derProbe = -0.0124834500456541542 ;
480
480
double gf_derProbe = gf.derivativeProbe (probeNormal1, source1, probe1);
481
481
INFO (" ref_derProbe = "
482
482
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -488,7 +488,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
488
488
}
489
489
AND_THEN (" the value of the Green's function directional derivative wrt the "
490
490
" probe point outside the droplet is" ) {
491
- double derProbe = -0.27916006373779334 ;
491
+ double derProbe = -0.279220523679291066 ;
492
492
double gf_derProbe = gf.derivativeProbe (probeNormal2, source2, probe2);
493
493
INFO (" ref_derProbe = "
494
494
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -500,7 +500,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
500
500
}
501
501
THEN (" the value of the Green's function directional derivative wrt the source "
502
502
" point inside the droplet is" ) {
503
- double derSource = 0.0124835526376168571 ;
503
+ double derSource = 0.0124835526869870872 ;
504
504
double gf_derSource = gf.derivativeSource (sourceNormal1, source1, probe1);
505
505
INFO (" ref_derSource = "
506
506
<< std::setprecision (std::numeric_limits<long double >::digits10)
@@ -512,7 +512,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
512
512
}
513
513
AND_THEN (" the value of the Green's function directional derivative wrt the "
514
514
" source point outside the droplet is" ) {
515
- double derSource = 3.4616694327527231 ;
515
+ double derSource = 4.79631817304426722 ;
516
516
double gf_derSource = gf.derivativeSource (sourceNormal2, source2, probe2);
517
517
INFO (" ref_derSource = "
518
518
<< std::setprecision (std::numeric_limits<long double >::digits10)
0 commit comments