m_hyperelastic.fpp.f90

Go to the documentation of this file.

# 1 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!>
!! @file
!! @brief Contains module m_hyperelastic
 
# 1 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp" 1
# 1 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 1
# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 46 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 58 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 68 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 98 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 110 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 120 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
# 1 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp" 1
# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 46 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 58 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 68 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 98 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 110 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 120 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp" 2
 
# 4 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 7 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 8 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 20 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 43 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 48 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 53 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 58 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 63 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 68 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 76 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 81 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 86 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 91 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 96 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 101 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 106 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 111 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 116 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 121 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 151 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 192 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 206 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 231 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 242 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 244 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 255 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 284 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 294 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 304 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 313 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 330 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 340 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 347 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 353 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 359 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 365 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 371 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 377 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
! New line at end of file is required for FYPP
# 3 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
# 1 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp" 1
# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 46 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 58 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 68 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 98 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 110 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 120 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp" 2
 
# 7 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 17 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 22 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 27 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 32 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 37 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 42 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 47 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 52 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 57 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 62 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 73 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 78 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 83 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 88 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 103 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 131 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 160 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 175 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 193 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 215 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 244 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 259 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 269 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 278 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 294 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 304 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 311 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
! New line at end of file is required for FYPP
# 4 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
 
! GPU parallel region (scalar reductions, maxval/minval)
# 23 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! GPU parallel loop over threads (most common GPU macro)
# 43 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Required closing for GPU_PARALLEL_LOOP
# 55 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Mark routine for device compilation
# 112 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Declare device-resident data
# 130 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Inner loop within a GPU parallel region
# 145 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Scoped GPU data region
# 164 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Host code with device pointers (for MPI with GPU buffers)
# 193 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Allocate device memory (unscoped)
# 207 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Free device memory
# 219 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Atomic operation on device
# 231 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! End atomic capture block
# 242 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Copy data between host and device
# 254 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Synchronization barrier
# 266 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Import GPU library module (openacc or omp_lib)
# 275 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code only for AMD compiler
# 282 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code for non-Cray compilers
# 289 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code only for Cray compiler
# 296 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code for non-NVIDIA compilers
# 303 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 305 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
# 306 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp" 2
 
# 14 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Caution: This macro requires the use of a binding script to set CUDA_VISIBLE_DEVICES, such that we have one GPU device per MPI
! rank. That's because for both cudaMemAdvise (preferred location) and cudaMemPrefetchAsync we use location = device_id = 0. For an
! example see misc/nvidia_uvm/bind.sh. NVIDIA unified memory page placement hint
# 57 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Allocate and create GPU device memory
# 77 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Free GPU device memory and deallocate
# 85 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for vector fields
# 109 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for scalar fields
# 125 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for acoustic source spatials
# 150 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 156 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 163 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
! New line at end of file is required for FYPP
# 6 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp" 2
 
!> @brief Computes the left Cauchy--Green deformation tensor and hyperelastic stress source terms
 
module m_hyperelastic
 
    use m_derived_types
    use m_global_parameters
    use m_variables_conversion
    use m_finite_differences
 
    implicit none
 
    private; public :: s_hyperelastic_rmt_stress_update, s_initialize_hyperelastic_module, s_finalize_hyperelastic_module
 
    ! The btensor at the cell-interior Gaussian quadrature points. These tensor is needed to be calculated once and make the code
    ! DRY.
    type(vector_field) :: btensor
 
# 23 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 23 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc declare create(btensor) 
# 23 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 23 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp declare target (btensor) 
# 23 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
 
    real(wp), allocatable, dimension(:,:) :: fd_coeff_x_hyper
    real(wp), allocatable, dimension(:,:) :: fd_coeff_y_hyper
    real(wp), allocatable, dimension(:,:) :: fd_coeff_z_hyper
 
# 28 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 28 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc declare create(fd_coeff_x_hyper, fd_coeff_y_hyper, fd_coeff_z_hyper) 
# 28 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 28 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp declare target (fd_coeff_x_hyper, fd_coeff_y_hyper, fd_coeff_z_hyper) 
# 28 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
    real(wp), allocatable, dimension(:) :: gs_hyper
 
