Users July 2022

users@lists.einsteintoolkit.org

9 participants
16 discussions

Meeting minutes for 2022-07-28
by Bing-Jyun Tsao 28 Jul '22

28 Jul '22

Present: Peter, Beyhan, Bruno, Helvi, Keith, Zach, local people from the workshop (Leo, Roland, Liwei, Federico, Alexandru, Steve, Johnny, Cheng-Hsin) Chair: Leo Minutes: Johnny - modules suggested for inclusion for next release [RH] - Canuda scalar and complex fields ticket <https://bitbucket.org/einsteintoolkit/tickets/issues/2485/include-complex-a…> (Miguel Zilhao): Currently adapted two punctures, scalar field, evolution. Will need review, possibly someone who planned to use Canuda [RH]. Compile without error in a few weeks, and will need to be done by Sep 15. ScalarBase, ScalarInit, ScalarEvolve will all be needed to include in ETK, will need documentation for the whole arrangement, refer the variables to the equations in the papers. - FLRW solver ticket <https://bitbucket.org/einsteintoolkit/tickets/issues/2549/inlcude-flrwsolve…> (Hayley Macpherson): Helvi was originally suggested but will be busy. Some students planning to use the code were expressing interests. *Zach mentioned he can review the code. - SelfForce-1D ticket <https://bitbucket.org/einsteintoolkit/tickets/issues/2602/include-regge-whe…> (Samuel Cupp): - NRPyEllipticET ticket <https://bitbucket.org/einsteintoolkit/tickets/issues/2616/add-nrpyelliptice…> (Leo Werneck, Thiago Assumpcao, Zach Etienne): Students in Helvi’s group will be reviewers. - Elliptica ticket <https://bitbucket.org/einsteintoolkit/tickets/issues/2619/include-ellipitca…> (Alireza Rashti): Beyhan volunteered to review, contact Alireza directly. Should ask if reader exists. - timeline <https://docs.einsteintoolkit.org/et-docs/Release_Details> for upcoming release [RH] - Proposed thorn in the master branch, all thorns and selfForce-1D [PD]. In two weeks, start testing codes on different clusters. Release in November. - missing updates for gallery examples [RH] - BBH <https://www.einsteintoolkit.org/gallery/bbh/index.html> [SB]: will be updated on the website - TOV <https://www.einsteintoolkit.org/gallery/ns/index.html> [PD]: the run is done, will need to update webpage - Tickets: - ticket #2624 <https://bitbucket.org/einsteintoolkit/tickets/issues/2624/not-possible-to-r…>: Run with extra environment variables . - ticket #2620 <https://bitbucket.org/einsteintoolkit/tickets/issues/2620/wvuthorns-baikal-…>: latest commit of Baikal fixes. Different orders are compiled at different files instead one file [ZE], would need review or more tests*. *Helvi can review Baikal, Zach can review FLRW - For the working workshop next year: If the workshop becomes only one week, maybe a one-week virtual workshop to get people to set up the environment and compile the code first [RH]. CarpetX will need review in the future (next year) [HV], maybe a review workshop [SB] - ET meeting in Dublin <https://sites.google.com/view/eetm2022> in a month Next call: Chair: Peter Minutes: Steve

1 0

Einstein Toolkit Meeting Reminder
by reminders＠einsteintoolkit.org 27 Jul '22

27 Jul '22

Hello, Please consider joining the weekly Einstein Toolkit phone call at 8:30 am US central time on Thursdays. For details on how to connect and what agenda items are to be discussed, use the link below. https://docs.einsteintoolkit.org/et-docs/Main_Page#Weekly_Users_Call --The Maintainers

1 0

Re: [Users] Users Digest, Vol 148, Issue 12
by Luciano Combi 26 Jul '22

26 Jul '22

