src/common/m_boundary_common.fpp

• Changed boundary condition array indexing from (-1:1) to (1:2) to
  support case optimization
  
• Updated all references from negative indices (-1, 1) to positive
  indices (1, 2) throughout the module
  
• Added conditional compilation directives #:if not
MFC_CASE_OPTIMIZATION or num_dims > 2 to skip 3D boundary operations
  in optimized 2D cases
  
• Changed data type from real(wp) to real(stp) for mixed precision
  support in buffer parameters
  
• Updated MPI I/O operations to handle mixed precision with MPI_BYTE and
  element count calculations

src/simulation/m_weno.fpp

• Added #:include 'case.fpp' directive for case optimization support
  
• Wrapped WENO order 5 and 7 reconstruction code with conditional
  compilation #:if not MFC_CASE_OPTIMIZATION or weno_num_stencils > N
  
• Nested TENO stencil calculations within additional optimization guards
  for higher-order schemes
  
• Maintained existing WENO algorithm logic while enabling compile-time
  optimization for reduced stencil cases

src/simulation/m_bubbles_EL.fpp

• Changed q_beta from type(vector_field) to type(scalar_field),
dimension(:), allocatable for better memory layout
  
• Updated allocation from q_beta%vf(1:q_beta_idx) to
  q_beta(1:q_beta_idx) with individual scalar field setup
  
• Changed data type from real(wp) to real(stp) for mixed precision
  support in gradient calculations
  
• Updated s_gradient_dir subroutine signature to accept raw arrays
  instead of scalar field types
  
• Updated all references from q_beta%vf(i)%sf to q_beta(i)%sf throughout
  the module

src/pre_process/m_boundary_conditions.fpp

• Updated boundary condition array indexing from (-1:1) to (1:2) in
  function signatures
  
• Added conditional logic to map old negative indices (-1) to new
  positive index (1) and (1) to (2)
  
• Updated all bc_type array accesses to use new indexing scheme with
  explicit if-statements for location mapping

src/simulation/m_viscous.fpp

• Added case.fpp include directive for case optimization support
  
• Wrapped viscous stress tensor computation blocks with
  MFC_CASE_OPTIMIZATION conditional to skip 3D-specific calculations in
  2D cases
  
• Reorganized gradient computation loops with conditional compilation
  for dimensional optimization
  
• Improved code indentation and structure for better readability

src/simulation/m_riemann_solvers.fpp

• Added case.fpp include for case optimization support
  
• Wrapped MHD-related 3D vector operations with MFC_CASE_OPTIMIZATION
  conditionals to optimize 2D cases
  
• Expanded private variable lists in GPU parallel loops for better
  memory management
  
• Fixed duplicate variable assignments and improved variable
  initialization in HLLC solver
  
• Added explicit type conversions for mixed precision support

src/common/m_mpi_common.fpp

• Added case.fpp include directive for case optimization
  
• Changed halo_size from integer to integer(kind=8) for 64-bit support
  
• Fixed dimension indexing to use num_dims instead of hardcoded value 3
  for 2D/3D compatibility
  
• Added explicit type conversions between wp and stp precision kinds in
  MPI buffer operations

src/simulation/m_cbc.fpp

• Added case.fpp include for case optimization support
  
• Moved dpres_ds variable declaration from module level to local scope
  within subroutine
  
• Added dpres_ds to private variable list in GPU parallel loop
  
• Wrapped dpi_inf_dt assignment with MFC_CASE_OPTIMIZATION conditional
  for single-fluid optimization
  
• Improved GPU update directive formatting

src/pre_process/m_global_parameters.fpp

• Added new simplex_perturb logical flag for simplex noise perturbation
  support
  
• Added simplex_params derived type variable to store simplex noise
  parameters
  
• Initialized simplex perturbation parameters in module initialization
  routine

src/simulation/m_qbmm.fpp

• Changed pb parameter type from real(wp) to real(stp) for mixed
  precision support
  
• Wrapped bubble model coefficient calculations with preprocessor
  conditionals for case optimization
  
• Added missing loop variables i1, i2, j to GPU parallel loop private
  list
  
• Added TODO comment regarding rhs_pb and rhs_mv precision types

src/pre_process/m_simplex_noise.fpp

• New module implementing 3D and 2D Perlin simplex noise functions
  
• Includes gradient lookup tables and permutation vectors for noise
  generation
  
• Provides f_simplex3d and f_simplex2d functions for procedural noise

src/simulation/m_fftw.fpp

• Refactored GPU FFT filter implementation to use direct device
  addresses instead of pointers
  
• Removed #:block UNDEF_CCE wrapper and simplified GPU data management
  
• Restructured loop logic to process initial ring separately before main
  loop
  
• Improved compatibility with both NVIDIA and AMD GPU compilers

src/simulation/m_start_up.fpp

• Changed WP_MOK from 8 bytes to 4 bytes for mixed precision I/O
  
• Updated MPI file read calls to use mpi_io_type multiplier for data
  size
  
• Added logic to copy q_cons_ts(1) to q_cons_ts(2) for multi-stage time
  steppers
  
• Fixed NaN checking to use explicit real() conversion for mixed
  precision
  
• Added GPU updates for bubble and chemistry parameters
  
• Removed unused GPU declare directive for q_cons_temp

src/simulation/m_rhs.fpp

• Added conditional GPU declare directives for OpenACC-specific
  declarations
  
• Changed loop variable types to integer(kind=8) for large array
  indexing
  
• Wrapped flux allocation logic with if (.not. igr) condition
  
• Updated bc_type parameter dimension from 1:3, -1:1 to 1:3, 1:2
  
• Changed pb_in, mv_in parameters to real(stp) type for mixed precision
  
• Added TODO comments about precision type consistency

src/pre_process/m_icpp_patches.fpp

• Added conditional compilation for patch_id_fp array type based on
  MFC_MIXED_PRECISION flag
  
• Uses integer(kind=1) for mixed precision mode, standard integer
  otherwise
  
• Applied consistently across all patch subroutines for memory
  efficiency

src/simulation/m_ibm.fpp

• Changed pb_in and mv_in parameters to real(stp) type for mixed
  precision
  
• Added missing variables to GPU parallel loop private list
  
• Added type casting for pb_in and mv_in in interpolation calls
  
• Changed levelset_in access to use explicit real() conversion
  
• Updated pb_in, mv_in intent from INOUT to IN in interpolation
  subroutine

src/simulation/include/inline_capillary.fpp

• Wrapped 3D capillary tensor calculations with preprocessor conditional
  
• Only computes 3D components when not using case optimization or when
  num_dims > 2
  
• Reduces unnecessary computations for 2D simulations

src/simulation/m_muscl.fpp

• Renamed MUSCL workspace arrays from v_rs_ws_x/y/z to
  v_rs_ws_x/y/z_muscl to avoid naming conflicts
  
• Fixed GPU parallel loop syntax and indentation consistency
  
• Corrected #:endcall directives to include full macro name for clarity
  
• Improved code formatting and nested subroutine indentation

src/common/m_chemistry.fpp

• Removed sequential GPU loop directives that were causing compilation
  issues
  
• Added temporary variable T_in for type conversion in temperature
  calculations
  
• Wrapped large GPU parallel loop with #:block UNDEF_AMD to handle AMD
  compiler differences
  
• Added missing variables to GPU parallel loop private and copyin lists