Post-processing: plot and analyze the results

This documentation is autogenerated from ./tests/exaparty/planet-orbit/plots-time-integration-test.ipynb. Once you clone the Peano repository, please work against this notebook or study Peano's generic description how to generate web docu from Jupyter notebooks.

import meshio
import numpy as np
import glob, os
 
import matplotlib.pyplot as plt
from matplotlib import gridspec
from matplotlib.colors import LogNorm
from matplotlib import mathtext
import math
from scipy import stats
from tkinter import Tcl
 
#Plot Setup
font = {'size':12}
plt.rc('font', **font)
plt.rc('font', family='serif')
plt.rc('text', usetex=False)

Measure total Energy

# Import snapshots
snapshot_dir = ''
#snapshot_dir = 'run_const_dt/'
files = glob.glob(snapshot_dir + "particles-*.vtu")
# Sort the snapshots
files = Tcl().call('lsort', '-dict', files)
 
# Simulation parameters
E_ini = 1.11362
time_step_sim = 1e-3
snapshot_rate = 25
SUN_X = 0.5
SUN_Y = 0.55
 
# Initialize arrays
E_dif_all = np.zeros( len(files) )
t_all = np.zeros( len(files) )
 
# Read snapshots and gather data
for N_snapshot, file in enumerate(files, start=0):
    
    print("N_snapshot: ", N_snapshot)
    time = N_snapshot * snapshot_rate * time_step_sim
    t_all[N_snapshot] = time
    
    print('Reading snapshot: ' + file )
    print("time of simulation = ", time)
    print("********************************************")
 
    mesh = meshio.read( file )
    print(mesh.point_data_to_sets)
    print("********************************************")
 
    # Read coordinates
    pos = mesh.point_data["x"]
    x = pos[:,0]
    y = pos[:,1]
 
    # Read energy
    E_tot = mesh.point_data["energyTot"][:,0]
    E_kin = mesh.point_data["energyKin"][:,0]
    E_pot = mesh.point_data["energyPot"][:,0]
    E_tot = np.sum(E_tot)
    E_dif = (E_tot - E_ini) / E_ini
    E_dif_all[N_snapshot] = E_dif
 
    print("Measure energy conservation:")
    print("Initial E = ", E_ini)
    print("Kinetic E = ", E_kin)
    print("Potential E = ", E_pot)
    print("Total E = ", E_tot)
    print("Relative Energy difference = ", E_dif )
    print("********************************************")
 
print("Done reading all snapshots. Total= ", N_snapshot )

N_snapshot:  0
Reading snapshot: particles-0.vtu
time of simulation =  0.0

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.2341]
Potential E =  [2.34772]
Total E =  1.11362
Relative Energy difference =  3.8434028902163246e-08

N_snapshot:  1
Reading snapshot: particles-1.vtu
time of simulation =  0.025

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.23276]
Potential E =  [2.34631]
Total E =  1.11355
Relative Energy difference =  -6.290499573695996e-05

N_snapshot:  2
Reading snapshot: particles-2.vtu
time of simulation =  0.05

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.22863]
Potential E =  [2.34213]
Total E =  1.1135
Relative Energy difference =  -0.0001077575417775998

N_snapshot:  3
Reading snapshot: particles-3.vtu
time of simulation =  0.075

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.22186]
Potential E =  [2.33531]
Total E =  1.11345
Relative Energy difference =  -0.00015261008781823965

N_snapshot:  4
Reading snapshot: particles-4.vtu
time of simulation =  0.1

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.21268]
Potential E =  [2.32609]
Total E =  1.1134
Relative Energy difference =  -0.00019756968050814116

N_snapshot:  5
Reading snapshot: particles-5.vtu
time of simulation =  0.125

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.20138]
Potential E =  [2.31475]
Total E =  1.11337
Relative Energy difference =  -0.00022454543612208206

N_snapshot:  6
Reading snapshot: particles-6.vtu
time of simulation =  0.15

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.1883]
Potential E =  [2.30164]
Total E =  1.11333
Relative Energy difference =  -0.0002604060636247416

N_snapshot:  7
Reading snapshot: particles-7.vtu
time of simulation =  0.17500000000000002

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.17382]
Potential E =  [2.28713]
Total E =  1.11331
Relative Energy difference =  -0.0002783899007007022

N_snapshot:  8
Reading snapshot: particles-8.vtu
time of simulation =  0.2

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.15833]
Potential E =  [2.27163]
Total E =  1.11329
Relative Energy difference =  -0.0002963737377766628

N_snapshot:  9
Reading snapshot: particles-9.vtu
time of simulation =  0.225

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.14222]
Potential E =  [2.25551]
Total E =  1.11329
Relative Energy difference =  -0.0002963737377766628

N_snapshot:  10
Reading snapshot: particles-10.vtu
time of simulation =  0.25

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.12586]
Potential E =  [2.23915]
Total E =  1.11328
Relative Energy difference =  -0.0003052586096653814

N_snapshot:  11
Reading snapshot: particles-11.vtu
time of simulation =  0.275

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.10961]
Potential E =  [2.2229]
Total E =  1.11329
Relative Energy difference =  -0.0002963737377766628

N_snapshot:  12
Reading snapshot: particles-12.vtu
time of simulation =  0.3

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.09379]
Potential E =  [2.20709]
Total E =  1.1133
Relative Energy difference =  -0.0002873818192386825

N_snapshot:  13
Reading snapshot: particles-13.vtu
time of simulation =  0.325

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.07868]
Potential E =  [2.19199]
Total E =  1.11332
Relative Energy difference =  -0.0002693979821627219

N_snapshot:  14
Reading snapshot: particles-14.vtu
time of simulation =  0.35000000000000003

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.06452]
Potential E =  [2.17786]
Total E =  1.11334
Relative Energy difference =  -0.0002514141450867613

N_snapshot:  15
Reading snapshot: particles-15.vtu
time of simulation =  0.375

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.05154]
Potential E =  [2.16491]
Total E =  1.11337
Relative Energy difference =  -0.00022454543612208206

N_snapshot:  16
Reading snapshot: particles-16.vtu
time of simulation =  0.4

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.03991]
Potential E =  [2.1533]
Total E =  1.1134
Relative Energy difference =  -0.00019756968050814116

