<div dir="auto"><div>Thank you! <div dir="auto"><br></div><div dir="auto">Of course, any help would be greatly appreciated. What I need to change in the code of the LeanBSSNMoL thorn?</div><div dir="auto"><br></div><div dir="auto">Erik</div><br><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Oct 24, 2019, 10:34 Miguel Zilhão <<a href="mailto:miguel.zilhao.nogueira@tecnico.ulisboa.pt" rel="noreferrer noreferrer" target="_blank">miguel.zilhao.nogueira@tecnico.ulisboa.pt</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">hi Erik,<br>
<br>
you can also try with the LeanBSSNMoL thorn, which is not Kranc-generated. adding another gauge <br>
condition should be fairly simple; i can point you exactly to what you'd need to change in the code <br>
if you need.<br>
<br>
cheers,<br>
Miguel<br>
<br>
On 24/10/19 16:21, Peter Diener wrote:<br>
> Hi Erik,<br>
> <br>
> McLachlan is a thorn generated from a tool named Kranc from high level tensorial expressions in <br>
> Mathematica. So if you have<br>
> access to Mathematica you can add additional lapse conditions in the<br>
> mathematica source file located in:<br>
> <br>
> arrangements/McLachlan/m/McLachlan_BSSN.m<br>
> <br>
> and then regenerate the McLachlan thorn using<br>
> <br>
> make McLachlan_BSSN.out<br>
> <br>
> in the arrangements/McLachlan/m directory.<br>
> <br>
> I admit the section on gauge conditions is rather messy as the selection of particular gauges are <br>
> done based on conditionals on thorn parameters, but hopefully you can figure this out. Otherwise let <br>
> me know and I'll be happy to help.<br>
> <br>
> Cheers,<br>
> <br>
> Peter<br>
> <br>
> On Thursday 2019-10-24 00:14, Erik Rodrígo Jiménez Vázquez wrote:<br>
> <br>
>> Date: Thu, 24 Oct 2019 00:14:05<br>
>> From: Erik Rodrígo Jiménez Vázquez <<a href="mailto:erjive@ciencias.unam.mx" rel="noreferrer noreferrer noreferrer" target="_blank">erjive@ciencias.unam.mx</a>><br>
>> To: <a href="mailto:users@einsteintoolkit.org" rel="noreferrer noreferrer noreferrer" target="_blank">users@einsteintoolkit.org</a><br>
>> Subject: [Users] Add a lapse condition<br>
>><br>
>> Dear all<br>
>><br>
>> Is there a way to add different gauge conditions for the lapse? I'm trying<br>
>> to add a gauge shock avoid lapse for the Bona-Masso family of slicing<br>
>> conditions<br>
>><br>
>> dot alpha = -alpha^2 f(alpha) K<br>
>><br>
>> with f(alpha)=1+constant/alpha^2<br>
>><br>
>> I don't know if I should add a thorn that complements the McLachlan thorn,<br>
>> or what is the correct procedure to implement the gauge condition?<br>
>> Furthermore, this kind of gauge condition allows to the lapse become<br>
>> negative which I know is a serious drawback when the lapse collapses to<br>
>> zero, however I want to do some numerical experiments with this kind of<br>
>> lapse.<br>
>><br>
>> Any help is very welcome. Thanks<br>
>> Erik<br>
>><br>
>><br>
>><br>
> <br>
> _______________________________________________<br>
> Users mailing list<br>
> <a href="mailto:Users@einsteintoolkit.org" rel="noreferrer noreferrer noreferrer" target="_blank">Users@einsteintoolkit.org</a><br>
> <a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" rel="noreferrer noreferrer noreferrer noreferrer" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
> <br>
</blockquote></div>
</div></div>