Exodus output skipped timestep when sync only

[BoZeng1997]

Hi,

I have an issue with exodus output when running a multiapp problem with output on sync time only. There are timesteps in my sync list that were skipped in exodus output. There is no error message.

To reproduce

I will provide the input file main.i and sub.i in the following comment. Simply run it with {moose-executable} -i main.i (either parallel or serial) then go check the mech.e. You will see step

0.0004
0.003 0.004 0.005 0.009
0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09
0.1 0.2 0.3
2.5 3.0 3.5 4.0 4.5

are missing in the file.

My moose version

Here is the info about my MOOSE. I built it on an Ubuntu machine with Conda environment.

Framework Information:
MOOSE Version:           git commit e7a20f90 on 2024-03-08
LibMesh Version:         
PETSc Version:           3.20.3
SLEPc Version:           3.20.1

[BoZeng1997]

file main.i

dt=2e-5
refine=2
filebase=mech
[Problem]
  type = FEProblem
  solve = False
  kernel_coverage_check = false
  material_coverage_check = false
[]

[MultiApps]
  [sub]
    type = TransientMultiApp
    input_files = sub.i
    cli_args = 'dt=${dt};refine=${refine}'
    execute_on = 'TIMESTEP_END'
  []
[]

[Mesh]
uniform_refine=${fparse max(0,refine)}
  [gen]
    type = GeneratedMeshGenerator
    dim = 2
    xmax = 25e-3
    ymin = 0
    ymax = 5e-3
    nx = 25
    ny = 5
  []
[]


[Variables]
  [disp_x]
  []
[]


[Executioner]
  type = Transient
  solve_type = NEWTON
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package'
  petsc_options_value = 'lu      superlu_dist    '

  start_time = 0
  end_time = 5.0
  # dt = ${dt}
    [TimeSteppers]
      [time_step]
        type = FunctionDT
        function = set_timestep
      []
    []
  fixed_point_max_its = 1
  accept_on_max_fixed_point_iteration = true
[]

[Functions]
  [set_timestep]
    type = PiecewiseLinear
    x = '0     0.00198  0.002   0.0089   0.0090  0.0095  0.01  0.098  0.1   0.99   1     5'
    y = '${dt} ${dt}    1e-4     1e-4     5e-4    5e-4    2e-3   2e-3  0.01  0.01  0.05   0.05'
  []
[]

[Outputs]
  [nemesis]
    type = Exodus
    additional_execute_on = 'FINAL'
    sync_times =  '0.0001 0.0002 0.0003 0.0004 0.0005 0.0006 0.0007 0.0008 0.0009 '
                  '0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009 '
                  '0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 '
                  '0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 '
                  '1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0'
    sync_only=true
  []
  file_base = '${filebase}'
  print_linear_residuals = false
[]

[BoZeng1997]

file sub.i

[Problem]
  type = FEProblem
  solve = False
  kernel_coverage_check = false
[]

[Mesh]
  uniform_refine=${fparse max(0,refine)}
  [gen]
    type = GeneratedMeshGenerator
    dim = 2
    xmax = 25e-3
    ymin = 0
    ymax = 5e-3
    nx = 25
    ny = 5
  []
[]

[Executioner]
  type = Transient
  solve_type = NEWTON
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package -snes_type'
  petsc_options_value = 'lu      superlu_dist     vinewtonrsls'

  start_time = 0
  end_time = 600
  # dt = ${dt}
    [TimeSteppers]
      [time_step]
        type = FunctionDT
        function = set_timestep
      []
    []
[]

[Functions]
  [set_timestep]
    type = PiecewiseLinear
    x = '0     0.00198  0.002   0.0089   0.0090  0.0095  0.01  0.098  0.1   0.99   1     5'
    y = '${dt} ${dt}    1e-4     1e-4     5e-4    5e-4    2e-3   2e-3  0.01  0.01  0.05   0.05'
  []
[]

[Outputs]
  print_linear_residuals = false
[]

[GiudGiud]

I actually get different time steps in my output. But not the full list.

I ll create an issue and we ll address it when we can

[GiudGiud]

in the mean time just output on all times?