N_snapshot:  17
Reading snapshot: particles-17.vtu
time of simulation =  0.425

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.02978]
Potential E =  [2.1432]
Total E =  1.11343
Relative Energy difference =  -0.00017059392489420023

N_snapshot:  18
Reading snapshot: particles-18.vtu
time of simulation =  0.45

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.02126]
Potential E =  [2.13472]
Total E =  1.11346
Relative Energy difference =  -0.00014372521592952102

N_snapshot:  19
Reading snapshot: particles-19.vtu
time of simulation =  0.47500000000000003

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.01446]
Potential E =  [2.12796]
Total E =  1.1135
Relative Energy difference =  -0.0001077575417775998

N_snapshot:  20
Reading snapshot: particles-20.vtu
time of simulation =  0.5

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00944]
Potential E =  [2.12297]
Total E =  1.11353
Relative Energy difference =  -8.078178616365888e-05

N_snapshot:  21
Reading snapshot: particles-21.vtu
time of simulation =  0.525

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00625]
Potential E =  [2.11982]
Total E =  1.11357
Relative Energy difference =  -4.492115866099935e-05

N_snapshot:  22
Reading snapshot: particles-22.vtu
time of simulation =  0.55

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00492]
Potential E =  [2.11852]
Total E =  1.11361
Relative Energy difference =  -8.953484509078139e-06

N_snapshot:  23
Reading snapshot: particles-23.vtu
time of simulation =  0.5750000000000001

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00546]
Potential E =  [2.11911]
Total E =  1.11364
Relative Energy difference =  1.79152244556011e-05

N_snapshot:  24
Reading snapshot: particles-24.vtu
time of simulation =  0.6

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00788]
Potential E =  [2.12156]
Total E =  1.11368
Relative Energy difference =  5.3882898607522306e-05

N_snapshot:  25
Reading snapshot: particles-25.vtu
time of simulation =  0.625

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.01214]
Potential E =  [2.12585]
Total E =  1.11372
Relative Energy difference =  8.985057275944351e-05

N_snapshot:  26
Reading snapshot: particles-26.vtu
time of simulation =  0.65

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.0182]
Potential E =  [2.13195]
Total E =  1.11375
Relative Energy difference =  0.00011671928172412275

N_snapshot:  27
Reading snapshot: particles-27.vtu
time of simulation =  0.675

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.02601]
Potential E =  [2.13979]
Total E =  1.11379
Relative Energy difference =  0.00015268695587604398

N_snapshot:  28
Reading snapshot: particles-28.vtu
time of simulation =  0.7000000000000001

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.03547]
Potential E =  [2.14929]
Total E =  1.11382
Relative Energy difference =  0.0001795556648407232

N_snapshot:  29
Reading snapshot: particles-29.vtu
time of simulation =  0.725

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.04649]
Potential E =  [2.16034]
Total E =  1.11385
Relative Energy difference =  0.00020653142045466412

N_snapshot:  30
Reading snapshot: particles-30.vtu
time of simulation =  0.75

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.05893]
Potential E =  [2.1728]
Total E =  1.11387
Relative Energy difference =  0.00022451525753062473

N_snapshot:  31
Reading snapshot: particles-31.vtu
time of simulation =  0.775

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.07262]
Potential E =  [2.18651]
Total E =  1.1139
Relative Energy difference =  0.00025138396649530397

N_snapshot:  32
Reading snapshot: particles-32.vtu
time of simulation =  0.8

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.08736]
Potential E =  [2.20127]
Total E =  1.11392
Relative Energy difference =  0.00026936780357126455

N_snapshot:  33
Reading snapshot: particles-33.vtu
time of simulation =  0.8250000000000001

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.10292]
Potential E =  [2.21685]
Total E =  1.11393
Relative Energy difference =  0.00027835972210924485

N_snapshot:  34
Reading snapshot: particles-34.vtu
time of simulation =  0.85

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.11902]
Potential E =  [2.23296]
Total E =  1.11394
Relative Energy difference =  0.00028735164064722514

N_snapshot:  35
Reading snapshot: particles-35.vtu
time of simulation =  0.875

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.13538]
Potential E =  [2.24932]
Total E =  1.11394
Relative Energy difference =  0.00028735164064722514

N_snapshot:  36
Reading snapshot: particles-36.vtu
time of simulation =  0.9

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.15164]
Potential E =  [2.26557]
Total E =  1.11393
Relative Energy difference =  0.00027835972210924485

N_snapshot:  37
Reading snapshot: particles-37.vtu
time of simulation =  0.925

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.16743]
Potential E =  [2.28135]
Total E =  1.11392
Relative Energy difference =  0.00026936780357126455

N_snapshot:  38
Reading snapshot: particles-38.vtu
time of simulation =  0.9500000000000001

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.18239]
Potential E =  [2.29629]
Total E =  1.1139
Relative Energy difference =  0.00025138396649530397

N_snapshot:  39
Reading snapshot: particles-39.vtu
time of simulation =  0.975

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.1961]
Potential E =  [2.30997]
Total E =  1.11387
Relative Energy difference =  0.00022451525753062473

N_snapshot:  40
Reading snapshot: particles-40.vtu
time of simulation =  1.0

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.20819]
Potential E =  [2.32202]
Total E =  1.11383
Relative Energy difference =  0.0001885475833787035

N_snapshot:  41
Reading snapshot: particles-41.vtu
time of simulation =  1.025

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.21866]
Potential E =  [2.33244]
Total E =  1.11379
Relative Energy difference =  0.00015268695587604398

N_snapshot:  42
Reading snapshot: particles-42.vtu
time of simulation =  1.05

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.22636]
Potential E =  [2.3401]
Total E =  1.11374
Relative Energy difference =  0.00010772736318614245

N_snapshot:  43
Reading snapshot: particles-43.vtu
time of simulation =  1.075

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.23151]
Potential E =  [2.3452]
Total E =  1.11369
Relative Energy difference =  6.287481714550261e-05

N_snapshot:  44
Reading snapshot: particles-44.vtu
time of simulation =  1.1

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.23393]
Potential E =  [2.34756]
Total E =  1.11363
Relative Energy difference =  9.030352566882466e-06

N_snapshot:  45
Reading snapshot: particles-45.vtu
time of simulation =  1.125

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.23353]
Potential E =  [2.34711]
Total E =  1.11357
Relative Energy difference =  -4.492115866099935e-05

N_snapshot:  46
Reading snapshot: particles-46.vtu
time of simulation =  1.1500000000000001

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.23033]
Potential E =  [2.34385]
Total E =  1.11352
Relative Energy difference =  -8.977370470163919e-05

N_snapshot:  47
Reading snapshot: particles-47.vtu
time of simulation =  1.175

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.22444]
Potential E =  [2.33791]
Total E =  1.11347
Relative Energy difference =  -0.0001347332973915407

N_snapshot:  48
Reading snapshot: particles-48.vtu
time of simulation =  1.2

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.21605]
Potential E =  [2.32947]
Total E =  1.11342
Relative Energy difference =  -0.00017958584343218055

N_snapshot:  49
Reading snapshot: particles-49.vtu
time of simulation =  1.225

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.20543]
Potential E =  [2.31881]
Total E =  1.11338
Relative Energy difference =  -0.00021555351758410175

N_snapshot:  50
Reading snapshot: particles-50.vtu
time of simulation =  1.25

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.19291]
Potential E =  [2.30626]
Total E =  1.11334
Relative Energy difference =  -0.0002514141450867613

N_snapshot:  51
Reading snapshot: particles-51.vtu
time of simulation =  1.2750000000000001

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.17886]
Potential E =  [2.29218]
Total E =  1.11332
Relative Energy difference =  -0.0002693979821627219

N_snapshot:  52
Reading snapshot: particles-52.vtu
time of simulation =  1.3

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.16367]
Potential E =  [2.27697]
Total E =  1.1133
Relative Energy difference =  -0.0002873818192386825

N_snapshot:  53
Reading snapshot: particles-53.vtu
time of simulation =  1.325

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.14773]
Potential E =  [2.26101]
Total E =  1.11329
Relative Energy difference =  -0.0002963737377766628

N_snapshot:  54
Reading snapshot: particles-54.vtu
time of simulation =  1.35

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.13141]
Potential E =  [2.2447]
Total E =  1.11328
Relative Energy difference =  -0.0003052586096653814

N_snapshot:  55
Reading snapshot: particles-55.vtu
time of simulation =  1.375

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.11509]
Potential E =  [2.22838]
Total E =  1.11329
Relative Energy difference =  -0.0002963737377766628

N_snapshot:  56
Reading snapshot: particles-56.vtu
time of simulation =  1.4000000000000001

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.09909]
Potential E =  [2.21238]
Total E =  1.1133
Relative Energy difference =  -0.0002873818192386825

N_snapshot:  57
Reading snapshot: particles-57.vtu
time of simulation =  1.425

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.0837]
Potential E =  [2.19701]
Total E =  1.11331
Relative Energy difference =  -0.0002783899007007022

N_snapshot:  58
Reading snapshot: particles-58.vtu
time of simulation =  1.45

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.0692]
Potential E =  [2.18253]
Total E =  1.11333
Relative Energy difference =  -0.0002604060636247416

N_snapshot:  59
Reading snapshot: particles-59.vtu
time of simulation =  1.475

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.0558]
Potential E =  [2.16915]
Total E =  1.11336
Relative Energy difference =  -0.0002334303080108007

N_snapshot:  60
Reading snapshot: particles-60.vtu
time of simulation =  1.5

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.04369]
Potential E =  [2.15707]
Total E =  1.11339
Relative Energy difference =  -0.00020656159904612145

N_snapshot:  61
Reading snapshot: particles-61.vtu
time of simulation =  1.5250000000000001

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.03304]
Potential E =  [2.14645]
Total E =  1.11342
Relative Energy difference =  -0.00017958584343218055

N_snapshot:  62
Reading snapshot: particles-62.vtu
time of simulation =  1.55

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.02396]
Potential E =  [2.13741]
Total E =  1.11345
Relative Energy difference =  -0.00015261008781823965

N_snapshot:  63
Reading snapshot: particles-63.vtu
time of simulation =  1.575

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.01657]
Potential E =  [2.13005]
Total E =  1.11348
Relative Energy difference =  -0.0001257413788535604

N_snapshot:  64
Reading snapshot: particles-64.vtu
time of simulation =  1.6

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.01094]
Potential E =  [2.12446]
Total E =  1.11352
Relative Energy difference =  -8.977370470163919e-05

N_snapshot:  65
Reading snapshot: particles-65.vtu
time of simulation =  1.625

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00712]
Potential E =  [2.12068]
Total E =  1.11356
Relative Energy difference =  -5.3913077198979655e-05

N_snapshot:  66
Reading snapshot: particles-66.vtu
time of simulation =  1.6500000000000001

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00516]
Potential E =  [2.11875]
Total E =  1.11359
Relative Energy difference =  -2.6937321585038743e-05

N_snapshot:  67
Reading snapshot: particles-67.vtu
time of simulation =  1.675

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00507]
Potential E =  [2.1187]
Total E =  1.11363
Relative Energy difference =  9.030352566882466e-06

N_snapshot:  68
Reading snapshot: particles-68.vtu
time of simulation =  1.7

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00685]
Potential E =  [2.12052]
Total E =  1.11367
Relative Energy difference =  4.489098006954201e-05

N_snapshot:  69
Reading snapshot: particles-69.vtu
time of simulation =  1.725

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.01049]
Potential E =  [2.12419]
Total E =  1.1137
Relative Energy difference =  7.186673568348292e-05

N_snapshot:  70
Reading snapshot: particles-70.vtu
time of simulation =  1.75

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.01595]
Potential E =  [2.12969]
Total E =  1.11374
Relative Energy difference =  0.00010772736318614245

N_snapshot:  71
Reading snapshot: particles-71.vtu
time of simulation =  1.7750000000000001

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.02317]
Potential E =  [2.13695]
Total E =  1.11378
Relative Energy difference =  0.00014369503733806366

N_snapshot:  72
Reading snapshot: particles-72.vtu
time of simulation =  1.8

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.03209]
Potential E =  [2.14589]
Total E =  1.11381
Relative Energy difference =  0.0001705637463027429

