[Users] Problem with CarpetRegrid2 for changing num_levels via trigger

[Shamim Haque]

Hi Erik,

I indeed failed to check that. The Carpet::max_refinement_levels is set to
7. I should change it to 10.

Thank you.

Regards
Shamim Haque
Senior Research Fellow (SRF)
Department of Physics
IISER Bhopal

ᐧ

On Wed, May 25, 2022 at 2:11 AM Erik Schnetter <schnetter at gmail.com> wrote:

> 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/
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/users/attachments/20220525/18c13912/attachment-0001.html