1064 Commits

Author SHA1 Message Date
1e6b6f71e2 IO: Rewrite subroutine read_attributes_h5().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-23 14:39:01 -03:00
5180d03f2f IO, MPITOOLS: Spell checking.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-22 09:35:28 -03:00
140daaea04 IO: Rewrite subroutine write_attribute_vector_integer_h5().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-21 17:03:26 -03:00
a169a55085 IO: Rewrite subroutine read_attribute_vector_integer_h5().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-21 16:55:14 -03:00
07ca675497 IO: Rewrite subroutine write_attribute_double_h5().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-21 16:18:11 -03:00
3cd9575685 IO: Rewrite subroutine read_attribute_double_h5().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-21 16:13:36 -03:00
6e6e4141e0 IO: Rewrite write_attribute_integer_h5().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-21 16:06:42 -03:00
624ee8be66 IO: Spell corrections in read_attribute_integer_h5().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-21 15:49:12 -03:00
740fedede9 IO: Rewrite read_attribute_integer_h5().
The attribute, provided by the group identifier to which it is linked
and its name, is opened inside the subroutine
read_attribute_integer_h5().

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-21 15:46:19 -03:00
29c842155f IO: Quit write_attributes_h5() if cannot create group.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-21 13:13:44 -03:00
b746cd1c7e IO: Rewrite write_attributes_h5().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-21 12:58:12 -03:00
3132fb8da8 IO: Correct the description of function next_tout().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-21 11:56:48 -03:00
67a9beaf65 IO: Make the snapshot number continuously increasing.
If we restart a job with a different snapshot interval, the snapshot
number will continue from a different value corresponding to the
simulation time.

The snapshot number shouldn't be related to the current simulation time.
It should be always increasing.  The simulation time is stored as an
attribute, anyway.  This patch makes the snapshot number always
increasing continuously, without any gaps for restarted jobs.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-21 11:49:55 -03:00
f012eefe9f IO: Slightly rewrite write_restart_snapshot().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-19 13:50:44 -03:00
0a81a01242 MESH: Micro optimization in prolong_block().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-17 11:52:06 -03:00
d3a5055cf5 BOUNDARIES: Micro optimization in boundary_prolong().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-17 11:44:47 -03:00
b226959097 MESH: Interpolation symmetry in prolong_block().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-17 11:35:15 -03:00
a4805761db BOUNDARIES: Interpolation symmetry in boundary_prolong().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-17 11:31:33 -03:00
f85db888ff BOUNDARIES: Replace string boundary types with integer ones.
Now, the boundary type is stored in an integer array instead of string
variables.  This should reduce the comparison time in
boundary_specific().

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-09 13:48:05 -03:00
9b7f0232cd PROBLEMS: Make implosion problem symmetric about axes.
Since we use the reflective boundary conditions for this problem, the
initial profile should respect that and also apply reflection at the
boundaries.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-09 11:01:20 -03:00
b528322c63 BOUNDARIES: Apply specific boundaries to primitive variables.
Convert the specific boundaries update to use primitive variables
instead of conservative ones.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-09 10:05:39 -03:00
842f6bdd51 BOUNDARIES: Prolong primitive variables instead of conservative.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-09 09:54:12 -03:00
0876badc16 BOUNDARIES: Restrict primitive variables instead of conservative.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-09 09:49:21 -03:00
10072f5c04 BOUNDARIES: Copy primitive variables instead of conservative.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-09 09:45:59 -03:00
722d623e0b BOUNDARIES: Convert update_corners() to update primitive variables.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-09 09:40:33 -03:00
125e267f97 BOUNDARIES: Add subroutine to update conservative variables.
This subroutine is called after updating the boundaries of primitive
variables.  It converts those primitive variables to their conservative
representation in all ghost cells.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-09 09:35:15 -03:00
f874511673 EVOLUTION: Update primitive variables before boundaries.
The new order is to first, convert the updated conserved variables to
their primitive representation in the block interiors, update boundaries
of the primitive variables, and finally, convert the primitive variables
at the ghost cells to conserved ones.

