81cdde8899
EQUATIONS: Check if the state is physical in nr_iterate_srhd_adi_1dw().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-17 13:19:56 -02:00
989495ef97
EQUATIONS: Return correct info flag in nr_initial_brackets_srmhd_adi().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-17 11:41:02 -02:00
1a27369635
EQUATIONS: Slightly rewrite nr_initial_brackets_srmhd_adi().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-17 11:31:43 -02:00
d4ecfcccf2
EQUATIONS: Slightly rewrite nr_iterate_srmhd_adi_1dw().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-17 11:10:34 -02:00
45b703ad34
EQUATIONS: Adjust strings to the left in correct_unphysical_states().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-17 10:57:43 -02:00
a5400dbfb8
EQUATIONS: Improve the warning messages in correct_unphysical_states().
...
Also, set the minimum number of physical cells for averaging to 3.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-17 08:41:57 -02:00
7da950f652
EQUATIONS: Fix a bug in correct_unphysical_states().
...
The condition for taking the surrounding region was supposed to be based
on the number of physical cells available. Unfortunately, the counter
for unphysical cell was checked instead. Fix it.
Also increase the maximum region to four cells away.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-17 08:24:53 -02:00
fe117327d1
EQUATIONS: Rewrite nr_iterate_srhd_adi_1dw().
...
Make sure all cases are considered. Also, in the case of profiling, the
subroutine counts timing correctly.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-16 11:44:08 -02:00
c7a4159af7
EQUATIONS: Fix alignment in cons2prim_srmhd_adi().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-16 10:34:55 -02:00
096b7bc25e
EQUATIONS: Do not quit if solution not found in cons2prim_srmhd_adi().
...
Just skip the cell which failed to be converted from conservative to
primitive variables and hope it can be fixed later by averaging.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-16 10:17:05 -02:00
6047622bd4
EQUATIONS: Do not quit if solution not found in cons2prim_srhd_adi().
...
Just skip the cell which failed to be converted from conservative to
primitive variables and hope it can be fixed later by averaging.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-16 10:10:10 -02:00
f638d53064
EQUATIONS: Remove positivity flag since it is not generic.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-16 10:01:40 -02:00
fedba83d17
EQUATIONS: Implement correction of unphysical cells.
...
The unphysical cells (those with negative density or pressure) will be
corrected by averaging their states from the values of physical
neighbors.
The option is off by default, but can be turned on by setting:
fix_unphysical_cells = "on"
in the parameter file.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-16 09:57:23 -02:00
ded1fd55b1
ALGEBRA: Invert roots order in quadratic().
...
Fixes Special Relativity HLLC solver.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-09 11:07:45 -02:00
253ec8d77d
INTERPOLATIONS: Implement PPM with new limiters by Colella & Sekora.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-08 21:58:41 -02:00
63de17cd21
PROBLEMS: Make perturbation in RT problem perfectly symmetric.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-08 12:27:53 -02:00
68e5aea50e
DRIVER: Print the progress info no more than once per second.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-05 11:36:36 -02:00
38c857c851
Remove all references to PATHSCALE compiler.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-05 11:15:35 -02:00
0e6c338fc3
PROBLEMS: Update Rayleigh-Taylor problem.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-05 11:11:27 -02:00
1f8ce476d4
IO: Store code version in HDF5 files.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-05 09:39:10 -02:00
2c32139d8c
Update copyright info.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2018-01-04 12:13:09 -02:00
198e6db77c
IO: Store an attribute indicating the domain periodicity.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-12-30 14:25:49 -02:00
88434201bf
IO: Store 'gamma' for adiabatic, and 'csnd' for isothermal EoS.
...
Depending on the equation of state, adiabatic or isothermal, we
store the adiabatic index 'gamma' or the isothermal sound speed
'csnd' in the attributes of HDF5 files.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-12-30 13:17:25 -02:00
745faa3742
IO: Store attribute 'code' in the root of HDF5 file.
...
This attribute should help to determine the format of HDF5 files more
easy.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-12-30 13:11:34 -02:00
2feb69e6b4
PROBLEMS: Implement Sedov-Taylor blast problem.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-12-30 13:01:44 -02:00
1a04260167
IO: Adjust default compression levels.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-09-05 12:13:39 -03:00
8696f15696
IO: Clean up old compression code.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-09-05 12:09:03 -03:00
bb2e4c9baa
IO: Initialize ZSTD compression.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-09-05 12:02:07 -03:00
5208df0c60
IO: Set deflate compression and its level only once.
...
Also, remove SZIP compression.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-09-05 11:49:22 -03:00
9ccd050794
IO: Introduce flags for supported compression types.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-09-05 11:33:56 -03:00
f3649f10c1
IO: Move compression property creation to initialize_io().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-09-05 11:25:31 -03:00
bb09798398
IO: One local pid definition skipped.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-09-05 11:14:39 -03:00
1d070a9132
IO: Move the property ID to module variables.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-09-05 11:13:36 -03:00
060e9b397c
IO: Initialize/close HDF5 Fortran interface only once.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-09-05 11:08:54 -03:00
a3ce8c409e
IO: Add finalize_io() subroutine.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-09-05 09:53:47 -03:00
32129f4948
IO: Print information about the used HDF5 compression.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-09-05 08:58:13 -03:00
44aec4a280
IO: Add compression type detection.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-09-05 08:52:33 -03:00
9d75c49059
SCHEMES: Cosmetic changes in the adiabatic HLLD solver.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-05-16 17:03:01 -03:00
4fd495bf80
USER_PROBLEM: Remove unnecessary printed line.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-05-11 10:02:37 -03:00
9ed6f2ad1c
USER_PROBLEM: Remove printing of empty line.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-05-11 09:57:25 -03:00
5a10b9384d
DRIVER: Print problem header.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-05-11 09:54:12 -03:00
f6bb80b32d
INTERPOLATIONS: Apply MLP to all ghost cells.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-05-10 08:45:50 -03:00
ba82978b5b
INTERPOLATIONS: Rewrite high order 1D GP method.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-05-05 08:41:26 -03:00
97b10008f0
INTERPOLATIONS: Implement 7th order MP reconstruction.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-05-05 07:40:35 -03:00
b43dbe50ee
INTERPOLATIONS: Rewrite 5th order MP method.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-05-05 07:31:57 -03:00
c507f5b2ad
INTERPOLATIONS: Rewrite 5th order Compact MP method.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-05-03 11:57:36 -03:00
39981f5667
INTERPOLATIONS: Fix di coefficients in CRMP5LD method.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-05-03 11:53:38 -03:00
b862c8f1af
INTERPOLATIONS: Rewrite 5th order Compact Low Dissipation MP method.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-05-03 11:50:27 -03:00
4321e40657
INTERPOLATIONS: Implement 7th order Compact MP reconstruction.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-05-03 11:35:31 -03:00
17000df285
INTERPOLATIONS: Implement subroutine for MP limiting.
...
The monotonicity preserving (MP) limiter can be applied to various
reconstruction methods, therefore it is good to put it in a separate
subroutine.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2017-05-03 11:32:28 -03:00
