#945: add higher order restriction parameter to cell-centerd Carpet
-----------------------------------+----------------------------------------
Reporter: rhaas | Owner:
Type: enhancement | Status: new
Priority: minor | Milestone:
Component: EinsteinToolkit thorn | Version:
Keywords: CarpetLib |
-----------------------------------+----------------------------------------
the attached patch adds a 3rd order accurate restriction operator to
Carpet.
It is used only for cell-centered runs, for grid functions whose
prolongation operator is *not* ENO or WENO (ie. non matter variables
only).
It should be completely invisible if CarpetLib:use_cc_o3 is false (which
of course is the default).
It offers third-order accurate restriction operators for cell centered
grids when use_cc_o3 is set. This interpolation is done for samples at
the cell centres (so this is not a ppm scheme or anything like that).
It really only fits a polynomial of the form
\sum_{i,j,k=0}^3 a_{i,j,k} x^i y^j z^k
to the fine cells and evaluates at x=y=z=0. So it is good for the
metric, but bad for matter (since it will destroy the conservation).
Because of this it is not used for grid functions whose transport
operator is not WENO or ENO which hopefully excludes all matter
variables.
I also attach a patch for a modified WaveToyMoL thorn that I used for
testing. When test_restriction is set, then it puts a 3rd order
polynomial (with some random coefficients) into phi and the
differences between the restricted value and what should be there into
psi.
The wavetoymol thorn right now inherits from CarpetEvolutionMask since I
am also working on making this functional again (seems to be ok now).
Finally attached is a parameter file to test it. You should find a
region inside of reflevel 0 of psi where psi is exactly zero. This is
where the restriction works as expected.
Further improvements would be to not have the selection for 3rd order done
by a global carpetlib parameter but instead by a grid function tag (or a
global carpet paarameter). This change could also be made to the eno/weno
operators which contain special logic to use ENO for
prolongation_order_space = 3 and 5 both.
Ok to apply (or are there large scale changes to Carpet that are currently
private)?
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/945>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#434: Keep track of masked-out volume in CarpetMask
-------------------------+--------------------------------------------------
Reporter: eschnett | Owner: eschnett
Type: enhancement | Status: new
Priority: minor | Milestone:
Component: Carpet | Version:
Keywords: |
-------------------------+--------------------------------------------------
Keep track of the volume that is masked out by CarpetMask, and take this
volume into account when checking in CarpetReduce that the integral over
the simulation domain equals the domain volume.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/434>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#663: thornlists within Cactus should get cleaned up
-------------------------+--------------------------------------------------
Reporter: knarf | Owner:
Type: enhancement | Status: new
Priority: minor | Milestone:
Component: Cactus | Version: development version
Keywords: |
-------------------------+--------------------------------------------------
Most of the thornlists in
https://svn.cactuscode.org/Utilities/trunk/Thornlist*
are out of date, and should get cleaned up.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/663>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#878: Tests using Exact should use EinsteinExact
-----------------------------------+----------------------------------------
Reporter: hinder | Owner:
Type: enhancement | Status: new
Priority: major | Milestone:
Component: EinsteinToolkit thorn | Version:
Keywords: |
-----------------------------------+----------------------------------------
The EinsteinExact arrangement provides initial data accurate to roundoff,
in contrast to the Exact thorn which does not. The initial data from the
Exact thorn is also highly sensitive to roundoff level differences.
Hence, tests should be converted to using the EinsteinExact arrangement
where possible. In some cases, EinsteinExact will not support the
required metrics or parameters. These could either be added to
EinsteinExact, or the Exact thorn could continue to be used.
Milestone: ET_2012_11.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/878>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#361: delta_time when setting up ID with explicit time dependence
------------------------------------------+---------------------------------
Reporter: eloisa.bentivegna@… | Owner: eschnett
Type: defect | Status: new
Priority: minor | Milestone:
Component: Carpet | Version: ET_2010_11
Keywords: |
------------------------------------------+---------------------------------
In CallInitial, delta_time is used to calculate the time corresponding to
the different timelevels (line 399 of Carpet/src/Initialise.cc). At this
stage, though, delta_time is always equal to 1, leading to potentially
very separated initial-data slices when using init_each_timelevel and an
initial-data thorn that uses cctk_time explicitly. Should
cctkGH->cctk_delta_time be used here instead?
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/361>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#323: GetComponents overwrites thorn list
---------------------------+------------------------------------------------
Reporter: eschnett | Owner: eric9
Type: defect | Status: new
Priority: critical | Milestone:
Component: GetComponents | Version:
Keywords: |
---------------------------+------------------------------------------------
I specified a thorn list located in the main Cactus directory.
GetComponents overwrote this thorn list with its own thorn list (the one
it generates automatically).
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/323>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#631: Documentation fails to build
---------------------+------------------------------------------------------
Reporter: hinder | Owner:
Type: defect | Status: new
Priority: blocker | Milestone:
Component: Cactus | Version:
Keywords: |
---------------------+------------------------------------------------------
With the current ET, I cannot build the documentation:
[ianhin@login-damiana EinsteinToolkit]$ make AllDoc
Creating user documentation UsersGuide.pdf
Running pdflatex....
! Argument of \CA@getsvndate has an extra }.
! Paragraph ended before \CA@getsvndate was complete.
! Argument of \CA@getsvndate has an extra }.
! Paragraph ended before \CA@getsvndate was complete.
! Paragraph ended before \@parse@version was complete.
! Missing number, treated as zero.
! Missing = inserted for \ifnum.
! Missing number, treated as zero.
! LaTeX Error: Missing \begin{document}.
! Undefined control sequence.
Problem in . See doc/UsersGuide/LATEX_MESSAGES.
make: *** [UsersGuide] Error 1
I am attaching LATEX_MESSAGES. I have tried this both on damiana and on
my laptop. I consider this a blocker issue which could delay the release.
We must have working documentation.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/631>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#831: Cactus FAQ is out of date
----------------------------+-----------------------------------------------
Reporter: hinder | Owner:
Type: defect | Status: new
Priority: major | Milestone: ET_2012_05
Component: Cactus website | Version:
Keywords: |
----------------------------+-----------------------------------------------
The FAQ at http://cactuscode.org/documentation/faq is ancient and some of
its content is irrelevant now. It would be nice for this to be revised
before the next release.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/831>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#797: N13CarpetRegrid211snap_coarseE does not hold after being enforced in cell
centering
--------------------+-------------------------------------------------------
Reporter: rhaas | Owner: eschnett
Type: defect | Status: new
Priority: major | Milestone:
Component: Carpet | Version:
Keywords: |
--------------------+-------------------------------------------------------
I have a run (on Caltech's shc and on lonestar) which runs fine for about
300M or so and then fails with a property enforcement error inside of
CarpetRegrid.
Attached please find stdout with Carpet::veryverbose=yes,
CarpetRegrid2::veryverbose=yes as well as the parameter file.
It takes about 1hr to reproduce this starting from the checkpoints on shc.
Most of this time is IO (the checkpoint is iteration 133120 failure occurs
at preregrid in iteration 133185).
Turning off cell centering makes the error go away. Turning on refluxing
makes the run fails (much) earlier with a boxes not contained error (re-
running this on shc right now, the XSEDE clusters are just too slow).
In case you have access to shc (or zwicky) and would look at the actual
output files, the path is:
/panfs/ds06/sxs/rhaas/cactus/simulations/noreflux/output-0003
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/797>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#517: archive CactusExternal
-------------------------+--------------------------------------------------
Reporter: knarf | Owner:
Type: enhancement | Status: new
Priority: minor | Milestone:
Component: Cactus | Version:
Keywords: |
-------------------------+--------------------------------------------------
We should move the CactusExternal thorns to archive. They cause trouble
and confusion where they are right now. This should happen for the next
Cactus release.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/517>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit