Clean up stale IGR integration comments

Author: StevenXue98

src/simulation/m_rhs.fpp

@@ -520,7 +520,7 @@ contains
             $:END_GPU_PARALLEL_LOOP()
         end if
 
-        ! TODO: Native IGR reconstruction still reads conservative ghost cells directly.
+        ! IGR reconstruction still reads conservative ghost-cell data directly.
         if (igr) then
             call nvtxStartRange("RHS-COMMUNICATION")
             call s_populate_variables_buffers(bc_type, q_cons_vf, pb_in, mv_in, q_T_sf)
@@ -548,7 +548,7 @@ contains
         if (qbmm) call s_mom_inv(q_cons_qp%vf, q_prim_qp%vf, mom_sp, mom_3d, pb_in, rhs_pb, mv_in, rhs_mv, idwbuff(1), &
             & idwbuff(2), idwbuff(3))
 
-        ! TODO: IGR viscous terms are still evaluated by legacy IGR logic.
+        ! IGR fills its viscous source fluxes during custom reconstruction.
         if ((viscous .and. .not. igr)) then
             call nvtxStartRange("RHS-VISCOUS")
             call s_get_viscous(qL_rsx_vf, dqL_prim_dx_n, dqL_prim_dy_n, dqL_prim_dz_n, qL_prim, qR_rsx_vf, dqR_prim_dx_n, &

src/simulation/m_start_up.fpp

@@ -827,12 +827,10 @@ contains
             call s_initialize_acoustic_src()
         end if
 
-        ! TODO: first igr divergence point
         if (viscous .and. (.not. igr)) then
             call s_initialize_viscous_module()
         end if
 
-        ! TODO: second igr divergence point (inside module)
         call s_initialize_rhs_module()
 
         if (surface_tension) call s_initialize_surface_tension_module()
@@ -843,7 +841,6 @@ contains
         call s_initialize_derived_variables_module()
         call s_initialize_time_steppers_module()
 
-        ! TODO: may need to alter for igr integration
         call s_initialize_boundary_common_module()
 
         if (down_sample) then
@@ -906,15 +903,10 @@ contains
 
         ! Computation of parameters, allocation of memory, association of pointers, and/or execution of any other tasks that are
         ! needed to properly configure the modules. The preparations below DO DEPEND on the grid being complete.
-        ! TODO: main igr divergence point
         if (igr) then
             call s_initialize_igr_module()
-            ! initialize riemann solvers module for inviscid igr branch
             call s_initialize_riemann_solvers_module()
         end if
-        ! non-igr implementation uses WENO, MUSCL reconstruction modules and riemann solvers module
-        ! goal is to initialize these modules alongside igr and make igr a flux modifier
-        ! To be determined
         if (.not. igr) then
             if (recon_type == recon_type_weno) then
                 call s_initialize_weno_module()
@@ -1070,7 +1062,6 @@ contains
         if (ib) then
             $:GPU_UPDATE(device='[ib_markers%sf]')
         end if
-        ! TODO: may need to copy data for igr integration
         #:if not MFC_CASE_OPTIMIZATION
             $:GPU_UPDATE(device='[igr, nb, igr_order]')
         #:endif
@@ -1085,12 +1076,9 @@ contains
         if (hyperelasticity) call s_finalize_hyperelastic_module()
         call s_finalize_derived_variables_module()
         call s_finalize_data_output_module()
-        ! TODO: potential igr divergence point
         call s_finalize_rhs_module()
-        ! TODO: igr divergence point, the same way they diverged in s_initialize_modules
         if (igr) then
             call s_finalize_igr_module()
-            ! finalize riemann solvers module for inviscid igr branch
             call s_finalize_riemann_solvers_module()
         else
             call s_finalize_cbc_module()
@@ -1109,7 +1097,6 @@ contains
         call s_finalize_boundary_common_module()
         if (relax) call s_finalize_relaxation_solver_module()
         if (bubbles_lagrange) call s_finalize_lagrangian_solver()
-        ! TODO: igr divergence point, mirroring s_initialize_modules
         if (viscous .and. (.not. igr)) then
             call s_finalize_viscous_module()
         end if

src/simulation/m_time_steppers.fpp

@@ -462,7 +462,7 @@ contains
 
             if (s == 1) then
                 if (run_time_info) then
-                    ! TODO: route IGR runtime information through q_prim_vf once primitive workspace semantics are fully integrated.
+                    ! IGR runtime information still uses conservative variables until primitive workspaces are fully integrated
                     if (igr) then
                         call s_write_run_time_information(q_cons_ts(1)%vf, t_step)
                     end if
@@ -625,8 +625,7 @@ contains
         integer                :: j, k, l  !< Generic loop iterators
         integer                :: fl       !< Fluid loop iterator
 
-        ! TODO: use the normal primitive workspace for IGR dt calculation once
-        ! q_prim_vf is authoritative for IGR diagnostics and stability checks.
+        ! IGR dt calculation still uses its custom conservative-variable path.
         if (.not. igr) then
             call s_convert_conservative_to_primitive_variables(q_cons_ts(1)%vf, q_T_sf, q_prim_vf, idwint)
         end if