aboutsummaryrefslogtreecommitdiff
path: root/src
Commit message (Collapse)AuthorAge
...
* This is a major reworking of this interpolator.jthorn2002-09-02
| | | | | | | | | The interpolation code in template.c is now split up into subfunctions (in common/, Lagrange/, and Hermite/), so it should now compile in finite cpu/memory even with optimization. git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@114 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* This is a major reworking of this interpolator.jthorn2002-09-02
| | | | | | | | | The interpolation code in template.c is now split up into subfunctions (in common/, Lagrange/, and Hermite/), so it should now compile in finite cpu/memory even with optimization. git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@113 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* more coefficients being renamedjthorn2002-09-01
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@112 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* coefficients being redone for finite compilation size/timejthorn2002-09-01
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@111 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* there are now .c routines in this directory!jthorn2002-09-01
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@110 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* reorganize common to Lagrange and Hermite interpolatorsjthorn2002-09-01
| | | | | | | | | (i.e. code depending only on molecule size/shape) --> now there are separate *functions* rather than just code fragments --> ../template.c should hopefully compile in finite time..... git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@109 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* these are all being replaced as part of the "finite compilation"jthorn2002-09-01
| | | | | | | restructuring git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@108 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* extend naming convention to include molecule shapejthorn2002-09-01
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@107 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* data and coefficient structures for use in splitting up the interpolatorjthorn2002-09-01
| | | | | | | into smaller functions (which I hope will be compilable in finite cpu/memory) git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@106 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* add Hermite interpolation directory to the list of subdirsjthorn2002-08-28
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@105 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* describe the files in this directoryjthorn2002-08-28
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@104 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* add Hermite interpolation functions to dispatch tablejthorn2002-08-28
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@103 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* add Maple code and C skeletons for Hermite polynomial interpolationjthorn2002-08-28
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@101 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* In Hermite_polynomial_interpolant(),jthorn2002-08-28
| | | | | | | | add explicit check for the case when the given conditions don't uniquely determine the interpolant git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@100 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* redo Hermite interpolant computation to hopefully be more efficientjthorn2002-08-28
| | | | | | | | (in practice it doesn't seem to make much difference, though :( :(), and to be a bit easier to read/understand/modify the code git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@99 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* add some more new coefficients needed for Hermite interpolationjthorn2002-08-20
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@98 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* refresh 3d coeffs again (I forgot some last time)jthorn2002-08-20
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@97 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* refresh machine-generated coefficient filesjthorn2002-08-20
| | | | | | | with latest tweaks to Maple code which generates them git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@96 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* switch from $Id:$ to $Header:$jthorn2002-08-20
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@95 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* add new files which will be (are) needed by Hermite interpolatorjthorn2002-08-20
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@94 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* switch from $Id:$ to $Header:$jthorn2002-08-20
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@93 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* add $Header:$ lines, switch a few $Id:$ --> $Header:$jthorn2002-08-20
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@92 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* fix grdoc comment spelling typojthorn2002-08-19
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@91 df1f8a13-aa1d-4dd4-9681-27ded5b42416
* 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
* 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
* 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
* 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
* 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
* 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