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"
 
# 207 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 232 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 243 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 245 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 255 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 283 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 293 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 303 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 312 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 329 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 339 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 346 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 352 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 358 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 364 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 370 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 376 "/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"
 
# 192 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 213 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 241 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 256 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 266 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 275 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 291 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 301 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 308 "/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
 
# 21 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 37 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 50 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 76 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 91 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 102 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 115 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 143 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 154 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 165 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 176 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 187 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 198 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 208 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 214 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 220 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 226 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 232 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 234 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
# 235 "/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.
# 63 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 81 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 88 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 111 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 127 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 153 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 159 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 167 "/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        !< definitions of the derived types
 
    use m_global_parameters    !< definitions of the global parameters
 
    use m_variables_conversion !< state variables type conversion procedures
 
    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 !<
 
# 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(btensor) 
# 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 (btensor) 
# 28 "/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
 
# 33 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 33 "/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) 
# 33 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 33 "/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) 
# 33 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
    real(wp), allocatable, dimension(:) :: gs_hyper
 
# 35 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 35 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc declare create(Gs_hyper) 
# 35 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 35 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp declare target (Gs_hyper) 
# 35 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
 
contains
 
    !>  The following subroutine handles the calculation of the btensor.
        !!   The calculation of the btensor takes qprimvf.
        !! calculate the grad_xi, grad_xi is a nxn tensor
        !! calculate the inverse of grad_xi to obtain F, F is a nxn tensor
        !! calculate the FFtranspose to obtain the btensor, btensor is nxn tensor
        !! btensor is symmetric, save the data space
    impure subroutine s_initialize_hyperelastic_module
        integer :: i !< generic iterator
 
#ifdef MFC_DEBUG
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:48: ', '@:ALLOCATE(btensor%vf(1:b_size))'
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    allocate (btensor%vf(1:b_size))
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc enter data create(btensor%vf) 
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target enter data map(always,alloc:btensor%vf) 
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
        do i = 1, b_size
#ifdef MFC_DEBUG
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:50: ', '@:ALLOCATE(btensor%vf(i)%sf(0:m, 0:n, 0:p))'
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    allocate (btensor%vf(i)%sf(0:m, 0:n, 0:p))
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 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 enter data create(btensor%vf(i)%sf) 
# 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 enter data map(always,alloc:btensor%vf(i)%sf) 
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
        end do
#ifdef _CRAYFTN
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        integer :: macros_setup_vfs_i
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#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: ', '@:ACC_SETUP_VFs(btensor)'
# 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"
 
# 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 copyin(btensor) 
# 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(to:btensor) 
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 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 copyin(btensor%vf) 
# 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(to:btensor%vf) 
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
            if (allocated(btensor%vf)) then
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                do macros_setup_vfs_i = lbound(btensor%vf, 1), ubound(btensor%vf, 1)
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                    if (associated(btensor%vf(macros_setup_vfs_i)%sf)) then
# 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 copyin(btensor%vf(macros_setup_vfs_i)) 
# 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(to:btensor%vf(macros_setup_vfs_i)) 
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 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 copyin(btensor%vf(macros_setup_vfs_i)%sf) 
# 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(to:btensor%vf(macros_setup_vfs_i)%sf) 
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                    end if
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                end do
# 52 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
            end if
# 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
 
#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(Gs_hyper(1:num_fluids))'
# 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 (gs_hyper(1:num_fluids))
# 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(Gs_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:Gs_hyper) 
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
 
# 55 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 55 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc loop seq 
# 55 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 55 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 55 "/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
 
# 59 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 59 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc update device(Gs_hyper) 
# 59 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 59 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target update to(Gs_hyper) 
# 59 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
 
#ifdef MFC_DEBUG
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:61: ', '@:ALLOCATE(fd_coeff_x_hyper(-fd_number:fd_number, 0:m))'
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    allocate (fd_coeff_x_hyper(-fd_number:fd_number, 0:m))
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc enter data create(fd_coeff_x_hyper) 
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target enter data map(always,alloc:fd_coeff_x_hyper) 
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
        if (n > 0) then