# 30 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 30 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc declare create(Gs_hyper) 
# 30 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 30 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp declare target (Gs_hyper) 
# 30 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
 
contains
 
    !> Initialize the hyperelastic module
    impure subroutine s_initialize_hyperelastic_module
 
        integer :: i  !< generic iterator
 
#ifdef MFC_DEBUG
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:39: ', '@:ALLOCATE(btensor%vf(1:b_size))'
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    allocate (btensor%vf(1:b_size))
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc enter data create(btensor%vf) 
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target enter data map(always,alloc:btensor%vf) 
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
        do i = 1, b_size
#ifdef MFC_DEBUG
# 41 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 41 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 41 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 41 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:41: ', '@:ALLOCATE(btensor%vf(i)%sf(0:m, 0:n, 0:p))'
# 41 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 41 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 41 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 41 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 41 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    allocate (btensor%vf(i)%sf(0:m, 0:n, 0:p))
# 41 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 41 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 41 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 41 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc enter data create(btensor%vf(i)%sf) 
# 41 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 41 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target enter data map(always,alloc:btensor%vf(i)%sf) 
# 41 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
        end do
#ifdef _CRAYFTN
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        integer :: macros_setup_vfs_i
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#ifdef MFC_DEBUG
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:43: ', '@:ACC_SETUP_VFs(btensor)'
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc enter data copyin(btensor) 
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target enter data map(to:btensor) 
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc enter data copyin(btensor%vf) 
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target enter data map(to:btensor%vf) 
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
            if (allocated(btensor%vf)) then
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                do macros_setup_vfs_i = lbound(btensor%vf, 1), ubound(btensor%vf, 1)
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                    if (associated(btensor%vf(macros_setup_vfs_i)%sf)) then
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc enter data copyin(btensor%vf(macros_setup_vfs_i)) 
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target enter data map(to:btensor%vf(macros_setup_vfs_i)) 
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc enter data copyin(btensor%vf(macros_setup_vfs_i)%sf) 
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target enter data map(to:btensor%vf(macros_setup_vfs_i)%sf) 
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                    end if
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                end do
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
            end if
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
 
#ifdef MFC_DEBUG
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:45: ', '@:ALLOCATE(Gs_hyper(1:num_fluids))'
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    allocate (gs_hyper(1:num_fluids))
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc enter data create(Gs_hyper) 
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target enter data map(always,alloc:Gs_hyper) 
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
 
# 46 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 46 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc loop seq 
# 46 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 46 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 46 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
        do i = 1, num_fluids
            gs_hyper(i) = fluid_pp(i)%G
        end do
 
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc update device(Gs_hyper) 
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target update to(Gs_hyper) 
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
 
#ifdef MFC_DEBUG
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:52: ', '@:ALLOCATE(fd_coeff_x_hyper(-fd_number:fd_number, 0:m))'
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    allocate (fd_coeff_x_hyper(-fd_number:fd_number, 0:m))
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc enter data create(fd_coeff_x_hyper) 
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target enter data map(always,alloc:fd_coeff_x_hyper) 
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
        if (n > 0) then
#ifdef MFC_DEBUG
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:54: ', '@:ALLOCATE(fd_coeff_y_hyper(-fd_number:fd_number, 0:n))'
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    allocate (fd_coeff_y_hyper(-fd_number:fd_number, 0:n))
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc enter data create(fd_coeff_y_hyper) 
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target enter data map(always,alloc:fd_coeff_y_hyper) 
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
        end if
        if (p > 0) then
#ifdef MFC_DEBUG
# 57 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 57 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 57 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 57 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:57: ', '@:ALLOCATE(fd_coeff_z_hyper(-fd_number:fd_number, 0:p))'
# 57 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 57 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 57 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 57 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 57 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    allocate (fd_coeff_z_hyper(-fd_number:fd_number, 0:p))
# 57 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 57 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 57 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 57 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc enter data create(fd_coeff_z_hyper) 
# 57 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 57 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target enter data map(always,alloc:fd_coeff_z_hyper) 
# 57 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
        end if
 
        ! Computing centered finite difference coefficients
        call s_compute_finite_difference_coefficients(m, x_cc, fd_coeff_x_hyper, buff_size, fd_number, fd_order)
 
