f8cdeaeffd
EQUATIONS: Fix Jacobian calculation in SRHD 2D(W,u²) variable solver.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-19 22:23:05 -02:00
c9c4f88f61
EQUATIONS: Bisect if W is out of brackets in SRHD 1D(W) variable solver.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-19 20:14:03 -02:00
5240cad7b3
EQUATIONS: Add initial brackets and guess for all SRHD variable solvers.
...
Before we perform Newton-Raphson step to find the roots W, we can use
some estimates of the brackets, and from them we can calculate the
initial guess for W.
This change implements this step for all special relativity HD primitive
variable solvers.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-19 19:00:23 -02:00
017c98dca0
EQUATIONS: Add parameter for the maximum allowed Lorentz factor.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-19 18:57:42 -02:00
8c0a574b80
EQUATIONS: Add additional names to SRHD variable solvers.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-19 17:44:14 -02:00
8f59203036
EQUATIONS: Micro optimize SRHD 2D(W,v²) primitive variable solver.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-19 17:38:47 -02:00
a470a2a874
EQUATIONS: Fix and optimize SRHD 2D(W,u²) primitive variable solver.
...
One of the derivative was calculated in a wrong way. This patch also
incorporates the pressure derivatives directly in the derivatives of
functions F(W,u²) and G(W,u²), which allows to some microoptimizations.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-19 17:17:38 -02:00
28badcb066
EQUATIONS: Implement SRHD 2D(W,u²) primitive variable solver.
...
The solver approximates the root by searching for W and u² at the same
time.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-19 12:35:03 -02:00
01e6d799aa
EQUATIONS: Rename nr_iterate_srhd_adi_2d() to nr_iterate_srhd_adi_2dwv().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-19 11:38:16 -02:00
56da64d935
EQUATIONS: Improve the guess in SRHD primitive variables solver.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-19 11:35:34 -02:00
f46981fb03
EQUATIONS: Avoid division errors in nr_iterate_srhd_adi_2d().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-19 10:01:05 -02:00
020142137d
EQUATIONS: Get the W guess from previous step in SRHD and SRMHD.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-18 15:36:01 -02:00
6448ad40e8
SHAPES: Reorder variables in update_shapes_jet().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-17 18:34:48 -02:00
3010c36cbd
EQUATIONS: Implement 2D primitive variable solver for SRMHD.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-17 18:15:11 -02:00
49c141a2c5
EQUATIONS: Fix some comments.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-17 17:29:05 -02:00
dc58f6686d
PROBLEMS: Separate jet and ambient magnetic fields in jet problem.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-17 17:16:37 -02:00
0d62c46535
EQUATIONS: Simplify some expressions in nr_iterate_srmhd_adi_1dw().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-17 12:46:22 -02:00
cd6608de0e
EQUATIONS: Remove dependence on nr_function() in SRMHD equations.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-17 12:35:52 -02:00
e2406e10e7
EQUATIONS: Remove dependence on ww_to_v() in SRMHD subroutines.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2015-02-17 12:23:42 -02:00
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
