#939: NullEvol_InitialData.F90 triggers Fortran exception
-----------------------------------+----------------------------------------
Reporter: rhaas | Owner:
Type: defect | Status: new
Priority: major | Milestone:
Component: EinsteinToolkit thorn | Version:
Keywords: NullEvolve |
-----------------------------------+----------------------------------------
I ran the SphericalHarmonicReconn testsuite (regression_test.par) with a
Lovelace executable compiled via simfactory build --debug and get Fortran
out-of bounds exceptions from NullEvolve when running the test (2
processes 2 threads each):
{{{
INFO (NullInterp): the guard point shell is set for a max. stencil size
of 3
INFO (NullSHRExtract): CCE world-tube index range is 4
4 20.404040404040405 20.404040404040405
INFO (NullEvolve): Null Initial Data
At line 269 of file
/mnt/data/rhaas/postdoc/gr/ET_Lovelace/arrangements/PITTNullCode/NullEvolve/src/NullEvol_InitialData.F90
Fortran runtime error: Array bound mismatch for dimension 1 of array
'jcn_rad' (15/31)
At line 269 of file
/mnt/data/rhaas/postdoc/gr/ET_Lovelace/arrangements/PITTNullCode/NullEvolve/src/NullEvol_InitialData.F90
Fortran runtime error: Array bound mismatch for dimension 1 of array
'jcn_rad' (16/31)
}}}
This happens when compiling with gcc 4.6.3 but does not seem to happen
when using the intel compiler (though I am not quite sure how to turn on
all the bounds checks in it).
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/939>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#927: Enable Simfactory to use CACTUS_CONFIGS_DIR
------------------------+---------------------------------------------------
Reporter: sbrandt | Owner: eschnett
Type: defect | Status: new
Priority: major | Milestone:
Component: SimFactory | Version:
Keywords: |
------------------------+---------------------------------------------------
This simple patch enables Simfactory to use CACTUS_CONFIGS_DIR
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/927>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#894: add support for CarpetInterp2 to InterpToArray
-----------------------------------+----------------------------------------
Reporter: rhaas | Owner:
Type: enhancement | Status: new
Priority: minor | Milestone:
Component: EinsteinToolkit thorn | Version:
Keywords: |
-----------------------------------+----------------------------------------
the attached patch adds support for CarpetInterp2 to interptoarray. For
after the release.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/894>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#938: out of bounds array access in thorn ADM
-----------------------------------+----------------------------------------
Reporter: rhaas | Owner:
Type: defect | Status: new
Priority: minor | Milestone:
Component: EinsteinToolkit thorn | Version:
Keywords: ADM |
-----------------------------------+----------------------------------------
I just ran the Lovelace test with an executable build via sim build
--debug. I get failing testsuites that did not fails with a normal build.
One is ADM/test_ADM_2 which fails with (tow threads per process, two
processes):
{{{
INFO (IDAnalyticBH): setting up Schwarzschild initial
INFO (IOBasic): Periodic scalar output requested for 'ADMANALYSIS::grr',
'ADMBASE::gxx', 'ADMBASE::kxx', 'ADMBASE::alp', 'ADMCONSTRAINTS::ham',
'ADMCONSTRAINTS::momx', 'ADMCONSTRAINTS::momy', 'ADMCONSTRAINTS::momz'
INFO (ADM): 1+log slicing
INFO (ADM): in the kleing style
INFO (ADM): but not using the initial lapse coeff
INFO (ADM): flat space diffusion .2000000E-04
INFO (ADM): and no curved space diffusion
INFO (ADM): Initializing Leapfrog with FTCS Step
Using a corrector step
At line 23 of file
/mnt/data/rhaas/postdoc/gr/ET_Lovelace/arrangements/CactusArchive/ADM/src/DoubleLeap.F
Fortran runtime error: Index '0' of dimension 1 of array 'betax' below
lower bound of 1
}}}
The line number is completely wrong since line 23 is not even inside a
function yet but recompiling without F_LINE_DIRECTIVES gives line 4553 of
file /mnt/data/rhaas/postdoc/gr/ET_Lovelace/configs/sim-
debug/build/ADM/DoubleLeap.f (after preprocessing). The code is
{{{
do k=1,nz
do j=1,ny
do i=1,nx
c The right-hand-side for the (conformal) metric evolution equation
cdgdt_cdgxxdt = - 2*alp(i,j,k)*kxx(i,j,k)
cdgdt_cdgxydt = - 2*alp(i,j,k)*kxy(i,j,k)
cdgdt_cdgxzdt = - 2*alp(i,j,k)*kxz(i,j,k)
cdgdt_cdgyydt = - 2*alp(i,j,k)*kyy(i,j,k)
cdgdt_cdgyzdt = - 2*alp(i,j,k)*kyz(i,j,k)
cdgdt_cdgzzdt = - 2*alp(i,j,k)*kzz(i,j,k)
IF (conformal_state .gt. 0) THEN
cdgdt_ipsi4 = 1D0/(psi(i,j,k)**4)
cdgdt_cdgxxdt = cdgdt_cdgxxdt*cdgdt_ipsi4
cdgdt_cdgxydt = cdgdt_cdgxydt*cdgdt_ipsi4
cdgdt_cdgxzdt = cdgdt_cdgxzdt*cdgdt_ipsi4
cdgdt_cdgyydt = cdgdt_cdgyydt*cdgdt_ipsi4
cdgdt_cdgyzdt = cdgdt_cdgyzdt*cdgdt_ipsi4
cdgdt_cdgzzdt = cdgdt_cdgzzdt*cdgdt_ipsi4
END IF
IF (shift_state .ne. 0) THEN
if (local_spatial_order.eq.2) then
dxdb_dxdbx = I2DX*(betax(i+1,j,k) - betax(i-1,j,k))
dxdb_dxdby = I2DX*(betay(i+1,j,k) - betay(i-1,j,k))
dxdb_dxdbz = I2DX*(betaz(i+1,j,k) - betaz(i-1,j,k))
else
}}}
which comes from line 300 of the unprocessed file. Clearly the limits of
the loop should be
{{{
do k=2,nz-1
do j=2,ny-1
do i=2,nx-1
}}}
instead of {{{k=1,nz}}} etc.
ok to apply (given that we want to deprecate ADM anyway)? The change has
(surprisingly) no affect on the actual generated test output.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/938>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#930: ADMBase timelevel parameters should be steerable
-----------------------------------+----------------------------------------
Reporter: eschnett | Owner:
Type: defect | Status: new
Priority: major | Milestone:
Component: EinsteinToolkit thorn | Version:
Keywords: |
-----------------------------------+----------------------------------------
It makes e.g. sense to reduce the number of timelevels upon recovery, if
e.g. a thorn that requires more timelevels is not active any more, or if
that thorn's parameters have been changed and the thorn now requires fewer
timelevels, or if some of the timelevels were never needed to begin with.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/930>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#932: Test system should output details of the compiler used
-------------------------+--------------------------------------------------
Reporter: hinder | Owner:
Type: enhancement | Status: new
Priority: minor | Milestone:
Component: Cactus | Version:
Keywords: |
-------------------------+--------------------------------------------------
Certain compiler versions (builds) might cause certain tests to fail.
Keeping track of this would be easier if the test system output the build
of the compiler used to summary.log. Compiler options might also be
useful to have here.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/932>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#850: New thorns ML_WaveToy and ML_WaveToy_CL
-----------------------------------+----------------------------------------
Reporter: eschnett | Owner:
Type: enhancement | Status: new
Priority: major | Milestone:
Component: EinsteinToolkit thorn | Version:
Keywords: |
-----------------------------------+----------------------------------------
I suggest to add two new thorns to the Einstein Toolkit: ML_WaveToy and
ML_WaveToy_CL. Both are generated by Kranc and distributed as part of
McLachlan. The former is a standard WaveToy example, showing how a Kranc
script looks like that is much simpler than that for the BSSN equations.
The latter is essentially the same script, but generates OpenCL code,
demonstrating the Cactus OpenCL capabilities.
The Kranc script is WaveToy.m, and is already distributed as part of the
standard McLachlan distribution.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/850>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#805: Add thorn OpenCLRunTime to Einstein Toolkit
----------------------+-----------------------------------------------------
Reporter: eschnett | Owner:
Type: task | Status: new
Priority: major | Milestone:
Component: Other | Version:
Keywords: |
----------------------+-----------------------------------------------------
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/805>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#926: ADMBase variable initialization not thread-safe
--------------------+-------------------------------------------------------
Reporter: knarf | Owner:
Type: defect | Status: new
Priority: major | Milestone: ET_2012_05
Component: Other | Version: development version
Keywords: |
--------------------+-------------------------------------------------------
The attached patch makes the loop variable loop-private, making sure this
works when using OpenMP.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/926>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#924: Improve sync performance
-------------------------+--------------------------------------------------
Reporter: hinder | Owner: eschnett
Type: enhancement | Status: new
Priority: minor | Milestone:
Component: SimFactory | Version:
Keywords: |
-------------------------+--------------------------------------------------
Running "sim sync" to a slow remote filesystem can take a long time due to
having to check the time stamp and size of all the remote files. It would
be good to speed this up.
One option would be to assume that files are only edited on the local
machine, which is a common situation. Each invocation of sim sync to a
given machine would record the time at which it completed, and subsequent
syncs would only transfer those files which had changed on the local
system since then.
A further optimisation would be to avoid statting all the local files by
using a filesystem notification API such as fsevents or inotify to
determine the required information. An advanced case would also use such
a system on the remote side, so you wouldn't have to assume that the
remote files were unchanged.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/924>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit