How to make execution sequence of a basic MultiApp become: 1. principal parent.i 2; postprocessors parent.i; 3. principal sub.i; 4. postprocessors sub.i #26661

Wang-Yihu · 2024-01-28T07:37:42Z

Wang-Yihu
Jan 28, 2024

It seems that in MOOSE basic MultiApp, all the postprocessors objects cannot be executed between parent.i and sub.i..., different from transfers objects.
I have written a parent.i MOOSE input file and a sub.i MOOSE input file. In parent.i input file, I set

[MultiApps]
  [sub]
    type = TransientMultiApp
    execute_on = TIMESTEP_END
    input_files = sub.i
  []
[]

However, in postprocessors execute_on settings, there are these options:
ADJOINT_TIMESTEP_BEGIN;ADJOINT_TIMESTEP_END;ADJOINT;CUSTOM;FINAL;FORWARD;HOMOGENEOUS_FORWARD;INITIAL;LINEAR;MULTIAPP_FIXED_POINT_BEGIN;MULTIAPP_FIXED_POINT_END;NONE;NONLINEAR;TIMESTEP_BEGIN;TIMESTEP_END
None of them can realize my desire execution sequence: 1. principal parent.i 2; postprocessors parent.i; 3. principal sub.i; 4. postprocessors sub.i
So how can I realize my desire execution sequence? Make postprocessors objects be executed between parent.i and sub.i like transfers objects...

Answered by Wang-Yihu

Feb 1, 2024

Guillaume
I try another time. The true answer is to add execute_on = TIMESTEP_END in the [Postprocessors] of parent.i and add execute_on = TIMESTEP_END in the [MultiApps] of parent.i. The execution sequence will be 1. principal parent.i; 2. postprocessors parent.i; 3. principal sub.i, though the return sequence in the shell is different: 1. principal parent.i; 2. principal sub.i; 3. postprocessors parent.i.

Your guess MULTIAPP_FIXED_POINT_BEGIN is wrong, it makes the execution sequence: 1. principal parent.i; 2. principal sub.i; 3. postprocessors parent.i. , not my desire.

Your suggestion is quite valuable. I find the execution sequence and return sequence in the shell might are different…

View full answer

GiudGiud · 2024-01-28T16:44:30Z

GiudGiud
Jan 28, 2024
Collaborator

Hello

Use Debug/show_execution_order and Executioner/verbose to see the actual execution order.

I think the order you want is the natural one but it is good to check

Guillaume

5 replies

Wang-Yihu Jan 29, 2024
Author

Hello Guillaume, I've tried MULTIAPP_FIXED_POINT_BEGIN for the execute_on of the postprocessor, but I still cannot realize my goal. The execution sequence is still: 1. principal parent.i ;2. principal sub.i; 3; postprocessors parent.i; 4. postprocessors sub.i. The postprocessors objects are still not be executed between parent.i and sub.i like transfers objects...
Here are my input files and terminal returns. Is this problem caused by different versions between my MOOSE code and yours? You can try my input files in your computer...
parent.i

[Mesh]
  [rectangle_mesh]
    type = GeneratedMeshGenerator
    dim = 2
    nx = 10
    ny = 10
    xmax = 10.0
    ymax = 1.0
  []
[]
[Variables]
  [./u]
    initial_condition = 0.01
  [../]
[]
[Kernels]
  [./diff]
    type = HeatConduction
    variable = u
  [../]
  [./td]
    type = TimeDerivative
    variable = u
  [../]
[]
[BCs]
  [./left]
    type = DirichletBC
    variable = u
    boundary = left
    value = 10
  [../]
  [./right]
    type = DirichletBC
    variable = u
    boundary = right
    value = 0
  [../]
[]
[Materials]
  [coef]
    type = HeatConductionMaterial
    thermal_conductivity = 1
  []
[]
[Executioner]
  type = Transient
  num_steps = 5
  dt = 5
  solve_type = 'PJFNK'
  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]
[Outputs]
  exodus = true
[]
[MultiApps]
  [sub]
    type = TransientMultiApp
    execute_on = timestep_end
    input_files = sub.i
  []
[]
[Postprocessors]
  [pps_MULTIAPP_FIXED_POINT_BEGIN]
    type = SideDiffusiveFluxIntegral
    variable = u
    boundary = right
    diffusivity = thermal_conductivity
    execute_on = MULTIAPP_FIXED_POINT_BEGIN
  []
