| Commit message (Collapse) | Author | Age |
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
| |
(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
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@93 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@83 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
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
* Add code to print C assignments of the form
COEFF(-1,1) = coeff_dx_m1_p1;
etc
* change all functions which print to name files to now explicit close
the files when they're done -- previously we just left the files open,
and relied on Maple's implicit close-on-exit :(
* reformat whitespace in [123]d.maple
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@36 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
|
|
| |
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@35 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|
|
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@2 df1f8a13-aa1d-4dd4-9681-27ded5b42416
|