# 62 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 62 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc update device(fd_coeff_x_hyper) 
# 62 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 62 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target update to(fd_coeff_x_hyper) 
# 62 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
        if (n > 0) then
            call s_compute_finite_difference_coefficients(n, y_cc, fd_coeff_y_hyper, buff_size, fd_number, fd_order)
 
# 65 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 65 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc update device(fd_coeff_y_hyper) 
# 65 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 65 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target update to(fd_coeff_y_hyper) 
# 65 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
        end if
        if (p > 0) then
            call s_compute_finite_difference_coefficients(p, z_cc, fd_coeff_z_hyper, buff_size, fd_number, fd_order)
 
# 69 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 69 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc update device(fd_coeff_z_hyper) 
# 69 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 69 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target update to(fd_coeff_z_hyper) 
# 69 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
        end if
 
    end subroutine s_initialize_hyperelastic_module
 
    !> Compute the left Cauchy-Green deformation tensor and update the hyperelastic stress
    subroutine s_hyperelastic_rmt_stress_update(q_cons_vf, q_prim_vf)
 
        type(scalar_field), dimension(sys_size), intent(inout) :: q_cons_vf
        type(scalar_field), dimension(sys_size), intent(inout) :: q_prim_vf
 
# 83 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
            real(wp), dimension(tensor_size) :: tensora, tensorb
# 85 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 89 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
            real(wp), dimension(num_fluids) :: alpha_k, alpha_rho_k
# 91 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        real(wp), dimension(2) :: re
        real(wp)               :: rho, gamma, pi_inf, qv
        real(wp)               :: g_local
        integer                :: j, k, l, i, r
 
 
# 96 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 96 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 96 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc parallel loop collapse(3) gang vector default(present) private(i, j, k, l, alpha_K, alpha_rho_K, rho, gamma, pi_inf, qv, G_local, Re, tensora, tensorb) 
# 96 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 96 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 96 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 96 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 96 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(3) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(i, j, k, l, alpha_K, alpha_rho_K, rho, gamma, pi_inf, qv, G_local, Re, tensora, tensorb) 
# 96 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 98 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        do l = 0, p
            do k = 0, n
                do j = 0, m
                    call s_compute_species_fraction(q_cons_vf, j, k, l, alpha_rho_k, alpha_k)
 
                    ! If in simulation, use acc mixture subroutines
                    call s_convert_species_to_mixture_variables_acc(rho, gamma, pi_inf, qv, alpha_k, alpha_rho_k, re, g_local, &
                        & gs_hyper)
                    rho = max(rho, sgm_eps)
                    g_local = max(g_local, sgm_eps)
 
                    if (g_local > verysmall) then
 
# 110 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 110 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc loop seq 
# 110 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 110 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 110 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
                        do i = 1, tensor_size
                            tensora(i) = 0._wp
                        end do
                        ! STEP 1: computing the grad_xi tensor using finite differences grad_xi definition / organization number for
                        ! the tensor 1-3: dxix_dx, dxiy_dx, dxiz_dx 4-6 : dxix_dy, dxiy_dy, dxiz_dy 7-9 : dxix_dz, dxiy_dz, dxiz_dz
 