[]

sub.i

[Mesh]
  [rectangle_mesh]
    type = GeneratedMeshGenerator
    dim = 2
    nx = 10
    ny = 10
    xmax = 10.0
    ymax = 1.0
  []
[]
[Variables]
  [./diffused]
    order = FIRST
    family = LAGRANGE
  [../]
[]
[Kernels]
  [./diff]
    type = Diffusion
    variable = diffused
  [../]
[]
[BCs]
  [leftbc]
    type = DirichletBC
    variable = diffused
    boundary = left
    value = 10
  []
  [rightbc]
    type = DirichletBC
    variable = diffused
    boundary = right
    value = 0
  []
[]
[Executioner]
  type = Transient
  num_steps = 5
  dt = 5
  nl_abs_tol = 1e-6
  solve_type = 'PJFNK'
  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]
[Outputs]
  exodus = true
[]

The terminal will return

Framework Information:
MOOSE Version:           git commit 96aeccd391 on 2023-10-11
LibMesh Version:         
PETSc Version:           3.19.3
SLEPc Version:           3.19.1
Current Time:            Mon Jan 29 23:20:31 2024
Executable Timestamp:    Sat Jan 20 13:06:24 2024

Parallelism:
  Num Processors:          1
  Num Threads:             1

Mesh: 
  Parallel Type:           replicated
  Mesh Dimension:          2
  Spatial Dimension:       2
  Nodes:                   121
  Elems:                   100
  Num Subdomains:          1

Nonlinear System:
  Num DOFs:                121
  Num Local DOFs:          121
  Variables:               "u" 
  Finite Element Types:    "LAGRANGE" 
  Approximation Orders:    "FIRST" 

Execution Information:
  Executioner:             Transient
  TimeStepper:             ConstantDT
  TimeIntegrator:          ImplicitEuler
  Solver Mode:             Preconditioned JFNK
  PETSc Preconditioner:    hypre boomeramg 

sub0: Parallelism:
sub0:   Num Processors:          1
sub0:   Num Threads:             1
sub0: 
sub0: Mesh: 
sub0:   Parallel Type:           replicated
sub0:   Mesh Dimension:          2
sub0:   Spatial Dimension:       2
sub0:   Nodes:                   121
sub0:   Elems:                   100
sub0:   Num Subdomains:          1
sub0: 
sub0: Nonlinear System:
sub0:   Num DOFs:                121
sub0:   Num Local DOFs:          121
sub0:   Variables:               "diffused" 
sub0:   Finite Element Types:    "LAGRANGE" 
sub0:   Approximation Orders:    "FIRST" 
sub0: 
sub0: Execution Information:
sub0:   Executioner:             Transient
sub0:   TimeStepper:             ConstantDT
sub0:   TimeIntegrator:          ImplicitEuler
sub0:   Solver Mode:             Preconditioned JFNK
sub0:   PETSc Preconditioner:    hypre boomeramg 
sub0: 
sub0: 
sub0: Time Step 0, time = 0

Time Step 0, time = 0

Postprocessor Values:
+----------------+--------------------------------+
| time           | pps_MULTIAPP_FIXED_POINT_BEGIN |
+----------------+--------------------------------+
|   0.000000e+00 |                   0.000000e+00 |
+----------------+--------------------------------+


Time Step 1, time = 5, dt = 5
 0 Nonlinear |R| = 2.976489e+00
      0 Linear |R| = 2.976489e+00
      1 Linear |R| = 8.678395e-01
      2 Linear |R| = 2.276220e-01
      3 Linear |R| = 1.041954e-01
      4 Linear |R| = 2.323557e-02
      5 Linear |R| = 8.932714e-03
      6 Linear |R| = 1.469299e-03
      7 Linear |R| = 5.112804e-04
      8 Linear |R| = 9.046868e-05
      9 Linear |R| = 2.524437e-05
 1 Nonlinear |R| = 2.515752e-05
      0 Linear |R| = 2.515752e-05
      1 Linear |R| = 1.134176e-05
      2 Linear |R| = 4.153676e-06
      3 Linear |R| = 1.059784e-06
      4 Linear |R| = 2.775753e-07
      5 Linear |R| = 1.250038e-07
      6 Linear |R| = 4.036809e-08
      7 Linear |R| = 9.373929e-09
      8 Linear |R| = 2.731261e-09
      9 Linear |R| = 9.684379e-10
     10 Linear |R| = 3.027080e-10
     11 Linear |R| = 1.172366e-10
 2 Nonlinear |R| = 1.192954e-10
 Solve Converged!
