<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 Gabriele Bozzola):</p>
<p>What I had in mind is that one can put <code>power.py</code> in their PATH and use it for all their simulations without having to touch the code. This requires the parameters to be adjustable via command-line. But if that is not the designed way of using the tool, that would be fine.</p>
<p>Also, for clarity, instead of defining a string <code>FROM_TWOPUNCTURES</code> that is going to be compared to something that is morally a float (<code>f0</code> and <code>ADMMass</code>), I would recommend setting the default of two parameters to <code>None</code> in the function definitions, and then checking if value is <code>None</code>. If not, one can assume that the values from TwoPunctures are used (which is what is happening now). So, something along the lines of</p>
<div class="codehilite"><pre><span></span><span class="k">if</span> <span class="n">f0</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">:</span>
<span class="n">f0</span> <span class="o">=</span> <span class="n">getCutoffFrequencyFromTwoPuncturesBBH</span><span class="p">(</span><span class="n">meta_name</span><span class="p">)</span>
</pre></div>
<p>(As someone that is not familiar with the code, I found the <code>FROM_TWOPUNCTURES</code> variable very confusing, but maybe its role will be clarified by the documentation.)</p>
<p>Two additional comments are:</p>
<ul>
<li>When the spin is computed from QLM a specific column is read. In my simulation, column 59 corresponds to <code>59:qlm_coordspinx[1]</code>, which is not the mass. Some assumptions are being made here and no checks are performed. This might lead to subtly wrong results.</li>
<li>It would be nice, and probably easy, to extend POWER to be able to compute the strain at a fixed radius (essentially, only doing the FFI). My use case is this. I tried POWER on my data and the output doesn’t look good at all. I don’t know if the problem is the interpolation to infinity or the data. Maybe the data becomes unreliable at a specific distance. If I can plot the strain
<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>