aboutsummaryrefslogtreecommitdiff
path: root/src/findnormals.F90
Commit message (Collapse)AuthorAge
* Add parameters to set the origin of the excision region for "fixedschnetter2004-02-14
| | | | | | | | | | | | | | | excision". When checking for (real-valued) mask values, do not check for equality, but instead allow for a fudge factor. Check that the thorn is activated when any of its routines are called. When there is only a single point excised, do not abort with an internal error, but pretend the normal extends into the x direction. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/LegoExcision/trunk@42 f75ba9e5-694f-0410-ac2c-87ea7ce7132b
* Include cctk_Functions for the DECLARE_CCTK_FUNCTIONS macros.hawke2003-10-28
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/LegoExcision/trunk@39 f75ba9e5-694f-0410-ac2c-87ea7ce7132b
* Use the generic routines sqrt(), abs(), log(), etc. instead of double-precisionpollney2002-06-18
| | | | | | | | | ones dsqrt(), dabs(), dlog(), detc. This makes the code compile in single precision. (Fix courtesy of ThomasR) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/LegoExcision/trunk@33 f75ba9e5-694f-0410-ac2c-87ea7ce7132b
* Changed the way in which the normals where found. Before, I just chosemiguel2001-04-27
| | | | | | | | | | that grid point closest to the optimal normal. But it turned out that in some cases this grid point was in fact masked. I now choose the closest unmasked grid point to the optimal normal. This means a second loop over nearest neighbours, but it can't be avoided. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/LegoExcision/trunk@20 f75ba9e5-694f-0410-ac2c-87ea7ce7132b
* Oooops! bug. Fixed it (divisions by zero).miguel2001-04-03
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/LegoExcision/trunk@18 f75ba9e5-694f-0410-ac2c-87ea7ce7132b
* Changed the way in which the normal was found so that when the "real"miguel2001-04-03
| | | | | | | | normal vector is close to a diagonal, we get in fact the diagnal and not a vector along the coordinate lines. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/LegoExcision/trunk@17 f75ba9e5-694f-0410-ac2c-87ea7ce7132b
* New findboundary routine without templates. It seems to work, but requiresmiguel2001-04-02
| | | | | | | more testing yet. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/LegoExcision/trunk@14 f75ba9e5-694f-0410-ac2c-87ea7ce7132b
* Debugging.miguel2001-03-30
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/LegoExcision/trunk@13 f75ba9e5-694f-0410-ac2c-87ea7ce7132b
* Fixing small bug.miguel2001-03-30
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/LegoExcision/trunk@12 f75ba9e5-694f-0410-ac2c-87ea7ce7132b
* Taking out ccthghmiguel2001-03-29
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/LegoExcision/trunk@9 f75ba9e5-694f-0410-ac2c-87ea7ce7132b
* Removed the original file that handled extrapolation directly.miguel2001-03-29
There are now three files for three purposes: findboundary takes a mask with only 0s and 1s and marks the boundary with 0.5s. findnormals takes a mask where the boundary has been marked and returns the normal directions to the mask for these locations. extrapolate takes a mask the normals and extrapolates a grid function. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/LegoExcision/trunk@6 f75ba9e5-694f-0410-ac2c-87ea7ce7132b