sub0: 
sub0: Time Step 1, time = 5, dt = 5
sub0:  0 Nonlinear |R| = 3.082207e+00
sub0:       0 Linear |R| = 3.082207e+00
sub0:       1 Linear |R| = 4.587452e-01
sub0:       2 Linear |R| = 1.336525e-01
sub0:       3 Linear |R| = 2.379502e-02
sub0:       4 Linear |R| = 1.470892e-02
sub0:       5 Linear |R| = 2.342500e-03
sub0:       6 Linear |R| = 7.673881e-04
sub0:       7 Linear |R| = 2.421086e-04
sub0:       8 Linear |R| = 6.234554e-05
sub0:       9 Linear |R| = 2.425057e-05
sub0:  1 Nonlinear |R| = 2.348556e-05
sub0:       0 Linear |R| = 2.348556e-05
sub0:       1 Linear |R| = 1.738484e-05
sub0:       2 Linear |R| = 4.770847e-06
sub0:       3 Linear |R| = 2.143374e-06
sub0:       4 Linear |R| = 4.556795e-07
sub0:       5 Linear |R| = 1.460167e-07
sub0:       6 Linear |R| = 6.627679e-08
sub0:       7 Linear |R| = 1.759526e-08
sub0:       8 Linear |R| = 4.461972e-09
sub0:       9 Linear |R| = 1.544348e-09
sub0:      10 Linear |R| = 7.239115e-10
sub0:      11 Linear |R| = 1.230347e-10
sub0:  2 Nonlinear |R| = 6.356498e-10
sub0:  Solve Converged!

Postprocessor Values:
+----------------+--------------------------------+
| time           | pps_MULTIAPP_FIXED_POINT_BEGIN |
+----------------+--------------------------------+
|   0.000000e+00 |                   0.000000e+00 |
|   5.000000e+00 |                   0.000000e+00 |
+----------------+--------------------------------+


Time Step 2, time = 10, dt = 5
 0 Nonlinear |R| = 5.255087e-01
      0 Linear |R| = 5.255087e-01
      1 Linear |R| = 2.134116e-01
      2 Linear |R| = 8.702496e-02
      3 Linear |R| = 2.748490e-02
      4 Linear |R| = 6.583767e-03
      5 Linear |R| = 2.346716e-03
      6 Linear |R| = 5.119259e-04
      7 Linear |R| = 9.869132e-05
      8 Linear |R| = 2.162195e-05
      9 Linear |R| = 9.498853e-06
     10 Linear |R| = 2.851269e-06
 1 Nonlinear |R| = 4.479992e-06
      0 Linear |R| = 4.479992e-06
      1 Linear |R| = 1.257410e-06
      2 Linear |R| = 4.677313e-07
      3 Linear |R| = 1.580319e-07
      4 Linear |R| = 4.340009e-08
      5 Linear |R| = 1.218637e-08
      6 Linear |R| = 5.880566e-09
      7 Linear |R| = 1.758683e-09
      8 Linear |R| = 4.543855e-10
      9 Linear |R| = 1.443835e-10
     10 Linear |R| = 4.588931e-11
     11 Linear |R| = 1.386184e-11
 2 Nonlinear |R| = 1.428141e-11
 Solve Converged!
sub0: 
sub0: Time Step 2, time = 10, dt = 5
sub0:  0 Nonlinear |R| = 6.356498e-10
sub0:  Solve Converged!

Postprocessor Values:
+----------------+--------------------------------+
| time           | pps_MULTIAPP_FIXED_POINT_BEGIN |
+----------------+--------------------------------+
|   0.000000e+00 |                   0.000000e+00 |
|   5.000000e+00 |                   0.000000e+00 |
|   1.000000e+01 |                   1.061243e-01 |
+----------------+--------------------------------+


Time Step 3, time = 15, dt = 5
 0 Nonlinear |R| = 2.185691e-01
      0 Linear |R| = 2.185691e-01
      1 Linear |R| = 1.485067e-01
      2 Linear |R| = 4.815993e-02
      3 Linear |R| = 1.735658e-02
      4 Linear |R| = 3.797050e-03
      5 Linear |R| = 1.118219e-03
      6 Linear |R| = 2.356297e-04
      7 Linear |R| = 5.823421e-05
      8 Linear |R| = 1.533279e-05
      9 Linear |R| = 5.411445e-06
     10 Linear |R| = 2.099659e-06
 1 Nonlinear |R| = 3.708459e-06
      0 Linear |R| = 3.708459e-06
      1 Linear |R| = 1.211572e-06
      2 Linear |R| = 3.690877e-07
      3 Linear |R| = 1.193817e-07
      4 Linear |R| = 2.856420e-08
      5 Linear |R| = 1.237382e-08
      6 Linear |R| = 4.792531e-09
      7 Linear |R| = 1.133186e-09
      8 Linear |R| = 4.216877e-10
      9 Linear |R| = 1.318355e-10
     10 Linear |R| = 3.946277e-11
     11 Linear |R| = 1.256721e-11
 2 Nonlinear |R| = 1.256378e-11
 Solve Converged!
sub0: 
sub0: Time Step 3, time = 15, dt = 5
sub0:  0 Nonlinear |R| = 6.356498e-10
sub0:  Solve Converged!

Postprocessor Values:
+----------------+--------------------------------+
| time           | pps_MULTIAPP_FIXED_POINT_BEGIN |
+----------------+--------------------------------+
|   0.000000e+00 |                   0.000000e+00 |
|   5.000000e+00 |                   0.000000e+00 |
|   1.000000e+01 |                   1.061243e-01 |
|   1.500000e+01 |                   2.842361e-01 |
+----------------+--------------------------------+


Time Step 4, time = 20, dt = 5
 0 Nonlinear |R| = 1.324911e-01
      0 Linear |R| = 1.324911e-01
      1 Linear |R| = 9.412384e-02
      2 Linear |R| = 2.641949e-02
      3 Linear |R| = 8.728063e-03
      4 Linear |R| = 1.611333e-03
      5 Linear |R| = 4.976305e-04
      6 Linear |R| = 1.063236e-04
      7 Linear |R| = 3.647481e-05
      8 Linear |R| = 9.732385e-06
      9 Linear |R| = 3.124005e-06
     10 Linear |R| = 1.222690e-06
 1 Nonlinear |R| = 2.502088e-06
      0 Linear |R| = 2.502088e-06
      1 Linear |R| = 7.185495e-07
      2 Linear |R| = 2.005761e-07
      3 Linear |R| = 6.191835e-08
      4 Linear |R| = 1.774005e-08
      5 Linear |R| = 7.944339e-09
      6 Linear |R| = 2.595579e-09
      7 Linear |R| = 7.611138e-10
      8 Linear |R| = 2.793350e-10
      9 Linear |R| = 7.237345e-11
     10 Linear |R| = 1.797982e-11
 2 Nonlinear |R| = 1.804603e-11
 Solve Converged!
sub0: 
sub0: Time Step 4, time = 20, dt = 5
sub0:  0 Nonlinear |R| = 6.356498e-10
sub0:  Solve Converged!

Postprocessor Values:
+----------------+--------------------------------+
| time           | pps_MULTIAPP_FIXED_POINT_BEGIN |
+----------------+--------------------------------+
|   0.000000e+00 |                   0.000000e+00 |
|   5.000000e+00 |                   0.000000e+00 |
|   1.000000e+01 |                   1.061243e-01 |
|   1.500000e+01 |                   2.842361e-01 |
|   2.000000e+01 |                   4.713031e-01 |
+----------------+--------------------------------+


Time Step 5, time = 25, dt = 5
 0 Nonlinear |R| = 8.652591e-02
      0 Linear |R| = 8.652591e-02
      1 Linear |R| = 6.177158e-02
      2 Linear |R| = 1.632813e-02
      3 Linear |R| = 5.117503e-03
      4 Linear |R| = 6.913475e-04
      5 Linear |R| = 2.601342e-04
      6 Linear |R| = 6.040914e-05
      7 Linear |R| = 2.230330e-05
      8 Linear |R| = 6.151167e-06
      9 Linear |R| = 1.981397e-06
     10 Linear |R| = 7.485812e-07
 1 Nonlinear |R| = 1.458158e-06
      0 Linear |R| = 1.458158e-06
      1 Linear |R| = 4.481580e-07
      2 Linear |R| = 1.200042e-07
      3 Linear |R| = 4.291735e-08
      4 Linear |R| = 1.301063e-08
      5 Linear |R| = 5.220129e-09
      6 Linear |R| = 1.690922e-09
      7 Linear |R| = 6.195129e-10
      8 Linear |R| = 2.129039e-10
      9 Linear |R| = 4.227435e-11
     10 Linear |R| = 1.178166e-11
 2 Nonlinear |R| = 1.182717e-11
 Solve Converged!
