aboutsummaryrefslogtreecommitdiff
path: root/src
Commit message (Collapse)AuthorAge
* implement UIUC's speed-up in "evaluation" of spectral solutionrhaas2013-06-11
| | | | | | | | | | | | | | | | | | | | | | | | | | The accompanying svn patch applies the recently released refinement of TwoPunctures by Vasileios Paschalidis and Zach Etienne at UIUC, greatly reducing the time taken to properly apply the solution of the spectral solve to all grid points. From the abstract of the notes in arXiv:1304.0457: TwoPunctures is perhaps the most widely-adopted code for generating binary black hole "puncture" initial data and interpolating these (spectral) data onto evolution grids. In typical usage, the bulk of this code's run time is spent in its spectral interpolation routine. We announce a new publicly-available spectral interpolation routine that improves the performance of the original interpolation routine by a factor of ~100, yielding results consistent with the original spectral interpolation routine to roundoff precision. This note serves as a guide for installing this routine both in the original standalone TwoPunctures code and the Einstein Toolkit supported version of this code. Patch kindly provided by Bernard Kelly Original code by Vasileios Paschalidis and Zach Etienne git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@133 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Declare all loop variables privateeschnett2012-05-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@127 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Add output of metadata file TwoPunctures.bbhhinder2011-09-25
| | | | | | | | The format and keys are defined in http://arxiv.org/abs/0709.0093 with additional draft keys defined for the NR-AR project (to be committed to the arXiv version when finalized). Also introduced grid scalars for the ADM energy and angular momentum and the puncture ADM masses. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@124 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* remove compiler warnings (declare functions, rename x1,y1,z1 (because of ↵knarf2011-01-10
| | | | | | bessel functions in math libraries), initialize variables) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@119 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* patch from Bernard Kelly:knarf2011-01-10
| | | | | | | Output the individual puncture ADM mass post-solve in TwoPunctures git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@118 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* remove (leftover) permanent change of the grid functions x y and zknarf2010-06-24
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@110 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* TwoPunctures: fix bug in d3tensorrhaas2010-05-25
| | | | | | | | | Had hard coded the lower bound to be zero and forgotten to initialize the first column of pointers git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@108 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* make PunctEvalAtArbitPosition work with more than one variablerhaas2010-05-24
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@105 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* TwoPunctures: improve search for ADM masses of black holesrhaas2010-05-16
| | | | | | | | | | | | | | | | | | | | | | | | | Incorporates the additions by Barry Wardell: 1) When the bare masses are not given as parameters (give_bare_mass=no), they are calculated at double precision, but the par_m_plus and par_m_minus parameters are then only set at single precision. 2) When the bare masses have been calculated, they are printed to stdout. In some cases, the initial guess for the bare mass is printed rather than the final calculated bare mass. 3) After looking at this section of code further, it seems that the termination criterion for the bare mass search could be improved. Attached is an updated patch which implements and improved check. It now continues trying to find the bare masses until the ADM masses are within a specified tolerance of their target. Other than the change in termination criterion, the algorithm is unchanged, although I have simplified the calculation of a new bare mass guess somewhat. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@103 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* TwoPunctures: remove NR memory management and complex arithmetic coderhaas2010-04-29
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@101 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* remove some code from numerical recipes and replace by code using the gnurhaas2010-03-24
| | | | | | | scientific library git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@100 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* a bit more output (b), make par_b steerableknarf2010-03-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@99 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* partial openMP support, one loop remains serialized, because it would have ↵knarf2009-10-20
| | | | | | to be rewritten to make use of openMP git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@97 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* output residuum and initial guess debug information only on processor 0tradke2009-02-26
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@91 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Added missing type declarations.pollney2008-11-27
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@89 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Erik's patch to improve parallelisation via OpenMP.pollney2008-11-26
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@88 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* This is much slower, but much more accurate and does not give diffrent ↵jese2008-05-07
| | | | | | values after the last post-mass-solving call to the Newton solver. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@86 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* yay! it works now. thanks, erik.jese2008-03-27
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@84 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Version which works for runs and calculates bare_mass, but may not properly ↵jese2008-03-26
| | | | | | set calculated bare masses. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@83 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Set the puncture masses on the refined levels as well. This shouldschnetter2008-02-29
| | | | | | | remove the nans that appeared recently. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@82 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Use snprintf instead of sprintfschnetter2008-02-29
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@81 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Added parameters boolean "give_bare_mass" and plus and minus target ADM ↵jese2008-02-19
| | | | | | masses for bare mass solver. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@80 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Added center_offset parameter to allow the initial data BHs to bepollney2006-10-18
| | | | | | | shifted around the grid. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@70 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* remove compiler warning about oo many arguments in printfknarf2006-08-29
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@69 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Added initial lapse profile "brownsville" accordint to the specificationpollney2006-08-28
| | | | | | | in Phys. Rev. D 74, 041501 (2006). git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@67 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* make swap() void because it does not return something, the return value is ↵knarf2006-08-28
| | | | | | never used and the compiler warns about this git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@66 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Prevent name clash for the macro SWAP. SWAP was already defined inschnetter2006-08-13
| | | | | | | TP_utilities.h, but unused. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@65 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Add a parameter "swap_xz" which swaps the x and z coordinates duringschnetter2006-08-12
| | | | | | | | interpolation. This makes it possible to set up binary systems which are separated along the z axis. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@64 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* start to prepare to use with more than 1 equationknarf2006-08-11
| | | | | | | nothing changes for the old case git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@63 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* again also interpolate into ghost zones until someone finds why it does not ↵knarf2006-08-10
| | | | | | work always git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@61 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Correct an error in the loop boundaries of the main loop. Instead ofschnetter2006-08-10
| | | | | | | cctk_nghostzones[d], the code contained cctk_nghostzones[2]. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@60 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* - add a parameter to choose when to schedule the IDknarf2006-08-09
| | | | | | | | | | | | | | - add a parameter for multiplying an already computed lapse to the one computed here This does not change the default behaviour. However, the testsuite fails - but it also failed before this change in the same way. Does anybody know anything about this? -> It fails at exactly one point in the 1D-slices, but only in y and z directions git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@59 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Use CCTK_REAL instead of double. This allows using higher precisionsschnetter2006-07-27
| | | | | | | | | | | that double. Do not initialise the ghost zones; synchronise instead. Since interpolating the initial data is expensive this should save some time. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@57 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Print info about the solver in the initial bin, not in the paramcheckschnetter2006-05-07
| | | | | | | | bin. The paramcheck bin is executed while recovering, and it is confusing to see solver messages while recovering. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@56 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Use the correct power of TP_epsilon for the extrinsic curvatureschnetter2006-05-07
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@55 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Added new intial lapse profile, 'psi^n' corresponding to the initialpollney2006-04-03
| | | | | | | | conformal factor raised to some exponent, given by the new parameter initial_lapse_psi_exponent. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@54 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Smooth out Aij as wellschnetter2006-03-22
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@53 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Introduce a parameter TP_epsilon that can be used to smooth outschnetter2006-03-01
| | | | | | | | | | | | initial data. It has the same meaning as for the Brill-Lindquist initial data in thorn AEIThorns/Exact, namely by r -> (r^4 + epsilon^4)^(1/4) By default, epsilon=0, which introduces no change. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@51 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* a few smaller fixesknarf2005-12-02
| | | | | | | | allow an extended region blended instead of the puncture % done - info (locally, but who cares ...) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@48 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* Do not define asinh, acosh, and atanh. The compiler usually alreadyschnetter2005-11-21
| | | | | | | | | | knows these, and defining them again can lead to errors. If these are now undefined on some systems, then we have to re-introduce them with modified names. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@47 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* new parameter to be able to avoid nans at or near the puncture locationsknarf2005-05-13
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@44 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* initialise variables correctlyknarf2005-05-11
| | | | | | | This did only show up in valgrind, but should be fixed anyway git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@43 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* pass NULL in case there is no conformal factorknarf2005-05-02
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@42 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* C89 compliantknarf2005-04-28
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@41 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* transform // style comments into /* */ style comments because some compilersknarf2005-02-26
| | | | | | | | | do not understand C99 yet, unfortunately :( no other change done git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@40 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* fix a compile bugknarf2005-02-23
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@39 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* insert a commented out (#if 0) section to cut very small spectral coefficientsknarf2005-02-23
| | | | | | | to 0. Does nothing if not enabled git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@38 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* do a lot more initial debug output - does nothing if not enabledknarf2005-02-23
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@37 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* improve/change debug output formatknarf2005-02-23
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@36 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
* bugfixknarf2005-02-16
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@35 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
generated by cgit v1.2.3 (git 2.39.1) at 2024-05-17 04:22:15 +0000