N_snapshot:  73
Reading snapshot: particles-73.vtu
time of simulation =  1.825

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.04259]
Potential E =  [2.15643]
Total E =  1.11384
Relative Energy difference =  0.0001975395019166838

N_snapshot:  74
Reading snapshot: particles-74.vtu
time of simulation =  1.85

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.05457]
Potential E =  [2.16843]
Total E =  1.11387
Relative Energy difference =  0.00022451525753062473

N_snapshot:  75
Reading snapshot: particles-75.vtu
time of simulation =  1.875

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.06785]
Potential E =  [2.18174]
Total E =  1.11389
Relative Energy difference =  0.0002424990946065853

N_snapshot:  76
Reading snapshot: particles-76.vtu
time of simulation =  1.9000000000000001

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.08226]
Potential E =  [2.19617]
Total E =  1.11391
Relative Energy difference =  0.00026037588503328426

N_snapshot:  77
Reading snapshot: particles-77.vtu
time of simulation =  1.925

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.09757]
Potential E =  [2.21149]
Total E =  1.11393
Relative Energy difference =  0.00027835972210924485

N_snapshot:  78
Reading snapshot: particles-78.vtu
time of simulation =  1.95

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.11352]
Potential E =  [2.22746]
Total E =  1.11394
Relative Energy difference =  0.00028735164064722514

N_snapshot:  79
Reading snapshot: particles-79.vtu
time of simulation =  1.975

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.12983]
Potential E =  [2.24377]
Total E =  1.11394
Relative Energy difference =  0.00028735164064722514

N_snapshot:  80
Reading snapshot: particles-80.vtu
time of simulation =  2.0

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.14616]
Potential E =  [2.26009]
Total E =  1.11394
Relative Energy difference =  0.00028735164064722514

N_snapshot:  81
Reading snapshot: particles-81.vtu
time of simulation =  2.025

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.16216]
Potential E =  [2.27608]
Total E =  1.11393
Relative Energy difference =  0.00027835972210924485

N_snapshot:  82
Reading snapshot: particles-82.vtu
time of simulation =  2.05

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.17744]
Potential E =  [2.29135]
Total E =  1.11391
Relative Energy difference =  0.00026037588503328426

N_snapshot:  83
Reading snapshot: particles-83.vtu
time of simulation =  2.075

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.19162]
Potential E =  [2.3055]
Total E =  1.11388
Relative Energy difference =  0.00023350717606860502

N_snapshot:  84
Reading snapshot: particles-84.vtu
time of simulation =  2.1

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.2043]
Potential E =  [2.31814]
Total E =  1.11385
Relative Energy difference =  0.00020653142045466412

N_snapshot:  85
Reading snapshot: particles-85.vtu
time of simulation =  2.125

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.21511]
Potential E =  [2.32892]
Total E =  1.11381
Relative Energy difference =  0.0001705637463027429

N_snapshot:  86
Reading snapshot: particles-86.vtu
time of simulation =  2.15

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.22373]
Potential E =  [2.33749]
Total E =  1.11376
Relative Energy difference =  0.00012571120026210305

N_snapshot:  87
Reading snapshot: particles-87.vtu
time of simulation =  2.1750000000000003

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.22988]
Potential E =  [2.34358]
Total E =  1.11371
Relative Energy difference =  8.085865422146322e-05

N_snapshot:  88
Reading snapshot: particles-88.vtu
time of simulation =  2.2

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.23334]
Potential E =  [2.347]
Total E =  1.11365
Relative Energy difference =  2.69071429935814e-05

N_snapshot:  89
Reading snapshot: particles-89.vtu
time of simulation =  2.225

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.23402]
Potential E =  [2.34761]
Total E =  1.1136
Relative Energy difference =  -1.794540304705844e-05

N_snapshot:  90
Reading snapshot: particles-90.vtu
time of simulation =  2.25

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.23187]
Potential E =  [2.34541]
Total E =  1.11354
Relative Energy difference =  -7.178986762567859e-05

N_snapshot:  91
Reading snapshot: particles-91.vtu
time of simulation =  2.275

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.22698]
Potential E =  [2.34046]
Total E =  1.11349
Relative Energy difference =  -0.0001167494603155801

N_snapshot:  92
Reading snapshot: particles-92.vtu
time of simulation =  2.3000000000000003

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.21951]
Potential E =  [2.33294]
Total E =  1.11344
Relative Energy difference =  -0.00016160200635621994

N_snapshot:  93
Reading snapshot: particles-93.vtu
time of simulation =  2.325

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.2097]
Potential E =  [2.32309]
Total E =  1.11339
Relative Energy difference =  -0.00020656159904612145

N_snapshot:  94
Reading snapshot: particles-94.vtu
time of simulation =  2.35

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.19786]
Potential E =  [2.31121]
Total E =  1.11336
Relative Energy difference =  -0.0002334303080108007

N_snapshot:  95
Reading snapshot: particles-95.vtu
time of simulation =  2.375

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.18434]
Potential E =  [2.29767]
Total E =  1.11333
Relative Energy difference =  -0.0002604060636247416

N_snapshot:  96
Reading snapshot: particles-96.vtu
time of simulation =  2.4

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.16954]
Potential E =  [2.28284]
Total E =  1.1133
Relative Energy difference =  -0.0002873818192386825

N_snapshot:  97
Reading snapshot: particles-97.vtu
time of simulation =  2.4250000000000003

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.15383]
Potential E =  [2.26712]
Total E =  1.11329
Relative Energy difference =  -0.0002963737377766628

N_snapshot:  98
Reading snapshot: particles-98.vtu
time of simulation =  2.45

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.13761]
Potential E =  [2.2509]
Total E =  1.11328
Relative Energy difference =  -0.0003052586096653814

N_snapshot:  99
Reading snapshot: particles-99.vtu
time of simulation =  2.475

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.12125]
Potential E =  [2.23454]
Total E =  1.11329
Relative Energy difference =  -0.0002963737377766628

N_snapshot:  100
Reading snapshot: particles-100.vtu
time of simulation =  2.5

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.10509]
Potential E =  [2.21839]
Total E =  1.11329
Relative Energy difference =  -0.0002963737377766628

N_snapshot:  101
Reading snapshot: particles-101.vtu
time of simulation =  2.525

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.08944]
Potential E =  [2.20275]
Total E =  1.11331
Relative Energy difference =  -0.0002783899007007022

