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

1650 Commits

This Branch
This Branch
All Branches
Author SHA1 Message Date
Grzegorz Kowal
704ec5754d SCHEMES: Make state variable argument 'inout' in Riemann solvers. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2016-08-14 19:15:42 -03:00
Grzegorz Kowal
bf38b8ad04 SCHEMES: Remove reconstruct_states_1d(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2016-08-14 19:04:36 -03:00
Grzegorz Kowal
87ac23663e INTERPOLATIONS: Implement interfaces separately along each direction. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2016-08-12 23:39:08 -03:00
Grzegorz Kowal
b2830f7979 INTERPOLATIONS: Limit derivatives instead in interfaces_tvd(). 
Try better for positive variable interpolation. Limit the derivatives
instead of zeroing them.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2016-08-12 20:37:53 -03:00
Grzegorz Kowal
a74240e3b5 Merge branch 'master' into full_reconstruction 2016-08-12 20:27:24 -03:00
Grzegorz Kowal
c58117d7ed BOUNDARIES: Guarantee positivity in block_face_prolong(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2016-08-12 20:25:28 -03:00
Grzegorz Kowal
68087041db BOUNDARIES: Guarantee positivity in block_edge_prolong(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2016-08-12 20:21:55 -03:00
Grzegorz Kowal
17f46cf1d6 BOUNDARIES: Guarantee positivity in block_corner_prolong(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2016-08-12 20:17:10 -03:00
Grzegorz Kowal
753b288b71 MESH: Make sure the negativity is detected in the full cell. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2016-08-12 20:11:11 -03:00
Grzegorz Kowal
25ba130da3 MESH: Guarantee the positivity of positive variables in prolongation. 
We use 1D TVD limiter to prolongate blocks, but in 2D or 3D they can
produce negative values of positive variables in the corners. Detect
such situation and limit the derivatives until the whole interpolated
region is positive.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2016-08-12 19:24:31 -03:00
Grzegorz Kowal
cbf9c081ec SCHEMES: Move interface interpolation to module INTERPOLATIONS. 
This will allow for implementations of truly multidimentional
interpolations.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2016-08-09 16:00:41 -03:00
Grzegorz Kowal
1ee7e0b402 SCHEMES: Adopt update_flux_hd_iso() to use separate reconstruction. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2016-08-01 23:44:09 -03:00
Grzegorz Kowal
e5cfae408a SCHEMES: Move reconstruction part to new subroutine. 
The reconstruction in all three directions has been moved to a separate
subroutine in order to make it generic. It will allow for introducing a
full two or three dimensional reconstruction. This is for adiabatic
hydrodynamics so far, but other set of equations will be supported in
the following changes.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2016-08-01 23:36:15 -03:00
Grzegorz Kowal
954a287129 SCHEMES: Separate reconstruction and solver in update_flux_hd_adi(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2016-07-30 00:22:20 -03:00
Grzegorz Kowal
677284b003 EVOLUTION, SCHEMES: Calculate all directional fluxes at once. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2016-07-26 22:47:55 -03:00
Grzegorz Kowal
744c213648 INTERPOLATIONS: Implement reconstruction based on Gaussian Processes. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2016-07-23 19:18:41 -03:00
Grzegorz Kowal
84b20d6d6a ALGEBRA: Add subroutine to invert arbitrary matrix. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2016-07-23 14:43:20 -03:00
Grzegorz Kowal
bc1a7b9c98 OPERATORS: Make derivative subroutines public as well. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2016-05-13 13:53:08 -03:00
Grzegorz Kowal
a8afe2d3d7 BOUNDARIES: Handle magnetic field in the reflective boundaries. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2016-03-10 07:23:09 -03:00
Grzegorz Kowal
1b8add50eb MESH: Fix block statistics calculation after job restart. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-18 06:55:32 -02:00
Grzegorz Kowal
14d2b82a02 INTERPOLATIONS: Implement low-dissipation version of CRMP5 method. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-14 18:59:11 -02:00
Grzegorz Kowal
6412e00d03 INTERPOLATIONS: Make two extrema points monotonic in WENO5NS. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-14 09:48:11 -02:00
Grzegorz Kowal
4ce1b6574b INTERPOLATIONS: Make two extrema points monotonic in WENO5YC. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-14 09:46:25 -02:00
Grzegorz Kowal
a66cb723e4 INTERPOLATIONS: Make two extrema points monotonic in WENO5Z. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-14 09:44:45 -02:00
Grzegorz Kowal
c7e0a3e17c INTERPOLATIONS: Make the extrema points monotonic in CRWENO5NS. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-13 12:55:34 -02:00
Grzegorz Kowal
e49a186f22 INTERPOLATIONS: Make the extrema points monotonic in CRWENO5YC. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-13 12:48:54 -02:00
Grzegorz Kowal
1185ee2072 INTERPOLATIONS: Make the extrema points monotonic in CRWENO5Z. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-13 12:42:18 -02:00
Grzegorz Kowal
b60840bc53 INTERPOLATIONS: Make the extrema points monotonic in WENO5NS. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-13 12:34:53 -02:00
Grzegorz Kowal
2f0d9f056b INTERPOLATIONS: Make the extrema points monotonic in WENO5YC. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-13 12:31:23 -02:00
Grzegorz Kowal
2f26b79e79 INTERPOLATIONS: Make the extrema points monotonic in CRMP5. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-13 11:34:03 -02:00
Grzegorz Kowal
c22a5437b8 INTERPOLATIONS: Make the extrema points monotonic in MP5. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-13 11:10:31 -02:00
Grzegorz Kowal
28324447c4 INTERPOLATIONS: Parameter eps should be at least 1.0d-12 for LIMO3. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-13 10:59:31 -02:00
Grzegorz Kowal
5b5efea431 INTERPOLATIONS: Make LIMO3 reconstruction symmetric. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-13 10:56:28 -02:00
Grzegorz Kowal
f96c0d0297 INTERPOLATIONS: Make the extrema points monotonic in WENO5Z. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-13 09:14:34 -02:00
Grzegorz Kowal
dbe5fe71a0 INTERPOLATIONS: Make the extrema points monotonic in LIMO3. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-13 08:59:21 -02:00
Grzegorz Kowal
a876581745 INTERPOLATIONS: Make the extrema points monotonic in WENO3. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-13 08:53:15 -02:00
Grzegorz Kowal
d2ca2551dd INTERPOLATIONS: Make the extrema points monotonic in TVD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-13 08:48:20 -02:00
Grzegorz Kowal
164f809ee5 Revert "INTERPOLATIONS: Make sure MC, VL, and SB limiters do not overshoot." 
This reverts commit 99cc1ca2abcf3e52ecfcf6bc66ee051cf034bbce.
 2015-12-10 11:22:23 -02:00
Grzegorz Kowal
9cf87ed959 EQUATIONS: Make optional arguments really optional in fluxspeed(). 
This commit makes the optional arguments working only in relativistic
systems.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-10 07:25:17 -02:00
Grzegorz Kowal
5056ac452b EQUATIONS: Make optional arguments really optional in fluxspeed(). 
This commit makes the optional arguments working only in classical systems.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-10 07:07:10 -02:00
Grzegorz Kowal
5957461d92 REFINEMENT: Add vorticity based criterion. 
Also, calculate vorticity and current density assuming dx=dy=dz=1, so
the criterion value is level independent.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-09 07:41:46 -02:00
Grzegorz Kowal
224cc37b53 INTERPOLATIONS: Remove unused variables. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-08 09:54:35 -02:00
Grzegorz Kowal
c382b714ad INTERPOLATIONS: Use tridiag() solver in CRWENO5NS reconstruction. 
Also fix the matrix coefficient calculation.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-08 09:53:34 -02:00
Grzegorz Kowal
e3d91d48f1 INTERPOLATIONS: Use tridiag() solver in CRWENO5YC reconstruction. 
Also fix the matrix coefficient calculation.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-08 09:23:13 -02:00
Grzegorz Kowal
52f520f215 INTERPOLATIONS: Use tridiag() solver in CRWENO5Z reconstruction. 
Also fix the matrix coefficient calculation.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-08 09:14:48 -02:00
Grzegorz Kowal
a00dc32c71 INTERPOLATIONS: Use tridiag() solver in CRMP5 reconstruction. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-08 08:21:39 -02:00
Grzegorz Kowal
eb118745a1 ALGEBRA: Add tridiagonal linear system solver. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-08 08:03:08 -02:00
Grzegorz Kowal
68b56b08ce INTERPOLATIONS: Fix weigth coeffcients in compact WENO methods. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-06 12:28:26 -02:00
Grzegorz Kowal
5794ab4b13 INTERPOLATIONS: Fix weigth coeffcients in explicit WENO methods. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2015-12-06 12:03:59 -02:00
Grzegorz Kowal
132cd3326f Revert "INTERPOLATIONS: Correct weights in WENO5Z reconstruction." 
This reverts commit 9cf1eb846a3eec90f30f9fe6499adefc3ecc5519.
 2015-12-06 11:44:45 -02:00