#ifdef MFC_DEBUG
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:63: ', '@:ALLOCATE(fd_coeff_y_hyper(-fd_number:fd_number, 0:n))'
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    allocate (fd_coeff_y_hyper(-fd_number:fd_number, 0:n))
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc enter data create(fd_coeff_y_hyper) 
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target enter data map(always,alloc:fd_coeff_y_hyper) 
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
        end if
        if (p > 0) then
#ifdef MFC_DEBUG
# 66 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 66 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 66 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 66 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:66: ', '@:ALLOCATE(fd_coeff_z_hyper(-fd_number:fd_number, 0:p))'
# 66 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 66 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 66 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 66 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 66 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    allocate (fd_coeff_z_hyper(-fd_number:fd_number, 0:p))
# 66 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 66 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 66 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 66 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc enter data create(fd_coeff_z_hyper) 
# 66 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 66 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target enter data map(always,alloc:fd_coeff_z_hyper) 
# 66 "/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)
 
# 72 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 72 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc update device(fd_coeff_x_hyper) 
# 72 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 72 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target update to(fd_coeff_x_hyper) 
# 72 "/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)
 
# 76 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 76 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc update device(fd_coeff_y_hyper) 
# 76 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 76 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target update to(fd_coeff_y_hyper) 
# 76 "/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)
 
# 81 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 81 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc update device(fd_coeff_z_hyper) 
# 81 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 81 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target update to(fd_coeff_z_hyper) 
# 81 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
        end if
 
    end subroutine s_initialize_hyperelastic_module
 
    !>  The following subroutine handles the calculation of the btensor.
        !!   The calculation of the btensor takes qprimvf.
        !! @param q_cons_vf Conservative variables
        !! @param q_prim_vf Primitive variables
        !! calculate the grad_xi, grad_xi is a nxn tensor
        !! calculate the inverse of grad_xi to obtain F, F is a nxn tensor
        !! calculate the FFtranspose to obtain the btensor, btensor is nxn tensor
        !! btensor is symmetric, save the data space
    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
# 101 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
            real(wp), dimension(tensor_size) :: tensora, tensorb
# 103 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 107 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
            real(wp), dimension(num_fluids) :: alpha_k, alpha_rho_k
# 109 "/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
 
 
# 114 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 114 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 114 "/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) 
# 114 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 114 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    
# 114 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 114 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 114 "/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) 
# 114 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 114 "/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 ) G_K = 0._wp
 
                    if (g_local > verysmall) then
 
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc loop seq 
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 129 "/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
 
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc loop seq 
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 138 "/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
 
# 172 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 172 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc loop seq 
# 172 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 172 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 172 "/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
# 189 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                                btensor%vf(1)%sf(j, k, l) = tensorb(1)
# 189 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                                btensor%vf(2)%sf(j, k, l) = tensorb(2)
# 189 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                                btensor%vf(3)%sf(j, k, l) = tensorb(5)
# 189 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                                btensor%vf(4)%sf(j, k, l) = tensorb(3)
# 189 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                                btensor%vf(5)%sf(j, k, l) = tensorb(6)
# 189 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
                                btensor%vf(6)%sf(j, k, l) = tensorb(9)
