arma-thesis

git clone https://git.igankevich.com/arma-thesis.git
Log | Files | Refs | LICENSE

commit a19c6c376f0e859bfb77dc445af7d945f01c8b00
parent 39d17aa6a15a95301edfa1fd31fd7047cc037673
Author: Ivan Gankevich <igankevich@ya.ru>
Date:   Wed,  5 Apr 2017 18:56:05 +0300

Update ARMA configuration files. Add configuration files for vel.

Diffstat:
arma-thesis.org | 364++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++-----------
config/parallel_mt.dat | 0
config/plain_wave_high_amplitude_solver.arma | 29+++++++++++++++++++++++++++++
config/plain_wave_linear_solver.arma | 29+++++++++++++++++++++++++++++
config/propagating_wave.arma | 16+++++++++++-----
config/standing_wave.arma | 12++++++++----
6 files changed, 393 insertions(+), 57 deletions(-)

diff --git a/arma-thesis.org b/arma-thesis.org @@ -27,32 +27,32 @@ #+PROPERTY: header-args:R :results graphics :exports results * Config :noexport: -** Produce data for Q-Q and ACF plots +** Produce data for Q-Q, ACF and velocity potential plots #+begin_src sh :exports none :results verbatim root=$(pwd) -for testname in propagating_wave standing_wave +for testname in \ + propagating_wave \ + standing_wave \ + plain_wave_linear_solver \ + plain_wave_high_amplitude_solver do wd=$root/build/$testname rm -rf $wd mkdir -p $wd cd $wd + cp $root/config/parallel_mt.dat . arma -c $root/config/$testname.arma 2>&1 done #+end_src #+RESULTS: #+begin_example -Input file = /home/igankevich/workspace/phd-diss/config/propagating_wave.arma -ACF grid size = (20,10,10) -ACF grid patch size = (0.526316,0.555556,0.555556) +Input file = /home/igankevich/workspace/arma-thesis/config/propagating_wave.arma Output grid size = (200,40,40) Output grid patch size = (1,1,1) -AR order = (10,10,10) -Do least squares = 0 -ACF function = propagating_wave Model = MA -MA algorithm = fixed_point_iteration Verification scheme = manual +MA model = order=(20,10,10),acf.shape=(20,10,10),algorithm=fixed_point_iteration ACF variance = 5 fixed_point_iteration:Iteration=0, var_wn=2.70831 fixed_point_iteration:Iteration=1, var_wn=1.93791 @@ -80,7 +80,7 @@ fixed_point_iteration:Iteration=22, var_wn=0.518666 fixed_point_iteration:Iteration=23, var_wn=0.509497 fixed_point_iteration:Iteration=24, var_wn=0.50103 fixed_point_iteration:Iteration=25, var_wn=0.493191 -fixed_point_iteration:Iteration=26, var_wn=0.485916 +fixed_point_iteration:Iteration=26, var_wn=0.485915 fixed_point_iteration:Iteration=27, var_wn=0.479148 fixed_point_iteration:Iteration=28, var_wn=0.472841 fixed_point_iteration:Iteration=29, var_wn=0.466951 @@ -93,23 +93,23 @@ fixed_point_iteration:Iteration=35, var_wn=0.43856 fixed_point_iteration:Iteration=36, var_wn=0.434752 fixed_point_iteration:Iteration=37, var_wn=0.431155 fixed_point_iteration:Iteration=38, var_wn=0.427755 -fixed_point_iteration:Iteration=39, var_wn=0.424538 +fixed_point_iteration:Iteration=39, var_wn=0.424537 fixed_point_iteration:Iteration=40, var_wn=0.42149 -fixed_point_iteration:Iteration=41, var_wn=0.418601 +fixed_point_iteration:Iteration=41, var_wn=0.4186 fixed_point_iteration:Iteration=42, var_wn=0.415859 fixed_point_iteration:Iteration=43, var_wn=0.413256 fixed_point_iteration:Iteration=44, var_wn=0.410782 -fixed_point_iteration:Iteration=45, var_wn=0.40843 +fixed_point_iteration:Iteration=45, var_wn=0.408429 fixed_point_iteration:Iteration=46, var_wn=0.406191 fixed_point_iteration:Iteration=47, var_wn=0.404059 -fixed_point_iteration:Iteration=48, var_wn=0.402029 +fixed_point_iteration:Iteration=48, var_wn=0.402028 fixed_point_iteration:Iteration=49, var_wn=0.400092 fixed_point_iteration:Iteration=50, var_wn=0.398246 fixed_point_iteration:Iteration=51, var_wn=0.396483 fixed_point_iteration:Iteration=52, var_wn=0.3948 fixed_point_iteration:Iteration=53, var_wn=0.393193 fixed_point_iteration:Iteration=54, var_wn=0.391656 -fixed_point_iteration:Iteration=55, var_wn=0.390188 +fixed_point_iteration:Iteration=55, var_wn=0.390187 fixed_point_iteration:Iteration=56, var_wn=0.388782 fixed_point_iteration:Iteration=57, var_wn=0.387438 fixed_point_iteration:Iteration=58, var_wn=0.386151 @@ -121,7 +121,7 @@ fixed_point_iteration:Iteration=63, var_wn=0.380482 fixed_point_iteration:Iteration=64, var_wn=0.379485 fixed_point_iteration:Iteration=65, var_wn=0.378528 fixed_point_iteration:Iteration=66, var_wn=0.37761 -fixed_point_iteration:Iteration=67, var_wn=0.376729 +fixed_point_iteration:Iteration=67, var_wn=0.376728 fixed_point_iteration:Iteration=68, var_wn=0.375882 fixed_point_iteration:Iteration=69, var_wn=0.37507 fixed_point_iteration:Iteration=70, var_wn=0.374289 @@ -129,17 +129,17 @@ fixed_point_iteration:Iteration=71, var_wn=0.373539 fixed_point_iteration:Iteration=72, var_wn=0.372818 fixed_point_iteration:Iteration=73, var_wn=0.372126 fixed_point_iteration:Iteration=74, var_wn=0.37146 -fixed_point_iteration:Iteration=75, var_wn=0.37082 +fixed_point_iteration:Iteration=75, var_wn=0.370819 fixed_point_iteration:Iteration=76, var_wn=0.370204 -fixed_point_iteration:Iteration=77, var_wn=0.369612 +fixed_point_iteration:Iteration=77, var_wn=0.369611 fixed_point_iteration:Iteration=78, var_wn=0.369042 -fixed_point_iteration:Iteration=79, var_wn=0.368494 +fixed_point_iteration:Iteration=79, var_wn=0.368493 fixed_point_iteration:Iteration=80, var_wn=0.367966 fixed_point_iteration:Iteration=81, var_wn=0.367458 fixed_point_iteration:Iteration=82, var_wn=0.366969 fixed_point_iteration:Iteration=83, var_wn=0.366499 fixed_point_iteration:Iteration=84, var_wn=0.366046 -fixed_point_iteration:Iteration=85, var_wn=0.36561 +fixed_point_iteration:Iteration=85, var_wn=0.365609 fixed_point_iteration:Iteration=86, var_wn=0.365189 fixed_point_iteration:Iteration=87, var_wn=0.364785 fixed_point_iteration:Iteration=88, var_wn=0.364395 @@ -152,8 +152,8 @@ fixed_point_iteration:Iteration=94, var_wn=0.362337 fixed_point_iteration:Iteration=95, var_wn=0.362036 fixed_point_iteration:Iteration=96, var_wn=0.361746 fixed_point_iteration:Iteration=97, var_wn=0.361466 -fixed_point_iteration:Iteration=98, var_wn=0.361197 -fixed_point_iteration:Iteration=99, var_wn=0.360937 +fixed_point_iteration:Iteration=98, var_wn=0.361196 +fixed_point_iteration:Iteration=99, var_wn=0.360936 fixed_point_iteration:Iteration=100, var_wn=0.360686 fixed_point_iteration:Iteration=101, var_wn=0.360444 fixed_point_iteration:Iteration=102, var_wn=0.360211 @@ -166,7 +166,7 @@ fixed_point_iteration:Iteration=108, var_wn=0.358975 fixed_point_iteration:Iteration=109, var_wn=0.358794 fixed_point_iteration:Iteration=110, var_wn=0.358619 fixed_point_iteration:Iteration=111, var_wn=0.35845 -fixed_point_iteration:Iteration=112, var_wn=0.358288 +fixed_point_iteration:Iteration=112, var_wn=0.358287 fixed_point_iteration:Iteration=113, var_wn=0.35813 fixed_point_iteration:Iteration=114, var_wn=0.357979 fixed_point_iteration:Iteration=115, var_wn=0.357832 @@ -177,7 +177,7 @@ fixed_point_iteration:Iteration=119, var_wn=0.357296 fixed_point_iteration:Iteration=120, var_wn=0.357173 fixed_point_iteration:Iteration=121, var_wn=0.357055 fixed_point_iteration:Iteration=122, var_wn=0.356941 -fixed_point_iteration:Iteration=123, var_wn=0.356831 +fixed_point_iteration:Iteration=123, var_wn=0.35683 fixed_point_iteration:Iteration=124, var_wn=0.356724 fixed_point_iteration:Iteration=125, var_wn=0.356621 fixed_point_iteration:Iteration=126, var_wn=0.356522 @@ -192,20 +192,20 @@ fixed_point_iteration:Iteration=134, var_wn=0.355842 fixed_point_iteration:Iteration=135, var_wn=0.355769 fixed_point_iteration:Iteration=136, var_wn=0.355699 fixed_point_iteration:Iteration=137, var_wn=0.355632 -fixed_point_iteration:Iteration=138, var_wn=0.355567 +fixed_point_iteration:Iteration=138, var_wn=0.355566 fixed_point_iteration:Iteration=139, var_wn=0.355504 fixed_point_iteration:Iteration=140, var_wn=0.355443 fixed_point_iteration:Iteration=141, var_wn=0.355384 fixed_point_iteration:Iteration=142, var_wn=0.355327 -fixed_point_iteration:Iteration=143, var_wn=0.355273 +fixed_point_iteration:Iteration=143, var_wn=0.355272 fixed_point_iteration:Iteration=144, var_wn=0.35522 -fixed_point_iteration:Iteration=145, var_wn=0.355169 +fixed_point_iteration:Iteration=145, var_wn=0.355168 fixed_point_iteration:Iteration=146, var_wn=0.355119 -fixed_point_iteration:Iteration=147, var_wn=0.355072 -fixed_point_iteration:Iteration=148, var_wn=0.355026 +fixed_point_iteration:Iteration=147, var_wn=0.355071 +fixed_point_iteration:Iteration=148, var_wn=0.355025 fixed_point_iteration:Iteration=149, var_wn=0.354981 fixed_point_iteration:Iteration=150, var_wn=0.354938 -fixed_point_iteration:Iteration=151, var_wn=0.354897 +fixed_point_iteration:Iteration=151, var_wn=0.354896 fixed_point_iteration:Iteration=152, var_wn=0.354856 fixed_point_iteration:Iteration=153, var_wn=0.354818 fixed_point_iteration:Iteration=154, var_wn=0.35478 @@ -216,11 +216,11 @@ fixed_point_iteration:Iteration=158, var_wn=0.354643 fixed_point_iteration:Iteration=159, var_wn=0.354612 fixed_point_iteration:Iteration=160, var_wn=0.354581 fixed_point_iteration:Iteration=161, var_wn=0.354552 -fixed_point_iteration:Iteration=162, var_wn=0.354524 +fixed_point_iteration:Iteration=162, var_wn=0.354523 fixed_point_iteration:Iteration=163, var_wn=0.354496 fixed_point_iteration:Iteration=164, var_wn=0.35447 fixed_point_iteration:Iteration=165, var_wn=0.354444 -fixed_point_iteration:Iteration=166, var_wn=0.35442 +fixed_point_iteration:Iteration=166, var_wn=0.354419 fixed_point_iteration:Iteration=167, var_wn=0.354396 fixed_point_iteration:Iteration=168, var_wn=0.354373 fixed_point_iteration:Iteration=169, var_wn=0.35435 @@ -231,15 +231,15 @@ fixed_point_iteration:Iteration=173, var_wn=0.354269 fixed_point_iteration:Iteration=174, var_wn=0.35425 fixed_point_iteration:Iteration=175, var_wn=0.354232 fixed_point_iteration:Iteration=176, var_wn=0.354214 -fixed_point_iteration:Iteration=177, var_wn=0.354198 +fixed_point_iteration:Iteration=177, var_wn=0.354197 fixed_point_iteration:Iteration=178, var_wn=0.354181 fixed_point_iteration:Iteration=179, var_wn=0.354165 fixed_point_iteration:Iteration=180, var_wn=0.35415 -fixed_point_iteration:Iteration=181, var_wn=0.354136 +fixed_point_iteration:Iteration=181, var_wn=0.354135 fixed_point_iteration:Iteration=182, var_wn=0.354121 -fixed_point_iteration:Iteration=183, var_wn=0.354108 +fixed_point_iteration:Iteration=183, var_wn=0.354107 fixed_point_iteration:Iteration=184, var_wn=0.354094 -fixed_point_iteration:Iteration=185, var_wn=0.354082 +fixed_point_iteration:Iteration=185, var_wn=0.354081 fixed_point_iteration:Iteration=186, var_wn=0.354069 fixed_point_iteration:Iteration=187, var_wn=0.354057 fixed_point_iteration:Iteration=188, var_wn=0.354046 @@ -247,23 +247,291 @@ fixed_point_iteration:Iteration=189, var_wn=0.354034 fixed_point_iteration:Iteration=190, var_wn=0.354024 fixed_point_iteration:Iteration=191, var_wn=0.354013 fixed_point_iteration:Iteration=192, var_wn=0.354003 -fixed_point_iteration:Iteration=193, var_wn=0.353994 -WN variance = 0.353994 -Input file = /home/igankevich/workspace/phd-diss/config/standing_wave.arma -ACF grid size = (10,10,10) -ACF grid patch size = (0.277778,0.555556,0.555556) +fixed_point_iteration:Iteration=193, var_wn=0.353993 +WN variance = 0.353993 +Partition size = (34,13,13) +Finished part [1/96] +Finished part [2/96] +Finished part [3/96] +Finished part [4/96] +Finished part [5/96] +Finished part [6/96] +Finished part [7/96] +Finished part [8/96] +Finished part [9/96] +Finished part [10/96] +Finished part [11/96] +Finished part [12/96] +Finished part [13/96] +Finished part [14/96] +Finished part [15/96] +Finished part [16/96] +Finished part [17/96] +Finished part [18/96] +Finished part [19/96] +Finished part [20/96] +Finished part [21/96] +Finished part [22/96] +Finished part [23/96] +Finished part [24/96] +Finished part [25/96] +Finished part [26/96] +Finished part [27/96] +Finished part [28/96] +Finished part [29/96] +Finished part [30/96] +Finished part [31/96] +Finished part [32/96] +Finished part [33/96] +Finished part [34/96] +Finished part [35/96] +Finished part [36/96] +Finished part [37/96] +Finished part [38/96] +Finished part [39/96] +Finished part [40/96] +Finished part [41/96] +Finished part [42/96] +Finished part [43/96] +Finished part [44/96] +Finished part [45/96] +Finished part [46/96] +Finished part [47/96] +Finished part [48/96] +Finished part [49/96] +Finished part [50/96] +Finished part [51/96] +Finished part [52/96] +Finished part [53/96] +Finished part [54/96] +Finished part [55/96] +Finished part [56/96] +Finished part [57/96] +Finished part [58/96] +Finished part [59/96] +Finished part [60/96] +Finished part [61/96] +Finished part [62/96] +Finished part [63/96] +Finished part [64/96] +Finished part [65/96] +Finished part [66/96] +Finished part [67/96] +Finished part [68/96] +Finished part [69/96] +Finished part [70/96] +Finished part [71/96] +Finished part [72/96] +Finished part [73/96] +Finished part [74/96] +Finished part [75/96] +Finished part [76/96] +Finished part [77/96] +Finished part [78/96] +Finished part [79/96] +Finished part [80/96] +Finished part [81/96] +Finished part [82/96] +Finished part [83/96] +Finished part [84/96] +Finished part [85/96] +Finished part [86/96] +Finished part [87/96] +Finished part [88/96] +Finished part [89/96] +Finished part [90/96] +Finished part [91/96] +Finished part [92/96] +Finished part [93/96] +Finished part [94/96] +Finished part [95/96] +Finished part [96/96] +Input file = /home/igankevich/workspace/arma-thesis/config/standing_wave.arma Output grid size = (200,40,40) Output grid patch size = (1,1,1) -AR order = (7,7,7) -Do least squares = 0 -ACF function = standing_wave Model = AR -MA algorithm = fixed_point_iteration Verification scheme = manual +AR model = order=(7,7,7),acf.shape=(10,10,10) ACF variance = 5 -WN variance = 0.00261323 -Zeta size = (193,33,33) -NaN: 29, -nan, 1.798e+36, -1.04284e+38, inf, -1.798e+36, -1.798e+36 +WN variance = 0.00259511 +Partition size = (21,10,10) +Finished part [1/160] +Finished part [2/160] +Finished part [3/160] +Finished part [4/160] +Finished part [5/160] +Finished part [6/160] +Finished part [7/160] +Finished part [8/160] +Finished part [9/160] +Finished part [10/160] +Finished part [11/160] +Finished part [12/160] +Finished part [13/160] +Finished part [14/160] +Finished part [15/160] +Finished part [16/160] +Finished part [17/160] +Finished part [18/160] +Finished part [19/160] +Finished part [20/160] +Finished part [21/160] +Finished part [22/160] +Finished part [23/160] +Finished part [24/160] +Finished part [25/160] +Finished part [26/160] +Finished part [27/160] +Finished part [28/160] +Finished part [29/160] +Finished part [30/160] +Finished part [31/160] +Finished part [32/160] +Finished part [33/160] +Finished part [34/160] +Finished part [35/160] +Finished part [36/160] +Finished part [37/160] +Finished part [38/160] +Finished part [39/160] +Finished part [40/160] +Finished part [41/160] +Finished part [42/160] +Finished part [43/160] +Finished part [44/160] +Finished part [45/160] +Finished part [46/160] +Finished part [47/160] +Finished part [48/160] +Finished part [49/160] +Finished part [50/160] +Finished part [51/160] +Finished part [52/160] +Finished part [53/160] +Finished part [54/160] +Finished part [55/160] +Finished part [56/160] +Finished part [57/160] +Finished part [58/160] +Finished part [59/160] +Finished part [60/160] +Finished part [61/160] +Finished part [62/160] +Finished part [63/160] +Finished part [64/160] +Finished part [65/160] +Finished part [66/160] +Finished part [67/160] +Finished part [68/160] +Finished part [69/160] +Finished part [70/160] +Finished part [71/160] +Finished part [72/160] +Finished part [73/160] +Finished part [74/160] +Finished part [75/160] +Finished part [76/160] +Finished part [77/160] +Finished part [78/160] +Finished part [79/160] +Finished part [80/160] +Finished part [81/160] +Finished part [82/160] +Finished part [83/160] +Finished part [84/160] +Finished part [85/160] +Finished part [86/160] +Finished part [87/160] +Finished part [88/160] +Finished part [89/160] +Finished part [90/160] +Finished part [91/160] +Finished part [92/160] +Finished part [93/160] +Finished part [94/160] +Finished part [95/160] +Finished part [96/160] +Finished part [97/160] +Finished part [98/160] +Finished part [99/160] +Finished part [100/160] +Finished part [101/160] +Finished part [102/160] +Finished part [103/160] +Finished part [104/160] +Finished part [105/160] +Finished part [106/160] +Finished part [107/160] +Finished part [108/160] +Finished part [109/160] +Finished part [110/160] +Finished part [111/160] +Finished part [112/160] +Finished part [113/160] +Finished part [114/160] +Finished part [115/160] +Finished part [116/160] +Finished part [117/160] +Finished part [118/160] +Finished part [119/160] +Finished part [120/160] +Finished part [121/160] +Finished part [122/160] +Finished part [123/160] +Finished part [124/160] +Finished part [125/160] +Finished part [126/160] +Finished part [127/160] +Finished part [128/160] +Finished part [129/160] +Finished part [130/160] +Finished part [131/160] +Finished part [132/160] +Finished part [133/160] +Finished part [134/160] +Finished part [135/160] +Finished part [136/160] +Finished part [137/160] +Finished part [138/160] +Finished part [139/160] +Finished part [140/160] +Finished part [141/160] +Finished part [142/160] +Finished part [143/160] +Finished part [144/160] +Finished part [145/160] +Finished part [146/160] +Finished part [147/160] +Finished part [148/160] +Finished part [149/160] +Finished part [150/160] +Finished part [151/160] +Finished part [152/160] +Finished part [153/160] +Finished part [154/160] +Finished part [155/160] +Finished part [156/160] +Finished part [157/160] +Finished part [158/160] +Finished part [159/160] +Finished part [160/160] +Zeta size = (192,32,32) +Input file = /home/igankevich/workspace/arma-thesis/config/plain_wave_linear_solver.arma +Output grid size = (200,128,40) +Output grid patch size = (1,0.0629921,0.205128) +Model = plain_wave +Verification scheme = manual +Plain wave model = func=cos,amplitudes=[4],wavenumbers=[0.125],phases=[0],velocities=[0.5] +Velocity potential solver name = N4arma8velocity13Linear_solverIdEE +Velocity potential solver = wnmax=(0,0.25),depth=12,domain=from (10,-12) to (10,4) npoints (1,128) +Input file = /home/igankevich/workspace/arma-thesis/config/plain_wave_high_amplitude_solver.arma +Output grid size = (200,128,40) +Output grid patch size = (1,0.0629921,0.205128) +Model = plain_wave +Verification scheme = manual +Plain wave model = func=cos,amplitudes=[4],wavenumbers=[0.125],phases=[0],velocities=[0.5] +Velocity potential solver name = N4arma8velocity21High_amplitude_solverIdEE +Velocity potential solver = wnmax=(0,0.25),depth=12,domain=from (10,-12) to (10,4) npoints (1,128) #+end_example * Introduction diff --git a/config/parallel_mt.dat b/config/parallel_mt.dat Binary files differ. diff --git a/config/plain_wave_high_amplitude_solver.arma b/config/plain_wave_high_amplitude_solver.arma @@ -0,0 +1,29 @@ +out_grid = (200,128,40) : (199,8,8) + +model = plain_wave +plain_wave = { + func = cos + amplitudes = [4] + wavenumbers = [0.125] + phases = [0] + velocities = [0.5] +} +verification = manual + +velocity_potential_solver = high_amplitude { + wnmax = (0,0.25) + depth = 12 + domain = from (10,-12) to (10,4) npoints (1,128) +} + +#velocity_potential_solver = plain { +# waves = { +# func = cos +# amplitudes = [4] +# wavenumbers = [0.125] +# phases = [0] +# velocities = [0.5] +# } +# depth = 12 +# domain = from (10,-12) to (10,4) npoints (1,128) +#} diff --git a/config/plain_wave_linear_solver.arma b/config/plain_wave_linear_solver.arma @@ -0,0 +1,29 @@ +out_grid = (200,128,40) : (199,8,8) + +model = plain_wave +plain_wave = { + func = cos + amplitudes = [4] + wavenumbers = [0.125] + phases = [0] + velocities = [0.5] +} +verification = manual + +velocity_potential_solver = linear { + wnmax = (0,0.25) + depth = 12 + domain = from (10,-12) to (10,4) npoints (1,128) +} + +#velocity_potential_solver = plain { +# waves = { +# func = cos +# amplitudes = [4] +# wavenumbers = [0.125] +# phases = [0] +# velocities = [0.5] +# } +# depth = 12 +# domain = from (10,-12) to (10,4) npoints (1,128) +#} diff --git a/config/propagating_wave.arma b/config/propagating_wave.arma @@ -1,9 +1,15 @@ out_grid = (200,40,40) -acf_grid = (20,10,10) : (10,5,5) -ma_order = (20,10,10) -least_squares = 0 -acf = propagating_wave model = MA -ma_algorithm = fixed_point_iteration +ma_model = { + acf = { + func = propagating_wave + grid = (20,10,10) : (10,5,5) + } + algorithm = fixed_point_iteration + order = (20,10,10) + max_iterations = 1000 + eps = 1e-5 + min_var_wn = 1e-6 +} verification = manual diff --git a/config/standing_wave.arma b/config/standing_wave.arma @@ -1,8 +1,12 @@ out_grid = (200,40,40) -acf_grid = (10,10,10) : (2.5,5,5) -ar_order = (7,7,7) -least_squares = 0 -acf = standing_wave model = AR verification = manual +ar_model = { + acf = { + func = standing_wave + grid = (10,10,10) : (2.5,5,5) + } + least_squares = 0 + order = (7,7,7) +}