[Users] Problem with CarpetRegrid2 for changing num_levels via trigger

[Erik Schnetter]

Shamim

There is a parameter Carpet::max_refinement_levels. What is its value?
The error message reads as if you set it to 7, and then try to
activate 10 levels. This won't work. Carpet::max_refinement_levels is
an upper bound for the number of levels that can be active. Much of
Carpet's internal bookkeeping depends on this value.

-erik


On Tue, May 24, 2022 at 3:58 PM Shamim Haque <shamims at iiserb.ac.in> wrote:
>
> Hi all,
>
> I am trying to change the num_levels for CarpetRegrid2 using the trigger, inspired from the gallery example at ET website. The trigger seems to kick in at the right place, as the simulation stopped while regridding for new num_levels.
>
> The dedicated lines in parfile are as follows:
>
> CarpetRegrid2::freeze_unaligned_levels = "yes"
> CarpetRegrid2::movement_threshold_1    = 0.10
> CarpetRegrid2::movement_threshold_2    = 0.10
> CarpetRegrid2::num_centres             = 2
> CarpetRegrid2::num_levels_1            = 7 #10
> CarpetRegrid2::num_levels_2            = 7
> CarpetRegrid2::position_x_1            = -13.543564
> CarpetRegrid2::position_x_2            = 13.543564
>
> CarpetRegrid2::radius_1[1]             =  270.0
> CarpetRegrid2::radius_1[2]             =  160.0
> CarpetRegrid2::radius_1[3]             =   90.0
> CarpetRegrid2::radius_1[4]             =   40.0
> CarpetRegrid2::radius_1[5]             =   25.0
> CarpetRegrid2::radius_1[6]             =   9.0
> CarpetRegrid2::radius_1[7]             =   6.0
> CarpetRegrid2::radius_1[8]             =   3.0
> CarpetRegrid2::radius_1[9]             =   1.5
>
> CarpetRegrid2::radius_2[1]             =  270.0
> CarpetRegrid2::radius_2[2]             =  160.0
> CarpetRegrid2::radius_2[3]             =   90.0
> CarpetRegrid2::radius_2[4]             =   40.0
> CarpetRegrid2::radius_2[5]             =   25.0
> CarpetRegrid2::radius_2[6]             =   9.0
>
> CarpetRegrid2::regrid_every            = 32
> CarpetRegrid2::symmetry_rotating180    = "no"
> CarpetRegrid2::verbose                 = "yes"
>
> Trigger::Trigger_Number = 1
> Trigger::Trigger_Checked_Variable[0]="ADMBase::alp"
> Trigger::Trigger_Reduction       [0]="minimum"
> Trigger::Trigger_Relation        [0]="<"
> Trigger::Trigger_Checked_Value   [0]=0.1
> Trigger::Trigger_Reaction        [0]="steerscalar"
> Trigger::Trigger_Steered_Scalar      [0] = "CarpetRegrid2::num_levels[0]" # == num_levels_1
> Trigger::Trigger_Steered_Scalar_Value[0] = "10"
>
> The outfile reports the following error:
>
> INFO (CarpetRegrid2): Centre 1 is at position [0.000191809,0.000291781,-0.00018496] with 10 levels
> INFO (CarpetRegrid2): Centre 2 is at position [0.000191809,0.000291781,-0.00018496] with 7 levels
> INFO (CarpetRegrid2): Regridding level 0 map 0 at iteration 4129 time 806.25
> INFO (CarpetRegrid2): Regridding levels 1 and up
> [1mERROR from host astro-Z390-AORUS-ELITE process 0
>   in thorn CarpetRegrid2, file /home/astro/Documents/ET/Cactus/configs/sim/build/CarpetRegrid2/regrid.cc:89:
>   -> [0m Region 1 has 10 levels active, which is larger than the maximum number of refinement levels 7
>
> I could not figure out where it went wrong. The minimum lapse at iteration 4128 is 0.095, so the condition in trigger seems to be working fine. Some insight into it would be helpful. Thanks in advance.
>
> Regards
> Shamim Haque
> Senior Research Fellow (SRF)
> Department of Physics
> IISER Bhopal
> _______________________________________________
> Users mailing list
> Users at einsteintoolkit.org
> http://lists.einsteintoolkit.org/mailman/listinfo/users



-- 
Erik Schnetter <schnetter at gmail.com>
http://www.perimeterinstitute.ca/personal/eschnetter/