@@ -68,7 +68,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
6868 Eigen::Vector3d sphereCenter = Eigen::Vector3d::Zero ();
6969 SphericalDiffuse<> gf (eps1, eps2, width, sphereRadius, sphereCenter, maxL);
7070 THEN (" the value of the Green's function inside the droplet is" ) {
71- double value = 0.0204749331147813275 ;
71+ double value = 0.0204749331147992819 ;
7272 double gf_value = gf.kernelS (source1, probe1);
7373 INFO (" ref_value = " << std::setprecision (
7474 std::numeric_limits<long double >::digits10)
@@ -79,7 +79,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
7979 REQUIRE (value == Approx (gf_value));
8080 }
8181 AND_THEN (" the value of the Green's function outside the droplet is" ) {
82- double value = 0.488060362514614432 ;
82+ double value = 0.487473605103865837 ;
8383 double gf_value = gf.kernelS (source2, probe2);
8484 INFO (" ref_value = " << std::setprecision (
8585 std::numeric_limits<long double >::digits10)
@@ -91,7 +91,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
9191 }
9292 THEN (" the value of the Green's function directional derivative wrt the probe "
9393 " point inside the droplet is" ) {
94- double derProbe = -0.0124999765635167015 ;
94+ double derProbe = -0.0124999765637075211 ;
9595 double gf_derProbe = gf.derivativeProbe (probeNormal1, source1, probe1);
9696 INFO (" ref_derProbe = "
9797 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -103,7 +103,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
103103 }
104104 AND_THEN (" the value of the Green's function directional derivative wrt the "
105105 " probe point outside the droplet is" ) {
106- double derProbe = -0.290662079075743041 ;
106+ double derProbe = -0.290662955233778053 ;
107107 double gf_derProbe = gf.derivativeProbe (probeNormal2, source2, probe2);
108108 INFO (" ref_derProbe = "
109109 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -115,7 +115,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
115115 }
116116 THEN (" the value of the Green's function directional derivative wrt the source "
117117 " point inside the droplet is" ) {
118- double derSource = 0.00937504723103055326 ;
118+ double derSource = 0.00937504723101320603 ;
119119 double gf_derSource = gf.derivativeSource (sourceNormal1, source1, probe1);
120120 INFO (" ref_derSource = "
121121 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -127,7 +127,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
127127 }
128128 AND_THEN (" the value of the Green's function directional derivative wrt the "
129129 " source point outside the droplet is" ) {
130- double derSource = -0.938830970128590181 ;
130+ double derSource = -0.901954472843635724 ;
131131 double gf_derSource = gf.derivativeSource (sourceNormal2, source2, probe2);
132132 INFO (" ref_derSource = "
133133 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -144,7 +144,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
144144 (Eigen::Vector3d () << 25.0 , 0.0 , 0.0 ).finished ();
145145 SphericalDiffuse<> gf (eps1, eps2, width, sphereRadius, sphereCenter, maxL);
146146 THEN (" the value of the Green's function inside the droplet is" ) {
147- double value = 0.0173845118865096904 ;
147+ double value = 0.0173845118865097424 ;
148148 double gf_value = gf.kernelS (source1, probe1);
149149 INFO (" ref_value = " << std::setprecision (
150150 std::numeric_limits<long double >::digits10)
@@ -155,7 +155,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
155155 REQUIRE (value == Approx (gf_value));
156156 }
157157 AND_THEN (" the value of the Green's function outside the droplet is" ) {
158- double value = 0.501844654162973303 ;
158+ double value = 0.501385603574595162 ;
159159 double gf_value = gf.kernelS (source2, probe2);
160160 INFO (" ref_value = " << std::setprecision (
161161 std::numeric_limits<long double >::digits10)
@@ -167,7 +167,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
167167 }
168168 THEN (" the value of the Green's function directional derivative wrt the probe "
169169 " point inside the droplet is" ) {
170- double derProbe = -0.012483438527558649 ;
170+ double derProbe = -0.0124834385275933435 ;
171171 double gf_derProbe = gf.derivativeProbe (probeNormal1, source1, probe1);
172172 INFO (" ref_derProbe = "
173173 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -179,7 +179,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
179179 }
180180 AND_THEN (" the value of the Green's function directional derivative wrt the "
181181 " probe point outside the droplet is" ) {
182- double derProbe = -0.291257922496734878 ;
182+ double derProbe = -0.291262383059764929 ;
183183 double gf_derProbe = gf.derivativeProbe (probeNormal2, source2, probe2);
184184 INFO (" ref_derProbe = "
185185 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -191,7 +191,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
191191 }
192192 THEN (" the value of the Green's function directional derivative wrt the source "
193193 " point inside the droplet is" ) {
194- double derSource = 0.0124835646475064677 ;
194+ double derSource = 0.012483564647454426 ;
195195 double gf_derSource = gf.derivativeSource (sourceNormal1, source1, probe1);
196196 INFO (" ref_derSource = "
197197 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -203,7 +203,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
203203 }
204204 AND_THEN (" the value of the Green's function directional derivative wrt the "
205205 " source point outside the droplet is" ) {
206- double derSource = 5.01504540067754245 ;
206+ double derSource = 2.720604738008503 ;
207207 double gf_derSource = gf.derivativeSource (sourceNormal2, source2, probe2);
208208 INFO (" ref_derSource = "
209209 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -221,7 +221,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
221221 SphericalDiffuse<OneLayerTanh> gf (
222222 eps1, eps2, width, sphereRadius, sphereCenter, maxL);
223223 THEN (" the value of the Green's function inside the droplet is" ) {
224- double value = 0.0204265162808782985 ;
224+ double value = 0.0204265162808963327 ;
225225 double gf_value = gf.kernelS (source1, probe1);
226226 INFO (" ref_value = " << std::setprecision (
227227 std::numeric_limits<long double >::digits10)
@@ -232,7 +232,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
232232 REQUIRE (value == Approx (gf_value));
233233 }
234234 AND_THEN (" the value of the Green's function outside the droplet is" ) {
235- double value = 0.428395676083591359 ;
235+ double value = 0.42831756921770181 ;
236236 double gf_value = gf.kernelS (source2, probe2);
237237 INFO (" ref_value = " << std::setprecision (
238238 std::numeric_limits<long double >::digits10)
@@ -244,7 +244,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
244244 }
245245 THEN (" the value of the Green's function directional derivative wrt the probe "
246246 " point inside the droplet is" ) {
247- double derProbe = -0.0124999769449823939 ;
247+ double derProbe = -0.0124999769448609632 ;
248248 double gf_derProbe = gf.derivativeProbe (probeNormal1, source1, probe1);
249249 INFO (" ref_derProbe = "
250250 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -256,7 +256,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
256256 }
257257 AND_THEN (" the value of the Green's function directional derivative wrt the "
258258 " probe point outside the droplet is" ) {
259- double derProbe = -0.288385047571004804 ;
259+ double derProbe = -0.288392889131738883 ;
260260 double gf_derProbe = gf.derivativeProbe (probeNormal2, source2, probe2);
261261 INFO (" ref_derProbe = "
262262 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -280,7 +280,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
280280 }
281281 AND_THEN (" the value of the Green's function directional derivative wrt the "
282282 " source point outside the droplet is" ) {
283- double derSource = 3.22957321418654297 ;
283+ double derSource = - 1.51576208806569745 ;
284284 double gf_derSource = gf.derivativeSource (sourceNormal2, source2, probe2);
285285 INFO (" ref_derSource = "
286286 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -298,7 +298,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
298298 SphericalDiffuse<OneLayerTanh> gf (
299299 eps1, eps2, width, sphereRadius, sphereCenter, maxL);
300300 THEN (" the value of the Green's function inside the droplet is" ) {
301- double value = 0.0173358023608628994 ;
301+ double value = 0.017335802360862948 ;
302302 double gf_value = gf.kernelS (source1, probe1);
303303 INFO (" ref_value = " << std::setprecision (
304304 std::numeric_limits<long double >::digits10)
@@ -309,7 +309,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
309309 REQUIRE (value == Approx (gf_value));
310310 }
311311 AND_THEN (" the value of the Green's function outside the droplet is" ) {
312- double value = 0.471193749836353093 ;
312+ double value = 0.470725113611886292 ;
313313 double gf_value = gf.kernelS (source2, probe2);
314314 INFO (" ref_value = " << std::setprecision (
315315 std::numeric_limits<long double >::digits10)
@@ -321,7 +321,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
321321 }
322322 THEN (" the value of the Green's function directional derivative wrt the probe "
323323 " point inside the droplet is" ) {
324- double derProbe = -0.0124834504454732209 ;
324+ double derProbe = -0.0124834504460630269 ;
325325 double gf_derProbe = gf.derivativeProbe (probeNormal1, source1, probe1);
326326 INFO (" ref_derProbe = "
327327 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -333,7 +333,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
333333 }
334334 AND_THEN (" the value of the Green's function directional derivative wrt the "
335335 " probe point outside the droplet is" ) {
336- double derProbe = -0.290867341330436346 ;
336+ double derProbe = -0.290870450536662162 ;
337337 double gf_derProbe = gf.derivativeProbe (probeNormal2, source2, probe2);
338338 INFO (" ref_derProbe = "
339339 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -345,7 +345,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
345345 }
346346 THEN (" the value of the Green's function directional derivative wrt the source "
347347 " point inside the droplet is" ) {
348- double derSource = 0.0124835522706361057 ;
348+ double derSource = 0.012483552270584064 ;
349349 double gf_derSource = gf.derivativeSource (sourceNormal1, source1, probe1);
350350 INFO (" ref_derSource = "
351351 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -357,7 +357,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
357357 }
358358 AND_THEN (" the value of the Green's function directional derivative wrt the "
359359 " source point outside the droplet is" ) {
360- double derSource = 4.9830666203676266 ;
360+ double derSource = 3.75908425086779463 ;
361361 double gf_derSource = gf.derivativeSource (sourceNormal2, source2, probe2);
362362 INFO (" ref_derSource = "
363363 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -376,7 +376,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
376376 SphericalDiffuse<OneLayerErf> gf (
377377 eps1, eps2, width, sphereRadius, sphereCenter, maxL);
378378 THEN (" the value of the Green's function inside the droplet is" ) {
379- double value = 0.0204465601345453427 ;
379+ double value = 0.0204465601345278428 ;
380380 double gf_value = gf.kernelS (source1, probe1);
381381 INFO (" ref_value = " << std::setprecision (
382382 std::numeric_limits<long double >::digits10)
@@ -387,7 +387,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
387387 REQUIRE (value == Approx (gf_value));
388388 }
389389 AND_THEN (" the value of the Green's function outside the droplet is" ) {
390- double value = 0.546459075558889285 ;
390+ double value = 0.546440841068159378 ;
391391 double gf_value = gf.kernelS (source2, probe2);
392392 INFO (" ref_value = " << std::setprecision (
393393 std::numeric_limits<long double >::digits10)
@@ -399,7 +399,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
399399 }
400400 THEN (" the value of the Green's function directional derivative wrt the probe "
401401 " point inside the droplet is" ) {
402- double derProbe = -0.0124999769317464537 ;
402+ double derProbe = -0.0124999769314168563 ;
403403 double gf_derProbe = gf.derivativeProbe (probeNormal1, source1, probe1);
404404 INFO (" ref_derProbe = "
405405 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -411,7 +411,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
411411 }
412412 AND_THEN (" the value of the Green's function directional derivative wrt the "
413413 " probe point outside the droplet is" ) {
414- double derProbe = -0.292047321249766512 ;
414+ double derProbe = -0.292045737975143993 ;
415415 double gf_derProbe = gf.derivativeProbe (probeNormal2, source2, probe2);
416416 INFO (" ref_derProbe = "
417417 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -423,7 +423,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
423423 }
424424 THEN (" the value of the Green's function directional derivative wrt the source "
425425 " point inside the droplet is" ) {
426- double derSource = 0.00937504649521289646 ;
426+ double derSource = 0.00937504649523024369 ;
427427 double gf_derSource = gf.derivativeSource (sourceNormal1, source1, probe1);
428428 INFO (" ref_derSource = "
429429 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -435,7 +435,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
435435 }
436436 AND_THEN (" the value of the Green's function directional derivative wrt the "
437437 " source point outside the droplet is" ) {
438- double derSource = -2.1857233870498094 ;
438+ double derSource = -1.60173403890240262 ;
439439 double gf_derSource = gf.derivativeSource (sourceNormal2, source2, probe2);
440440 INFO (" ref_derSource = "
441441 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -453,7 +453,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
453453 SphericalDiffuse<OneLayerErf> gf (
454454 eps1, eps2, width, sphereRadius, sphereCenter, maxL);
455455 THEN (" the value of the Green's function inside the droplet is" ) {
456- double value = 0.0173558529449408631 ;
456+ double value = 0.0173558563009676217 ;
457457 double gf_value = gf.kernelS (source1, probe1);
458458 INFO (" ref_value = " << std::setprecision (
459459 std::numeric_limits<long double >::digits10)
@@ -464,7 +464,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
464464 REQUIRE (value == Approx (gf_value));
465465 }
466466 AND_THEN (" the value of the Green's function outside the droplet is" ) {
467- double value = 0.372262125530776311 ;
467+ double value = 0.37238664692877399 ;
468468 double gf_value = gf.kernelS (source2, probe2);
469469 INFO (" ref_value = " << std::setprecision (
470470 std::numeric_limits<long double >::digits10)
@@ -476,7 +476,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
476476 }
477477 THEN (" the value of the Green's function directional derivative wrt the probe "
478478 " point inside the droplet is" ) {
479- double derProbe = -0.0124834509424368023 ;
479+ double derProbe = -0.0124834500456541542 ;
480480 double gf_derProbe = gf.derivativeProbe (probeNormal1, source1, probe1);
481481 INFO (" ref_derProbe = "
482482 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -488,7 +488,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
488488 }
489489 AND_THEN (" the value of the Green's function directional derivative wrt the "
490490 " probe point outside the droplet is" ) {
491- double derProbe = -0.27916006373779334 ;
491+ double derProbe = -0.279220523679291066 ;
492492 double gf_derProbe = gf.derivativeProbe (probeNormal2, source2, probe2);
493493 INFO (" ref_derProbe = "
494494 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -500,7 +500,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
500500 }
501501 THEN (" the value of the Green's function directional derivative wrt the source "
502502 " point inside the droplet is" ) {
503- double derSource = 0.0124835526376168571 ;
503+ double derSource = 0.0124835526869870872 ;
504504 double gf_derSource = gf.derivativeSource (sourceNormal1, source1, probe1);
505505 INFO (" ref_derSource = "
506506 << std::setprecision (std::numeric_limits<long double >::digits10)
@@ -512,7 +512,7 @@ SCENARIO("Evaluation of the spherical diffuse Green's function and its derivativ
512512 }
513513 AND_THEN (" the value of the Green's function directional derivative wrt the "
514514 " source point outside the droplet is" ) {
515- double derSource = 3.4616694327527231 ;
515+ double derSource = 4.79631817304426722 ;
516516 double gf_derSource = gf.derivativeSource (sourceNormal2, source2, probe2);
517517 INFO (" ref_derSource = "
518518 << std::setprecision (std::numeric_limits<long double >::digits10)
0 commit comments