[Users] Refinement levels 3 to 5.

[Bilal Hameed]

Awesome, thanks!

On Wed, Aug 26, 2020, 5:37 PM Roland Haas <rhaas at illinois.edu> wrote:

> Hello Bilal,
>
> > When I change, Carpet::max_refinement_levels = 5 (instead of 3). and
> > I also changed,
> >
> > CarpetRegrid2::num_centres  = 1
> > #CarpetRegrid2::num_levels_1 = 3
> > #CarpetRegrid2::radius_1[1]  = 1
> > #CarpetRegrid2::radius_1[2]  = 1
> >
> > CarpetRegrid2::num_levels_1 = 5
> > CarpetRegrid2::radius_1[1]  = 5
> > CarpetRegrid2::radius_1[2]  = 5
> > CarpetRegrid2::radius_1[3]  = 5
> > CarpetRegrid2::radius_1[4]  = 5
> You are correct in changing both Carpet::max_refinement_levels which
> controls the maximum possible number of refinement levels as well as
> CarpetRegrid2::num_levels_1 which controls that actual number of
> refinement levels used.
>
> > My simulation runs successfully. but fourth and fifth levels shown
> > only nan as follows,
> >
> > INFO (CT_MultiLevel): === Level 2 ===
> > INFO (CT_MultiLevel):  * Equation #0: 10 iterations, norm of final
> residual = 2.5741620337e-05
> > INFO (CT_MultiLevel):  * Equation #1: 10 iterations, norm of final
> residual = 8.7340004508e-06
> > INFO (CT_MultiLevel):  * Equation #2: 10 iterations, norm of final
> residual = 8.7156879010e-06
> > INFO (CT_MultiLevel):  * Equation #3: 10 iterations, norm of final
> residual = 8.6982806224e-06
> > INFO (CT_MultiLevel): === Level 3 ===
> > INFO (CT_MultiLevel):  * Equation #0: 1 iterations, norm of final
> residual = nan
> > INFO (CT_MultiLevel):  * Equation #1: 1 iterations, norm of final
> residual = nan
> > INFO (CT_MultiLevel):  * Equation #2: 1 iterations, norm of final
> residual = nan
> > INFO (CT_MultiLevel):  * Equation #3: 1 iterations, norm of final
> residual = nan
>
> > What else should I fix?
>
> I am not really familiar (at all, never having used it) with
> CT_MultiLevel itself. However looking at the parfile
> repos/ctthorns/CT_MultiLevel/par/constraints_spherical.par you can see
> that it sets:
>
> CT_MultiLevel::topMGlevel          = 2
>
> and that parameter is described in the thorn's documentation (in the
> doc directory or online)
>
> http://einsteintoolkit.org/thornguide/CTThorns/CT_MultiLevel/documentation.html
> as:
>
> --8<--
> The only other parameter which must be set to ensure correct operation
> is CT_MultiLevel::topMGlevel, which tells the solver which is the finest
> refinement level that covers the entire domain in which the equation is
> to be solved. All levels below and including this will be used in the
> classical multigrid sense (e.g., in a V- or FMG-cycle). The ones above
> will be treated as local AMR boxes, and be solved via a multilevel
> prescription; presently, this involves simply interpolating the
> solution from the topMGlevel refinement level at the end of the
> multigrid part and relaxing it progressively, from coarser to finer, on
> all the local grids.
> --8<--
>
> so since you changed max_refinement_levels from 3 to 5 I would also
> change topMGlevel from 2 to 4.
>
> Yours,
> Roland
>
> --
> My email is as private as my paper mail. I therefore support encrypting
> and signing email messages. Get my PGP key from http://pgp.mit.edu .
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/users/attachments/20200826/de8c5c45/attachment.html