# 116 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 116 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc loop seq 
# 116 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 116 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 116 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
                        do r = -fd_number, fd_number
                            ! derivatives in the x-direction
                            tensora(1) = tensora(1) + q_prim_vf(xibeg)%sf(j + r, k, l)*fd_coeff_x_hyper(r, j)
                            tensora(2) = tensora(2) + q_prim_vf(xibeg + 1)%sf(j + r, k, l)*fd_coeff_x_hyper(r, j)
                            tensora(3) = tensora(3) + q_prim_vf(xiend)%sf(j + r, k, l)*fd_coeff_x_hyper(r, j)
                            ! derivatives in the y-direction
                            tensora(4) = tensora(4) + q_prim_vf(xibeg)%sf(j, k + r, l)*fd_coeff_y_hyper(r, k)
                            tensora(5) = tensora(5) + q_prim_vf(xibeg + 1)%sf(j, k + r, l)*fd_coeff_y_hyper(r, k)
                            tensora(6) = tensora(6) + q_prim_vf(xiend)%sf(j, k + r, l)*fd_coeff_y_hyper(r, k)
                            ! derivatives in the z-direction
                            tensora(7) = tensora(7) + q_prim_vf(xibeg)%sf(j, k, l + r)*fd_coeff_z_hyper(r, l)
                            tensora(8) = tensora(8) + q_prim_vf(xibeg + 1)%sf(j, k, l + r)*fd_coeff_z_hyper(r, l)
                            tensora(9) = tensora(9) + q_prim_vf(xiend)%sf(j, k, l + r)*fd_coeff_z_hyper(r, l)
                        end do
                        ! STEP 2a: computing the adjoint of the grad_xi tensor for the inverse
                        tensorb(1) = tensora(5)*tensora(9) - tensora(6)*tensora(8)
                        tensorb(2) = -(tensora(2)*tensora(9) - tensora(3)*tensora(8))
                        tensorb(3) = tensora(2)*tensora(6) - tensora(3)*tensora(5)
                        tensorb(4) = -(tensora(4)*tensora(9) - tensora(6)*tensora(7))
                        tensorb(5) = tensora(1)*tensora(9) - tensora(3)*tensora(7)
                        tensorb(6) = -(tensora(1)*tensora(6) - tensora(4)*tensora(3))
                        tensorb(7) = tensora(4)*tensora(8) - tensora(5)*tensora(7)
                        tensorb(8) = -(tensora(1)*tensora(8) - tensora(2)*tensora(7))
                        tensorb(9) = tensora(1)*tensora(5) - tensora(2)*tensora(4)
 
                        ! STEP 2b: computing the determinant of the grad_xi tensor
                        tensorb(tensor_size) = tensora(1)*(tensora(5)*tensora(9) - tensora(6)*tensora(8)) - tensora(2)*(tensora(4) &
                                & *tensora(9) - tensora(6)*tensora(7)) + tensora(3)*(tensora(4)*tensora(8) - tensora(5)*tensora(7))
 
                        if (tensorb(tensor_size) > verysmall) then
                            ! STEP 2c: computing the inverse of grad_xi tensor = F tensorb is the adjoint, tensora becomes F
 
# 148 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 148 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc loop seq 
# 148 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 148 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 148 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
                            do i = 1, tensor_size - 1
                                tensora(i) = tensorb(i)/tensorb(tensor_size)
                            end do
 
                            ! STEP 2d: computing the J = det(F) = 1/det(\grad{\xi})
                            tensorb(tensor_size) = 1._wp/tensorb(tensor_size)
 
                            ! STEP 3: computing F transpose F
                            tensorb(1) = tensora(1)**2 + tensora(2)**2 + tensora(3)**2
                            tensorb(5) = tensora(4)**2 + tensora(5)**2 + tensora(6)**2
                            tensorb(9) = tensora(7)**2 + tensora(8)**2 + tensora(9)**2
                            tensorb(2) = tensora(1)*tensora(4) + tensora(2)*tensora(5) + tensora(3)*tensora(6)
                            tensorb(3) = tensora(1)*tensora(7) + tensora(2)*tensora(8) + tensora(3)*tensora(9)
                            tensorb(6) = tensora(4)*tensora(7) + tensora(5)*tensora(8) + tensora(6)*tensora(9)
                            ! STEP 4: update the btensor, this is consistent with Riemann solvers