Excellent, Miguel, that worked! Cheers. El dom, 24 jul 2022 a las 14:19, <users-request(a)einsteintoolkit.org> escribió: > Send Users mailing list submissions to > users(a)einsteintoolkit.org > > To subscribe or unsubscribe via the World Wide Web, visit > http://lists.einsteintoolkit.org/mailman/listinfo/users > or, via email, send a message with subject or body 'help' to > users-request(a)einsteintoolkit.org > > You can reach the person managing the list at > users-owner(a)einsteintoolkit.org > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of Users digest..." > > > Today's Topics: > > 1. Re: AHFinderDirect test in cosmology (Miguel Zilhão) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sat, 23 Jul 2022 23:26:28 +0100 > From: Miguel Zilhão <mzilhao(a)ua.pt> > Subject: Re: [Users] AHFinderDirect test in cosmology > To: Luciano Combi <combi.luciano(a)gmail.com>, > <users(a)einsteintoolkit.org> > Message-ID: <4aea2a6a-777e-eeb7-6613-59489eeb999f(a)ua.pt> > Content-Type: text/plain; charset="UTF-8"; format=flowed > > hi Luciano, > > we've done something like this in https://arxiv.org/abs/1204.2019 > you can indeed find cosmological apparent horizons using AHFinderDirect, > the trick is to add the > following parameter: > > AHFinderDirect::surface_definition[1] = "inner expansion" > > hope this helps. > > cheers, > Miguel > > On 23/07/2022 06:48, Luciano Combi wrote: > > Hi all, > > > > I'm trying to use AHFinderDirect to find the apparent horizon in the de > Sitter metric, implemented > > in the Exact thorn in comoving coordinates. I tried a few things but > cannot get it to work. Have > > someone tried something similar, i.e. using the AHFD in cosmological > spacetimes instead of BHs? > > Below I attach a typical parfile I used > > > > Thanks! > > ------------------------------------------------------- > > # flesh > > cactus::cctk_itlast = 0 > > > > ActiveThorns = "PUGH" > > driver::ghost_size = 3 > > driver::global_nx = 120 > > driver::global_ny = 120 > > driver::global_nz = 120 ##57 > > > > ActiveThorns = "CoordBase SymBase CartGrid3D" > > grid::avoid_origin = "false" > > grid::type = "byspacing" > > grid::dxyz = 0.2 > > > > ActiveThorns = "ADMBase ADMCoupling StaticConformal CoordGauge SpaceMask > Exact" > > ADMBase::initial_lapse = "exact" > > ADMBase::initial_shift = "exact" > > ADMBase::initial_data = "exact" > > ADMBase::lapse_evolution_method = "exact" > > ADMBase::shift_evolution_method = "exact" > > ADMBase::metric_type = "physical" > > Exact::exact_model = "de Sitter+Lambda" > > Exact::de_Sitter_Lambda__scale = 0.50 > > > > > > ######################################## > > > > ActiveThorns = "IOUtil" > > IOUtil::parfile_write = "no" > > > > ######################################## > > > > ActiveThorns = "AEILocalInterp PUGHInterp PUGHReduce SphericalSurface > AHFinderDirect" > > > > AHFinderDirect::print_timing_stats = "true" > > > > AHFinderDirect::output_Theta_every = 1 > > AHFinderDirect::h_base_file_name = "dS.h" > > AHFinderDirect::Theta_base_file_name = "dS.Theta" > > > > AHFinderDirect::N_horizons = 1 > > AHFinderDirect::origin_x[1] = 0.0 > > AHFinderDirect::origin_y[1] = 0.0 > > AHFinderDirect::origin_z[1] = 0.0 > > > > AHFinderDirect::initial_guess_method[1] = "coordinate sphere" > > AHFinderDirect::initial_guess__coord_sphere__x_center[1] = 0.0 > > AHFinderDirect::initial_guess__coord_sphere__y_center[1] = 0.0 > > AHFinderDirect::initial_guess__coord_sphere__z_center[1] = 0.0 > > AHFinderDirect::initial_guess__coord_sphere__radius[1] = 1.0 > > ------------------------ > > -- > > *Luc Combi* > > CONICET PhD fellow (Buenos Aires, Argentina) > > Instituto Argentino de Radioastronomía > > > > _______________________________________________ > > Users mailing list > > Users(a)einsteintoolkit.org > > http://lists.einsteintoolkit.org/mailman/listinfo/users > > > ------------------------------ > > _______________________________________________ > Users mailing list > Users(a)einsteintoolkit.org > http://lists.einsteintoolkit.org/mailman/listinfo/users > > > End of Users Digest, Vol 148, Issue 12 > ************************************** > -- *Luc Combi* CONICET PhD fellow (Buenos Aires, Argentina) Instituto Argentino de Radioastronomía

1 0

Agenda for Thursday's Meeting
by reminders＠einsteintoolkit.org 25 Jul '22

25 Jul '22

Please update the Wiki with agenda items for Thursday's meeting. Thanks! https://docs.einsteintoolkit.org/et-docs/meeting_agenda --The Maintainers

1 0

AHFinderDirect test in cosmology
by Luciano Combi 23 Jul '22

23 Jul '22

