diff --git a/kalman/compute_errors_kalman.py b/kalman/compute_errors_kalman.py index 17b4d62..5e48288 100644 --- a/kalman/compute_errors_kalman.py +++ b/kalman/compute_errors_kalman.py @@ -6,8 +6,7 @@ def compute_relative_error(vcur,vref): for i in range(len(vref)): eps += ((vcur[i] - vref[i])/vref[i])**2 - return np.sqrt(eps) - + return np.sqrt(eps)/len(vcur) @@ -15,57 +14,117 @@ def compute_relative_error(vcur,vref): vref = [10, 250, 250, 250, 30] #case 1 -v_no_case1 = [ 10.05 , 248.34, 250.52 , 254.43 , 29.82] +v_120_case1 = [ 10.22 , 249.1 , 254.37 , 259.15 , 29.81 ] v_70_case1 = [20.38, 211.19, 189.01, 207.09, 28.78] v_50_case1 = [90.45, 240.98, 241.03, 316.89, 21.27] +v_25_case1 = [2012.53, 1004.3, 1172.01, 1755.02, 12.69] + +vnew_120_case1 = [ 10.26 , 249.84, 254.6, 259.69 , 29.81] vnew_70_case1 = [ 10.26 , 250.73, 252.78, 256.62 , 29.77] vnew_50_case1 = [ 10.27 , 248.41 ,254.71 ,253.01 , 29.82] +vnew_25_case1 = [ 10.18 , 251.13 ,253.78 ,257.4 , 29.82] + #case 2 -v_no_case2 = [ 9.77, 249.7 , 251.94, 254.91 , 29.82] -v_70_case2 = [ 20.7 , 214.77 ,196.47 ,209.31 , 28.76] -v_50_case2 = [ 88.4, 236.73 ,240.39 ,312.85, 21.22] +v_120_case2 = [ 9.95, 248.5, 254.4, 253.4, 29.80] +v_70_case2 = [ 20.7 , 214.77 , 196.47 , 209.31 , 28.76] +v_50_case2 = [ 88.40, 236.7, 240.4, 312.9, 21.22] +v_25_case2 = [ 880.4, 747.2, 822.7, 1380, 9.37] -vnew_70_case2 = [ 10 , 250.24 ,252.05 ,254.98 , 29.76] -vnew_50_case2 = [ 9.92 ,247.61 ,253.35 ,250.64 ,29.81] +vnew_120_case2 = [ 10, 246.6, 250.5, 254.5, 29.80] +vnew_70_case2 = [ 10.00, 250.2, 252.1, 255.0, 29.76] +vnew_50_case2 = [ 9.92, 247.6, 253.4, 250.6, 29.81] +vnew_25_case2 = [ 10.04 , 247.95 , 250.51 , 253.49 , 29.8 ] -eps_ref_case1 = compute_relative_error(v_no_case1,vref) + +#case 3 +v_120_case3 = [ 9.91 , 236.69 , 236.44 , 224.36 ,30.15] +v_70_case3 = [ 16.48 , 192.8 , 155.69 , 148.34 , 29.01] +v_50_case3 = [231.07 , 369.12 , 521.28 , 458.69 , 23.66] +v_25_case3 = [ 0, 0, 0, 0, 0] + +vnew_120_case3 = [ 10.17 , 247.95, 249.76 , 248.29 ,29.91] +vnew_70_case3 = [ 9.98 , 249.29 , 248.97 , 251.16 ,29.8 ] +vnew_50_case3 = [ 9.98 , 246.89 , 251.51 , 245.14 , 29.89] +vnew_25_case3 = [ 9.6 , 244.93 , 244.42 , 245.22 ,29.84] + + +####################################################################### + +eps_120_case1 = compute_relative_error(v_120_case1,vref) eps_70_case1 = compute_relative_error(v_70_case1,vref) -eps_new70_case1 = compute_relative_error(vnew_70_case1,vref) eps_50_case1 = compute_relative_error(v_50_case1,vref) +eps_25_case1 = compute_relative_error(v_25_case1,vref) + +eps_new120_case1 = compute_relative_error(vnew_120_case1,vref) +eps_new70_case1 = compute_relative_error(vnew_70_case1,vref) eps_new50_case1 = compute_relative_error(vnew_50_case1,vref) +eps_new25_case1 = compute_relative_error(vnew_25_case1,vref) -eps_ref_case2 = compute_relative_error(v_no_case2,vref) +eps_120_case2 = compute_relative_error(v_120_case2,vref) eps_70_case2 = compute_relative_error(v_70_case2,vref) -eps_new70_case2 = compute_relative_error(vnew_70_case2,vref) eps_50_case2 = compute_relative_error(v_50_case2,vref) +eps_25_case2 = compute_relative_error(v_25_case2,vref) + +eps_new120_case2 = compute_relative_error(vnew_120_case2,vref) +eps_new70_case2 = compute_relative_error(vnew_70_case2,vref) eps_new50_case2 = compute_relative_error(vnew_50_case2,vref) +eps_new25_case2 = compute_relative_error(vnew_25_case2,vref) +eps_120_case3 = compute_relative_error(v_120_case3,vref) +eps_70_case3 = compute_relative_error(v_70_case3,vref) +eps_50_case3 = compute_relative_error(v_50_case3,vref) +eps_25_case3 = compute_relative_error(v_25_case3,vref) +eps_new120_case3 = compute_relative_error(vnew_120_case3,vref) +eps_new70_case3 = compute_relative_error(vnew_70_case3,vref) +eps_new50_case3 = compute_relative_error(vnew_50_case3,vref) +eps_new25_case3 = compute_relative_error(vnew_25_case3,vref) print('--------- CASE 1 ---------') -print('reference error = {e}'.format(e=np.round(100*eps_ref_case1,2))) +print('venc 120% error = {e}'.format(e=np.round(100*eps_120_case1,2))) print('venc 70% error = {e}'.format(e=np.round(100*eps_70_case1,2))) -print('odv venc 70% error = {e}'.format(e=np.round(100*eps_new70_case1,2))) print('venc 50% error = {e}'.format(e=np.round(100*eps_50_case1,2))) +print('venc 25% error = {e}'.format(e=np.round(100*eps_25_case1,2))) + +print('odv venc 120% error = {e}'.format(e=np.round(100*eps_new120_case1,2))) +print('odv venc 70% error = {e}'.format(e=np.round(100*eps_new70_case1,2))) print('odv venc 50% error = {e}'.format(e=np.round(100*eps_new50_case1,2))) +print('odv venc 25% error = {e}'.format(e=np.round(100*eps_new25_case1,2))) print('--------- CASE 2 ---------') -print('reference error = {e}'.format(e=np.round(100*eps_ref_case2,2))) -print('venc 70% error = {e}'.format(e=np.round(100*eps_70_case2,2))) -print('odv venc 70% error = {e}'.format(e=np.round(100*eps_new70_case2,2))) -print('venc 50% error = {e}'.format(e=np.round(100*eps_50_case2,2))) -print('odv venc 50% error = {e}'.format(e=np.round(100*eps_new50_case2,2))) +print('venc 120% error = {e}'.format(e=np.round(100*eps_120_case2,2))) +print('venc 70% error = {e}'.format(e=np.round(100*eps_70_case2,2))) +print('venc 50% error = {e}'.format(e=np.round(100*eps_50_case2,2))) +print('venc 25% error = {e}'.format(e=np.round(100*eps_25_case2,2))) + +print('odv venc 120% error = {e}'.format(e=np.round(100*eps_new120_case2,2))) +print('odv venc 70% error = {e}'.format(e=np.round(100*eps_new70_case2,2))) +print('odv venc 50% error = {e}'.format(e=np.round(100*eps_new50_case2,2))) +print('odv venc 25% error = {e}'.format(e=np.round(100*eps_new25_case2,2))) + + +print('--------- CASE 3 ---------') + +print('venc 120% error = {e}'.format(e=np.round(100*eps_120_case3,2))) +print('venc 70% error = {e}'.format(e=np.round(100*eps_70_case3,2))) +print('venc 50% error = {e}'.format(e=np.round(100*eps_50_case3,2))) +print('venc 25% error = {e}'.format(e=np.round(100*eps_25_case3,2))) + +print('odv venc 120% error = {e}'.format(e=np.round(100*eps_new120_case3,2))) +print('odv venc 70% error = {e}'.format(e=np.round(100*eps_new70_case3,2))) +print('odv venc 50% error = {e}'.format(e=np.round(100*eps_new50_case3,2))) +print('odv venc 25% error = {e}'.format(e=np.round(100*eps_new25_case3,2))) \ No newline at end of file diff --git a/kalman/input_files/aorta.yaml b/kalman/input_files/aorta.yaml index 1a162a2..6e36f84 100755 --- a/kalman/input_files/aorta.yaml +++ b/kalman/input_files/aorta.yaml @@ -25,28 +25,28 @@ boundary_conditions: type: 'dirichlet' value: ['0','0','-U*sin(DOLFIN_PI*t/Th)*(t<=Th) + (Th + + + + + + + u.h5:/Mesh/0/mesh/topology + + + u.h5:/Mesh/0/mesh/geometry + + + + + u.h5:/Mesh/1/mesh/topology + + + u.h5:/Mesh/1/mesh/geometry + + + + + u.h5:/Mesh/2/mesh/topology + + + u.h5:/Mesh/2/mesh/geometry + + + + + u.h5:/Mesh/3/mesh/topology + + + u.h5:/Mesh/3/mesh/geometry + + + + + u.h5:/Mesh/4/mesh/topology + + + u.h5:/Mesh/4/mesh/geometry + + + + + u.h5:/Mesh/5/mesh/topology + + + u.h5:/Mesh/5/mesh/geometry + + + + + u.h5:/Mesh/6/mesh/topology + + + u.h5:/Mesh/6/mesh/geometry + + + + + u.h5:/Mesh/7/mesh/topology + + + u.h5:/Mesh/7/mesh/geometry + + + + + u.h5:/Mesh/8/mesh/topology + + + u.h5:/Mesh/8/mesh/geometry + + + + + u.h5:/Mesh/9/mesh/topology + + + u.h5:/Mesh/9/mesh/geometry + + + + + u.h5:/Mesh/10/mesh/topology + + + u.h5:/Mesh/10/mesh/geometry + + + + + u.h5:/Mesh/11/mesh/topology + + + u.h5:/Mesh/11/mesh/geometry + + + + + u.h5:/Mesh/12/mesh/topology + + + u.h5:/Mesh/12/mesh/geometry + + + + + u.h5:/Mesh/13/mesh/topology + + + u.h5:/Mesh/13/mesh/geometry + + + + + u.h5:/Mesh/14/mesh/topology + + + u.h5:/Mesh/14/mesh/geometry + + + + + u.h5:/Mesh/15/mesh/topology + + + u.h5:/Mesh/15/mesh/geometry + + + + + u.h5:/Mesh/16/mesh/topology + + + u.h5:/Mesh/16/mesh/geometry + + + + + u.h5:/Mesh/17/mesh/topology + + + u.h5:/Mesh/17/mesh/geometry + + + + + u.h5:/Mesh/18/mesh/topology + + + u.h5:/Mesh/18/mesh/geometry + + + + + u.h5:/Mesh/19/mesh/topology + + + u.h5:/Mesh/19/mesh/geometry + + + + + u.h5:/Mesh/20/mesh/topology + + + u.h5:/Mesh/20/mesh/geometry + + + + + diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u0.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u0.h5 new file mode 100644 index 0000000..42da1e1 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u0.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u12.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u12.h5 new file mode 100644 index 0000000..d2c5664 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u12.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u16.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u16.h5 new file mode 100644 index 0000000..64141aa Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u16.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u20.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u20.h5 new file mode 100644 index 0000000..605daa2 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u20.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u24.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u24.h5 new file mode 100644 index 0000000..9176246 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u24.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u28.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u28.h5 new file mode 100644 index 0000000..865e0b3 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u28.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u32.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u32.h5 new file mode 100644 index 0000000..fb0beb2 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u32.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u36.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u36.h5 new file mode 100644 index 0000000..67cfc08 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u36.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u4.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u4.h5 new file mode 100644 index 0000000..61cc8b2 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u4.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u40.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u40.h5 new file mode 100644 index 0000000..f7b8c84 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u40.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u44.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u44.h5 new file mode 100644 index 0000000..2e74cf5 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u44.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u48.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u48.h5 new file mode 100644 index 0000000..e64604a Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u48.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u52.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u52.h5 new file mode 100644 index 0000000..2694ae2 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u52.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u56.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u56.h5 new file mode 100644 index 0000000..a962377 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u56.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u60.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u60.h5 new file mode 100644 index 0000000..2fa6b87 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u60.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u64.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u64.h5 new file mode 100644 index 0000000..0309b95 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u64.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u68.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u68.h5 new file mode 100644 index 0000000..0f471d5 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u68.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u72.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u72.h5 new file mode 100644 index 0000000..24e3040 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u72.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u76.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u76.h5 new file mode 100644 index 0000000..b965112 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u76.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u8.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u8.h5 new file mode 100644 index 0000000..1d7b0f3 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u8.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V120/u80.h5 b/kalman/measurements/aorta/Perturbation/Ks12V120/u80.h5 new file mode 100644 index 0000000..6f05976 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V120/u80.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u.h5 new file mode 100644 index 0000000..d8723db Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u.xdmf b/kalman/measurements/aorta/Perturbation/Ks12V25/u.xdmf new file mode 100644 index 0000000..e936dfb --- /dev/null +++ b/kalman/measurements/aorta/Perturbation/Ks12V25/u.xdmf @@ -0,0 +1,260 @@ + + + + + + + + u.h5:/Mesh/0/mesh/topology + + + u.h5:/Mesh/0/mesh/geometry + + + + + u.h5:/Mesh/1/mesh/topology + + + u.h5:/Mesh/1/mesh/geometry + + + + + u.h5:/Mesh/2/mesh/topology + + + u.h5:/Mesh/2/mesh/geometry + + + + + u.h5:/Mesh/3/mesh/topology + + + u.h5:/Mesh/3/mesh/geometry + + + + + u.h5:/Mesh/4/mesh/topology + + + u.h5:/Mesh/4/mesh/geometry + + + + + u.h5:/Mesh/5/mesh/topology + + + u.h5:/Mesh/5/mesh/geometry + + + + + u.h5:/Mesh/6/mesh/topology + + + u.h5:/Mesh/6/mesh/geometry + + + + + u.h5:/Mesh/7/mesh/topology + + + u.h5:/Mesh/7/mesh/geometry + + + + + u.h5:/Mesh/8/mesh/topology + + + u.h5:/Mesh/8/mesh/geometry + + + + + u.h5:/Mesh/9/mesh/topology + + + u.h5:/Mesh/9/mesh/geometry + + + + + u.h5:/Mesh/10/mesh/topology + + + u.h5:/Mesh/10/mesh/geometry + + + + + u.h5:/Mesh/11/mesh/topology + + + u.h5:/Mesh/11/mesh/geometry + + + + + u.h5:/Mesh/12/mesh/topology + + + u.h5:/Mesh/12/mesh/geometry + + + + + u.h5:/Mesh/13/mesh/topology + + + u.h5:/Mesh/13/mesh/geometry + + + + + u.h5:/Mesh/14/mesh/topology + + + u.h5:/Mesh/14/mesh/geometry + + + + + u.h5:/Mesh/15/mesh/topology + + + u.h5:/Mesh/15/mesh/geometry + + + + + u.h5:/Mesh/16/mesh/topology + + + u.h5:/Mesh/16/mesh/geometry + + + + + u.h5:/Mesh/17/mesh/topology + + + u.h5:/Mesh/17/mesh/geometry + + + + + u.h5:/Mesh/18/mesh/topology + + + u.h5:/Mesh/18/mesh/geometry + + + + + u.h5:/Mesh/19/mesh/topology + + + u.h5:/Mesh/19/mesh/geometry + + + + + u.h5:/Mesh/20/mesh/topology + + + u.h5:/Mesh/20/mesh/geometry + + + + + diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u0.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u0.h5 new file mode 100644 index 0000000..5152edf Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u0.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u12.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u12.h5 new file mode 100644 index 0000000..9a3d1de Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u12.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u16.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u16.h5 new file mode 100644 index 0000000..d76dfb5 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u16.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u20.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u20.h5 new file mode 100644 index 0000000..ec911e6 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u20.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u24.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u24.h5 new file mode 100644 index 0000000..12e619c Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u24.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u28.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u28.h5 new file mode 100644 index 0000000..f410a55 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u28.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u32.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u32.h5 new file mode 100644 index 0000000..b367620 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u32.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u36.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u36.h5 new file mode 100644 index 0000000..070eb57 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u36.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u4.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u4.h5 new file mode 100644 index 0000000..d1e42d7 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u4.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u40.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u40.h5 new file mode 100644 index 0000000..1ac19ae Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u40.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u44.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u44.h5 new file mode 100644 index 0000000..c4d7be6 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u44.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u48.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u48.h5 new file mode 100644 index 0000000..5f93652 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u48.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u52.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u52.h5 new file mode 100644 index 0000000..349716c Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u52.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u56.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u56.h5 new file mode 100644 index 0000000..b1a17fb Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u56.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u60.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u60.h5 new file mode 100644 index 0000000..2a6520b Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u60.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u64.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u64.h5 new file mode 100644 index 0000000..0f5e1d9 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u64.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u68.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u68.h5 new file mode 100644 index 0000000..00079c4 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u68.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u72.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u72.h5 new file mode 100644 index 0000000..a8c67cc Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u72.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u76.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u76.h5 new file mode 100644 index 0000000..e5eff87 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u76.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u8.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u8.h5 new file mode 100644 index 0000000..6e5392c Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u8.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V25/u80.h5 b/kalman/measurements/aorta/Perturbation/Ks12V25/u80.h5 new file mode 100644 index 0000000..179c50d Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V25/u80.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u.h5 new file mode 100644 index 0000000..c5e8491 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u.xdmf b/kalman/measurements/aorta/Perturbation/Ks12V50/u.xdmf new file mode 100644 index 0000000..e936dfb --- /dev/null +++ b/kalman/measurements/aorta/Perturbation/Ks12V50/u.xdmf @@ -0,0 +1,260 @@ + + + + + + + + u.h5:/Mesh/0/mesh/topology + + + u.h5:/Mesh/0/mesh/geometry + + + + + u.h5:/Mesh/1/mesh/topology + + + u.h5:/Mesh/1/mesh/geometry + + + + + u.h5:/Mesh/2/mesh/topology + + + u.h5:/Mesh/2/mesh/geometry + + + + + u.h5:/Mesh/3/mesh/topology + + + u.h5:/Mesh/3/mesh/geometry + + + + + u.h5:/Mesh/4/mesh/topology + + + u.h5:/Mesh/4/mesh/geometry + + + + + u.h5:/Mesh/5/mesh/topology + + + u.h5:/Mesh/5/mesh/geometry + + + + + u.h5:/Mesh/6/mesh/topology + + + u.h5:/Mesh/6/mesh/geometry + + + + + u.h5:/Mesh/7/mesh/topology + + + u.h5:/Mesh/7/mesh/geometry + + + + + u.h5:/Mesh/8/mesh/topology + + + u.h5:/Mesh/8/mesh/geometry + + + + + u.h5:/Mesh/9/mesh/topology + + + u.h5:/Mesh/9/mesh/geometry + + + + + u.h5:/Mesh/10/mesh/topology + + + u.h5:/Mesh/10/mesh/geometry + + + + + u.h5:/Mesh/11/mesh/topology + + + u.h5:/Mesh/11/mesh/geometry + + + + + u.h5:/Mesh/12/mesh/topology + + + u.h5:/Mesh/12/mesh/geometry + + + + + u.h5:/Mesh/13/mesh/topology + + + u.h5:/Mesh/13/mesh/geometry + + + + + u.h5:/Mesh/14/mesh/topology + + + u.h5:/Mesh/14/mesh/geometry + + + + + u.h5:/Mesh/15/mesh/topology + + + u.h5:/Mesh/15/mesh/geometry + + + + + u.h5:/Mesh/16/mesh/topology + + + u.h5:/Mesh/16/mesh/geometry + + + + + u.h5:/Mesh/17/mesh/topology + + + u.h5:/Mesh/17/mesh/geometry + + + + + u.h5:/Mesh/18/mesh/topology + + + u.h5:/Mesh/18/mesh/geometry + + + + + u.h5:/Mesh/19/mesh/topology + + + u.h5:/Mesh/19/mesh/geometry + + + + + u.h5:/Mesh/20/mesh/topology + + + u.h5:/Mesh/20/mesh/geometry + + + + + diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u0.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u0.h5 new file mode 100644 index 0000000..16f08e0 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u0.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u12.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u12.h5 new file mode 100644 index 0000000..97f8f9c Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u12.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u16.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u16.h5 new file mode 100644 index 0000000..db39a3f Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u16.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u20.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u20.h5 new file mode 100644 index 0000000..74f26c7 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u20.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u24.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u24.h5 new file mode 100644 index 0000000..32971e9 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u24.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u28.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u28.h5 new file mode 100644 index 0000000..158a0c9 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u28.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u32.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u32.h5 new file mode 100644 index 0000000..e483edd Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u32.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u36.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u36.h5 new file mode 100644 index 0000000..79271e9 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u36.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u4.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u4.h5 new file mode 100644 index 0000000..b86e3a7 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u4.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u40.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u40.h5 new file mode 100644 index 0000000..1dce6a2 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u40.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u44.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u44.h5 new file mode 100644 index 0000000..5479afb Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u44.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u48.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u48.h5 new file mode 100644 index 0000000..1177ab5 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u48.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u52.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u52.h5 new file mode 100644 index 0000000..f343235 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u52.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u56.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u56.h5 new file mode 100644 index 0000000..67f882d Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u56.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u60.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u60.h5 new file mode 100644 index 0000000..9b669e4 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u60.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u64.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u64.h5 new file mode 100644 index 0000000..8574a74 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u64.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u68.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u68.h5 new file mode 100644 index 0000000..fea206d Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u68.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u72.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u72.h5 new file mode 100644 index 0000000..ce33ef6 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u72.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u76.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u76.h5 new file mode 100644 index 0000000..e8a6bfe Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u76.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u8.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u8.h5 new file mode 100644 index 0000000..b5358ba Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u8.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V50/u80.h5 b/kalman/measurements/aorta/Perturbation/Ks12V50/u80.h5 new file mode 100644 index 0000000..f534bb2 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V50/u80.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u.h5 new file mode 100644 index 0000000..33a52fd Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u.xdmf b/kalman/measurements/aorta/Perturbation/Ks12V70/u.xdmf new file mode 100644 index 0000000..e936dfb --- /dev/null +++ b/kalman/measurements/aorta/Perturbation/Ks12V70/u.xdmf @@ -0,0 +1,260 @@ + + + + + + + + u.h5:/Mesh/0/mesh/topology + + + u.h5:/Mesh/0/mesh/geometry + + + + + u.h5:/Mesh/1/mesh/topology + + + u.h5:/Mesh/1/mesh/geometry + + + + + u.h5:/Mesh/2/mesh/topology + + + u.h5:/Mesh/2/mesh/geometry + + + + + u.h5:/Mesh/3/mesh/topology + + + u.h5:/Mesh/3/mesh/geometry + + + + + u.h5:/Mesh/4/mesh/topology + + + u.h5:/Mesh/4/mesh/geometry + + + + + u.h5:/Mesh/5/mesh/topology + + + u.h5:/Mesh/5/mesh/geometry + + + + + u.h5:/Mesh/6/mesh/topology + + + u.h5:/Mesh/6/mesh/geometry + + + + + u.h5:/Mesh/7/mesh/topology + + + u.h5:/Mesh/7/mesh/geometry + + + + + u.h5:/Mesh/8/mesh/topology + + + u.h5:/Mesh/8/mesh/geometry + + + + + u.h5:/Mesh/9/mesh/topology + + + u.h5:/Mesh/9/mesh/geometry + + + + + u.h5:/Mesh/10/mesh/topology + + + u.h5:/Mesh/10/mesh/geometry + + + + + u.h5:/Mesh/11/mesh/topology + + + u.h5:/Mesh/11/mesh/geometry + + + + + u.h5:/Mesh/12/mesh/topology + + + u.h5:/Mesh/12/mesh/geometry + + + + + u.h5:/Mesh/13/mesh/topology + + + u.h5:/Mesh/13/mesh/geometry + + + + + u.h5:/Mesh/14/mesh/topology + + + u.h5:/Mesh/14/mesh/geometry + + + + + u.h5:/Mesh/15/mesh/topology + + + u.h5:/Mesh/15/mesh/geometry + + + + + u.h5:/Mesh/16/mesh/topology + + + u.h5:/Mesh/16/mesh/geometry + + + + + u.h5:/Mesh/17/mesh/topology + + + u.h5:/Mesh/17/mesh/geometry + + + + + u.h5:/Mesh/18/mesh/topology + + + u.h5:/Mesh/18/mesh/geometry + + + + + u.h5:/Mesh/19/mesh/topology + + + u.h5:/Mesh/19/mesh/geometry + + + + + u.h5:/Mesh/20/mesh/topology + + + u.h5:/Mesh/20/mesh/geometry + + + + + diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u0.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u0.h5 new file mode 100644 index 0000000..83a4a9a Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u0.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u12.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u12.h5 new file mode 100644 index 0000000..aafed89 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u12.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u16.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u16.h5 new file mode 100644 index 0000000..c407fd5 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u16.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u20.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u20.h5 new file mode 100644 index 0000000..4ae6545 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u20.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u24.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u24.h5 new file mode 100644 index 0000000..2c8f2e2 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u24.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u28.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u28.h5 new file mode 100644 index 0000000..dfdc774 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u28.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u32.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u32.h5 new file mode 100644 index 0000000..361346f Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u32.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u36.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u36.h5 new file mode 100644 index 0000000..ec738ec Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u36.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u4.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u4.h5 new file mode 100644 index 0000000..152561e Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u4.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u40.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u40.h5 new file mode 100644 index 0000000..0cce170 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u40.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u44.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u44.h5 new file mode 100644 index 0000000..c7b1f21 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u44.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u48.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u48.h5 new file mode 100644 index 0000000..401a0fb Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u48.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u52.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u52.h5 new file mode 100644 index 0000000..713dc4c Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u52.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u56.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u56.h5 new file mode 100644 index 0000000..eb00b2e Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u56.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u60.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u60.h5 new file mode 100644 index 0000000..07cbf69 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u60.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u64.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u64.h5 new file mode 100644 index 0000000..54917b0 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u64.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u68.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u68.h5 new file mode 100644 index 0000000..0dfa150 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u68.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u72.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u72.h5 new file mode 100644 index 0000000..fa09d28 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u72.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u76.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u76.h5 new file mode 100644 index 0000000..08bf8d1 Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u76.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u8.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u8.h5 new file mode 100644 index 0000000..92cbd0f Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u8.h5 differ diff --git a/kalman/measurements/aorta/Perturbation/Ks12V70/u80.h5 b/kalman/measurements/aorta/Perturbation/Ks12V70/u80.h5 new file mode 100644 index 0000000..d6c42ab Binary files /dev/null and b/kalman/measurements/aorta/Perturbation/Ks12V70/u80.h5 differ diff --git a/kalman/measurements/aorta/input.yaml b/kalman/measurements/aorta/input.yaml new file mode 100755 index 0000000..e3002bb --- /dev/null +++ b/kalman/measurements/aorta/input.yaml @@ -0,0 +1,155 @@ + +mesh: './meshes/coaortaH1.h5' +# Physical parameters of the fluid +fluid: + density: 1.2 + dynamic_viscosity: 0.035 + +io: + write_path: 'results/aorta' + restart: + path: '' # './projects/nse_coa3d/results/test_restart2/' + time: 0 + write_xdmf: True + write_checkpoints: True + write_hdf5_timeseries: False + write_velocity: 'update' # update or tentative + +boundary_conditions: + - + id: 1 + type: 'dirichlet' + value: ['0','0','0'] + - + id: 2 + type: 'dirichlet' + value: ['0','0','-U*sin(DOLFIN_PI*t/Th)*(t<=Th) + (Th + + + + + + + u_all.h5:/Mesh/0/mesh/topology + + + u_all.h5:/Mesh/0/mesh/geometry + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +