Différences

Ci-dessous, les différences entre deux révisions de la page.

--- sunfluidh:homedata_poissonsolver_namelist [2017/09/25 11:22] – [Full data set of the namelist] yann
+++ sunfluidh:homedata_poissonsolver_namelist [2021/04/12 16:30] (Version actuelle) – yann
@@ Ligne 1: / Ligne 1: @@
 ===== Namelist "HomeData_PoissonSolver" =====
 <WRAP info>
-This namelist contains all parameters on the iterative solvers and the nV-cycle multigrid method directly implemented in the code for solving the Poisson's equation ("homemade" development).
+This namelist contains all parameters on the iterative solvers and the nV-cycle multigrid method directly implemented in the code for solving the Poisson's equation (in-house development). The iterative solvers are based on the Successive Over Relaxation methods (SOR).
 </WRAP>
 <note important>
-   * These data are taken into account only when the variable "Poisson_NumericalMethod" in the Namelist "Numerical_Methods" is correctly set (see ...).
+   * These data are only associated with the SOR-type solvers directly implemented in the code ("homemade" development)
+   * These data are taken into account only when the variable "Poisson_NumericalMethod" in the Namelist "Numerical_Methods" is correctly set (see [[New_Numerical_Methods_Setup_Namelist |Numerical_Methods (new version)]]).
      * Poisson_NumericalMethod= Home-Multigrid-ConstantMatrixCoef
      * Poisson_NumericalMethod= Home-Multigrid-VariableMatrixCoef
    * __** Not for the release SUNFLUIDH_EDU**__
 </note>
@@ Ligne 16: / Ligne 16: @@
    &HomeData_PoissonSolver  SolverName="SOR-Redblack-MPI",
-                            Convergence_Criterion= 1.D-15,
                             Number_max_Grid= 5,
                             Number_max_Cycle= 10,
@@ Ligne 24: / Ligne 23: @@
                             Number_Iteration_CoarsestGrid= 10 ,
                             Relaxation_Coefficient= 1.7 ,
+                            Convergence_Criterion= 1.D-15,
+                            Prolongation_Type="Uniform",
+                            Coarsening_Id(1:5)= 7,7,7,7,7 ,
                             Off_Set_Poisson_Source_term= .false. /
@@ Ligne 30: / Ligne 32: @@
 <note important>
    * Following the problem treated, some numerical methods are better adapted than others, this point is precised for each available option.
-   * The partial diagonalization method (Poisson_NumericalMethod = "Home-PartialDiagonalization") used for solving the Poisson equation is only permitted for separable problems.
-   * The HYPRE library solvers for solving the Poisson's equation are only available if the HYPRE library has been installed.
 </note>
 ===== Definition of the data set =====
+==== SolverName ====
+   * Type: Character string
+   * Selection of the iterative solver coupled with the multigrid method.
+       * __Jacobi__           : Successive Over Relaxed solver based on the Jacobi's method
+       * __SOR__              : Successive Over Relaxed solver based on the Gauss-Seidel method
+       * __SOR-Redblack__     : Successive Over Relaxed solver based on a "Red-Black" algorithm (useful with OpenMP-parallel computations or MPI domain decompositions because it can ensure the replicability of the solution whatever the number of MPI Processus and OpenMP threads).
+       * __SOR-Redblack-MPI__ : Successive Over Relaxed solver based on a "Red-Black" algorithm optimized for the domain decomposition approach (MPI-parallel computations). ** DO NOT USE IT IN Hydrid parallelisation (MPI+ OpenMP)**.
+       * __SOR-Redblack-MPIOMP__ : Successive Over Relaxed solver based on a "Red-Black" algorithm optimized for the domain decomposition approach (MPI-parallel computations) and the multi-threading parallelism (OpenMP).
+       * __SOR-RRB__ : Successive Over Relaxed solver based on a "Red-Black" algorithm with optimization based on data reordering. This method can be used for MPI-parallel computations but it is not yet ready for OpenMP and hydrid parallel ones.
 ==== Number_max_Grid (multigrid method only) ====
-   * Type : Integer value
+   * Type : Integer value (maximum value 16)
    * Maximum number of grid levels used in the multigrid method.
 ==== Number_max_Cycle (multigrid method only) ====
@@ Ligne 61: / Ligne 70: @@
    * Type : real value
    * The nV-cycle procedure of the multigrid method is stopped if the L2-norm of the Poisson equation's residu is lower than the criterion chosen.
+==== Coarsening_Id(:) ====
+   * Type : array of integer (max. size = 16)
+   * Define how the grid is coarsened during the restriction stage for each grid level "ng". The array size must be equal to the number of grid levels. For the grid level "ng", the coarsening is defined in regard to the previous (finner) grid level "ng-1". The value in Coarsening_Id(ng) "defines the direction(s) in which the coarsening occurs. Values are defined as :
+         * along the I-direction : 1
+         * along the J-direction : 2
+         * along the K-direction : 4
+         * along I and J directions: 3
+         * along I and K directions: 5
+         * along J and K directions: 6
+         * along I,J and K directions: 7
+   * 7 is the default value considered by the code when Coarsening_Id is not explicitely set (coarsening in all directions).
+   * The value of Coarsening_Id(1) is not considered since ng= 1 is related to the main grid, the finnest one, and it does not result from a restriction stage.
+   * For example, in order to define a coarsening process over 5 grid levels by considering the I-direction on the 2nd grid level, the J and K directions on the 3rd level, and all directions for any other levels : Coarsening_Id(1:5)= 7,1,6,7,7 or Coarsening_Id(2:5)= 1,6,7,7, or Coarsening_Id(2:3)= 1,6 are equivalent.
+====  Prolongation_Type ====
+   * Type : character string
+   * "Uniform " : The prolongation operator is of zero order. By considering two successive grid levels, errors on the coarser grid are directly reported on the finner one.This is the default value when "Prolongation_Type" is not explicitly set.
+   * "Linear" : The prolongation operator is of first-order. By considering two successive grid levels, errors on the coarsest grid are reported on the finest one by means of a tri-linear interpolation (in 3D). __BEWARE__ when the numerical simulation contains immersed bodies which can induce small inaccuracies.
 ==== Off_Set_Poisson_Source_term ====
    * Type : Boolean value. __** Not for the release SUNFLUIDH_EDU**__.
    * The source term is shifted from its averaged value such as its integral value over the domain is zero. This could facilitate the convergence in some particuliar cases. this option must be used with caution. For expert users only.