| Commit message (Collapse) | Author | Age |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
| |
this file yesterday with my other commits. :( :( :(
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@87 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@86 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@85 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@84 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@83 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@79 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@78 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@77 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@76 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@75 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@71 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@70 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@68 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@66 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@65 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
* 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
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@63 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@62 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@61 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
* 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
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@58 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@57 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@56 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@55 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
|
|
|
| |
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
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@53 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
|
|
|
|
| |
--> 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
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@51 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@50 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@49 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@48 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@47 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@46 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@45 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
| |
--> 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
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@41 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
|
|
|
| |
-- no actual Jacobian query code yet, though
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@40 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
|
|
|
|
| |
(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
|
|
|
|
|
|
|
|
| |
--> 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
|
|
|
|
|
|
|
|
|
|
|
| |
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
|