Hi all, I'm trying to use AHFinderDirect to find the apparent horizon in the de Sitter metric, implemented in the Exact thorn in comoving coordinates. I tried a few things but cannot get it to work. Have someone tried something similar, i.e. using the AHFD in cosmological spacetimes instead of BHs? Below I attach a typical parfile I used Thanks! ------------------------------------------------------- # flesh cactus::cctk_itlast = 0 ActiveThorns = "PUGH" driver::ghost_size = 3 driver::global_nx = 120 driver::global_ny = 120 driver::global_nz = 120 ##57 ActiveThorns = "CoordBase SymBase CartGrid3D" grid::avoid_origin = "false" grid::type = "byspacing" grid::dxyz = 0.2 ActiveThorns = "ADMBase ADMCoupling StaticConformal CoordGauge SpaceMask Exact" ADMBase::initial_lapse = "exact" ADMBase::initial_shift = "exact" ADMBase::initial_data = "exact" ADMBase::lapse_evolution_method = "exact" ADMBase::shift_evolution_method = "exact" ADMBase::metric_type = "physical" Exact::exact_model = "de Sitter+Lambda" Exact::de_Sitter_Lambda__scale = 0.50 ######################################## ActiveThorns = "IOUtil" IOUtil::parfile_write = "no" ######################################## ActiveThorns = "AEILocalInterp PUGHInterp PUGHReduce SphericalSurface AHFinderDirect" AHFinderDirect::print_timing_stats = "true" AHFinderDirect::output_Theta_every = 1 AHFinderDirect::h_base_file_name = "dS.h" AHFinderDirect::Theta_base_file_name = "dS.Theta" AHFinderDirect::N_horizons = 1 AHFinderDirect::origin_x[1] = 0.0 AHFinderDirect::origin_y[1] = 0.0 AHFinderDirect::origin_z[1] = 0.0 AHFinderDirect::initial_guess_method[1] = "coordinate sphere" AHFinderDirect::initial_guess__coord_sphere__x_center[1] = 0.0 AHFinderDirect::initial_guess__coord_sphere__y_center[1] = 0.0 AHFinderDirect::initial_guess__coord_sphere__z_center[1] = 0.0 AHFinderDirect::initial_guess__coord_sphere__radius[1] = 1.0 ------------------------ -- *Luc Combi* CONICET PhD fellow (Buenos Aires, Argentina) Instituto Argentino de Radioastronomía

2 1

Meeting minutes for 2022-07-21
by Leonardo Rosa Werneck 21 Jul '22

21 Jul '22

Dear all, Here are the meetings for today's meeting. Participants: Peter, Steve, Keith, Beyhan, Leo, Roland, Alexandru, Fede, Liwei, Cheng, Johnny • Reviewers • Canuda: Miguel Zilhão will be contacted by Roland • FLRW solver: Helvi might be able to do it, but we might want to find a backup reviewer. No volunteers in the call. • Self-force 1D: already have reviewers. • NRPyElliptic: Leo will talk to someone in Helvi’s group about it ASAP. • Elliptica initial data reader: Beyhan volunteered. Roland will put Beyhan in touch with Alireza, the code’s “champion.” • Gallery examples • BBH: Steve ran the simulation half-way through completion, but a crash (unrelated with physical issues) prevented the run from finishing. Steve is confident enough that the results are still the same as in previous runs. Will update the website soon. • TOV: Peter already finished running and will update the website soon. • BNS: Atul and Yosef are not present, so no comment. • Open tickets • #2621: The star is likely very under-resolved, which is why the NaNs are appearing at merger. Roland will provide suggestions on how to modify the parfile to improve this behavior. • Other discussions • Peter: Github bots are sending emails to mailing lists, which are then getting stuck for approval by the maintainers. Are these intended? Roland: Yes, had to create a few robot accounts. • Peter: Brief discussion about a presentation poster, which has been cleared out by Roland and Steve. • Next week’s meeting • Chair: Leo • Minute taker: Johnny Cheers, Leo ------ Leonardo R. Werneck, Ph.D. Postdoctoral researcher Office EP 314 | Department of Physics | University of Idaho 875 Perimeter Dr. MS 0903 Moscow, ID 83844-0903, USA leonardo(a)uidaho.edu <mailto:leonardo@uidaho.edu> https://leowerneck.github.io <https://leowerneck.github.io/>

1 0

Einstein Toolkit Meeting Reminder
by reminders＠einsteintoolkit.org 20 Jul '22

20 Jul '22

1 0

Agenda for Thursday's Meeting
by reminders＠einsteintoolkit.org 18 Jul '22

18 Jul '22

Please update the Wiki with agenda items for Thursday's meeting. Thanks! https://docs.einsteintoolkit.org/et-docs/meeting_agenda --The Maintainers