The change updates the primitive variables in the block interiors before
the boundary update.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-09 08:51:04 -03:00
7255561d24 EQUATIONS: Conserved to primitive variables only in interioris.
Conversion of the conserved to primitive variables can be expensive for
some sets of equations, such as the relativistic MHD.  Therefore, we
will convert only the interior parts of blocks and update the ghost
zones later in the boundary update.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-09 08:45:44 -03:00
f9529def8c PROBLEMS: Implement implosion test problem.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-02 16:17:30 -03:00
22f71ac94e DRIVER: Set the maximum number of iterations to a huge number.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-04-02 15:25:08 -03:00
e18dfc5977 BOUNDARIES: Reintroduce edge/corner update.
This guarantees that the very edges or corners are always physical and
minimizes their influence on the variable evolution in the ghost zones.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-03-30 12:23:38 -03:00
be40a6113c IO: Store ghost cells in the snapshots by default.
This is just in case we would like to use interpolation to reconstruct
the full domain resolution.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-03-19 19:04:42 -03:00
88edaec837 EQUATIONS, SCHEMES: Implement Roe's solver for adiabatic MHD.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-03-06 14:19:31 -03:00
42c738c78d EQUATIONS, SCHEMES: Implement Roe's solver for isothermal MHD.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-03-06 13:55:09 -03:00
78b15d5596 EQUATIONS, SCHEMES: Implement Roe's solver for adiabatic HD.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-03-06 13:16:57 -03:00
a8d575e471 EQUATIONS, SCHEMES: Implement Roe's solver for isothermal HD.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-03-06 12:59:51 -03:00
a49dad5ad0 SCHEMES: Separate Riemann state reconstruction from Riemann solvers.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-03-06 12:15:44 -03:00
f1378c6e4e SCHEMES: Group subroutines according to the equations.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-03-05 19:04:23 -03:00
73cbed8192 PROBLEMS: Rewrite module.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-02-07 13:57:28 -02:00
55de38609d SHAPES: Update module dependencies and spell check.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-02-07 13:50:17 -02:00
e0ca4b0392 SHAPES: Add new module to handle shapes embedded in domain.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-02-07 12:12:27 -02:00
4bbd84147f BLOCKS: Rewrite neighbor selection for refinement.
In this patch we make a general subroutine iterating over all neighbors
iterate_over_neighbors() which takes as arguments a pointer to the block
which neighbors are evaluated and pointer to subroutine which should be
called with each neighbor. So far two such subroutines are used, one to
select the proper refinement flag, and another to select the neighbor
for update.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-02-07 10:25:53 -02:00
9efa568ddf PROBLEMS: Change default B angle to 0.0 in KH problem.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-01-31 11:08:44 -02:00
613ef775bb REFINEMENTS: Modify the way of estimating error.
We take the maximum error along each direction instead of calculating
the norm sqrt(ex*ex+ey*ey+ez*ez).  This solves the symmetry problems
with some of the reconstructions.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-01-31 10:50:18 -02:00
35cbbc9313 PROBLEMS: Add Kelvin-Helmholtz instability problem setup.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-01-31 10:47:05 -02:00
ffc6d0b538 TIMERS: Avoid usingg MAX() with INTEGER(KIND=8) variables.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-01-30 18:03:33 -02:00
32fb3177e0 REFINEMENT: Fix return value kind in check_refinement_criterion().
The meta block field %refine is 4 byte integer. Make the value returned
by check_refinement_criterion() of the same type.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-01-30 17:57:01 -02:00
29a1aa67d7 DRIVER: Add missing comma.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-01-30 17:55:41 -02:00
57e3e7ec42 DOMAINS: Rewrite module.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-01-29 16:20:42 -02:00
fa3ea8ca5d COORDINATES: Rewrite this module.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2014-01-29 16:05:45 -02:00