N_snapshot:  102
Reading snapshot: particles-102.vtu
time of simulation =  2.5500000000000003

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.07457]
Potential E =  [2.1879]
Total E =  1.11333
Relative Energy difference =  -0.0002604060636247416

N_snapshot:  103
Reading snapshot: particles-103.vtu
time of simulation =  2.575

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.06073]
Potential E =  [2.17408]
Total E =  1.11335
Relative Energy difference =  -0.000242422226548781

N_snapshot:  104
Reading snapshot: particles-104.vtu
time of simulation =  2.6

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.04811]
Potential E =  [2.16149]
Total E =  1.11337
Relative Energy difference =  -0.00022454543612208206

N_snapshot:  105
Reading snapshot: particles-105.vtu
time of simulation =  2.625

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.03689]
Potential E =  [2.1503]
Total E =  1.1134
Relative Energy difference =  -0.00019756968050814116

N_snapshot:  106
Reading snapshot: particles-106.vtu
time of simulation =  2.65

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.02721]
Potential E =  [2.14064]
Total E =  1.11344
Relative Energy difference =  -0.00016160200635621994

N_snapshot:  107
Reading snapshot: particles-107.vtu
time of simulation =  2.6750000000000003

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.01917]
Potential E =  [2.13264]
Total E =  1.11347
Relative Energy difference =  -0.0001347332973915407

N_snapshot:  108
Reading snapshot: particles-108.vtu
time of simulation =  2.7

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.01286]
Potential E =  [2.12637]
Total E =  1.11351
Relative Energy difference =  -9.876562323961949e-05

N_snapshot:  109
Reading snapshot: particles-109.vtu
time of simulation =  2.725

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00835]
Potential E =  [2.12189]
Total E =  1.11354
Relative Energy difference =  -7.178986762567859e-05

N_snapshot:  110
Reading snapshot: particles-110.vtu
time of simulation =  2.75

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00568]
Potential E =  [2.11926]
Total E =  1.11358
Relative Energy difference =  -3.5929240123019044e-05

N_snapshot:  111
Reading snapshot: particles-111.vtu
time of simulation =  2.775

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00488]
Potential E =  [2.1185]
Total E =  1.11362
Relative Energy difference =  3.8434028902163246e-08

N_snapshot:  112
Reading snapshot: particles-112.vtu
time of simulation =  2.8000000000000003

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00596]
Potential E =  [2.11961]
Total E =  1.11365
Relative Energy difference =  2.69071429935814e-05

N_snapshot:  113
Reading snapshot: particles-113.vtu
time of simulation =  2.825

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00889]
Potential E =  [2.12258]
Total E =  1.11369
Relative Energy difference =  6.287481714550261e-05

N_snapshot:  114
Reading snapshot: particles-114.vtu
time of simulation =  2.85

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.01367]
Potential E =  [2.12739]
Total E =  1.11373
Relative Energy difference =  9.873544464816214e-05

N_snapshot:  115
Reading snapshot: particles-115.vtu
time of simulation =  2.875

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.02023]
Potential E =  [2.13399]
Total E =  1.11376
Relative Energy difference =  0.00012571120026210305

N_snapshot:  116
Reading snapshot: particles-116.vtu
time of simulation =  2.9

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.02852]
Potential E =  [2.14231]
Total E =  1.1138
Relative Energy difference =  0.00016167887441402427

N_snapshot:  117
Reading snapshot: particles-117.vtu
time of simulation =  2.9250000000000003

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.03843]
Potential E =  [2.15226]
Total E =  1.11383
Relative Energy difference =  0.0001885475833787035

N_snapshot:  118
Reading snapshot: particles-118.vtu
time of simulation =  2.95

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.04987]
Potential E =  [2.16372]
Total E =  1.11386
Relative Energy difference =  0.0002155233389926444

N_snapshot:  119
Reading snapshot: particles-119.vtu
time of simulation =  2.975

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.06267]
Potential E =  [2.17655]
Total E =  1.11388
Relative Energy difference =  0.00023350717606860502

N_snapshot:  120
Reading snapshot: particles-120.vtu
time of simulation =  3.0

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.07668]
Potential E =  [2.19058]
Total E =  1.1139
Relative Energy difference =  0.00025138396649530397

N_snapshot:  121
Reading snapshot: particles-121.vtu
time of simulation =  3.025

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.09167]
Potential E =  [2.2056]
Total E =  1.11392
Relative Energy difference =  0.00026936780357126455

N_snapshot:  122
Reading snapshot: particles-122.vtu
time of simulation =  3.0500000000000003

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.10742]
Potential E =  [2.22135]
Total E =  1.11393
Relative Energy difference =  0.00027835972210924485

N_snapshot:  123
Reading snapshot: particles-123.vtu
time of simulation =  3.075

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.12363]
Potential E =  [2.23757]
Total E =  1.11394
Relative Energy difference =  0.00028735164064722514

N_snapshot:  124
Reading snapshot: particles-124.vtu
time of simulation =  3.1

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.13999]
Potential E =  [2.25393]
Total E =  1.11394
Relative Energy difference =  0.00028735164064722514

N_snapshot:  125
Reading snapshot: particles-125.vtu
time of simulation =  3.125

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.15616]
Potential E =  [2.27009]
Total E =  1.11393
Relative Energy difference =  0.00027835972210924485

N_snapshot:  126
Reading snapshot: particles-126.vtu
time of simulation =  3.15

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.17176]
Potential E =  [2.28567]
Total E =  1.11392
Relative Energy difference =  0.00026936780357126455

N_snapshot:  127
Reading snapshot: particles-127.vtu
time of simulation =  3.1750000000000003

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.1864]
Potential E =  [2.30029]
Total E =  1.11389
Relative Energy difference =  0.0002424990946065853

N_snapshot:  128
Reading snapshot: particles-128.vtu
time of simulation =  3.2

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.19969]
Potential E =  [2.31355]
Total E =  1.11386
Relative Energy difference =  0.0002155233389926444

N_snapshot:  129
Reading snapshot: particles-129.vtu
time of simulation =  3.225

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.21126]
Potential E =  [2.32508]
Total E =  1.11382
Relative Energy difference =  0.0001795556648407232

