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

227 Commits

Author SHA1 Message Date
Grzegorz Kowal
e6e73dd4f1 SCHEMES: Unify numerical_flux() subroutines. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:46:02 -03:00
Grzegorz Kowal
1ad1ca8890 SCHEMES: Introduce pointers to Riemann solvers. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:58 -03:00
Grzegorz Kowal
733dca391d SCHEMES: Separate ROE Riemann solver for adiabatic MHD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:52 -03:00
Grzegorz Kowal
45222ea63b SCHEMES: Separate ROE Riemann solver for isothermal MHD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:47 -03:00
Grzegorz Kowal
20a48c2543 SCHEMES: Separate ROE Riemann solver for adiabatic HD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:42 -03:00
Grzegorz Kowal
ea4b69050b SCHEMES: Separate ROE Riemann solver for isothermal HD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:36 -03:00
Grzegorz Kowal
a9483f320a SCHEMES: Separate HLLC Riemann solver for adiabatic SRMHD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:31 -03:00
Grzegorz Kowal
8f27c838df SCHEMES: Separate HLLC Riemann solver for adiabatic SRHD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:25 -03:00
Grzegorz Kowal
bae69d1a95 SCHEMES: Separate HLLD Riemann solver for adiabatic MHD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:20 -03:00
Grzegorz Kowal
366ad64f43 SCHEMES: Separate Mignone HLLD Riemann solver for isothermal MHD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:15 -03:00
Grzegorz Kowal
1bddacc2d8 SCHEMES: Separate HLLD Riemann solver for isothermal MHD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:10 -03:00
Grzegorz Kowal
5f366ab2eb EQUATIONS: Add indicator of positive variables. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:06 -03:00
Grzegorz Kowal
c2c2a65dbf SCHEMES: Unify subroutine update_flux(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:01 -03:00
Grzegorz Kowal
87216b1b73 EQUATIONS: Add the vector of variable indices for update_flux(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:56 -03:00
Grzegorz Kowal
60539e0558 SCHEMES: Move subroutines update_flux() up. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:52 -03:00
Grzegorz Kowal
be548fde57 SCHEMES: Separate HLLC Riemann solver for adiabatic MHD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:45 -03:00
Grzegorz Kowal
dc44a3ba41 SCHEMES: Separate HLLC Riemann solver for adiabatic HD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:40 -03:00
Grzegorz Kowal
f3aa475f63 SCHEMES: Unify HLL Riemann solvers. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:33 -03:00
Grzegorz Kowal
bf9ef52fb9 SCHEMES: Rename riemann pointer to numerical_flux. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:28 -03:00
Grzegorz Kowal
c54505917a SCHEMES: Use rather subroutines instead of pointers from equations. 
Use compatible (for the same set of equations and equation of state) in
Riemann solvers, not the procedure pointers.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:23 -03:00
Grzegorz Kowal
9801ac6d8b SCHEMES: Fix description of riemann_hll(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:19 -03:00
Grzegorz Kowal
fed3900810 SCHEMES: Separate universal HLL solver. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:14 -03:00
Grzegorz Kowal
1154cb5ec0 SCHEMES: Use higher_order_flux_correction(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:10 -03:00
Grzegorz Kowal
1acb7cca04 SCHEMES: Add a separate subroutine for higher order flux corrections. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:04 -03:00
Grzegorz Kowal
ad37de3a99 EQUATIONS, SCHEMES: Merge speeds cm and cp into one variable. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:43:55 -03:00
Grzegorz Kowal
f3670ddf37 EQUATIONS: Reorganize variable indices. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:43:45 -03:00
Grzegorz Kowal
ad0303d3ba Update copyright. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-02 13:18:53 -03:00
Grzegorz Kowal
f36979d74a SCHEMES: Let user control the high order flux correction - corrected. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-04 11:12:05 -03:00
Grzegorz Kowal
7578f1e9f2 Revert "SCHEMES: Let user control the high order flux correction." 
This reverts commit 0dcc3cc95c384c8ed63d6b62aa1cdbcaa58e6672.
 2019-10-04 11:11:05 -03:00
Grzegorz Kowal
0dcc3cc95c SCHEMES: Let user control the high order flux correction. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-04 11:02:05 -03:00
Grzegorz Kowal
0806bfc1e0 EQUATIONS: Do not specify variable range if not necessary. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-03 22:17:09 -03:00
Grzegorz Kowal
92b905cecd SCHEMES: Pass only flux variables to Riemann solvers. 
There is no need to interpolate or pass all variables to the Riemann
solvers. Just pass those related to flux calculation.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-03 21:48:04 -03:00
Grzegorz Kowal
e26d9b23e4 PROBLEMS: Implement passive scalar for KH test problem. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-03 21:14:05 -03:00
Grzegorz Kowal
6506083785 EQUATIONS: Update scalars in correct_unphysical_states() too. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-03 13:52:34 -03:00
Grzegorz Kowal
c1f8233613 BOUNDARIES: Update conservative scalars after boundary update too. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-03 13:47:33 -03:00
Grzegorz Kowal
fe4a712c76 PROBLEMS: Calculate conservative representation of scalars in RT. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-03 13:42:22 -03:00
Grzegorz Kowal
1037377414 PROBLEMS: Add passive scalar initialization for RT problem. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-03 13:39:16 -03:00
Grzegorz Kowal
31ac085d70 EQUATIONS: Take into account passive scalars in prim2cons(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-03 13:38:27 -03:00
Grzegorz Kowal
f22cc55b58 EQUATIONS: Take into account passive scalars in cons2prim(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-03 13:27:47 -03:00
Grzegorz Kowal
8e9e8f196c EVOLUTION: Update the increment for passive scalars. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-02 21:16:29 -03:00
Grzegorz Kowal
a80b8cf41f BOUNDARIES, DRIVER, EVOLUTION: Use nf for the number of block fluxes. 
Now, since the number of variables and fluxes may be different, it is
important to keep the size of the block flux array using nf instead of
nv.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-02 15:39:25 -03:00
Grzegorz Kowal
a5094d1e6b EQUATIONS: Print number of fluxes and passive scalars. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-02 14:03:12 -03:00
Grzegorz Kowal
5760b6c40c EQUATIONS: Add parameters for the number of fluxes and passive scalars. 
Initialize these parameters and generate the names of passive scalars.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-02 12:28:33 -03:00
Grzegorz Kowal
0133adfc53 MPITOOLS, BOUNDARIES: Simplify arguments of exchange_real_arrays(). 
Some arguments, like sproc or stag, are redundant, so remove them.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-01 11:15:31 -03:00
Grzegorz Kowal
05ec954e16 COORDINATES: Print zmn & zmx, not ymn & ymx if zmn > zmx. 
Reported-by: Lorenzo Nativi <lorenzo.nativi@astro.su.se>
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-09-23 22:32:45 -03:00
Grzegorz Kowal
5847a4f086 SCHEMES: Correct numerical fluxes for higher-orders. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-03-30 22:50:08 -03:00
Grzegorz Kowal
368894560b INTERPOLATIONS: Add order module parameter. 
This parameter keeps the order of the used reconstruction.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-03-30 21:12:28 -03:00
Grzegorz Kowal
4b7615e72f INTEGRALS: Use 'magnetized' instead of 'ibx > 0'. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-03-08 07:03:39 -03:00
Grzegorz Kowal
055f67c874 MESH: Replace allocatable array with module one in prolong_block(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-03-07 11:12:49 -03:00
Grzegorz Kowal
7b241c90a8 MESH: Add status argument to prolong_block(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-03-07 10:51:00 -03:00