gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity
2,748 Commits 5 Branches 0 Tags
Commit Graph

899 Commits

Author SHA1 Message Date
Grzegorz Kowal
9829505650 Update copyright year to 2022. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-02-02 09:51:41 -03:00
Grzegorz Kowal
a1c465a3a1 SCHEMES: Keep alpha_f,s in the limites for ROE and KEPES. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-29 18:24:33 -03:00
Grzegorz Kowal
4645652e41 EQUATIONS: Remove all functions maxspeed(). 
These functions were replaced by get_maximum_speeds().

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-28 13:06:53 -03:00
Grzegorz Kowal
57c6a7c6ec EQUATIONS: Remove all subroutines for Roe solver. 
Roe Riemann solver with their eigensystems are constructed in SCHEMES
now.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-28 13:01:38 -03:00
Grzegorz Kowal
68ca97ac6c SCHEME: Simplify speed calculation in riemann_mhd_iso_roe(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-28 12:58:31 -03:00
Grzegorz Kowal
415d2a8efb SCHEME: Simplify speed calculation in riemann_mhd_iso_kepes(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-28 11:59:55 -03:00
Grzegorz Kowal
85a2f6ed32 SCHEME: Rewrite speed calculation in riemann_mhd_adi_kepes(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-28 11:49:46 -03:00
Grzegorz Kowal
509e126583 SCHEMES: Rewrite riemann_mhd_adi_roe(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-28 10:36:05 -03:00
Grzegorz Kowal
13332a7377 SCHEMES: Simplify symmetry preserving in riemann_mhd_adi_kepes(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-26 17:31:57 -03:00
Grzegorz Kowal
2fdb52d972 SCHEMES: Rewrite riemann_mhd_iso_roe(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-26 13:10:55 -03:00
Grzegorz Kowal
5c96f3fd5a SOURCES: Add KEPES source terms for isothermal GLM-MHD. 
These source terms guarantee the KEPES solver to conserve the entropy.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-22 11:57:01 -03:00
Grzegorz Kowal
a343a400e6 SCHEMES: Simplify symmetry preserving in riemann_mhd_iso_kepes(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-21 23:31:54 -03:00
Grzegorz Kowal
1514fc0872 SCHEMES: Micro optimizations in riemann_mhd_iso_hlld(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-21 21:29:38 -03:00
Grzegorz Kowal
03f5606461 SCHEMES: Remove unnecessary parenthesis in riemann_hd_adi_kepes(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-21 19:40:55 -03:00
Grzegorz Kowal
51a9794a84 SCHEMES: Simplify symmetry preserving in riemann_hd_adi_kepes(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-21 19:35:58 -03:00
Grzegorz Kowal
889882c48f SCHEMES: Simplify symmetry preserving in riemann_hd_iso_kepes(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-21 19:06:28 -03:00
Grzegorz Kowal
64d194933d SCHEMES: Simplify symmetry preserving in riemann_hd_iso_roe(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-21 18:59:20 -03:00
Grzegorz Kowal
257b23078c SCHEMES: Rewrite riemann_hd_adi_roe(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-21 18:47:17 -03:00
Grzegorz Kowal
6363ee5323 SCHEMES: Rewrite riemann_hd_hllc(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-20 18:02:05 -03:00
Grzegorz Kowal
1180d22cda SCHEMES: Rewrite riemann_hd_iso_kepes(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-20 17:48:33 -03:00
Grzegorz Kowal
beb31d0fe6 SCHEMES: Rewrite riemann_hd_iso_roe(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-20 16:57:06 -03:00
Grzegorz Kowal
982637c745 SCHEMES: Polish riemann_hll(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-20 12:51:34 -03:00
Grzegorz Kowal
f52539e34a SCHEMES: Reorganize Riemann solvers. 
Move calculation of the conservative variables and fluxes at the states
to the solvers themselfes. Also group solvers by the system of
equations.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-19 11:53:56 -03:00
Grzegorz Kowal
ce1280d00c IO: Rebuild the leaf and datablock lists in read_restart_snapshot(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-18 15:23:19 -03:00
Grzegorz Kowal
2b21d29da7 IO: Store/restore cglm. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-11 12:45:26 -03:00
Grzegorz Kowal
b7f5d28a45 SCHEMES: Micro optimizations in riemann_mhd_iso_kepes(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-09 17:48:50 -03:00
Grzegorz Kowal
124ad01eca EVOLUTION: Fix time step estimation with OpenMP. 
The variable um should be also thread private in initialize_time_step()
and new_time_step().

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-09 17:24:56 -03:00
Grzegorz Kowal
a634d9bed5 VARIOUS: Make thread number thread private. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-08 11:52:16 -03:00
Grzegorz Kowal
af4b83cd7c BOUNDARIES: Make variable initialization OpenMP conformant. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-08 11:39:48 -03:00
Grzegorz Kowal
c56daa6991 EVOLUTION: Make variable initialization OpenMP conformant. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-08 11:12:15 -03:00
Grzegorz Kowal
977fd180ec FORCING: Make variable initialization OpenMP conformant. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-08 11:09:45 -03:00
Grzegorz Kowal
865a9e3013 MESH: Make variable initialization OpenMP conformant. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-08 11:06:19 -03:00
Grzegorz Kowal
55576b5211 PROBLEMS: Make variable initialization OpenMP conformant. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-08 11:05:08 -03:00
Grzegorz Kowal
301d4fb449 REFINEMENT: Make variable initialization OpenMP conformant. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-08 11:02:05 -03:00
Grzegorz Kowal
e8facf4647 SHAPES: Make variable initialization OpenMP conformant. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-08 11:00:08 -03:00
Grzegorz Kowal
c57194aef7 BLOCKS: Make variable initialization OpenMP conformant. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-08 10:54:34 -03:00
Grzegorz Kowal
d8e959b219 SOURCES: Make variable initialization OpenMP conformant. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-08 10:52:05 -03:00
Grzegorz Kowal
9b8b578931 STATISTICS: Make nt save in store_statistics(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-08 10:49:12 -03:00
Grzegorz Kowal
67a0934242 EQUATIONS: Make variable initialization OpenMP conformant. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-08 10:45:02 -03:00
Grzegorz Kowal
210689da08 SCHEMES: Make variable initialization OpenMP conformant. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-08 10:37:53 -03:00
Grzegorz Kowal
9913468203 INTERPOLATIONS: Make variable initialization OpenMP conformant. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-08 10:36:59 -03:00
Grzegorz Kowal
9b5198a5ff EQUATIONS, EVOLUTION: Rework get_maximum_speeds(). 
There is no need for the maximum characteristic speed. The divergence
correcting speed cglm is calculated using the maximum eigenvalue and
physical velocity.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-07 15:38:54 -03:00
Grzegorz Kowal
c8ed7fc07a SCHEMES: Rewrite slighly riemann_mhd_adi_kepes(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-07 11:43:25 -03:00
Grzegorz Kowal
b744ab7180 SCHEMES: Use cglm speed in KEPES MHD Riemann solvers. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-06 22:49:45 -03:00
Grzegorz Kowal
c32e4adcde EVOLUTION: Use get_maximum_speeds() in the time step estimation. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-06 22:48:35 -03:00
Grzegorz Kowal
0cf2481992 EQUATIONS: Add subroutines get_maximum_speeds(). 
These subroutines determine the maximum physical and characteristic
speeds, and the maximum eigenvalue in the system.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-06 22:44:21 -03:00
Grzegorz Kowal
226998e80e SOURCES: Add source terms related to the MHD KEPES solver. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-06 16:46:36 -03:00
Grzegorz Kowal
25ea3c3daf EQUATIONS: Take into account ½ψ² in total energy for adiabatic MHD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-06 16:32:58 -03:00
Grzegorz Kowal
4f3a05f9ba SCHEMES: Optimize slightly riemann_mhd_adi_kepes(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-04 12:51:05 -03:00
Grzegorz Kowal
9af9357032 SCHEMES: Optimize slightly riemann_hd_adi_kepes(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2022-01-04 12:44:10 -03:00