N_snapshot:  130
Reading snapshot: particles-130.vtu
time of simulation =  3.25

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.22074]
Potential E =  [2.33452]
Total E =  1.11378
Relative Energy difference =  0.00014369503733806366

N_snapshot:  131
Reading snapshot: particles-131.vtu
time of simulation =  3.275

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.22785]
Potential E =  [2.34158]
Total E =  1.11373
Relative Energy difference =  9.873544464816214e-05

N_snapshot:  132
Reading snapshot: particles-132.vtu
time of simulation =  3.3000000000000003

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.23235]
Potential E =  [2.34603]
Total E =  1.11367
Relative Energy difference =  4.489098006954201e-05

N_snapshot:  133
Reading snapshot: particles-133.vtu
time of simulation =  3.325

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.23409]
Potential E =  [2.34771]
Total E =  1.11362
Relative Energy difference =  3.8434028902163246e-08

N_snapshot:  134
Reading snapshot: particles-134.vtu
time of simulation =  3.35

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.23301]
Potential E =  [2.34657]
Total E =  1.11356
Relative Energy difference =  -5.3913077198979655e-05

N_snapshot:  135
Reading snapshot: particles-135.vtu
time of simulation =  3.375

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.22915]
Potential E =  [2.34265]
Total E =  1.11351
Relative Energy difference =  -9.876562323961949e-05

N_snapshot:  136
Reading snapshot: particles-136.vtu
time of simulation =  3.4

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.22263]
Potential E =  [2.33608]
Total E =  1.11345
Relative Energy difference =  -0.00015261008781823965

N_snapshot:  137
Reading snapshot: particles-137.vtu
time of simulation =  3.4250000000000003

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.21367]
Potential E =  [2.32708]
Total E =  1.11341
Relative Energy difference =  -0.00018857776197016084

N_snapshot:  138
Reading snapshot: particles-138.vtu
time of simulation =  3.45

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.20256]
Potential E =  [2.31593]
Total E =  1.11337
Relative Energy difference =  -0.00022454543612208206

N_snapshot:  139
Reading snapshot: particles-139.vtu
time of simulation =  3.475

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.18964]
Potential E =  [2.30297]
Total E =  1.11334
Relative Energy difference =  -0.0002514141450867613

N_snapshot:  140
Reading snapshot: particles-140.vtu
time of simulation =  3.5

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.17528]
Potential E =  [2.28859]
Total E =  1.11331
Relative Energy difference =  -0.0002783899007007022

N_snapshot:  141
Reading snapshot: particles-141.vtu
time of simulation =  3.525

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.15986]
Potential E =  [2.27316]
Total E =  1.1133
Relative Energy difference =  -0.0002873818192386825

N_snapshot:  142
Reading snapshot: particles-142.vtu
time of simulation =  3.5500000000000003

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.1438]
Potential E =  [2.25708]
Total E =  1.11329
Relative Energy difference =  -0.0002963737377766628

N_snapshot:  143
Reading snapshot: particles-143.vtu
time of simulation =  3.575

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.12745]
Potential E =  [2.24073]
Total E =  1.11328
Relative Energy difference =  -0.0003052586096653814

N_snapshot:  144
Reading snapshot: particles-144.vtu
time of simulation =  3.6

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.11117]
Potential E =  [2.22446]
Total E =  1.11329
Relative Energy difference =  -0.0002963737377766628

N_snapshot:  145
Reading snapshot: particles-145.vtu
time of simulation =  3.625

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.09529]
Potential E =  [2.20859]
Total E =  1.1133
Relative Energy difference =  -0.0002873818192386825

N_snapshot:  146
Reading snapshot: particles-146.vtu
time of simulation =  3.65

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.0801]
Potential E =  [2.19342]
Total E =  1.11332
Relative Energy difference =  -0.0002693979821627219

N_snapshot:  147
Reading snapshot: particles-147.vtu
time of simulation =  3.6750000000000003

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.06584]
Potential E =  [2.17918]
Total E =  1.11334
Relative Energy difference =  -0.0002514141450867613

N_snapshot:  148
Reading snapshot: particles-148.vtu
time of simulation =  3.7

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.05274]
Potential E =  [2.1661]
Total E =  1.11336
Relative Energy difference =  -0.0002334303080108007

N_snapshot:  149
Reading snapshot: particles-149.vtu
time of simulation =  3.725

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.04097]
Potential E =  [2.15436]
Total E =  1.11339
Relative Energy difference =  -0.00020656159904612145

N_snapshot:  150
Reading snapshot: particles-150.vtu
time of simulation =  3.75

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.03069]
Potential E =  [2.14411]
Total E =  1.11342
Relative Energy difference =  -0.00017958584343218055

N_snapshot:  151
Reading snapshot: particles-151.vtu
time of simulation =  3.775

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.02201]
Potential E =  [2.13547]
Total E =  1.11346
Relative Energy difference =  -0.00014372521592952102

N_snapshot:  152
Reading snapshot: particles-152.vtu
time of simulation =  3.8000000000000003

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.01504]
Potential E =  [2.12853]
Total E =  1.11349
Relative Energy difference =  -0.0001167494603155801

N_snapshot:  153
Reading snapshot: particles-153.vtu
time of simulation =  3.825

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00985]
Potential E =  [2.12337]
Total E =  1.11353
Relative Energy difference =  -8.078178616365888e-05

N_snapshot:  154
Reading snapshot: particles-154.vtu
time of simulation =  3.85

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00648]
Potential E =  [2.12004]
Total E =  1.11356
Relative Energy difference =  -5.3913077198979655e-05

N_snapshot:  155
Reading snapshot: particles-155.vtu
time of simulation =  3.875

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00497]
Potential E =  [2.11857]
Total E =  1.1136
Relative Energy difference =  -1.794540304705844e-05

N_snapshot:  156
Reading snapshot: particles-156.vtu
time of simulation =  3.9

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00533]
Potential E =  [2.11897]
Total E =  1.11364
Relative Energy difference =  1.79152244556011e-05

N_snapshot:  157
Reading snapshot: particles-157.vtu
time of simulation =  3.9250000000000003

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00756]
Potential E =  [2.12124]
Total E =  1.11368
Relative Energy difference =  5.3882898607522306e-05

