<div dir="auto"><div>Hi Krishiv,</div><div dir="auto"><br></div><div dir="auto">kuibit implements the fixed-frequency integration method to reconstruct the strain from Psi4. This algorithm requires a threshold frequency, or, equivalently, a threshold period. The error that you see tells you that the amount of data that you have is approximately smaller than this period, and the result would not be trustworthy.</div><div dir="auto"><br></div><div dir="auto">I typically choose the threshold period looking at the longest physical timescale in the system. In case of a binary, this is of the order of the first orbit. You can change that value, or run your simulation for longer. Both will allow you to extract a gravitational wave strain.</div><div dir="auto"><br></div><div dir="auto">Best,</div><div dir="auto">Gabriele<br><br><div class="gmail_quote" dir="auto"><div dir="ltr" class="gmail_attr">On Tue, Oct 12, 2021, 13:34 Erik Schnetter <<a href="mailto:schnetter@cct.lsu.edu" target="_blank" rel="noreferrer">schnetter@cct.lsu.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Krishiv<div><br></div><div>Unfortunately I don't know what causes this error message; I am not familiar with this tool. I'm sure someone else here on the mailing list knows!</div><div><br></div><div>-erik</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Oct 11, 2021 at 3:30 AM Krishiv Bhatia <<a href="mailto:flyingunicorns5724@gmail.com" rel="noreferrer noreferrer" target="_blank">flyingunicorns5724@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hello Dr. Schnetter,</div><div><br></div>I also ran trajectories and waveforms. While trajectories ran successfully and gave the snip below, waveform gave the following error:<div><br></div><div>
<p style="margin:0px;font-style:normal;font-variant-caps:normal;font-weight:normal;font-stretch:normal;font-size:12px;line-height:normal;font-family:"Andale Mono";color:rgb(47,255,18);background-color:rgba(0,0,0,0.9)"><span style="font-variant-ligatures:no-common-ligatures">$<span> </span>python3 waveform.py --datadir ~/simulations/GW150914_30/</span></p><p style="margin:0px;font-style:normal;font-variant-caps:normal;font-weight:normal;font-stretch:normal;font-size:12px;line-height:normal;font-family:"Andale Mono";color:rgb(47,255,18);background-color:rgba(0,0,0,0.9)"><span style="font-variant-ligatures:no-common-ligatures">Adm Mass: 0.989937</span></p><p style="margin:0px;font-style:normal;font-variant-caps:normal;font-weight:normal;font-stretch:normal;font-size:12px;line-height:normal;font-family:"Andale Mono";color:rgb(47,255,18);background-color:rgba(0,0,0,0.9)"><span style="font-variant-ligatures:no-common-ligatures">Radius: 500.0</span></p><p style="margin:0px;font-style:normal;font-variant-caps:normal;font-weight:normal;font-stretch:normal;font-size:12px;line-height:normal;font-family:"Andale Mono";color:rgb(47,255,18);background-color:rgba(0,0,0,0.9)"><span style="font-variant-ligatures:no-common-ligatures">Traceback (most recent call last):</span></p><p style="margin:0px;font-style:normal;font-variant-caps:normal;font-weight:normal;font-stretch:normal;font-size:12px;line-height:normal;font-family:"Andale Mono";color:rgb(47,255,18);background-color:rgba(0,0,0,0.9)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>File "/Users/<userdir>/simulations/Misc/waveform.py", line 144, in <module></span></p><p style="margin:0px;font-style:normal;font-variant-caps:normal;font-weight:normal;font-stretch:normal;font-size:12px;line-height:normal;font-family:"Andale Mono";color:rgb(47,255,18);background-color:rgba(0,0,0,0.9)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>strain = detector.get_strain_lm(</span></p><p style="margin:0px;font-style:normal;font-variant-caps:normal;font-weight:normal;font-stretch:normal;font-size:12px;line-height:normal;font-family:"Andale Mono";color:rgb(47,255,18);background-color:rgba(0,0,0,0.9)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>File "/Library/Frameworks/Python.framework/Versions/3.9/lib/python3.9/site-packages/kuibit/cactus_waves.py", line 247, in get_strain_lm</span></p><p style="margin:0px;font-style:normal;font-variant-caps:normal;font-weight:normal;font-stretch:normal;font-size:12px;line-height:normal;font-family:"Andale Mono";color:rgb(47,255,18);background-color:rgba(0,0,0,0.9)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>raise ValueError("pcut too large for timeseries")</span></p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:12px;line-height:normal;font-family:"Andale Mono";color:rgb(47,255,18);background-color:rgba(0,0,0,0.9)">
</p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:12px;line-height:normal;font-family:"Andale Mono";color:rgb(47,255,18);background-color:rgba(0,0,0,0.9)"><span style="font-variant-ligatures:no-common-ligatures">ValueError: pcut too large for timeseries</span></p></div><div><br></div><div><img src="https://mail.google.com/mail/?ui=2&ik=1caabb6aad&attid=0.2&th=17c765d91331fd4c&view=fimg&rm=17c765d91331fd4c&sz=w1600-h1000&attbid=ANGjdJ9mh2K5j1l03Gf9nbFOLLo2M4WBj-T_AwzHVS7jsLKfWfT7OrpW1-THSKw8Q5GapzRnvE6yCKCv8HEz8sDXQNsU_n3LaFV8xqtexiVdlVvCJTexF4MqOyGjuZo&disp=emb&realattid=ii_kumbyy4z6&zw" alt="image.png" width="562" height="458"><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Oct 10, 2021 at 11:54 PM Krishiv Bhatia <<a href="mailto:flyingunicorns5724@gmail.com" rel="noreferrer noreferrer" target="_blank">flyingunicorns5724@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div>Hello Dr. Schnetter,</div><div><br></div><div>My 12 hr simulation just now finished. Yes, as can be noted from the first snip, they did not go very far yet.</div><div>The second and third snips, which plot the trajectory of the 2 BHs individually, do show they going in a circular path. </div><div><br></div><div>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:12px;line-height:normal;font-family:"Andale Mono";color:rgb(47,255,18);background-color:rgba(0,0,0,0.9)"><span style="font-variant-ligatures:no-common-ligatures">gnuplot> plot 'puncturetracker-pt_loc..asc' u 23:33 title "BH 1", 'puncturetracker-pt_loc..asc' u 24:34 title "BH 2";</span></p></div><div><img src="https://mail.google.com/mail/?ui=2&ik=1caabb6aad&attid=0.1&th=17c765d91331fd4c&view=fimg&rm=17c765d91331fd4c&sz=w1600-h1000&attbid=ANGjdJ_LrGJUFpo4JJcNn9J4K3JzMq38qzi4Sfwfz_T2Kms1h1hchiFf5xiXHB40aofA5hfMBfAUgoPHBIxDVYtJb-ddFkcMDtWVSACBeiNhRyU81GaStDuspwdDTFQ&disp=emb&realattid=ii_kumao2ha3&zw" alt="image.png" width="562" height="551"><br></div><div><br></div><div>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:12px;line-height:normal;font-family:"Andale Mono";color:rgb(47,255,18);background-color:rgba(0,0,0,0.9)"><span style="font-variant-ligatures:no-common-ligatures">gnuplot> plot 'puncturetracker-pt_loc..asc' u 23:33 title "BH 1"</span></p></div><div><img src="https://mail.google.com/mail/?ui=2&ik=1caabb6aad&attid=0.6&th=17c765d91331fd4c&view=fimg&rm=17c765d91331fd4c&sz=w1600-h1000&attbid=ANGjdJ9XKrF-x13wEHZOsim5X_PC1zdI7rUAbNL7PEWbKutQS-lhjucXZF9RHgKWKpy5puXtw47ZyJIfdPTL20LUMClog0quCbGOcMBnKt09jmSe8HJbhkQkKbh2P8g&disp=emb&realattid=ii_kumap0ll4&zw" alt="image.png" width="562" height="541"><br></div><div><br></div><div>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:12px;line-height:normal;font-family:"Andale Mono";color:rgb(47,255,18);background-color:rgba(0,0,0,0.9)"><span style="font-variant-ligatures:no-common-ligatures">gnuplot> plot 'puncturetracker-pt_loc..asc' u 24:34 title "BH 2";</span></p></div><div><img src="https://mail.google.com/mail/?ui=2&ik=1caabb6aad&attid=0.4&th=17c765d91331fd4c&view=fimg&rm=17c765d91331fd4c&sz=w1600-h1000&attbid=ANGjdJ9ENobAe7kMhUFXMvGcHqCj57kYBRQbT99loidxXF8ZrCd68qmubVneJY00Bem9eIqO5CVsIj51RQ8SIfJRgc78lqnErY31RtonRsope4UtDkGeegaKaUzItbE&disp=emb&realattid=ii_kumaq2c55&zw" alt="image.png" width="562" height="540"><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Oct 10, 2021 at 6:14 PM Krishiv Bhatia <<a href="mailto:flyingunicorns5724@gmail.com" rel="noreferrer noreferrer" target="_blank">flyingunicorns5724@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr">Thanks Dr. Schnetter. I really value your time. <div><br></div><div>Good to know I am on the right path. I am not running on a cluster. It is my local mac. <div><br></div><div>I am also running a simulation for 10 hrs and then will run for more time. I will update you with the results. </div><div><br></div><div>- Krishiv</div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Oct 10, 2021 at 4:55 PM Erik Schnetter <<a href="mailto:schnetter@cct.lsu.edu" rel="noreferrer noreferrer" target="_blank">schnetter@cct.lsu.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr">Congratulations! These are the two black holes moving counterclockwise in the x-y plane. If you use an aspect ratio of 1:1, then you'll see that they haven't actually moved very far yet. You also see that the black holes start at rest (in terms of coordinate speed) because the gauge conditions in the initial conditions is suboptimal in puncture initial conditions. If you zoom in onto one of the black holes you'll see that they don't move in straight lines, but rather on a circle with a radius of about 6.<div><br></div><div>You should, of course, be able to run this setup on 128 cores instead of just 8, where it would run faster. This is likely a problem with either you *.ini, *.sub, or *.run file. If you tell us more about your cluster (maybe point to its documentation) and share the respective Simfactory files, then we might be able to help out. Unfortunately, many HPC systems are designed to be different in their details from other systems, and thus setting up correct Simfactory files for a new cluster is difficult. (As you've noticed.)</div><div><br></div><div>-erik</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Oct 10, 2021 at 3:58 PM Krishiv Bhatia <<a href="mailto:flyingunicorns5724@gmail.com" rel="noreferrer noreferrer" target="_blank">flyingunicorns5724@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr">In Jupyter, it crashed. But it ran on my mac. I ran it for 6 hrs. Got this:<div><br></div><div>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:12px;line-height:normal;font-family:"Andale Mono";color:rgb(47,255,18);background-color:rgba(0,0,0,0.9)"><span style="font-variant-ligatures:no-common-ligatures">> plot 'puncturetracker-pt_loc..asc' u 23:33 title "BH 1", 'puncturetracker-pt_loc..asc' u 24:34 title "BH 2";</span></p></div><div><br></div><div><img src="https://mail.google.com/mail/?ui=2&ik=1caabb6aad&attid=0.5&th=17c765d91331fd4c&view=fimg&rm=17c765d91331fd4c&sz=w1600-h1000&attbid=ANGjdJ-GE4ZGKnlznhXwN4yqmIQWWXWOXnuQVO4TXwocL5LXDEOivD7hRSzlL2t3kEyp0p-oTEC5wMuG7mVrJ7wXPSFcEGTKuCRgg0LllLCbt68dzrUK-bLujN89pec&disp=emb&realattid=ii_kulnf8712&zw" alt="image.png" width="562" height="499"><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Oct 9, 2021 at 3:46 PM Krishiv Bhatia <<a href="mailto:flyingunicorns5724@gmail.com" rel="noreferrer noreferrer" target="_blank">flyingunicorns5724@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr">Ah... it is working with procs = 8 instead of 128.<div>16 * 8 = 128. Seems to equate to node = 1<br><div><img src="https://mail.google.com/mail/?ui=2&ik=1caabb6aad&attid=0.3&th=17c765d91331fd4c&view=fimg&rm=17c765d91331fd4c&sz=w1600-h1000&attbid=ANGjdJ8AtwPQZsdSEbub7l59yNH2LdOdMj-F18Ajzkv0hoL-DGzrZ8mZeL7mhta620HeVNuHNDmvOe1cFB3r3K_tsAzQmHi4nDKKX56UIEJ5GFjUsy-jWwUtUeKLLBI&disp=emb&realattid=ii_kukdyask1&zw" alt="image.png" width="562" height="174"><br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Oct 9, 2021 at 1:15 PM Krishiv Bhatia <<a href="mailto:flyingunicorns5724@gmail.com" rel="noreferrer noreferrer" target="_blank">flyingunicorns5724@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>In the below mail, the first execution is in:</div><div> <a href="https://etk.cct.lsu.edu/user/flyingunicorns5724/notebooks/GW150914.ipynb" rel="noreferrer noreferrer" target="_blank">https://etk.cct.lsu.edu/user/flyingunicorns5724/notebooks/GW150914.ipynb</a></div><div><br></div><div>while the second execution is in my local mac.</div><div><br></div><div>- Krishiv</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Oct 9, 2021 at 1:08 PM Krishiv Bhatia <<a href="mailto:flyingunicorns5724@gmail.com" rel="noreferrer noreferrer" target="_blank">flyingunicorns5724@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hello,</div><div><br></div><div>I was running the commands in the following tutorial:</div><div> <a href="http://einsteintoolkit.org/gallery/bbh/cr.html" rel="noreferrer noreferrer" target="_blank">http://einsteintoolkit.org/gallery/bbh/cr.html</a><br></div><div><br></div><div>and got the following error message by running the below command.</div><div><br></div><div>Will appreciate your help.</div><div><br></div><div><img src="https://mail.google.com/mail/?ui=2&ik=1caabb6aad&attid=0.7&th=17c765d91331fd4c&view=fimg&rm=17c765d91331fd4c&sz=w1600-h1000&attbid=ANGjdJ9REr-e52EqBn4f_W4_amrIr-qtBsdUmifU1JY01y6xL4rUz9uQWnfVc_hx6eJwzl8SP60tSedFlyAPBNE5dq2bQFHSwdPShzso8TkeuFwbpyJ9DbCpyCCBgUw&disp=emb&realattid=ii_kuk8252c0&zw" alt="image.png" width="542" height="159"><br></div><div><br></div><div>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:12px;line-height:normal;font-family:"Andale Mono";color:rgb(47,255,18);background-color:rgba(0,0,0,0.9)"><span style="font-variant-ligatures:no-common-ligatures">$ simfactory/bin/sim create-submit GW150914_28 --define N 28 --parfile par/GW150914/GW150914.rpar --procs 128 --walltime 24:00:00</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:12px;line-height:normal;font-family:"Andale Mono";color:rgb(47,255,18);background-color:rgba(0,0,0,0.9)"><br></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:12px;line-height:normal;font-family:"Andale Mono";color:rgb(47,255,18);background-color:rgba(0,0,0,0.9)"><span style="font-variant-ligatures:no-common-ligatures">Assigned restart id: 0<span> </span></span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:12px;line-height:normal;font-family:"Andale Mono";color:rgb(47,255,18);background-color:rgba(0,0,0,0.9)"><span style="font-variant-ligatures:no-common-ligatures">Error: Too many nodes specified: nodes=16 (maxnodes is 1)</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:12px;line-height:normal;font-family:"Andale Mono";color:rgb(47,255,18);background-color:rgba(0,0,0,0.9)"><span style="font-variant-ligatures:no-common-ligatures">Aborting Simfactory.</span></p></div></div>
</blockquote></div></div>
</blockquote></div></div>
</blockquote></div>
_______________________________________________<br>
Users mailing list<br>
<a href="mailto:Users@einsteintoolkit.org" rel="noreferrer noreferrer" target="_blank">Users@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" rel="noreferrer noreferrer noreferrer" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div>Erik Schnetter <<a href="mailto:schnetter@cct.lsu.edu" rel="noreferrer noreferrer" target="_blank">schnetter@cct.lsu.edu</a>><br><a href="http://www.perimeterinstitute.ca/personal/eschnetter/" rel="noreferrer noreferrer" target="_blank">http://www.perimeterinstitute.ca/personal/eschnetter/</a></div><div><br></div></div></div>
</blockquote></div>
</blockquote></div></div>
</blockquote></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div>Erik Schnetter <<a href="mailto:schnetter@cct.lsu.edu" rel="noreferrer noreferrer" target="_blank">schnetter@cct.lsu.edu</a>><br><a href="http://www.perimeterinstitute.ca/personal/eschnetter/" rel="noreferrer noreferrer" target="_blank">http://www.perimeterinstitute.ca/personal/eschnetter/</a></div><div><br></div></div></div>
_______________________________________________<br>
Users mailing list<br>
<a href="mailto:Users@einsteintoolkit.org" rel="noreferrer noreferrer" target="_blank">Users@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" rel="noreferrer noreferrer noreferrer" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
</blockquote></div></div></div>