| Commit message (Collapse) | Author | Age |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
* support the same off-centering and boundary options as the new global
interpolator (which needs the support here to work): this means changing
the former parameter-table entry off_centering_tolerance to the
new entries
boundary_off_centering_tolerance
boundary_extrapolation_tolerance
* implement the N_boundary_points_to_omit parameter-table entry,
(previously this was documented but not implemented)
* completely rewrite the LocalInterp_molecule_posn() function
and its test driver to support the above changes
* change error return code from old CCTK_ERROR_INTERP_POINT_X_RANGE
to new CCTK_ERROR_INTERP_POINT_OUTSIDE to avoid confusion with
the old name suggesting that this was specific to the X coordinate
* fix documentation of interpolation operator names for older
(uniform Cartesian) CCTK_InterpLocal() interpolator -- this
closes Erik Schnetter's bug report CactusBase/1367
Documentation changes
* document the above code changes
* much expanded description of molecule centering and boundary handling,
including 2 new diagrams (done in latex picture mode :)
* update README files to list all code files
* move caching parameters (not implemented, and probably won't be
implemented any time soon) out of documentation.tex into a new file
future-ideas.tex
* general cleanup of documentation.tex (not quite done, but 80% or so...)
to more clearly describe what this thorn actually implements, rather
than a generic spec of which this thorn only implements some pieces
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@125 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@58 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
|
|
|
|
|
|
| |
1.0e4*DBL_EPSILON
--> this at least scales with machine precision,
abeit it does suffer from the serious problem that C double
might not be the right type to use (CCTK_REAL)
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@20 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@2 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|