N_snapshot:  158
Reading snapshot: particles-158.vtu
time of simulation =  3.95

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.01164]
Potential E =  [2.12536]
Total E =  1.11371
Relative Energy difference =  8.085865422146322e-05

N_snapshot:  159
Reading snapshot: particles-159.vtu
time of simulation =  3.975

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.01754]
Potential E =  [2.13129]
Total E =  1.11375
Relative Energy difference =  0.00011671928172412275

N_snapshot:  160
Reading snapshot: particles-160.vtu
time of simulation =  4.0

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.02518]
Potential E =  [2.13896]
Total E =  1.11378
Relative Energy difference =  0.00014369503733806366

N_snapshot:  161
Reading snapshot: particles-161.vtu
time of simulation =  4.025

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.03449]
Potential E =  [2.1483]
Total E =  1.11382
Relative Energy difference =  0.0001795556648407232

N_snapshot:  162
Reading snapshot: particles-162.vtu
time of simulation =  4.05

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.04536]
Potential E =  [2.15921]
Total E =  1.11385
Relative Energy difference =  0.00020653142045466412

N_snapshot:  163
Reading snapshot: particles-163.vtu
time of simulation =  4.075

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.05767]
Potential E =  [2.17154]
Total E =  1.11387
Relative Energy difference =  0.00022451525753062473

N_snapshot:  164
Reading snapshot: particles-164.vtu
time of simulation =  4.1

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.07124]
Potential E =  [2.18514]
Total E =  1.1139
Relative Energy difference =  0.00025138396649530397

N_snapshot:  165
Reading snapshot: particles-165.vtu
time of simulation =  4.125

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.08589]
Potential E =  [2.19981]
Total E =  1.11391
Relative Energy difference =  0.00026037588503328426

N_snapshot:  166
Reading snapshot: particles-166.vtu
time of simulation =  4.15

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.10138]
Potential E =  [2.21531]
Total E =  1.11393
Relative Energy difference =  0.00027835972210924485

N_snapshot:  167
Reading snapshot: particles-167.vtu
time of simulation =  4.175

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.11745]
Potential E =  [2.23139]
Total E =  1.11394
Relative Energy difference =  0.00028735164064722514

N_snapshot:  168
Reading snapshot: particles-168.vtu
time of simulation =  4.2

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.13379]
Potential E =  [2.24773]
Total E =  1.11394
Relative Energy difference =  0.00028735164064722514

N_snapshot:  169
Reading snapshot: particles-169.vtu
time of simulation =  4.225

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.15008]
Potential E =  [2.26401]
Total E =  1.11393
Relative Energy difference =  0.00027835972210924485

N_snapshot:  170
Reading snapshot: particles-170.vtu
time of simulation =  4.25

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.16594]
Potential E =  [2.27986]
Total E =  1.11392
Relative Energy difference =  0.00026936780357126455

N_snapshot:  171
Reading snapshot: particles-171.vtu
time of simulation =  4.275

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.18099]
Potential E =  [2.29489]
Total E =  1.1139
Relative Energy difference =  0.00025138396649530397

N_snapshot:  172
Reading snapshot: particles-172.vtu
time of simulation =  4.3

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.19484]
Potential E =  [2.30871]
Total E =  1.11387
Relative Energy difference =  0.00022451525753062473

N_snapshot:  173
Reading snapshot: particles-173.vtu
time of simulation =  4.325

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.2071]
Potential E =  [2.32094]
Total E =  1.11384
Relative Energy difference =  0.0001975395019166838

N_snapshot:  174
Reading snapshot: particles-174.vtu
time of simulation =  4.3500000000000005

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.21741]
Potential E =  [2.33121]
Total E =  1.11379
Relative Energy difference =  0.00015268695587604398

N_snapshot:  175
Reading snapshot: particles-175.vtu
time of simulation =  4.375

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.22546]
Potential E =  [2.3392]
Total E =  1.11375
Relative Energy difference =  0.00011671928172412275

N_snapshot:  176
Reading snapshot: particles-176.vtu
time of simulation =  4.4

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.23097]
Potential E =  [2.34466]
Total E =  1.11369
Relative Energy difference =  6.287481714550261e-05

N_snapshot:  177
Reading snapshot: particles-177.vtu
time of simulation =  4.425

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.23377]
Potential E =  [2.3474]
Total E =  1.11364
Relative Energy difference =  1.79152244556011e-05

N_snapshot:  178
Reading snapshot: particles-178.vtu
time of simulation =  4.45

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.23375]
Potential E =  [2.34733]
Total E =  1.11358
Relative Energy difference =  -3.5929240123019044e-05

N_snapshot:  179
Reading snapshot: particles-179.vtu
time of simulation =  4.4750000000000005

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.23093]
Potential E =  [2.34446]
Total E =  1.11353
Relative Energy difference =  -8.078178616365888e-05

N_snapshot:  180
Reading snapshot: particles-180.vtu
time of simulation =  4.5

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.2254]
Potential E =  [2.33887]
Total E =  1.11347
Relative Energy difference =  -0.0001347332973915407

N_snapshot:  181
Reading snapshot: particles-181.vtu
time of simulation =  4.525

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.21733]
Potential E =  [2.33076]
Total E =  1.11343
Relative Energy difference =  -0.00017059392489420023

N_snapshot:  182
Reading snapshot: particles-182.vtu
time of simulation =  4.55

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.207]
Potential E =  [2.32039]
Total E =  1.11338
Relative Energy difference =  -0.00021555351758410175

N_snapshot:  183
Reading snapshot: particles-183.vtu
time of simulation =  4.575

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.19472]
Potential E =  [2.30807]
Total E =  1.11335
Relative Energy difference =  -0.000242422226548781

N_snapshot:  184
Reading snapshot: particles-184.vtu
time of simulation =  4.6000000000000005



<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.18086]
Potential E =  [2.29418]
Total E =  1.11332
Relative Energy difference =  -0.0002693979821627219

N_snapshot:  185
Reading snapshot: particles-185.vtu
time of simulation =  4.625

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.1658]
Potential E =  [2.2791]
Total E =  1.1133
Relative Energy difference =  -0.0002873818192386825

