#66: IOHDF5::out3D_ghosts and friends doesn't work and corrupts 2D data slices
----------------------------------+-----------------------------------------
Reporter: bcmsma@… | Owner: eschnett
Type: defect | Status: new
Priority: minor | Milestone:
Component: Carpet | Version:
Keywords: |
----------------------------------+-----------------------------------------
While the assigned values for the parameters
IOASCII::output_symmetry_points = "no"
IOASCII::out3D_ghosts = "no"
IOASCII::out3D_outer_ghosts = "no"
work as intended when CarpetIOASCII is active, it doesn't work for
CarpetIOHDF5:
IOHDF5::out3D_ghosts = "no"
IOHDF5::out3D_outer_ghosts = "no"
IOHDF5::output_symmetry_points = "no"
It doesn't do anything to the 3D data but it corrupts the 2D data slices
while
chopping out those regions. It would be nice to have them working for
IOHDF5
method, specially when visualising Pi symmetric data.
Note also that this report is for git version of Carpet. These parameters
seem to become deprecated for the hg version. Does anyone use them
regularly? Any
substitute in mind?
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/66>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#223: Make declaration of CCTK_ARGUMENTS safer
-------------------------+--------------------------------------------------
Reporter: eschnett | Owner:
Type: enhancement | Status: new
Priority: minor | Milestone:
Component: Cactus | Version:
Keywords: |
-------------------------+--------------------------------------------------
I suggest to change the declaration of CCTK_ARGUMENTS in C from
cGH * cctkGH
to
cGH const * CCTK_RESTRICT const cctkGH
which should lead to safer code and may even enable some optimisations.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/223>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#425: Case sensitivity GSL_DIR=BUILD
----------------------------------------+-----------------------------------
Reporter: reisswig@… | Type: enhancement
Status: new | Priority: minor
Milestone: | Component: Cactus
Version: | Keywords: GSL_DIR BUILD
----------------------------------------+-----------------------------------
The thorn GSL would not compile if "GSL_DIR=build" in the option list. In
order to work, it requires "GSL_DIR=BUILD", i.e. it is really case
sensitive
although I think, this is not required and should be case insensitive.
Possibly, HDF5_DIR=build may also not work for the same reason (have not
checked; should be checked as well).
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/425>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#234: Carpet mailing list is not functioning
--------------------+-------------------------------------------------------
Reporter: hinder | Owner: eschnett
Type: defect | Status: new
Priority: major | Milestone:
Component: Carpet | Version:
Keywords: |
--------------------+-------------------------------------------------------
The Carpet mailing list has been unavailable since 22-Mar-2010. This list
provides a forum for discussion of Carpet features and issues, and would
be the first port-of-call for new users. As such I think it would be good
to resurrect the mailing list.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/234>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#145: EOS_Omni requires HDF5 library. What about if it only uses HDF5 instead?
---------------------------+------------------------------------------------
Reporter: bmundim | Owner:
Type: enhancement | Status: new
Priority: major | Milestone:
Component: Other | Version:
Keywords: EOS_Omni HDF5 |
---------------------------+------------------------------------------------
At the moment EOS_Omni requires the HDF5 library to read
the nuclear EOS values from a table. Since not all distros have
fortran support for the HDF5 library as a default, I propose
to change this requirement and provide a fall back in case this
support is not present. Perhaps not using the nuclear EOS, or
reading from an ASCII file instead. This could be just a temporary
solution until the HDF5 library build script becomes more robust.
Opinions?
Thanks,
Bruno.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/145>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#365: ADMBase: evolving with evolution_method="static" does not apply boundary
conditions
----------------------+-----------------------------------------------------
Reporter: eschnett | Owner:
Type: defect | Status: new
Priority: major | Milestone:
Component: Cactus | Version:
Keywords: |
----------------------+-----------------------------------------------------
When one uses the evolution method ADMBASE::evolution_method = "static"
(which is very tempting e.g. for Cowling), then no boundary conditions are
applied. This means that the ADMBase variables remain unset on the
boundaries after regridding, leading to all sorts of (rather unexpected)
problems.
The easy solution is to use thorn Exact instead, which applies an exact
solution at each time step, including on the boundary.
I suggest to remove the ADMBase functionality to provide a "static"
evolution method, since this evolution method does not know which boundary
condition to apply, and hence is bound to fail in non-trivial situations.
Alternatively, ADMBase needs to synchronise in postregrid (and probably a
few other bins as well). ADMBase should then also offer to apply an outer
boundary condition, in this case probably a Minkowski Dirichlet boundary
condition or a von Neumann condition.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/365>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#410: CarpetIOHDF5 may segfault when one_file_per_group is selected
--------------------------------+-------------------------------------------
Reporter: bmundim | Owner: eschnett
Type: defect | Status: new
Priority: major | Milestone:
Component: Carpet | Version:
Keywords: Carpet git version |
--------------------------------+-------------------------------------------
I was facing one of those cryptic MPI error messages on Ranger,
where your job quits without any useful error message. It turns out
that this parameter file had
CarpetIOHDF5::one_file_per_group = "yes"
CarpetIOHDF5::out2D_vars = " ADMBase::gxx"
and CarpetIOHDF5 would try to loop over the group variables when
just one of them was selected. As a result it has accessed memory
address that it was not supposed to, consequently killing the job.
I have attached the stack back trace.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/410>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#431: Update some external libraries
-----------------------------------+----------------------------------------
Reporter: eschnett | Owner:
Type: enhancement | Status: new
Priority: minor | Milestone:
Component: EinsteinToolkit thorn | Version:
Keywords: |
-----------------------------------+----------------------------------------
Several external libraries have new versions available:
LAPACK
PETSc
curl
git
These are only minor updates. I suggest to apply them.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/431>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#136: Don't rebuild external libraries so often
-------------------------+--------------------------------------------------
Reporter: eschnett | Owner:
Type: enhancement | Status: new
Priority: minor | Milestone:
Component: Cactus | Version:
Keywords: |
-------------------------+--------------------------------------------------
One way to keep the external libraries that have been built would be the
following. Create a dummy configuration "ext" where all external libraries
are built. When another configuration is built, it should be easy to
specify to look there for the external libraries, or maybe this should
even be the default.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/136>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#222: Dubious code in Hydro_InitExcision.c
----------------------+-----------------------------------------------------
Reporter: eschnett | Owner:
Type: defect | Status: new
Priority: minor | Milestone:
Component: Cactus | Version:
Keywords: |
----------------------+-----------------------------------------------------
Hydro_InitExcision.c contains the following code:
if ( (hydro_initexcision_coordinate_length <= 0.0) &&
( ( x_frac > 0.5 - hydro_initexcision_fraction) &&
( x_frac < 0.5 + hydro_initexcision_fraction) &&
( y_frac > 0.5 - hydro_initexcision_fraction) &&
( y_frac < 0.5 + hydro_initexcision_fraction) &&
( z_frac > 0.5 - hydro_initexcision_fraction) &&
( z_frac < 0.5 + hydro_initexcision_fraction)
) ||
( (hydro_initexcision_coordinate_length > 0.0) &&
( fabs(x[point]-hydro_initexcision_position_x) <=
hydro_initexcision_coordinate_length*0.5) &&
( fabs(y[point]-hydro_initexcision_position_y) <=
hydro_initexcision_coordinate_length*0.5) &&
( fabs(z[point]-hydro_initexcision_position_z) <=
hydro_initexcision_coordinate_length*0.5)
)
)
This code has an "and" (&&) and an "or" (||) operation at top level. Is
this intended? The code would be clearer with an additional set of
parenthesis, or by introducing a suitable set of temporaries for sub-
expressions.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/222>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit