diff --git a/kalman/graphics/figure2.py b/kalman/graphics/figure2.py
index 8d84cae..de42690 100644
--- a/kalman/graphics/figure2.py
+++ b/kalman/graphics/figure2.py
@@ -51,12 +51,12 @@ def plot_parameters(dat, input_file, deparameterize=False, ref=None):
with open(inputfile_path) as file:
inputfile = yaml.full_load(file)
- #true_val = [10,250,250,250,30]
+
true_values = {
- 3: 10,
- 4: 60,
- 5: 220,
- 6: 160,
+ 3: 2400,
+ 4: 4200,
+ 5: 11000,
+ 6: 7800,
2: 100
}
@@ -71,7 +71,7 @@ def plot_parameters(dat, input_file, deparameterize=False, ref=None):
for bnd_set in inputfile['boundary_conditions']:
if bnd_c['id'] == bnd_set['id']:
ids.append(bnd_c['id'])
- current_val.append(bnd_set['value'][0])
+ current_val.append(bnd_set['parameters']['R_d'])
labels.append('$R_' + str(bnd_c['id']))
elif 'dirichlet' in bnd_c['type']:
diff --git a/kalman/input_files/aorta.yaml b/kalman/input_files/aorta.yaml
index 8badcf1..e007e9c 100755
--- a/kalman/input_files/aorta.yaml
+++ b/kalman/input_files/aorta.yaml
@@ -4,10 +4,10 @@ mesh: './meshes/coaortaH1.h5'
fluid:
density: 1.2
dynamic_viscosity: 0.035
- stokes: False
+ stokes: True
io:
- write_path: 'results/aorta'
+ write_path: 'results/aorta_1ms'
restart:
path: '' # './projects/nse_coa3d/results/test_restart2/'
time: 0
@@ -32,34 +32,39 @@ boundary_conditions:
-
id: 3
type: 'windkessel'
- #value: [10,0,0]
- value: [10,0.0008,2400]
- p0: [80,1333.223874]
- #p0: [0,1333.223874]
+ parameters:
+ R_p: 10
+ C: 0.0008
+ R_d: 2400
+ p0: 80
+ conv: 1333.223874
-
id: 4
type: 'windkessel'
- #value: [60,0,0]
- #value: [250,0,0]
- value: [60,0.00034,4200]
- p0: [80,1333.223874]
- #p0: [0,1333.223874]
+ parameters:
+ R_p: 60
+ C: 0.00034
+ R_d: 4200
+ p0: 80
+ conv: 1333.223874
-
id: 5
type: 'windkessel'
- #value: [220,0,0]
- #value: [250,0,0]
- value: [220,0.00034,11000]
- p0: [80,1333.223874]
- #p0: [0,1333.223874]
+ parameters:
+ R_p: 220
+ C: 0.00034
+ R_d: 11000
+ p0: 80
+ conv: 1333.223874
-
id: 6
type: 'windkessel'
- #value: [160,0,0]
- #value: [250,0,0]
- value: [160,0.00034,7800]
- p0: [80,1333.223874]
- #p0: [0,1333.223874]
+ parameters:
+ R_p: 160
+ C: 0.00034
+ R_d: 7800
+ p0: 80
+ conv: 1333.223874
timemarching:
velocity_pressure_coupling: 'fractionalstep' # monolithic, fractionalstep
@@ -148,13 +153,14 @@ estimation:
measurements:
-
- mesh: '/home/yeye/NuMRI/kalman/meshes/coaortaH3_leo2.0.h5'
+ #mesh: '/home/yeye/NuMRI/kalman/meshes/coaortaH3_leo2.0.h5'
+ mesh: './meshes/coaortaH1.h5'
fe_degree: 1
xdmf_file: 'measurements/aorta/u_all.xdmf'
file_root: 'measurements/aorta/u{i}.h5'
indices: 0 # indices of checkpoints to be processed. 0 == all
velocity_direction: ~
- noise_stddev: 500 # standard deviation of Gaussian noise
+ noise_stddev: 0 # standard deviation of Gaussian noise
roukf:
particles: 'simplex' # unique or simplex
diff --git a/poi/checkpoint/0/p.h5 b/poi/checkpoint/0/p.h5
new file mode 100644
index 0000000..e8434a9
Binary files /dev/null and b/poi/checkpoint/0/p.h5 differ
diff --git a/poi/checkpoint/0/u.h5 b/poi/checkpoint/0/u.h5
new file mode 100644
index 0000000..17f9535
Binary files /dev/null and b/poi/checkpoint/0/u.h5 differ
diff --git a/poi/checkpoint/19/p.h5 b/poi/checkpoint/19/p.h5
new file mode 100644
index 0000000..f32cd18
Binary files /dev/null and b/poi/checkpoint/19/p.h5 differ
diff --git a/poi/checkpoint/19/u.h5 b/poi/checkpoint/19/u.h5
new file mode 100644
index 0000000..50d4be8
Binary files /dev/null and b/poi/checkpoint/19/u.h5 differ
diff --git a/poi/checkpoint/29/p.h5 b/poi/checkpoint/29/p.h5
new file mode 100644
index 0000000..f14587b
Binary files /dev/null and b/poi/checkpoint/29/p.h5 differ
diff --git a/poi/checkpoint/29/u.h5 b/poi/checkpoint/29/u.h5
new file mode 100644
index 0000000..0de0dc6
Binary files /dev/null and b/poi/checkpoint/29/u.h5 differ
diff --git a/poi/checkpoint/39/p.h5 b/poi/checkpoint/39/p.h5
new file mode 100644
index 0000000..63a628a
Binary files /dev/null and b/poi/checkpoint/39/p.h5 differ
diff --git a/poi/checkpoint/39/u.h5 b/poi/checkpoint/39/u.h5
new file mode 100644
index 0000000..8250e13
Binary files /dev/null and b/poi/checkpoint/39/u.h5 differ
diff --git a/poi/checkpoint/49/p.h5 b/poi/checkpoint/49/p.h5
new file mode 100644
index 0000000..5b96b90
Binary files /dev/null and b/poi/checkpoint/49/p.h5 differ
diff --git a/poi/checkpoint/49/u.h5 b/poi/checkpoint/49/u.h5
new file mode 100644
index 0000000..25d2258
Binary files /dev/null and b/poi/checkpoint/49/u.h5 differ
diff --git a/poi/checkpoint/59/p.h5 b/poi/checkpoint/59/p.h5
new file mode 100644
index 0000000..99b98e7
Binary files /dev/null and b/poi/checkpoint/59/p.h5 differ
diff --git a/poi/checkpoint/59/u.h5 b/poi/checkpoint/59/u.h5
new file mode 100644
index 0000000..13294d0
Binary files /dev/null and b/poi/checkpoint/59/u.h5 differ
diff --git a/poi/checkpoint/69/p.h5 b/poi/checkpoint/69/p.h5
new file mode 100644
index 0000000..444eb0b
Binary files /dev/null and b/poi/checkpoint/69/p.h5 differ
diff --git a/poi/checkpoint/69/u.h5 b/poi/checkpoint/69/u.h5
new file mode 100644
index 0000000..072cfbe
Binary files /dev/null and b/poi/checkpoint/69/u.h5 differ
diff --git a/poi/checkpoint/79/p.h5 b/poi/checkpoint/79/p.h5
new file mode 100644
index 0000000..fd05f8d
Binary files /dev/null and b/poi/checkpoint/79/p.h5 differ
diff --git a/poi/checkpoint/79/u.h5 b/poi/checkpoint/79/u.h5
new file mode 100644
index 0000000..7c22bc9
Binary files /dev/null and b/poi/checkpoint/79/u.h5 differ
diff --git a/poi/checkpoint/89/p.h5 b/poi/checkpoint/89/p.h5
new file mode 100644
index 0000000..5704218
Binary files /dev/null and b/poi/checkpoint/89/p.h5 differ
diff --git a/poi/checkpoint/89/u.h5 b/poi/checkpoint/89/u.h5
new file mode 100644
index 0000000..2087579
Binary files /dev/null and b/poi/checkpoint/89/u.h5 differ
diff --git a/poi/checkpoint/9/p.h5 b/poi/checkpoint/9/p.h5
new file mode 100644
index 0000000..2c643a4
Binary files /dev/null and b/poi/checkpoint/9/p.h5 differ
diff --git a/poi/checkpoint/9/u.h5 b/poi/checkpoint/9/u.h5
new file mode 100644
index 0000000..977bd2e
Binary files /dev/null and b/poi/checkpoint/9/u.h5 differ
diff --git a/poi/checkpoint/99/p.h5 b/poi/checkpoint/99/p.h5
new file mode 100644
index 0000000..4de5831
Binary files /dev/null and b/poi/checkpoint/99/p.h5 differ
diff --git a/poi/checkpoint/99/u.h5 b/poi/checkpoint/99/u.h5
new file mode 100644
index 0000000..524d40f
Binary files /dev/null and b/poi/checkpoint/99/u.h5 differ
diff --git a/poi/p.h5 b/poi/p.h5
new file mode 100644
index 0000000..c56d6e9
Binary files /dev/null and b/poi/p.h5 differ
diff --git a/poi/p.xdmf b/poi/p.xdmf
new file mode 100644
index 0000000..e6410ee
--- /dev/null
+++ b/poi/p.xdmf
@@ -0,0 +1,140 @@
+
+
+
+
+
+
+
+ p.h5:/Mesh/0/mesh/topology
+
+
+ p.h5:/Mesh/0/mesh/geometry
+
+
+
+ p.h5:/VisualisationVector/0
+
+
+
+
+ p.h5:/Mesh/1/mesh/topology
+
+
+ p.h5:/Mesh/1/mesh/geometry
+
+
+
+ p.h5:/VisualisationVector/1
+
+
+
+
+ p.h5:/Mesh/2/mesh/topology
+
+
+ p.h5:/Mesh/2/mesh/geometry
+
+
+
+ p.h5:/VisualisationVector/2
+
+
+
+
+ p.h5:/Mesh/3/mesh/topology
+
+
+ p.h5:/Mesh/3/mesh/geometry
+
+
+
+ p.h5:/VisualisationVector/3
+
+
+
+
+ p.h5:/Mesh/4/mesh/topology
+
+
+ p.h5:/Mesh/4/mesh/geometry
+
+
+
+ p.h5:/VisualisationVector/4
+
+
+
+
+ p.h5:/Mesh/5/mesh/topology
+
+
+ p.h5:/Mesh/5/mesh/geometry
+
+
+
+ p.h5:/VisualisationVector/5
+
+
+
+
+ p.h5:/Mesh/6/mesh/topology
+
+
+ p.h5:/Mesh/6/mesh/geometry
+
+
+
+ p.h5:/VisualisationVector/6
+
+
+
+
+ p.h5:/Mesh/7/mesh/topology
+
+
+ p.h5:/Mesh/7/mesh/geometry
+
+
+
+ p.h5:/VisualisationVector/7
+
+
+
+
+ p.h5:/Mesh/8/mesh/topology
+
+
+ p.h5:/Mesh/8/mesh/geometry
+
+
+
+ p.h5:/VisualisationVector/8
+
+
+
+
+ p.h5:/Mesh/9/mesh/topology
+
+
+ p.h5:/Mesh/9/mesh/geometry
+
+
+
+ p.h5:/VisualisationVector/9
+
+
+
+
+ p.h5:/Mesh/10/mesh/topology
+
+
+ p.h5:/Mesh/10/mesh/geometry
+
+
+
+ p.h5:/VisualisationVector/10
+
+
+
+
+
diff --git a/poi/poiseuille.h5 b/poi/poiseuille.h5
new file mode 100644
index 0000000..ef53a1f
Binary files /dev/null and b/poi/poiseuille.h5 differ
diff --git a/poi/seq.npz b/poi/seq.npz
new file mode 100644
index 0000000..59ff599
Binary files /dev/null and b/poi/seq.npz differ
diff --git a/poi/seq_new.npz b/poi/seq_new.npz
new file mode 100644
index 0000000..1b74d47
Binary files /dev/null and b/poi/seq_new.npz differ
diff --git a/poi/u.h5 b/poi/u.h5
new file mode 100644
index 0000000..c51f7e9
Binary files /dev/null and b/poi/u.h5 differ
diff --git a/poi/u.xdmf b/poi/u.xdmf
new file mode 100644
index 0000000..4d65a18
--- /dev/null
+++ b/poi/u.xdmf
@@ -0,0 +1,140 @@
+
+
+
+
+
+
+
+ u.h5:/Mesh/0/mesh/topology
+
+
+ u.h5:/Mesh/0/mesh/geometry
+
+
+
+ u.h5:/VisualisationVector/0
+
+
+
+
+ u.h5:/Mesh/1/mesh/topology
+
+
+ u.h5:/Mesh/1/mesh/geometry
+
+
+
+ u.h5:/VisualisationVector/1
+
+
+
+
+ u.h5:/Mesh/2/mesh/topology
+
+
+ u.h5:/Mesh/2/mesh/geometry
+
+
+
+ u.h5:/VisualisationVector/2
+
+
+
+
+ u.h5:/Mesh/3/mesh/topology
+
+
+ u.h5:/Mesh/3/mesh/geometry
+
+
+
+ u.h5:/VisualisationVector/3
+
+
+
+
+ u.h5:/Mesh/4/mesh/topology
+
+
+ u.h5:/Mesh/4/mesh/geometry
+
+
+
+ u.h5:/VisualisationVector/4
+
+
+
+
+ u.h5:/Mesh/5/mesh/topology
+
+
+ u.h5:/Mesh/5/mesh/geometry
+
+
+
+ u.h5:/VisualisationVector/5
+
+
+
+
+ u.h5:/Mesh/6/mesh/topology
+
+
+ u.h5:/Mesh/6/mesh/geometry
+
+
+
+ u.h5:/VisualisationVector/6
+
+
+
+
+ u.h5:/Mesh/7/mesh/topology
+
+
+ u.h5:/Mesh/7/mesh/geometry
+
+
+
+ u.h5:/VisualisationVector/7
+
+
+
+
+ u.h5:/Mesh/8/mesh/topology
+
+
+ u.h5:/Mesh/8/mesh/geometry
+
+
+
+ u.h5:/VisualisationVector/8
+
+
+
+
+ u.h5:/Mesh/9/mesh/topology
+
+
+ u.h5:/Mesh/9/mesh/geometry
+
+
+
+ u.h5:/VisualisationVector/9
+
+
+
+
+ u.h5:/Mesh/10/mesh/topology
+
+
+ u.h5:/Mesh/10/mesh/geometry
+
+
+
+ u.h5:/VisualisationVector/10
+
+
+
+
+