fmutil_base.F90

Go to the documentation of this file.

!############################################################################################
!     ________  _____  ______________ 
!    / ____/  |/  / / / /_  __/  _/ / 
!   / /_  / /|_/ / / / / / /  / // /  
!  / __/ / /  / / /_/ / / / _/ // /___
! /_/   /_/  /_/\____/ /_/ /___/_____/                                     
!
! Copyright 2020 Bharat Mahajan
!
! Licensed under the Apache License, Version 2.0 (the "License");
! you may not use this file except in compliance with the License.
! You may obtain a copy of the License at
!
!    http://www.apache.org/licenses/LICENSE-2.0
!
! Unless required by applicable law or agreed to in writing, software
! distributed under the License is distributed on an "AS IS" BASIS,
! WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
! See the License for the specific language governing permissions and
! limitations under the License.
!
!
!############################################################################################
 
    
module fmutilbase
 
    use, intrinsic :: iso_fortran_env, only: real64, real128
 
    implicit none
    
    integer, parameter, public :: dp = real64
    integer, parameter, public :: qp = real128
    integer, parameter, public :: wp = dp
 
    real(wp), parameter, public :: eps = epsilon(1.0_wp)
    
    real(wp), parameter, public :: inf = huge(1.0_wp)
    
    integer, private ::  ind
    logical, dimension(*), parameter :: nzind3 = [ (ind == (3*(ind/3)+(ind/3+1)), ind=1,8), .true.]
    real(wp), dimension(*,*), parameter, public :: i3 = &
                         reshape([merge([(1.0_wp, ind=1,9)], [(0.0_wp, ind=1,9)], nzind3)],[3,3])
    
    
    logical, dimension(*), parameter :: nzind6 = [ (ind == (6*(ind/6)+(ind/6+1)), ind=1,6**2-1), .true.]
    real(wp), dimension(*,*), parameter, public :: i6 = &
                         reshape([merge([(1.0_wp, ind=1,6**2)], [(0.0_wp, ind=1,6**2)], nzind6)],[6,6])
 
 
 
 
 
    contains
    
 
    
end module fmutilbase