2339 Commits

Author SHA1 Message Date
6362fe13ad EVOLUTION: Fix partial times and timesteps in 3S*+ methods.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-23 20:52:03 -03:00
8297482639 IO: Store 'cmax' so we can properly resume jobs.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-23 11:23:58 -03:00
44ea63bc3d EVOLUTION: Implement optimized embedded methods by Ranocha et al.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-20 18:08:54 -03:00
27041d0feb DRIVER, EVOLUTION, IO: Rewrite handling precise snapshots.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-19 14:03:57 -03:00
55b0d24b77 EVOLUTION: Rework slightly the way new time step is calculated.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-19 13:28:25 -03:00
c604467a5c DRIVER: Use the initialize_time_step() instead of new_time_step().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-19 13:15:04 -03:00
d501b19041 EVOLUTION: Implement estimation of the initial time step.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-19 13:11:41 -03:00
6ca36c93e3 DRIVER, EVOLUTION: Move dtnext to evolution module.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-18 11:33:14 -03:00
0c8d8fac9c EVOLUTION: Rewrite SSPRK4(2)M integration method.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-16 19:03:25 -03:00
e7612f5cfd PYTHON: It should be 'block x of y'.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-15 23:03:26 -03:00
418a580888 EVOLUTION: Print 'err/tol' not the error during the integration.
This is more meaninful, since it shows the ratio of the maximum error to
the tolerance atol + rtol * max(U).

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-15 15:29:42 -03:00
155b5d7575 EVOLUTION: Implement L2 and maximum norms for error calculation.
Alow the user to choose between these two norms.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-15 13:51:25 -03:00
271af020df EVOLUTION: Implement 3rd order embedded SSP3(2)4 with error control.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-15 12:15:14 -03:00
27b6daebca PYTHON: Version 0.9.5
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-06 18:06:43 -03:00
0d4e955eb0 PYTHON: Reduce dimensions for 2D datasets in AmunH5.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-06 18:06:07 -03:00
9a7235b72d PYTHON: Use None for fields in chunks with dblocks == 0.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-06 18:05:58 -03:00
51cc3ea83b PYTHON: Determine toplev only for chunks with dblocks > 0.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-06 18:05:46 -03:00
d217d092ab PYTHON: Version 0.9.4
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-06 13:31:05 -03:00
e13bd00a0e PYTHON: Replace numpy.empty with numpy.zeros in AmunXML.__fill_chunks().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-06 13:07:57 -03:00
07fc5c20b0 PYTHON: Create empty chunk fields for chunks with no dblocks.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-06 13:02:31 -03:00
5eec85e598 PYTHON: Version 0.9.3.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-06 10:06:16 -03:00
ec8ddac795 PYTHON: Add velocity and magnetic field components.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-06 10:05:18 -03:00
f4df543b2c PYTHON: Import numpy without abbreviation.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-06 09:44:49 -03:00
b7496e0b3d PYTHON: Make sure dataset() accepts maxlev > 9.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-06 09:37:15 -03:00
c2c9c37b45 PYTHON: Version 0.9.2
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-05 11:12:07 -03:00
6945b90b66 PYTHON: Export WriteVTK through the AmunPy interface.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-05 11:11:17 -03:00
5a8526ef82 PYTHON: Version 0.9.1.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-05 10:59:39 -03:00
1c4fa6035d PYTHON: Rewrite vtkio.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-05 10:58:44 -03:00
e153cd417e PYTHON: Rewrite interpolations.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-05 10:19:18 -03:00
d714974b33 PYTHON: Restore zmin/zmax/zblocks for 2D AmunXML snapshots.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-05 09:38:49 -03:00
b6edd6ca97 PYTHON: Version 0.9
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-05 09:13:24 -03:00
764c1596b5 PYTHON: Improve comments in AmunH5 structures.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-05 09:09:24 -03:00
919e0ce6c2 PYTHON: Separate deprecated AmunH5 functions.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-05 09:05:57 -03:00
dbf977de36 PYTHON: Rewrite AmunXML class to use Amun base class.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-05 09:03:23 -03:00
61d5eeb839 PYTHON: Convert returns to exceptions in Amun class.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-04 22:24:08 -03:00
8e125253f7 PYTHON: Make AmunH5 class public.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-04 22:19:43 -03:00
da339f1ac9 PYTHON: Implement AmunH5 class for Amun's HDF5 format.
Mark the previous subroutines deprecated.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-04 22:17:53 -03:00
d28194dc23 PYTHON: Add amun.py defining the base class for Amun snapshots.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-04 22:15:59 -03:00
abd1749945 PYTHON: Use None for no compression of VTK files.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-04 22:14:30 -03:00
6ffb3f5910 INTERPOLATIONS: Add a couple of limiters.
Use limiter_average() which is not a limiter, just returns average of
left and right derivatives, as the default limiter for prolongation.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-09-21 11:29:16 -03:00
d3c0737c8d BOUNDARIES, MESH: Allow to choose prolongation limiter.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-09-21 11:28:23 -03:00
61ed993821 EVOLUTION: Fix decay rate of ψ field.
This scalar potential tracks the development of div(B). Its decay
depends on the cell size, so make it calculated per refinement level.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-09-20 08:44:38 -03:00
7b79d88b6d CMAKE: Turn off HDF5 by default.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-09-12 12:02:04 -03:00
9d48cc9367 OPERATORS: Implement high-order differential operators.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-09-10 18:25:11 -03:00
7713df3ebf EVOLUTION: Add dth to track hydrodynamic timestep.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-08-17 11:23:36 -03:00
c766bb9059 MESH: Redistribute data blocks after each refinement.
There are situations when one process can suddenly have many data
blocks refined resulting in a significant memory stress before
the new blocks are redistributed among all processes. Just
redistribute data blocks each time the block refinement is
performed to avoid such a situation.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-07-20 17:58:44 -03:00
6dbda8bf14 IO: Use the same array to store each variable in write_snapshot_xml().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-07-20 17:56:51 -03:00
a0b6e60d7e EVOLUTION: Simplify conditions for the new time step.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-07-19 15:12:08 -03:00
b2a7bd08b6 EVOLUTION: Consider dte only for methods with error control.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-07-19 13:16:31 -03:00
8ec801228c EVOLUTION: Fix formatting once again.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-07-19 13:03:40 -03:00