N_snapshot:  186
Reading snapshot: particles-186.vtu
time of simulation =  4.65

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.14994]
Potential E =  [2.26323]
Total E =  1.11329
Relative Energy difference =  -0.0002963737377766628

N_snapshot:  187
Reading snapshot: particles-187.vtu
time of simulation =  4.675

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.13365]
Potential E =  [2.24694]
Total E =  1.11328
Relative Energy difference =  -0.0003052586096653814

N_snapshot:  188
Reading snapshot: particles-188.vtu
time of simulation =  4.7

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.11731]
Potential E =  [2.2306]
Total E =  1.11329
Relative Energy difference =  -0.0002963737377766628

N_snapshot:  189
Reading snapshot: particles-189.vtu
time of simulation =  4.7250000000000005

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.10124]
Potential E =  [2.21454]
Total E =  1.1133
Relative Energy difference =  -0.0002873818192386825

N_snapshot:  190
Reading snapshot: particles-190.vtu
time of simulation =  4.75

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.08576]
Potential E =  [2.19907]
Total E =  1.11331
Relative Energy difference =  -0.0002783899007007022

N_snapshot:  191
Reading snapshot: particles-191.vtu
time of simulation =  4.775

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.07112]
Potential E =  [2.18445]
Total E =  1.11333
Relative Energy difference =  -0.0002604060636247416

N_snapshot:  192
Reading snapshot: particles-192.vtu
time of simulation =  4.8

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.05756]
Potential E =  [2.17091]
Total E =  1.11335
Relative Energy difference =  -0.000242422226548781

N_snapshot:  193
Reading snapshot: particles-193.vtu
time of simulation =  4.825

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.04526]
Potential E =  [2.15864]
Total E =  1.11338
Relative Energy difference =  -0.00021555351758410175

N_snapshot:  194
Reading snapshot: particles-194.vtu
time of simulation =  4.8500000000000005

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.0344]
Potential E =  [2.14781]
Total E =  1.11341
Relative Energy difference =  -0.00018857776197016084

N_snapshot:  195
Reading snapshot: particles-195.vtu
time of simulation =  4.875

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.0251]
Potential E =  [2.13855]
Total E =  1.11344
Relative Energy difference =  -0.00016160200635621994

N_snapshot:  196
Reading snapshot: particles-196.vtu
time of simulation =  4.9

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.01748]
Potential E =  [2.13096]
Total E =  1.11348
Relative Energy difference =  -0.0001257413788535604

N_snapshot:  197
Reading snapshot: particles-197.vtu
time of simulation =  4.925

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.0116]
Potential E =  [2.12512]
Total E =  1.11351
Relative Energy difference =  -9.876562323961949e-05

N_snapshot:  198
Reading snapshot: particles-198.vtu
time of simulation =  4.95

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00754]
Potential E =  [2.12109]
Total E =  1.11355
Relative Energy difference =  -6.290499573695996e-05

N_snapshot:  199
Reading snapshot: particles-199.vtu
time of simulation =  4.9750000000000005

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00532]
Potential E =  [2.11891]
Total E =  1.11359
Relative Energy difference =  -2.6937321585038743e-05

N_snapshot:  200
Reading snapshot: particles-200.vtu
time of simulation =  5.0

<bound method Mesh.point_data_to_sets of <meshio mesh object>
  Number of points: 1
  Number of cells:
    vertex: 1
  Point data: x, cut-off-radius, association, v, energyKin, energyPot, energyTot>

Measure energy conservation:
Initial E =  1.11362
Kinetic E =  [1.00495]
Potential E =  [2.11857]
Total E =  1.11362
Relative Energy difference =  3.8434028902163246e-08

Done reading all snapshots. Total=  200

# Plot radial profile
fig = plt.figure(figsize=(6, 4)) 
gs = gridspec.GridSpec(1,1) 
 
ax0 = plt.subplot(gs[0])
plt.plot(t_all, E_dif_all, c='r', lw=0.8, zorder=1)
plt.xlabel(r'$t$')
plt.ylabel(r'$\Delta E/E_{\rm ini}$')
plt.margins(0.05, 0.1)
plt.tick_params(direction='in', top=True, right=True)
plt.axhline(0, c='gray', ls=':', lw=1.2, zorder=-49 )
plt.xlim(0, t_all.max())
plt.ylim(-0.005, 0.005 )
 
plt.show()
#fig.savefig('energy_conservation_profile.pdf')
plt.close()
 
# Check energy conservation
E_tolerance = 1e-2
 
E_min = E_dif_all.min()
E_max = E_dif_all.max()
E_mean= np.mean(E_dif_all)
 
assert abs(E_min) < E_tolerance, f"Energy not conserved. E_min: {E_min}"
assert abs(E_mean)< E_tolerance, f"Energy not conserved. E_mean: {E_min}"
assert abs(E_max) < E_tolerance, f"Energy not conserved. E_max: {E_max}"
 
# @TODO test if there is any long-term drift, e.g. of the mean.
 
print('Test passed!')

png

Test passed!

E_dif_all.mean()

-2.7212127908516467e-05

E_dif_all.min()

-0.0003052586096653814

E_dif_all.max()

0.00028735164064722514

mesh = meshio.read( "particles.pvd")

---

ReadError Traceback (most recent call last)

Input In [39], in <cell line: 1>() -—> 1 mesh = meshio.read( "particles.pvd")

File ~/.local/lib/python3.10/site-packages/meshio/_helpers.py:71, in read(filename, file_format)
     68 if is_buffer(filename, "r"):
     69     return _read_buffer(filename, file_format)
---> 71 return _read_file(Path(filename), file_format)


File ~/.local/lib/python3.10/site-packages/meshio/_helpers.py:96, in _read_file(path, file_format)
     93     possible_file_formats = [file_format]
     94 else:
     95     # deduce possible file formats from extension
---> 96     possible_file_formats = _filetypes_from_path(path)
     98 for file_format in possible_file_formats:
     99     if file_format not in reader_map:


File ~/.local/lib/python3.10/site-packages/meshio/_helpers.py:56, in _filetypes_from_path(path)
     53         pass
     55 if not out:
---> 56     raise ReadError(f"Could not deduce file format from path '{path}'.")
     57 return out


ReadError: Could not deduce file format from path 'particles.pvd'.