# 165 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                                btensor%vf(1)%sf(j, k, l) = tensorb(1)
# 165 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                                btensor%vf(2)%sf(j, k, l) = tensorb(2)
# 165 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                                btensor%vf(3)%sf(j, k, l) = tensorb(5)
# 165 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                                btensor%vf(4)%sf(j, k, l) = tensorb(3)
# 165 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                                btensor%vf(5)%sf(j, k, l) = tensorb(6)
# 165 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                                btensor%vf(6)%sf(j, k, l) = tensorb(9)
# 167 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                            ! store the determinant at the last entry of the btensor
                            btensor%vf(b_size)%sf(j, k, l) = tensorb(tensor_size)
                            ! STEP 5a: updating the Cauchy stress primitive scalar field
                            if (hyper_model == 1) then
                                call s_neohookean_cauchy_solver(btensor%vf, q_prim_vf, g_local, j, k, l)
                            else if (hyper_model == 2) then
                                call s_mooney_rivlin_cauchy_solver(btensor%vf, q_prim_vf, g_local, j, k, l)
                            end if
                            ! STEP 5b: updating the pressure field
                            q_prim_vf(e_idx)%sf(j, k, l) = q_prim_vf(e_idx)%sf(j, k, l) - g_local*q_prim_vf(xiend + 1)%sf(j, k, &
                                      & l)/gamma
                            ! STEP 5c: updating the Cauchy stress conservative scalar field
 
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc loop seq 
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
                            do i = 1, b_size - 1
                                q_cons_vf(strxb + i - 1)%sf(j, k, l) = rho*q_prim_vf(strxb + i - 1)%sf(j, k, l)
                            end do
                        end if
                    end if
                end do
            end do
        end do
 
# 188 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 188 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc end parallel loop
# 188 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 188 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 188 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp end target teams loop
# 188 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
 
    end subroutine s_hyperelastic_rmt_stress_update
 
    !> Compute the neo-Hookean Cauchy stress from the left Cauchy-Green tensor
    subroutine s_neohookean_cauchy_solver(btensor_in, q_prim_vf, G_param, j, k, l)
 
 
# 195 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if MFC_OpenACC
# 195 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc routine seq 
# 195 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif MFC_OpenMP
# 195 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 195 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 195 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp declare target device_type(any) 
# 195 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
        type(scalar_field), dimension(sys_size), intent(inout) :: q_prim_vf
        type(scalar_field), dimension(b_size), intent(inout)   :: btensor_in
        real(wp), intent(in)                                   :: G_param
        integer, intent(in)                                    :: j, k, l
        real(wp)                                               :: trace
        real(wp), parameter                                    :: f13 = 1._wp/3._wp
        integer                                                :: i  !< Generic loop iterators
        ! tensor is the symmetric tensor & calculate the trace of the tensor
        trace = btensor_in(1)%sf(j, k, l) + btensor_in(3)%sf(j, k, l) + btensor_in(6)%sf(j, k, l)
 
        ! Deviatoric left Cauchy-Green tensor: dev(b) = b - (tr(b)/3)*I
# 208 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
            btensor_in(1)%sf(j, k, l) = btensor_in(1)%sf(j, k, l) - f13*trace
# 208 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
            btensor_in(3)%sf(j, k, l) = btensor_in(3)%sf(j, k, l) - f13*trace
# 208 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
            btensor_in(6)%sf(j, k, l) = btensor_in(6)%sf(j, k, l) - f13*trace
# 210 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        ! dividing by the jacobian for neo-Hookean model setting the tensor to the stresses for riemann solver
 
# 211 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 211 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc loop seq 
# 211 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 211 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 211 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
        do i = 1, b_size - 1
            q_prim_vf(strxb + i - 1)%sf(j, k, l) = g_param*btensor_in(i)%sf(j, k, l)/btensor_in(b_size)%sf(j, k, l)
        end do
        ! First invariant strain energy: W = G/2 * (I1 - 3), neo-Hookean model
        q_prim_vf(xiend + 1)%sf(j, k, l) = 0.5_wp*(trace - 3.0_wp)/btensor_in(b_size)%sf(j, k, l)
 
    end subroutine s_neohookean_cauchy_solver
 
    !> Compute the Mooney-Rivlin Cauchy stress from the left Cauchy-Green tensor
    subroutine s_mooney_rivlin_cauchy_solver(btensor_in, q_prim_vf, G_param, j, k, l)
 
 
