@@ -1342,27 +1342,27 @@ validation_potential = validation_system.potential[1,:]
1342
1342
system3 = generate_gravitational (seed, n_bodies; radius_factor= 0.1 )
1343
1343
fmm. fmm! (system3; expansion_order= expansion_order, n_per_branch= n_per_branch, multipole_acceptance_criterion= multipole_acceptance_criterion, nearfield= true , farfield= true , unsort_bodies= true )
1344
1344
potential3 = system3. potential[1 ,:]
1345
- @test isapprox (maximum (abs .(potential3 - validation_potential)), 0.0 ; atol= 1e-10 )
1345
+ @test isapprox (maximum (abs .(potential3 - validation_potential)), 0.0 ; atol= 1e-9 )
1346
1346
1347
1347
system4 = generate_gravitational (seed, n_bodies; radius_factor= 0.1 )
1348
1348
fmm. fmm! ((system4,); expansion_order= expansion_order, n_per_branch= n_per_branch, multipole_acceptance_criterion= multipole_acceptance_criterion, nearfield= true , farfield= true , unsort_bodies= true )
1349
1349
potential4 = system4. potential[1 ,:]
1350
- @test isapprox (maximum (abs .(potential4 - validation_potential)), 0.0 ; atol= 1e-10 )
1350
+ @test isapprox (maximum (abs .(potential4 - validation_potential)), 0.0 ; atol= 1e-9 )
1351
1351
1352
1352
system5 = generate_gravitational (seed, n_bodies; radius_factor= 0.1 )
1353
1353
fmm. fmm! (system5, system5; expansion_order= expansion_order, n_per_branch_source= n_per_branch, n_per_branch_target= n_per_branch, multipole_acceptance_criterion= multipole_acceptance_criterion, nearfield= true , farfield= true , unsort_source_bodies= true , unsort_target_bodies= true )
1354
1354
potential5 = system5. potential[1 ,:]
1355
- @test isapprox (maximum (abs .(potential5 - validation_potential)), 0.0 ; atol= 1e-10 )
1355
+ @test isapprox (maximum (abs .(potential5 - validation_potential)), 0.0 ; atol= 1e-9 )
1356
1356
1357
1357
system6 = generate_gravitational (seed, n_bodies; radius_factor= 0.1 )
1358
1358
fmm. fmm! ((system6,), system6; expansion_order= expansion_order, n_per_branch_source= n_per_branch, n_per_branch_target= n_per_branch, multipole_acceptance_criterion= multipole_acceptance_criterion, nearfield= true , farfield= true , unsort_source_bodies= true , unsort_target_bodies= true )
1359
1359
potential6 = system6. potential[1 ,:]
1360
- @test isapprox (maximum (abs .(potential6 - validation_potential)), 0.0 ; atol= 1e-10 )
1360
+ @test isapprox (maximum (abs .(potential6 - validation_potential)), 0.0 ; atol= 1e-9 )
1361
1361
1362
1362
system7 = generate_gravitational (seed, n_bodies; radius_factor= 0.1 )
1363
1363
fmm. fmm! ((system7,), (system7,); expansion_order= expansion_order, n_per_branch_source= n_per_branch, n_per_branch_target= n_per_branch, multipole_acceptance_criterion= multipole_acceptance_criterion, nearfield= true , farfield= true , unsort_source_bodies= true , unsort_target_bodies= true )
1364
1364
potential7 = system7. potential[1 ,:]
1365
- @test isapprox (maximum (abs .(potential7 - validation_potential)), 0.0 ; atol= 1e-10 )
1365
+ @test isapprox (maximum (abs .(potential7 - validation_potential)), 0.0 ; atol= 1e-9 )
1366
1366
1367
1367
end
1368
1368
@@ -1379,27 +1379,27 @@ validation_potential = validation_system.potential[1,:]
1379
1379
system8 = fmm. SortWrapper (generate_gravitational (seed, n_bodies; radius_factor= 0.1 ))
1380
1380
fmm. fmm! (system8; expansion_order= expansion_order, n_per_branch= n_per_branch, multipole_acceptance_criterion= multipole_acceptance_criterion, nearfield= true , farfield= true , unsort_bodies= true )
1381
1381
potential8 = system8. system. potential[1 ,:]
1382
- @test isapprox (maximum (abs .(potential8 - validation_potential)), 0.0 ; atol= 1e-10 )
1382
+ @test isapprox (maximum (abs .(potential8 - validation_potential)), 0.0 ; atol= 1e-9 )
1383
1383
1384
1384
system9 = fmm. SortWrapper (generate_gravitational (seed, n_bodies; radius_factor= 0.1 ))
1385
1385
fmm. fmm! ((system9,); expansion_order= expansion_order, n_per_branch= n_per_branch, multipole_acceptance_criterion= multipole_acceptance_criterion, nearfield= true , farfield= true , unsort_bodies= true )
1386
1386
potential9 = system9. system. potential[1 ,:]
1387
- @test isapprox (maximum (abs .(potential9 - validation_potential)), 0.0 ; atol= 1e-10 )
1387
+ @test isapprox (maximum (abs .(potential9 - validation_potential)), 0.0 ; atol= 1e-9 )
1388
1388
1389
1389
system10 = fmm. SortWrapper (generate_gravitational (seed, n_bodies; radius_factor= 0.1 ))
1390
1390
fmm. fmm! (system10, system10; expansion_order= expansion_order, n_per_branch_source= n_per_branch, n_per_branch_target= n_per_branch, multipole_acceptance_criterion= multipole_acceptance_criterion, nearfield= true , farfield= true , unsort_source_bodies= true , unsort_target_bodies= true )
1391
1391
potential10 = system10. system. potential[1 ,:]
1392
- @test isapprox (maximum (abs .(potential10 - validation_potential)), 0.0 ; atol= 1e-10 )
1392
+ @test isapprox (maximum (abs .(potential10 - validation_potential)), 0.0 ; atol= 1e-9 )
1393
1393
1394
1394
system11 = fmm. SortWrapper (generate_gravitational (seed, n_bodies; radius_factor= 0.1 ))
1395
1395
fmm. fmm! ((system11,), system11; expansion_order= expansion_order, n_per_branch_source= n_per_branch, n_per_branch_target= n_per_branch, multipole_acceptance_criterion= multipole_acceptance_criterion, nearfield= true , farfield= true , unsort_source_bodies= true , unsort_target_bodies= true )
1396
1396
potential11 = system11. system. potential[1 ,:]
1397
- @test isapprox (maximum (abs .(potential11 - validation_potential)), 0.0 ; atol= 1e-10 )
1397
+ @test isapprox (maximum (abs .(potential11 - validation_potential)), 0.0 ; atol= 1e-9 )
1398
1398
1399
1399
system12 = fmm. SortWrapper (generate_gravitational (seed, n_bodies; radius_factor= 0.1 ))
1400
1400
fmm. fmm! ((system12,), (system12,); expansion_order= expansion_order, n_per_branch_source= n_per_branch, n_per_branch_target= n_per_branch, multipole_acceptance_criterion= multipole_acceptance_criterion, nearfield= true , farfield= true , unsort_source_bodies= true , unsort_target_bodies= true )
1401
1401
potential12 = system12. system. potential[1 ,:]
1402
- @test isapprox (maximum (abs .(potential12 - validation_potential)), 0.0 ; atol= 1e-10 )
1402
+ @test isapprox (maximum (abs .(potential12 - validation_potential)), 0.0 ; atol= 1e-9 )
1403
1403
1404
1404
end
1405
1405
0 commit comments