aboutsummaryrefslogtreecommitdiff
path: root/Carpet/CarpetInterp
Commit message (Collapse)AuthorAge
* CarpetInterp: Update waveinterp test casesErik Schnetter2005-08-11
| | | | | | | Update the waveinterp test cases after correcting three timelevel initialisation. darcs-hash:20050811113402-891bb-b40f5401b44d626bda9611de2622a77395c16295.gz
* CarpetInterp: Use explicit source map specifications even in single map modeErik Schnetter2005-08-10
| | | | darcs-hash:20050810170009-891bb-168facf88ab10cc5c1732d037527a66bc134453c.gz
* CarpetInterp: Make time_deriv_order per output array instead of per grid pointErik Schnetter2005-08-10
| | | | | | | Interpret the option table parameter time_deriv_order per output array instead of per grid point. darcs-hash:20050810165758-891bb-a4b84c90771da42c20878ebd59a797ef44a29b78.gz
* CarpetInterp: Remove commented out codeErik Schnetter2005-08-10
| | | | darcs-hash:20050810170107-891bb-bc5b8fa8a90d7289c0c7d76f3ad4a00aa3c1038f.gz
* CarpetInterp: Correct time interpolation stencilsErik Schnetter2005-08-11
| | | | darcs-hash:20050811170351-891bb-b7c586cb8494f23fac7254f169a1a2aeec523f97.gz
* CarpetInterp: fix waveinterp testsuiteThomas Radke2005-08-08
| | | | | | | | | | | | | 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
* CarpetInterp: Rename return valueErik Schnetter2005-08-01
| | | | | | | Do not resize source_map twice. Rename Util_TableGetIntArray return value from ierr to iret. darcs-hash:20050801124730-891bb-4b25ee01ef6a3a21d75a19bc1500eec4ccaf310f.gz
* CarpetInterp: modularize time interpolationSteve White2005-07-18
| | | | | | | | | | | | | | 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
* CarpetInterp: Plug memory leak in the recent function breakupErik Schnetter2005-07-16
| | | | darcs-hash:20050716140952-891bb-c7665e3f4b371287f54cfc167140719ce54120b8.gz
* CarpetInterp: break routine into piecesErik Schnetter2005-07-15
| | | | | | | | | | | | | | | | | | | | 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
* CarpetInterp: Correct error in calculating the past levels' current timesErik Schnetter2005-07-04
| | | | darcs-hash:20050704170346-891bb-9afcde47badcd4fc165a3ea9d4cd73d6975482f3.gz
* CarpetInterp: Correct error in calculating the past levels' current timesErik Schnetter2005-07-04
| | | | darcs-hash:20050704163535-891bb-e95030a4787fd6b65dc17d8c959e714d36ac575c.gz
* CarpetInterp: Handle multiple mapsErik Schnetter2005-07-04
| | | | | | | | | 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
* global: Add varying refinement factorsErik Schnetter2005-05-01
| | | | | | | | | | 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
* CarpetInterp: Correct error in allocating data objectsErik Schnetter2005-04-11
| | | | | | | | | 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
* CarpetInterp: Correct check whether a processor owns a componentErik Schnetter2005-03-07
| | | | darcs-hash:20050307170026-891bb-03754477692ad245563fda22ecdd4510da4549ab.gz
* CarpetInterp: Handle groups that have not all time levels allocatedErik Schnetter2005-03-03
| | | | darcs-hash:20050303102506-891bb-118df8e38d1ca7e981bb4f97b14d5a747c1b3928.gz
* global: Change the way in which the grid hierarchy is storedErik Schnetter2005-02-01
| | | | | | | | | | | | | | | | | | | | 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
* global: Turn CarpetLib templates into classesErik Schnetter2005-01-01
| | | | | | | | | | | | | | | | | | | 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
* global: Remove $Header$ keywords from Carpet arrangementErik Schnetter2005-01-01
| | | | darcs-hash:20050101162121-891bb-ac9d070faecc19f91b4b57389d3507bfc6c6e5ee.gz
* CarpetLib:gh exposed data clamp-downswhite2004-12-10
| | | | | | | | | | | | | | | | | | | | | | | | | | | 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
* Do not use CVS $Date$ for latex documentation any moreErik Schnetter2004-11-18
| | | | darcs-hash:20041118233344-891bb-44b65f16949b0d78bee0a66ad6f65c262652c671.gz
* CarpetInterp: return local_interpolator_statusErik Schnetter2004-10-20
| | | | | | | 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
* CarpetInterp: clean up time interpolationErik Schnetter2004-10-08
| | | | darcs-hash:20041008214254-891bb-222fb022121e1d3b6f18ca8176640d400039a3d4.gz
* CarpetInterp: Fix segfault when determining interpolation timesErik Schnetter2004-10-02
| | | | darcs-hash:20041002183708-891bb-96a5a3227becd0e40834fb9ab25ef5f9d75c1df5.gz
* Whitespace change in CarpetInterpErik Schnetter2004-09-22
| | | | darcs-hash:20040922203859-891bb-7e628f331129d9afeded4ea54732c37ca48cf26c.gz
* Implement time derivatives during interpolationErik Schnetter2004-09-21
| | | | | | | Add interpolation option key "time_interpolation_order" that allows taking time derivatives while interpolating. darcs-hash:20040921160917-891bb-236dc322c4bacdd2fa8d13b5c535541640e02537.gz
* Add true interpolation in timeErik Schnetter2004-09-18
| | | | | | | | | 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
* Remove CVS header tagsErik Schnetter2004-09-18
| | | | | | Replace all CVS header tags with the standard "$Header:$". darcs-hash:20040918132147-891bb-dea889bdd94a479ec412d14d08e9efca63e5c24d.gz
* Map output indices to input indices when interpolating in time.schnetter2004-09-13
| | | | darcs-hash:20040913123613-07bb3-060a47d1900cb82aa8260be38049e1e4b70f4e50.gz
* Add a new tags table entry to check: InterpNumTimelevels.hawke2004-08-20
| | | | | | | | | | | | | | | | | | | | | | | | 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 whenschnetter2004-08-02
| | | | | | | 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.schnetter2004-08-02
| | | | | | | Ensure that the array pointers, if given, are not null. Allow zero interpolation points. darcs-hash:20040802093953-07bb3-2af122bfdc454e59d31b4a7f73cc44f59c9bb9d7.gz
* A more useful error message if time interpolation can't be done.hawke2004-07-02
| | | | darcs-hash:20040702131619-58737-b02419fc8b4f4514ab1bb31f341054785bef0a21.gz
* Correct a bug in the time interpolation. We probably neverschnetter2004-06-27
| | | | | | | | 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 theschnetter2004-06-08
| | | | | | | Use "requires thorns" to ensure that the thorns are activated in the correct order. darcs-hash:20040608190925-07bb3-70b0ac8ee8e570c0b7fff50b02395a288482257d.gz
* Remove restrict qualifier.schnetter2004-05-28
| | | | darcs-hash:20040528115004-07bb3-dfa6908bd400b1096f8406de3b7e4cdd1283ecfc.gz
* Add support for symmetry interpolation.schnetter2004-05-21
| | | | darcs-hash:20040521161223-07bb3-6cc7ddc3777345de379eb54daf2cf6d532925c23.gz
* Undo interface/capability update, because now Carpet conflicts withschnetter2004-05-04
| | | | | | | | 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 forschnetter2004-05-04
| | | | | | | | 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
* Begin to clean up inheritance and requires relationshipsschnetter2004-05-04
| | | | darcs-hash:20040504200954-07bb3-cd9f17c981ce37247fa2122858b10b46425f678a.gz
* Use "requires" instead of "requires thorns".schnetter2004-04-18
| | | | darcs-hash:20040418120959-07bb3-05050a71e475171cac6b293de35118e4af8df5ba.gz
* When a grid point is not in any component, map it to an arbitraryschnetter2004-04-16
| | | | | | | When a grid point is not in any component, map it to an arbitrary component instead of aborting. darcs-hash:20040416185756-07bb3-9faf77dee8876c69e83780620702df6d2c660bec.gz
* Explicitly cast constness awayschnetter2004-03-31
| | | | darcs-hash:20040331143832-07bb3-7a1e263661a7034f7ad6d38b048ce79015d4f20d.gz
* Handle looping over refinement levels more elegantly.schnetter2004-03-31
| | | | darcs-hash:20040331105058-07bb3-0c0495c75cd248152699703e2e6d1520c9c80052.gz
* Change operator[] to at()schnetter2004-03-08
| | | | darcs-hash:20040308081244-07bb3-353fed99bc85bc4ba6c8aa6f6c9739a31831ffa3.gz
* Use the new vect::floor instead of vect::map(floor,)schnetter2004-02-27
| | | | darcs-hash:20040227152205-07bb3-33135e23f97ccb01f347189ca57b01003732863d.gz
* When interpolating in time, call the local interpolator many timesschnetter2004-02-15
| | | | | | | | | | 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
* Fix array size allocation.schnetter2004-02-15
| | | | darcs-hash:20040215093414-07bb3-1b9b2a4fd94a77be876fef35aa15a2eef7d55941.gz
* Use (rl-minrl) instead of rl to index into the count arrays (that areschnetter2004-02-14
| | | | | | | | 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