# 223 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if MFC_OpenACC
# 223 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc routine seq 
# 223 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif MFC_OpenMP
# 223 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 223 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 223 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp declare target device_type(any) 
# 223 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
        type(scalar_field), dimension(sys_size), intent(inout) :: q_prim_vf
        type(scalar_field), dimension(b_size), intent(inout)   :: btensor_in
        real(wp), intent(in)                                   :: G_param
        integer, intent(in)                                    :: j, k, l
        real(wp)                                               :: trace
        real(wp), parameter                                    :: f13 = 1._wp/3._wp
        integer                                                :: i  !< Generic loop iterators
        ! TODO: Make 1D and 2D capable
        trace = btensor_in(1)%sf(j, k, l) + btensor_in(3)%sf(j, k, l) + btensor_in(6)%sf(j, k, l)
 
        ! Deviatoric left Cauchy-Green tensor: dev(b) = b - (tr(b)/3)*I
        btensor_in(1)%sf(j, k, l) = btensor_in(1)%sf(j, k, l) - f13*trace
        btensor_in(3)%sf(j, k, l) = btensor_in(3)%sf(j, k, l) - f13*trace
        btensor_in(6)%sf(j, k, l) = btensor_in(6)%sf(j, k, l) - f13*trace
 
        ! dividing by the jacobian for neo-Hookean model setting the tensor to the stresses for riemann solver
 
# 240 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 240 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc loop seq 
# 240 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 240 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 240 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
        do i = 1, b_size - 1
            q_prim_vf(strxb + i - 1)%sf(j, k, l) = g_param*btensor_in(i)%sf(j, k, l)/btensor_in(b_size)%sf(j, k, l)
        end do
        ! First invariant strain energy: W = G/2 * (I1 - 3), neo-Hookean model
        q_prim_vf(xiend + 1)%sf(j, k, l) = 0.5_wp*(trace - 3.0_wp)/btensor_in(b_size)%sf(j, k, l)
 
    end subroutine s_mooney_rivlin_cauchy_solver
 
    !> Finalize the hyperelastic module
    impure subroutine s_finalize_hyperelastic_module()
 
        integer :: i  !< iterator
        ! Deallocating memory
 
        do i = 1, b_size
#ifdef MFC_DEBUG
# 256 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 256 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 256 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 256 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:256: ', '@:DEALLOCATE(btensor%vf(i)%sf)'
# 256 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 256 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 256 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 256 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 256 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 256 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 256 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc exit data delete(btensor%vf(i)%sf) 
# 256 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 256 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target exit data map(release:btensor%vf(i)%sf) 
# 256 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 256 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    deallocate (btensor%vf(i)%sf)
        end do
#ifdef MFC_DEBUG
# 258 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 258 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 258 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 258 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:258: ', '@:DEALLOCATE(fd_coeff_x_hyper)'
# 258 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 258 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 258 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 258 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 258 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 258 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 258 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc exit data delete(fd_coeff_x_hyper) 
# 258 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 258 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target exit data map(release:fd_coeff_x_hyper) 
# 258 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 258 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    deallocate (fd_coeff_x_hyper)
        if (n > 0) then
#ifdef MFC_DEBUG
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:260: ', '@:DEALLOCATE(fd_coeff_y_hyper)'
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc exit data delete(fd_coeff_y_hyper) 
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target exit data map(release:fd_coeff_y_hyper) 
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    deallocate (fd_coeff_y_hyper)
            if (p > 0) then
#ifdef MFC_DEBUG
# 262 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 262 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 262 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 262 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:262: ', '@:DEALLOCATE(fd_coeff_z_hyper)'
# 262 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 262 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 262 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 262 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 262 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 262 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 262 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc exit data delete(fd_coeff_z_hyper) 
# 262 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 262 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target exit data map(release:fd_coeff_z_hyper) 
# 262 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 262 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    deallocate (fd_coeff_z_hyper)
            end if
        end if
 
    end subroutine s_finalize_hyperelastic_module
 
end module m_hyperelastic