| Commit message (Collapse) | Author | Age |
|
|
|
|
|
|
| |
Update the waveinterp test cases after correcting three timelevel
initialisation.
darcs-hash:20050811113402-891bb-b40f5401b44d626bda9611de2622a77395c16295.gz
|
|
|
|
| |
darcs-hash:20050810170009-891bb-168facf88ab10cc5c1732d037527a66bc134453c.gz
|
|
|
|
|
|
|
| |
Interpret the option table parameter time_deriv_order per output array
instead of per grid point.
darcs-hash:20050810165758-891bb-a4b84c90771da42c20878ebd59a797ef44a29b78.gz
|
|
|
|
| |
darcs-hash:20050810170107-891bb-bc5b8fa8a90d7289c0c7d76f3ad4a00aa3c1038f.gz
|
|
|
|
| |
darcs-hash:20050811170351-891bb-b7c586cb8494f23fac7254f169a1a2aeec523f97.gz
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
The waveinterp testsuite was broken, many output files had different values.
Also, the old output files didn't contain the recently introduced header comment
lines saying which thorn did output them.
The waveinterp testsuite has now been split into two separate tests,
waveinterp-1p and waveinterp-2p to run on 1 and 2 processors respectively.
The only differences are in the CarpetIOASCII files which contain either
1 or 2 components in the output.
darcs-hash:20050808164941-776a0-fbecb78b32c48a16f030996afcf1259276c0a0a0.gz
|
|
|
|
|
|
|
| |
Do not resize source_map twice.
Rename Util_TableGetIntArray return value from ierr to iret.
darcs-hash:20050801124730-891bb-4b25ee01ef6a3a21d75a19bc1500eec4ccaf310f.gz
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
M ./Carpet/CarpetInterp/src/interp.cc -105 +148
Added a couple of small classes InterpolationTimes and InterpolationWeights
that embody arrays calculated in the remaining monster function
interpolate_within_component, with the result that the formulas are
easier to read, and the logic within the monster function is a little
easier to follow.
More stuff like this can easily be done.
darcs-hash:20050718170031-90671-3ab51afe56043f082e26c7bb6f73c845ee2a70b4.gz
|
|
|
|
| |
darcs-hash:20050716140952-891bb-c7665e3f4b371287f54cfc167140719ce54120b8.gz
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
This is the first stage of a breakup of the monster function
Carpet_DriverInterpolate
The routine is segmented into several functional subroutines, with
names given (mostly) by the one-line comments that were in the
original.
Todo:
1) One of these new subroutines, DoLocalInterpolation, is itself a
monster function, and needs to be dealt with.
2) All of the subroutines have very long argument lists. There are
repeated patterns of arguments, which suggest that a structure could
be created which embodies the pattern.
darcs-hash:20050715103538-891bb-c0757adc59d29b09c63b7e4dbe11b57cd46b6682.gz
|
|
|
|
| |
darcs-hash:20050704170346-891bb-9afcde47badcd4fc165a3ea9d4cd73d6975482f3.gz
|
|
|
|
| |
darcs-hash:20050704163535-891bb-e95030a4787fd6b65dc17d8c959e714d36ac575c.gz
|
|
|
|
|
|
|
|
|
| |
Handle multiple maps. In singlemap mode, interpolate from the current
map. If there is only one map, interpolate from this one. Otherwise,
require a table argument "CCTK_INT source_map[N_interp_points]" that
specifies the source map for each interpolation point.
darcs-hash:20050704163505-891bb-653681ef3a619a7fa3462809c2bd408aae60284d.gz
|
|
|
|
|
|
|
|
|
|
| |
Add support for varying refinement factors. The spatial refinement
factors can be different in different directions, can be different
from the time refinement factor, and can be different on each level.
(However, the underlying spatial transport operators do currently not
handle any factors except two.)
darcs-hash:20050501205010-891bb-8d3a74abaad55ee6c77ef18d51fca2a2b69740de.gz
|
|
|
|
|
|
|
|
|
| |
data<> objects cannot be implicitly copied. The standard template
library containers copy objects arbitrarily. That means that one has
to store pointers to data<> objects instead of the objects themselves,
and has to allocate and free them manually.
darcs-hash:20050411170907-891bb-406b9b6bb6b97df2f47c349f32d91398338439ae.gz
|
|
|
|
| |
darcs-hash:20050307170026-891bb-03754477692ad245563fda22ecdd4510da4549ab.gz
|
|
|
|
| |
darcs-hash:20050303102506-891bb-118df8e38d1ca7e981bb4f97b14d5a747c1b3928.gz
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
Change the way in which the grid hierarchy is stored. The new hierarchy is
map
mglevel
reflevel
component
timelevel
i.e., mglevel moved from the bottom to almost the top. This is
because mglevel used to be a true multigrid level, but is now meant to
be a convergence level.
Do not allocate all storage all the time. Allow storage to be
switched on an off per refinement level (and for a single mglevel,
which prompted the change above). Handle storage management with
CCTK_{In,De}creaseGroupStorage instead of
CCTK_{En,Dis}ableGroupStorage.
darcs-hash:20050201225827-891bb-eae3b6bd092ae8d6b5e49be84c6f09f0e882933e.gz
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
Turn most of the templates in CarpetLib, which used to have the form
template<int D> class XXX
into classes, i.e., into something like
class XXX
by setting D to the new global integer constant dim, which in turn is set to 3.
The templates gf and data, which used to be of the form
template<typename T, int D> class XXX
are now of the form
template<typename T> class XXX
The templates vect, bbox, and bboxset remain templates.
This change simplifies the code somewhat.
darcs-hash:20050101182234-891bb-c3063528841f0d078b12cc506309ea27d8ce730d.gz
|
|
|
|
| |
darcs-hash:20050101162121-891bb-ac9d070faecc19f91b4b57389d3507bfc6c6e5ee.gz
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
Made all variables in CarpetLib:gh to be const or non-public, to avoid
confusion and insanity.
Toward this end, created four member accessors for gh:
const rexts & extents() const
const rbnds & outer_boundaries() const
const rprocs & processors() const
const vector<vector<ibbox> > & bases() const
This involved a few changes in several files throughout the repository:
M ./Carpet/Carpet/src/Initialise.cc -2 +2
M ./Carpet/Carpet/src/Recompose.cc -8 +8
M ./Carpet/Carpet/src/SetupGH.cc -2 +2
M ./Carpet/Carpet/src/modes.cc -3 +3
M ./Carpet/CarpetIOHDF5/src/Recover.cc -4 +4
M ./Carpet/CarpetInterp/src/interp.cc -5 +5
M ./Carpet/CarpetLib/src/dh.cc -3 +3
M ./Carpet/CarpetLib/src/gh.cc -33 +33
M ./Carpet/CarpetLib/src/gh.hh -10 +28
M ./Carpet/CarpetReduce/src/mask_carpet.cc -4 +4
M ./Carpet/CarpetRegrid/src/automatic.cc -2 +2
darcs-hash:20041210131459-32473-ff0835ff0e57f6693fbe1ec23b350d1b1d18e3bc.gz
|
|
|
|
| |
darcs-hash:20041118233344-891bb-44b65f16949b0d78bee0a66ad6f65c262652c671.gz
|
|
|
|
|
|
|
| |
Return the local interpolators' statuses in a table entry
local_interpolator_status. This is e.g. required by AHFinderDirect.
darcs-hash:20041020092737-891bb-495add6874e22d32e23eacbc340493a1e4656da3.gz
|
|
|
|
| |
darcs-hash:20041008214254-891bb-222fb022121e1d3b6f18ca8176640d400039a3d4.gz
|
|
|
|
| |
darcs-hash:20041002183708-891bb-96a5a3227becd0e40834fb9ab25ef5f9d75c1df5.gz
|
|
|
|
| |
darcs-hash:20040922203859-891bb-7e628f331129d9afeded4ea54732c37ca48cf26c.gz
|
|
|
|
|
|
|
| |
Add interpolation option key "time_interpolation_order" that allows
taking time derivatives while interpolating.
darcs-hash:20040921160917-891bb-236dc322c4bacdd2fa8d13b5c535541640e02537.gz
|
|
|
|
|
|
|
|
|
| |
Add an interpolation option "interpolation_times" that specifies for
each grid point at what time the value should be calculated. This
allows for some interpolation in time when there are sufficient time
levels.
darcs-hash:20040918132542-891bb-887b348b0e510d7f17ac32b41cb78176079ff1fb.gz
|
|
|
|
|
|
| |
Replace all CVS header tags with the standard "$Header:$".
darcs-hash:20040918132147-891bb-dea889bdd94a479ec412d14d08e9efca63e5c24d.gz
|
|
|
|
| |
darcs-hash:20040913123613-07bb3-060a47d1900cb82aa8260be38049e1e4b70f4e50.gz
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
Add a new tags table entry to check: InterpNumTimelevels.
When this is set then the interpolator will take the given number of
timelevels as the number to use when interpolating in time. So if you
have grid functions with fewer than 3 timelevels mixed with grid
functions with 3 timelevels you can use this tag, set it to the number
of timelevels that you have available, and the interpolator will "do
the right thing".
Note that this means that GFs with only one timelevel are implicitly
treated as being constant in time.
The obvious examples, soon to be committed (I hope) are the static
conformal factor and the coordinates (although why you'd want to
interpolate the coordinates...). E.g.,
REAL confac TYPE = GF timelevels = 1 tags='tensortypealias="Scalar"
Prolongation="None" InterpNumTimelevels=1'
...
darcs-hash:20040820072138-58737-5ff68c487af14a1d72228aa6697cc9eea1ae3a89.gz
|
|
|
|
|
|
|
| |
Correct error in calculating the components' coordinate extents when
convergence_level is not zero.
darcs-hash:20040802163927-07bb3-639b6b816f790ab12935ccb4052c7c8808b7a6ff.gz
|
|
|
|
|
|
|
| |
Ensure that the array pointers, if given, are not null.
Allow zero interpolation points.
darcs-hash:20040802093953-07bb3-2af122bfdc454e59d31b4a7f73cc44f59c9bb9d7.gz
|
|
|
|
| |
darcs-hash:20040702131619-58737-b02419fc8b4f4514ab1bb31f341054785bef0a21.gz
|
|
|
|
|
|
|
|
| |
Correct a bug in the time interpolation. We probably never
interpolated in time after all, and if so, then with a sign error.
This correction is by analogy from CarpetReduce and not tested.
darcs-hash:20040627191732-07bb3-326b9ecadca08013a8a463f5ce1a7883b0e49bac.gz
|
|
|
|
|
|
|
| |
Use "requires thorns" to ensure that the thorns are activated in the
correct order.
darcs-hash:20040608190925-07bb3-70b0ac8ee8e570c0b7fff50b02395a288482257d.gz
|
|
|
|
| |
darcs-hash:20040528115004-07bb3-dfa6908bd400b1096f8406de3b7e4cdd1283ecfc.gz
|
|
|
|
| |
darcs-hash:20040521161223-07bb3-6cc7ddc3777345de379eb54daf2cf6d532925c23.gz
|
|
|
|
|
|
|
|
| |
Undo interface/capability update, because now Carpet conflicts with
PUGH. Need to sort things out, maybe switch temporarily back to the
"requires thorns" mechanism.
darcs-hash:20040504205239-07bb3-771c6ec96a2b454d63c87a6dc7a1f44f62cd8fd7.gz
|
|
|
|
|
|
|
|
| |
Use the local instead of a registered coordinate system for
interpolation. This will enable patch-local multi-patch
interpolation.
darcs-hash:20040504201254-07bb3-2c00ac4fb8a00d211289814abc892908d4c78a42.gz
|
|
|
|
| |
darcs-hash:20040504200954-07bb3-cd9f17c981ce37247fa2122858b10b46425f678a.gz
|
|
|
|
| |
darcs-hash:20040418120959-07bb3-05050a71e475171cac6b293de35118e4af8df5ba.gz
|
|
|
|
|
|
|
| |
When a grid point is not in any component, map it to an arbitrary
component instead of aborting.
darcs-hash:20040416185756-07bb3-9faf77dee8876c69e83780620702df6d2c660bec.gz
|
|
|
|
| |
darcs-hash:20040331143832-07bb3-7a1e263661a7034f7ad6d38b048ce79015d4f20d.gz
|
|
|
|
| |
darcs-hash:20040331105058-07bb3-0c0495c75cd248152699703e2e6d1520c9c80052.gz
|
|
|
|
| |
darcs-hash:20040308081244-07bb3-353fed99bc85bc4ba6c8aa6f6c9739a31831ffa3.gz
|
|
|
|
| |
darcs-hash:20040227152205-07bb3-33135e23f97ccb01f347189ca57b01003732863d.gz
|
|
|
|
|
|
|
|
|
|
| |
When interpolating in time, call the local interpolator many times
instead of only once. This is slower, but it avoids having to rewrite
all input arguments -- and there are many input arguments, most of
them hidden in the options table. In fact, Carpet can not even know
how many input arguments there are.
darcs-hash:20040215101516-07bb3-fac5be0b7f5557733f81b66cf70cea86e79f67ec.gz
|
|
|
|
| |
darcs-hash:20040215093414-07bb3-1b9b2a4fd94a77be876fef35aa15a2eef7d55941.gz
|
|
|
|
|
|
|
|
| |
Use (rl-minrl) instead of rl to index into the count arrays (that are
used to communicate which processor works on how many grid points on
what levels).
darcs-hash:20040214140929-07bb3-e0305bc753a075f79d1d476c22b5773e7e6edd88.gz
|