# 191 "/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)
                            elseif (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
 
# 203 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 203 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc loop seq 
# 203 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 203 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 203 "/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
 
# 213 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 213 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 213 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc end parallel loop
# 213 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 213 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    
# 213 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 213 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp end target teams loop
# 213 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 213 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
    end subroutine s_hyperelastic_rmt_stress_update
 
    !>  The following subroutine handles the calculation of the btensor.
        !!   The calculation of the btensor takes qprimvf.
        !! @param btensor_in Left Cauchy-Green deformation tensor
        !! @param q_prim_vf Primitive variables
        !! @param G_param Elastic shear modulus
        !! @param j x-direction cell index
        !! @param k y-direction cell index
        !! @param l z-direction cell index
        !! calculate the grad_xi, grad_xi is a nxn tensor
        !! calculate the inverse of grad_xi to obtain F, F is a nxn tensor
        !! calculate the FFtranspose to obtain the btensor, btensor is nxn tensor
        !! btensor is symmetric, save the data space
    subroutine s_neohookean_cauchy_solver(btensor_in, q_prim_vf, G_param, j, k, l)
 
# 229 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if MFC_OpenACC
# 229 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc routine seq 
# 229 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif MFC_OpenMP
# 229 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    
# 229 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 229 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp declare target device_type(any) 
# 229 "/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)
 
        ! calculate the deviatoric of the tensor
# 244 "/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
# 244 "/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
# 244 "/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
# 246 "/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
 
# 248 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 248 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc loop seq 
# 248 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 248 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 248 "/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
        ! compute the invariant without the elastic modulus
        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
 
    !>  The following subroutine handles the calculation of the btensor.
        !!   The calculation of the btensor takes qprimvf.
        !! @param btensor_in Left Cauchy-Green deformation tensor
        !! @param q_prim_vf Primitive variables
        !! @param G_param Elastic shear modulus
        !! @param j x-direction cell index
        !! @param k y-direction cell index
        !! @param l z-direction cell index
        !! calculate the grad_xi, grad_xi is a nxn tensor
        !! calculate the inverse of grad_xi to obtain F, F is a nxn tensor
        !! calculate the FFtranspose to obtain the btensor, btensor is nxn tensor
        !! btensor is symmetric, save the data space
    subroutine s_mooney_rivlin_cauchy_solver(btensor_in, q_prim_vf, G_param, j, k, l)
 
# 272 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if MFC_OpenACC
# 272 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc routine seq 
# 272 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif MFC_OpenMP
# 272 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    
# 272 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 272 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp declare target device_type(any) 
# 272 "/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 this 1D and 2D capable
        ! 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)
 
        ! calculate the deviatoric of the tensor
        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
 
# 293 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 293 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc loop seq 
# 293 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 293 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 293 "/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
        ! compute the invariant without the elastic modulus
        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
 
    !> @brief Deallocates memory for hyperelastic deformation tensor and finite-difference coefficients.
    impure subroutine s_finalize_hyperelastic_module()
 
        integer :: i !< iterator
 
        ! Deallocating memory
        do i = 1, b_size
#ifdef MFC_DEBUG
# 311 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 311 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 311 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 311 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:311: ', '@:DEALLOCATE(btensor%vf(i)%sf)'
# 311 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 311 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 311 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 311 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 311 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 311 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 311 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc exit data delete(btensor%vf(i)%sf) 
# 311 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 311 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target exit data map(release:btensor%vf(i)%sf) 
# 311 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 311 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    deallocate (btensor%vf(i)%sf)
        end do
#ifdef MFC_DEBUG
# 313 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 313 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 313 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 313 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:313: ', '@:DEALLOCATE(fd_coeff_x_hyper)'
# 313 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 313 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 313 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 313 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 313 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 313 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 313 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc exit data delete(fd_coeff_x_hyper) 
# 313 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 313 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target exit data map(release:fd_coeff_x_hyper) 
# 313 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 313 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    deallocate (fd_coeff_x_hyper)
        if (n > 0) then
#ifdef MFC_DEBUG
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:315: ', '@:DEALLOCATE(fd_coeff_y_hyper)'
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc exit data delete(fd_coeff_y_hyper) 
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target exit data map(release:fd_coeff_y_hyper) 
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    deallocate (fd_coeff_y_hyper)
            if (p > 0) then
#ifdef MFC_DEBUG
# 317 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    block
# 317 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        use iso_fortran_env, only: output_unit
# 317 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 317 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        print *, 'm_hyperelastic.fpp:317: ', '@:DEALLOCATE(fd_coeff_z_hyper)'
# 317 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 317 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
        call flush (output_unit)
# 317 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
    end block
# 317 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 317 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
 
# 317 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#if defined(MFC_OpenACC)
# 317 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$acc exit data delete(fd_coeff_z_hyper) 
# 317 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#elif defined(MFC_OpenMP)
# 317 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
!$omp target exit data map(release:fd_coeff_z_hyper) 
# 317 "/home/runner/work/MFC/MFC/src/simulation/m_hyperelastic.fpp"
#endif
# 317 "/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