1 0

Meeting minutes for 2022-07-14
by Bing-Jyun Tsao 14 Jul '22

14 Jul '22

Present: Leo, Roland, Yosef, Zach, Johnny, Hrishikesh, Beyhan Chair: Leo Minutes: Johnny 1. GitHub actions based CI framework <https://rhaas80.github.io/einsteintoolkit/> by Hrishikesh Kalyanaraman - add feature to plotting: hovering over shows the data - goal is to run ETK tests on GitHub actions 2. modules suggested for inclusion for next release <https://docs.einsteintoolkit.org/et-docs/Release_Details> [RH] - Canuda scalar and complex fields ticket <https://bitbucket.org/einsteintoolkit/tickets/issues/2485/include-complex-a…>-- proposed by Miguel Zilhao - Reviewed by Helvi Witek but might need help due to time constraint - FLRW (cosmological initial data) solver ticket <https://bitbucket.org/einsteintoolkit/tickets/issues/2549/inlcude-flrwsolve…> -- proposed by Hayley Macpherson - SelfForce-1D ticket <https://bitbucket.org/einsteintoolkit/tickets/issues/2602/include-regge-whe…> -- proposed by Samuel Cupp - Reviewed by Peter - NRPyEllipticET (elliptical solver for initial data) ticket <https://bitbucket.org/einsteintoolkit/tickets/issues/2616/add-nrpyelliptice…> -- proposed by Leo Werneck, Thiago Assumpcao, Zach Etienne - ID readers for Elliptica ticket <https://bitbucket.org/einsteintoolkit/tickets/issues/2619/include-ellipitca…> -- proposed by Alireza Rashti - issue: closed source, need to include library of initial data [ZE] help with writing the thorn and include data [RH] 3. timeline <https://docs.einsteintoolkit.org/et-docs/Release_Details> for upcoming release [RH] - Next week: choose features to include in ETK, finish review by 11/15, release is 11/10, name <https://docs.einsteintoolkit.org/et-docs/Release_Details#Release_names> will be chosen 10/06 (need suggestions) [RH] 4. missing updates for gallery examples [RH] - BNS <https://www.einsteintoolkit.org/gallery/bns/index.html> [YZ]: compiled with gcc 5. Tickets: - Ticket #2601 <https://bitbucket.org/einsteintoolkit/tickets/issues/2601/add-four-more-par…>: Add more parameters for carpet: Liwei provided a small parameter file as an example - Ticket #2620 <https://bitbucket.org/einsteintoolkit/tickets/issues/2620/wvuthorns-baikal-…>: issue with WVUThorns/Baikal: doesn’t compile on a cluster. Easy fix (missing include of CoordGauge) for the issue has been pushed to Master branch [ZE], if the fix will not affect other working thorns it should be fixed [YZ], suggestion: remove CoordGauge [YZ] - Ticket #181 <https://bitbucket.org/einsteintoolkit/tickets/issues/181/aeilocalinterp-sho…>: insufficient ghost points caused off-centering [RH] 6. Other: - ET working workshop/hackathon begins next week - ET summer school 2022 recording is available <https://www.youtube.com/playlist?list=PLRxi-yB7cTGfaVnPTkqg_Fw6Axy7nOBzK> [LW] Next call: Chair: Beyhan Minutes: Leo

2 1

Reduction operations with weights; 3 refinement centres
by Konrad Topolski 13 Jul '22

13 Jul '22

Dear ETK community, I'm writing with questions about two unrelated topics. 1. I would like to ask if there is a simple way of obtaining an L1 or L2 norm that comes from the integral with the sqrt(det(g)) measure, or a self-defined and stored weighing function? I realize it could be done most simply by defining an auxiliary variable with storage, whose value is the relevant field times "sqrt(det(g))", but perhaps this is already incorporated in some Cactus routines? Otherwise, I suppose the most efficient alternative is to assign storage for (and perfom reductions on) a new, detg-weighted variable, at only certain iterations, in the relevant thorn's schedule.ccl and only set values for it at these iterations? 2. I would like to learn, if there's a general guide for choosing the radii of 3 refinement centres, meaning a third, central, stationary one, in addition to the two moving ones. I want to prevent the edges of the 3rd refinement centre's levels from cutting into the objects, or the levels being unoptimally placed. For example, if I would like to have max_ref_lvl = L for the two moving centres, and L-1 (or L-2) for the third one, should the finest level of the third refinement centre (call it level N) have the radius R=(initial_separation) + r where r is the radius chosen for the other two refinement centres at level N? Best regards Konrad

2 1

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Users July 2022