aboutsummaryrefslogtreecommitdiff
Commit message (Collapse)AuthorAge
* This commit changes the semantics of the out_of_range_tolerance[]jthorn2002-08-19
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | parameter-table entry to allow the min/max ends of the grid to be handled differently. (We'll need this for the multiprocessor interpolation.) ***** IMPORTANT ***** This change is *NOT* backwards-compatible: if you are calling CCTK_InterpLocalUniform() from this thorn (CactusBase/LocalInterp), and you have set out_of_range_tolerance[] to a non-default value, you must (slightly) change your code, or else you'll get an error return from the interpolator. If you don't use out_of_range_tolerance[] , then you don't need to do anything. The semantics of out_of_range_tolerance[] are now as follows: @var out_of_range_tolerance @vdesc Specifies how out-of-range interpolation points should be handled. The array elements are matched up with the axes and minimum/maximum ends of the grid in the order [x_min, x_max, y_min, y_max, z_min, z_max, ...]. An array value TOL is interpreted as follows: If TOL >= 0.0, then an interpolation point is considered to be "out of range" if and only if the interpolation point is > TOL * coord_delta[axis] outside the grid in this coordinate direction. If TOL == -1.0, then an interpolation point is considered to be "out of range" if and only if a centered molecule (or more generally, a molecule whose centering is chosen pretending that the grid is of infinite extent), would require data from outside the grid in this direction. Other values of TOL are illegal (reserved for future use). @vtype const CCTK_REAL out_of_range_tolerance[2*N_dims] @endvar The change is that it used to have only N_dims elements, one for each axis -- now there are 2*N_dims elements, one for each min/max end of each axis. ---------------------------------------------------------------------- Modified Files: doc/documentation.tex src/GeneralizedPolynomial-Uniform/InterpLocalUniform.c src/GeneralizedPolynomial-Uniform/InterpLocalUniform.h src/GeneralizedPolynomial-Uniform/interpolate.maple src/GeneralizedPolynomial-Uniform/molecule_posn.c src/GeneralizedPolynomial-Uniform/template.c src/GeneralizedPolynomial-Uniform/test_molecule_posn.c git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@88 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* add new directory with Lagrange interp coeffs -- I forgot to commitjthorn2002-08-19
| | | | | | | this file yesterday with my other commits. :( :( :( git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@87 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* fix syntax botch which confused origin compilerjthorn2002-08-19
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@86 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* expand description of naming conventionsjthorn2002-08-18
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@85 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* describe this directory and what files are in itjthorn2002-08-18
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@84 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* fix typo in commentjthorn2002-08-18
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@83 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* Another change in preparation for adding Hermite interpolation:jthorn2002-08-18
| | | | | | | | | | Change the documented name of the Lagrange interpolation operator from "generalized polynomial interpolation" to "Lagrange polynomial interpolation" (& document that the old name is still accepted for backwards compatability). git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@82 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* This commit reorganizes the Maple scripts which generate the interpolationjthorn2002-08-18
| | | | | | | | | | coefficients, and the coefficient files themselves, to properly support multiple interpolation operators, and adds some skeleton support for Hermite interpolation. I will add full support for Hermite in a future checkin. git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@81 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* new directory for files common to multiple interpolation operatorsjthorn2002-08-18
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@79 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* expand description a bitjthorn2002-07-18
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@78 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* add Maintainer: line as suggested by Erik Schnetter"s recent E-mailjthorn2002-07-16
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@77 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* change C++- (and C99-) style // comments to C-style /* ... */jthorn2002-07-16
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@76 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* slight reformatting of whitespace; no change in codejthorn2002-07-02
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@75 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* Remove a const qualifier that shouldn't have been there -- if you wantjthorn2002-07-02
| | | | | | | | | | to declare an int and then change the value later, it has to be an int , not a const int . Oops! Thanks to Peter Diener (and icc :) for spotting this! git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@74 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* fix a bug where the default value for the Jacobian-query parameter-tablejthorn2002-07-01
| | | | | | | | | | entry Jacobian_part_stride (= stride in Jacobian array for real/complex part) wasn't handled: docs have always said this is optional, and defaults to 0 if omitted from the parameter table, but the code treated this as an error --> fix to supply the default value (0) git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@73 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* change the way we handle NULL pointers for the input/output/Jacobianjthorn2002-07-01
| | | | | | | | arrays to give more flexibility: now you can do just a Jacobian query by passing null input/output pointers, git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@72 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* update to point to CactusTest/TestInterp for testsuitejthorn2002-06-24
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@71 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* fix a couple of typosjthorn2002-06-24
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@70 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* Added email addresses of author(s)/maintainer(s) of this thorn.tradke2002-06-07
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@68 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* small revisions to section 5.2 (notes on the implementation)jthorn2002-06-04
| | | | | | | | to improve the latex formatting and make the comments on spurious compiler warnings a bit clearer git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@67 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* fix spacing after i.e. and e.g.jthorn2002-06-04
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@66 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* convert to use the new thornguide template (looks great!)jthorn2002-06-04
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@65 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* template.cjthorn2002-06-03
| | | | | | | | | | | | | | | | | | | | | | | | | | | | * change comments to note that we never return UTIL_ERROR_NO_MEMORY InterpLocalUniform.c * add this kludge to handle malloc(0): > /* > * Implementation Note: > * > * We malloc() several scratch arrays, some with sizes determined by > * N_{input,output}_arrays. Thus if N_{input,output}_arrays == 0, with > * the obvious code we would malloc(0). Alas, the C standard permits > * malloc(0) to return a NULL pointer, which the usual malloc() idiom > * CCTK_INT *const p = malloc(N * sizeof(CCTK_INT)); > * if (p == NULL) > * then return UTIL_ERROR_NO_MEMORY > * would falsely detect as an out-of-memory condition. > * > * As a workaround, we pad all our malloc request sizes, i.e. > * CCTK_INT *const p = malloc((N+1) * sizeof(CCTK_INT)); > * if (p == NULL) > * then return UTIL_ERROR_NO_MEMORY > * This is a kludge, but so are the other plausible solutions. :( :( > */ git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@64 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* switch from $Id:$ to $Header:$jthorn2002-06-02
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@63 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* these things have now been fixed!jthorn2002-06-02
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@62 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* add $Header:$ tag for cvsjthorn2002-06-02
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@61 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* three changes in this checkin:jthorn2002-06-02
| | | | | | | | | | | | | | * change wording in error messages to refer to CCTK_InterpLocalUniform(), not InterpLocalUniform() * fix a couple of bugs where we tried to pass a CCTK_INT to a %d format in CCTK_VWarn -- CCTK_INT may be a short/long rather than an int, so this checkin changes the code to explicitly cast to int * NULL pointers in input_arrays[] used to cause a core dump; this checkin changes the code so they now cause that interpolation to be skipped (this is a useful feature IMHO) git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@60 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* * fix a bug in Jacobian querying where the code read (roughly)jthorn2002-05-31
| | | | | | | | | | | | | | | | | | | status = ... if (status == ...) { ... } else if (status2 == ...) /* should be status, not status2 */ { ... } else ... * change wording in error messages to refer to CCTK_InterpLocalUniform(), not InterpLocalUniform() git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@59 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* add $Header$jthorn2002-05-29
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@58 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* change $Id$ to $Header$jthorn2002-05-29
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@57 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* fix CCTK_FILEVERSION to say that were in CactusBase now, not CactusPUGHjthorn2002-05-29
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@56 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* remove an obselete FIXMEjthorn2002-05-28
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@55 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* slight wording tweaks to clarify a footnotejthorn2002-05-28
| | | | | | | and to properly describe part of the Jacobian querying git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@54 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* remove all latex macros so we only use plain-vanilla latex2ejthorn2002-05-28
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@53 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* add a test that output_arrays[out] != NULL before doing the interpolationjthorn2002-05-19
| | | | | | | | --> we can now selectively skip interpolation (eg if we just want a Jacobian) by passing a null pointer in the appropriate place git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@52 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* add Jacobian_info argument to grdoc commentjthorn2002-05-19
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@51 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* add a "guided tour" to the source codejthorn2002-05-19
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@50 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* this function used to be in ../src/GeneralizedPolynomial-Uniform/util.maplejthorn2002-05-19
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@49 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* remove unused function C_str() (now in ../../archive/)jthorn2002-05-19
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@48 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* tweak comment wordingjthorn2002-05-19
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@47 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* tweak code formatting in examplesjthorn2002-05-16
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@46 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* finish support for Jacobian queries (not properly tested yet)jthorn2002-05-16
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@45 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* regenerate machine-generated coeffsjthorn2002-05-16
| | | | | | | these now all have -kb (binary) so CVS doesn't try to merge changes to them git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@44 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* generate coeff-*.store.c files with non-I coefficients multiplied by "factor"jthorn2002-05-14
| | | | | | | --> needed to handle 1/dx factors for derivatives git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@43 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* modify code to store interpolation coefficients in COEFF(...) expressions,jthorn2002-05-14
| | | | | | | | to optionally multiply by "factor" to handle the case when we're doing derivs git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@42 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* tweak a commentjthorn2002-05-14
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@41 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* add infrastructure for Jacobian queryjthorn2002-05-14
| | | | | | | -- no actual Jacobian query code yet, though git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@40 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* * delete file which shouldn't have been therejthorn2002-05-14
| | | | | | | | (can't do dxx with size=2 molecules) * regenerate coeffs git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@39 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* regenerate all machine-generated coefficient filesjthorn2002-05-14
| | | | | | | | --> now includes stuff to store interpolation coeffs in COEFF(...) experessions --> will be used to compute Jacobian coeffs git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@38 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* * Don't destroy the entire [123]d.coeffs/ directory when we regenerate thejthorn2002-05-14
| | | | | | | | | | | coefficients, because that would confuse CVS. Instead, just remove any old coefficients. * provide separate makefile targets to create empty directories for coefficients; these are only needed for adding a new interpolator order or suchlike git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@37 df1f8a13-aa1d-4dd4-9681-27ded5b42416