sub0: 
sub0: Time Step 5, time = 25, dt = 5
sub0:  0 Nonlinear |R| = 6.356498e-10
sub0:  Solve Converged!

Postprocessor Values:
+----------------+--------------------------------+
| time           | pps_MULTIAPP_FIXED_POINT_BEGIN |
+----------------+--------------------------------+
|   0.000000e+00 |                   0.000000e+00 |
|   5.000000e+00 |                   0.000000e+00 |
|   1.000000e+01 |                   1.061243e-01 |
|   1.500000e+01 |                   2.842361e-01 |
|   2.000000e+01 |                   4.713031e-01 |
|   2.500000e+01 |                   6.275704e-01 |
+----------------+--------------------------------+

GiudGiud Jan 29, 2024
Collaborator

Are you sure you are not figuring out when it is executed by the output of the table? because this is not right

the only thing that tells you when it is executed is the Debug/show_execution_order option.

The table is always output in this order (at the end of the parent app timestep)

Wang-Yihu Jan 29, 2024
Author

Oh! I do not know this before, I just think the execution order is the same as the solver output order!....

Wang-Yihu Feb 1, 2024
Author

Guillaume
I try another time. The true answer is to add execute_on = TIMESTEP_END in the [Postprocessors] of parent.i and add execute_on = TIMESTEP_END in the [MultiApps] of parent.i. The execution sequence will be 1. principal parent.i; 2. postprocessors parent.i; 3. principal sub.i, though the return sequence in the shell is different: 1. principal parent.i; 2. principal sub.i; 3. postprocessors parent.i.

Your guess MULTIAPP_FIXED_POINT_BEGIN is wrong, it makes the execution sequence: 1. principal parent.i; 2. principal sub.i; 3. postprocessors parent.i. , not my desire.

Your suggestion is quite valuable. I find the execution sequence and return sequence in the shell might are different for the first time. And I also find we can input

[Debug]
  show_execution_order = ALWAYS
[]

to track the execution sequence.

I also insert several printf in the MOOSE Postprocessors class code(my MOOSEApp name is lianxi_mesh, and rename my printf-inserted-Receiver.C, Receiver.h into ReceiverReport.C, ReceiverReport.h and rename my printf-inserted-SideDiffusiveFluxIntegral.C, SideDiffusiveFluxIntegral.h into SideDiffusiveFluxIntegralReport.C, SideDiffusiveFluxIntegralReport.h) to track the execution sequence. My code is here to give examples for learners.

ReceiverReport.C.txt
ReceiverReport.h.txt
SideDiffusiveFluxIntegralReport.C.txt
SideDiffusiveFluxIntegralReport.h.txt
sub.i.txt
parent.i.txt

Answer selected by Wang-Yihu

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

How to make execution sequence of a basic MultiApp become: 1. principal parent.i 2; postprocessors parent.i; 3. principal sub.i; 4. postprocessors sub.i #26661

{{title}}

Replies: 1 comment 5 replies

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

{{editor}}'s edit

{{editor}}'s edit

Select a reply

How to make execution sequence of a basic MultiApp become: 1. principal parent.i 2; postprocessors parent.i; 3. principal sub.i; 4. postprocessors sub.i #26661

Wang-Yihu Jan 28, 2024

Replies: 1 comment · 5 replies

GiudGiud Jan 28, 2024 Collaborator

Wang-Yihu Jan 29, 2024 Author

GiudGiud Jan 29, 2024 Collaborator

Wang-Yihu Jan 29, 2024 Author

Wang-Yihu Feb 1, 2024 Author

Wang-Yihu
Jan 28, 2024

Replies: 1 comment 5 replies

GiudGiud
Jan 28, 2024
Collaborator

Wang-Yihu Jan 29, 2024
Author

GiudGiud Jan 29, 2024
Collaborator

Wang-Yihu Jan 29, 2024
Author

Wang-Yihu Feb 1, 2024
Author