<div dir="ltr"><div>I have been trying a variety of initial guesses, but to no avail.</div><div><br></div><div>The setup is, in short:</div><div>par_b = 3.0</div><div>target_m_plus = 0.25</div><div>target_m_minus = 0.75</div><div><br></div><div>And the smaller black hole&#39;s apparent horizon gets detected easily, which cannot be said of the bigger one. My estimates for the AHFinderDirect were at -1.05 to -0.95 for the center of the ellipsoid/sphere, with initial guesses for the x_radius ranging from 0.25 to 0.4, with y_radius, z _radius about 10-20% smaller. <br></div><div><br></div><div>Am I still not getting enough resolution? The refinement levels in my finest setup around x=-1 are:</div><div><br></div><div>Carpet::max_refinement_levels    =8<br>CarpetRegrid2::movement_threshold_1 =   0.16<br>CarpetRegrid2::num_levels_2         =  8<br>CarpetRegrid2::position_x_2         = -1.0<br>CarpetRegrid2::radius_2[ 1]         =  64.0<br>CarpetRegrid2::radius_2[ 2]         =  16.0<br>CarpetRegrid2::radius_2[ 3]         =   8.0<br>CarpetRegrid2::radius_2[ 4]         =   4.0<br>CarpetRegrid2::radius_2[ 5]         =   3.0<br>CarpetRegrid2::radius_2[ 6]         =   2.5<br>CarpetRegrid2::radius_2[ 7]         =   2.0<br>CarpetRegrid2::movement_threshold_2 =   0.16<br></div><div><br></div><div>I do not know a-priori if this apparent horizon is best approximated by an ellipsoid (can one see this before performing the numerical simulation?) or not, but I don&#39;t know if I&#39;m on the right track. <br></div><div>It seems weird to me that I cannot find the x=-1 AH, even though the corresponding black hole it is ~3 times more massive. <br></div><div><br></div><div>In another attempt, I tried to copy the philosophy of the setup in the GW150914.rpar file, <br></div><div>setting the center_offset for the two punctures and other parameters, but that didn&#39;t work either.</div><div>I took</div><div> target_mp=0.75, target_mm=0.25, <br></div><div>separation D=6, so xp=D*mm = 1.5, xm=D*mp=-4.5, <br></div><div>par_b=D/2 = 3.0, center_offset=xp-half_D=-1.5<br></div><div>(all here: <a href="https://gist.github.com/konrad-topolski/1c208415d3edac2f428ded637153f639">https://gist.github.com/konrad-topolski/1c208415d3edac2f428ded637153f639</a>)</div><div>I&#39;m starting to doubt if any of these setups is correct. <br></div><div><br></div><div>By the way, can I somehow redirect the outputs of all MPI processes to one file at real time, better still, override the creation of CCTK_Proc files altogether?</div><div><br></div><div>Best regards</div><div>Konrad Topolski<br></div><div><br></div><div><br></div></div>