1909 Commits

Author SHA1 Message Date
3b352c89ee INTERPOLATIONS, DRIVER: Correct nghosts in initialize_interpolations().
The number of ghosts zone is important for the selected recontruction.
For higher order reconstruction methods, nghosts >= 4 usually.
Therefore, initialize interpolations before coordinates and correct the
number of ghost zone according to the selected interpolation method.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-11 11:27:22 -02:00
601adc8037 IO, DRIVER: Allow for restart with different number of ghost zones.
Now we can restart from restart snapshots with different number of ghost
zones. Also, remove the old restart snapshot format support from
read_datablocks_h5().

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-11 11:02:36 -02:00
6d421e6efb IO: Update the initialization/finalization status flag.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-11 10:01:35 -02:00
db1d3c6d4f SCHEMES: Make initialize_schemes() resistant to wrong parameters.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-11 09:51:02 -02:00
5e7b0962bc MESH: Update the initialization/finalization status flag.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-11 09:39:06 -02:00
cfe9614db2 MPITOOLS: Update the initialization/finalization status flag.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-11 09:36:02 -02:00
1b2251f7fc OPERATORS: Update the initialization/finalization status flag.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-11 09:21:55 -02:00
f84250d25a PROBLEMS: Update the initialization/finalization status flag.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-11 09:20:31 -02:00
2badc8c062 SHAPES: Update the initialization/finalization status flag.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-11 09:18:34 -02:00
f3229559e3 USER_PROBLEM: Update the initialization/finalization status flag.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-11 09:16:54 -02:00
ffdba19bb6 INTERPOLATIONS: Make initialize_interpolations() resistant to wrong parameters.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-11 09:14:51 -02:00
40f15e3117 INTEGRALS: Update the initialization/finalization status flag.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-11 08:38:15 -02:00
bfbc972e24 GRAVITY: Update the initialization/finalization status flag.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-11 08:35:38 -02:00
a772a9cf55 DOMAINS: Update the initialization/finalization status flag.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-11 08:34:01 -02:00
7eb75ba08c BOUNDARIES: Update the initialization/finalization status flag.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-11 08:32:27 -02:00
e6a1206919 BLOCKS: Update the initialization/finalization status flag.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-11 08:28:36 -02:00
1e3dae9add REFINEMENT: Make initialize_refinement() resistant to wrong parameters.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-10 21:47:45 -02:00
38d52cd6cd COORDINATES: Make initialize_coordinates() resistant to wrong parameters.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-08 16:30:53 -02:00
b2009f03be SOURCES: Make initialize_sources() resistant to wrong parameters.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-08 15:37:06 -02:00
2b131cf849 EQUATIONS: Make initialize_equations() resistant to wrong parameters.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-08 15:28:23 -02:00
603622886b EVOLUTION: Fix intermediate times in evolve_ssprk3_m().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-08 13:33:11 -02:00
3c649ec71f EVOLUTION: Fix intermediate time calculation in evolve_ssprk4_10().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-08 12:52:30 -02:00
d5c351e0ac EVOLUTION: dh has to be update each time in evolve_ssprk3_m().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-08 12:30:22 -02:00
7d50915b42 EVOLUTION: Add optimal m-stage 3rd order SSPRK method.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-08 12:03:10 -02:00
1273c3d0fd TIMERS: Print error messages to error unit.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-08 10:10:01 -02:00
dc8c1aaaa7 DRIVER: Do not print timings in the case of initialization errors.
If there was any problem with the initialization, there is no point in
printing the timings. Just leave the last error message and quit the
program.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-08 09:56:51 -02:00
a7fae8eb2d EVOLUTION: Move du(:,:,:,:) initialization to initialize_evolution().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-08 09:39:17 -02:00
d8dee86237 EVOLUTION: Make du(:,:,:,:) module allocatable array.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-08 09:28:44 -02:00
752a3600c6 EVOLUTION: Add optimal 10-stage 4th order SSPRK method.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-08 08:48:32 -02:00
bd2d46e6cf INTERPOLATIONS: Make arguments compatible for reconstruct_*().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-06 20:55:43 -02:00
dec44886a0 INTERPOLATIONS: Remove the vector length argument from reconstruc().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-06 12:00:04 -02:00
c91cc1bf51 INTERPOLATIONS: Rewrite mp_limiting().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-06 08:36:07 -02:00
fa1c934f13 INTERPOLATIONS: Rewrite fix_positivity().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-06 08:29:35 -02:00
1bcf69c4e3 INTERPOLATIONS: Rewrite clip_extrema().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-06 08:25:22 -02:00
5193ff3250 SCHEMES: Get rid of argument n in Riemann solvers.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-05 22:08:30 -02:00
8af6143f72 EQUATIONS: Remove explicit argument dimensions in esystem_roe().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-05 18:11:38 -02:00
e9d8bcc6cc EQUATIONS: Remove the vector length argument from prim2cons().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-05 18:07:59 -02:00
218ca06190 EQUATIONS: Remove the vector length argument from fluxspeed().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-05 17:26:29 -02:00
0c092cb19f EQUATIONS: Remove the vector length argument in cons2prim().
The vector length is determined inside the subroutine.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-05 15:48:42 -02:00
2c66f83d80 SCHEMES: Use qi and fi for directional states and fluxes.
Since block dimensions are the same in all directions, only one state
vector and flux vector is required in update_flux_*() subroutines.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-05 15:34:35 -02:00
6728086834 COORDINATES: Remove ib, ie, ubl, ieu, etc.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-05 15:15:29 -02:00
a53e65611f BOUNDARIES: Use nb, ne instead of ib, ie, jb, je, kb, ke.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-05 15:09:08 -02:00
ead582e451 INTERPOLATIONS: Use nb, ne instead of ib, ie, jb, je, kb, ke.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-05 14:53:07 -02:00
41709a49cb EQUATIONS: Use nb, ne instead of ib, ie, jb, je, kb, ke.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-05 14:40:45 -02:00
e65d9137d6 SCHEMES: Use nb, ne instead of ib, ie, jb, je, kb, ke.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-05 14:19:49 -02:00
46cae6fe6b INTEGRALS: Use nb, ne instead of ib, ie, jb, je, kb, ke.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-05 13:14:45 -02:00
05bbead94c EVOLUTION: Use nb, ne instead of ib, ie, jb, je, kb, ke.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-05 13:06:24 -02:00
ec9a5ee2bd IO: Use nb, ne instead of ib, ie, jb, je, kb, ke.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-05 13:02:27 -02:00
9b7f34d7ed REFINEMENT: Use nb, ne instead of ib, ie, jb, je, kb, ke.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-05 12:59:44 -02:00
b3404c186e MESH: Use nb, ne instead of ib, ie, jb, je, kb, ke.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-05 12:52:17 -02:00