<html>#2413: Include "POWER" code to extrapolate waveforms to scri+ into ET
<table style='border-spacing: 1ex 0pt; '>
<tr><td style='text-align:right'> Reporter:</td><td>Roland Haas</td></tr>
<tr><td style='text-align:right'> Status:</td><td>open</td></tr>
<tr><td style='text-align:right'>Milestone:</td><td></td></tr>
<tr><td style='text-align:right'> Version:</td><td></td></tr>
<tr><td style='text-align:right'> Type:</td><td>task</td></tr>
<tr><td style='text-align:right'> Priority:</td><td>major</td></tr>
<tr><td style='text-align:right'>Component:</td><td>EinsteinToolkit thorn</td></tr>
</table>
<p>Comment (by Roland Haas):</p>
<p>Right now using <code>FROM_TWOPUNCTURES</code> serves the purpose that is being suggested for <code>f0 = None</code> (<code>FROM_TWOPUNCTURES</code> is just a string and I only use the variable name to try and avoid typos as much as possible given that Python is fully dynamically typed and requires no variable declarations at all). I would rather not use <code>None</code> for this since there is only one <code>None</code> and one can get the total mass also in other ways which is easy to do by having another <code>FROM_XXX</code> but not if <code>None</code> is already taken.</p>
<p>You can use <code>power.py</code> as a Python module <code>import power</code> then call the Python functions directly. This means you have to write Python code rather than shell code.</p>
<p>I have code around somewhere that would parse the header and then picks the correct column. That should be easy to incorporate. </p>
<p>My first guesses for bad POWER output would be: </p>
<ul>
<li>make sure your psi4 starts and ends at 0 amplitude</li>
<li>make you may have to change the list of detectors being used, using detectors too far out or in can give bad results</li>
<li>TODO: probably there should be options to control the order of the polynomials that are being fit to the data</li>
</ul>
<p></p>
<p>--<br/>
Ticket URL: <a href='https://bitbucket.org/einsteintoolkit/tickets/issues/2413/include-power-code-to-extrapolate'>https://bitbucket.org/einsteintoolkit/tickets/issues/2413/include-power-code-to-extrapolate</a></p>
</html>