1505 Commits

Author SHA1 Message Date
c438f1f366 EQUATIONS: Remove lower limit for enthalpy in nr_function().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-17 11:57:47 -02:00
9cc25b790e EQUATIONS: Implement subroutines for special relativity MHD.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-17 11:28:40 -02:00
e5e0d293cf PROBLEMS: Add magnetic field to the jet problem.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-17 10:25:52 -02:00
1930bda0cc ALGEBRA: Add cubic and quartic solvers.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-17 09:53:19 -02:00
b8fed51178 SCHEMES: Implement subroutines for special relativity MHD.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-16 20:06:43 -02:00
a36dc4b638 SCHEMES: Cite the paper next to SRHD HLLC solver name.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-16 12:56:26 -02:00
8e7d63a1d2 SCHEMES: Rewrite the HLLC solver for SRHD.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-16 12:52:49 -02:00
18a8182c56 SCHEMES: Rename HLLC Riemann solver to HLLC-M.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-16 10:56:46 -02:00
ba470d5ec9 EQUATIONS: Simplify the W guess and pressure in cons2prim_srhd_adi().
The initia guess for enthalphy can be estimated from the total energy
and minimum pressure. In the similar way, the pressure can be obtained
by subtracting the energy from enthalphy.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-15 08:17:51 -02:00
d203702254 SCHEMES: Implement HLLC Riemann solver for SRHD.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-14 23:08:32 -02:00
fa5173547b ALGEBRA: Add new module which provides algebra subroutines.
The initial module version provides the subroutines to solve quadratic
equation.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-14 20:47:04 -02:00
46a0022a01 EQUATIONS: Add timers to primitive variable solvers.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-14 18:22:29 -02:00
e46eafa0b9 EQUATIONS: Rewrite 1Dw primitive variable solver for SRHD.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-14 18:12:51 -02:00
9c74e2ed70 EQUATIONS: Implement 2D primitive variable solver for SRHD.
This solver finds approximation of enthalpy W and squared velocity |V|²
at the same time.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-14 17:20:03 -02:00
188cbc6231 EQUATIONS: Reorder SRHD flux components in fluxspeed().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-08 17:38:28 -02:00
fdc6f9e3e8 SCHEMES: Add reconstruction on 4-velocities to SRHD.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-06 17:41:17 -02:00
b4253ef8a7 PROBLEMS: Reorganize jet problem parameters.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-06 12:53:34 -02:00
2ed5f74aa3 PROBLEMS: Fix wrong density profile for jet problem.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-06 10:51:01 -02:00
e4140f2058 PROBLEMS: Implement relativistic jet problem.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-06 10:13:23 -02:00
2f26993c7b SCHEMES: Implement subroutines for special relativity hydrodynamics.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-06 09:09:22 -02:00
f96e0a5b5b EQUATIONS: Implement subroutines for special relativity hydrodynamics.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-06 09:02:56 -02:00
7e9a79abcb EQUATIONS: Calculate left-, right-going speeds in fluxspeed().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-05 22:07:19 -02:00
f84abc7d76 SCHEMES: Remove preprocessor conditions on NDIMS.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-26 20:08:35 -02:00
1a2a312775 MESH: Correct calculation of the adaptive mesh efficiency.
The efficiency should be large, if the adaptive mesh uses few blocks
comparing to the corresponding uniform mesh. When the number of adaptive
blocks increase, the efficiency should drop.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-25 22:01:54 -02:00
cbac66e430 MESH: Use the list of leafs in store_mesh_stats().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 15:50:05 -02:00
2d1b97642f BOUNDARIES: Use the list of leafs in boundaries_corner_prolong().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 15:44:59 -02:00
e4acb7028e BOUNDARIES: Use the list of leafs in boundaries_corner_restrict().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 15:41:29 -02:00
b1c8bedbf2 BOUNDARIES: Use the list of leafs in boundaries_corner_copy().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 15:38:55 -02:00
b0e7489c79 BOUNDARIES: Use the list of leafs in boundaries_edge_prolong().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 15:34:13 -02:00
14b7b0f7ff BOUNDARIES: Use the list of leafs in boundaries_edge_restrict().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 15:30:29 -02:00
56a05fde98 BOUNDARIES: Use the list of leafs in boundaries_edge_copy().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 15:26:23 -02:00
3a314e0dd8 BOUNDARIES: Use the list of leafs in boundaries_face_prolong().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 15:21:46 -02:00
fb84013114 BOUNDARIES: Use the list of leafs in boundaries_face_restrict().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 15:18:15 -02:00
64dfbe901b BOUNDARIES: Use the list of leafs in boundaries_face_copy().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 15:14:27 -02:00
193b182022 BOUNDARIES: Reduce some comments.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 15:08:19 -02:00
f6e7456fa8 BOUNDARIES: Split timers for the boundary update subroutines.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 15:05:48 -02:00
1c75056007 BOUNDARIES: Use the list of leafs in boundary_specific().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 14:51:06 -02:00
4198d35b6e BOUNDARIES: Use the list of leafs in boundary_fluxes().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 14:41:12 -02:00
94e8b5ad93 BLOCKS: Make BLOCK_LEAF structure public.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 13:15:29 -02:00
fbedbb21a5 MESH: Regenerate the leaf list after each mesh change.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 13:13:41 -02:00
fa33402072 BLOCKS: Add subroutine to build the list of leafs.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 13:12:21 -02:00
ff68847858 BLOCKS: Add subroutine to clean up the leaf list.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 12:47:32 -02:00
d9666c113c BLOCKS: Add new list to keep track of blocks.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 12:22:03 -02:00
8294b09a33 BLOCKS: Introduce new BLOCK_LEAF structure.
This structure contains only two pointers, one to the next block leaf,
and another to the leaf meta block. It's gonna be use to keep the list
of all leaf meta blocks.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-13 12:18:35 -02:00
a29ec0bea5 INTERPOLATIONS: Remove unused parameters.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-12 17:21:19 -02:00
8424846f9b SCHEMES: Remove multiply by reciprocal in HD Riemann solvers.
Multiply by reciprocal is faster but value unsafe.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-11 16:09:05 -02:00
13d4d6fe71 Revert "EQUATIONS: Avoid multiple divisions. Replace by inverted value."
This reverts commit be7f200bf4da2077b9ed4f0c0d3f622b9bc164d9.

Multiplication by reciprocal decreases precision slightly, by increasing
the performance. This is a known problem. We prefer accurate results,
instead of fast.
2015-01-11 15:59:41 -02:00
b8e708d99f SCHEMES: Slightly optimize riemann_hd_adi_hllc().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-09 18:40:10 -02:00
82b8d51697 SCHEMES: Slightly optimize riemann_hd_adi_hll().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-09 17:32:26 -02:00
be7f200bf4 EQUATIONS: Avoid multiple divisions. Replace by inverted value.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-01-09 17:31:28 -02:00