SIMP_mapping_cpu.f90

Go to the documentation of this file.

 
!
module simp_mapping_cpu
  use num_types, only: rp
  implicit none
  private
 
  public :: simp_mapping_apply_cpu, simp_mapping_apply_backward_cpu
 
contains
 
  subroutine simp_mapping_apply_cpu(f_min, f_max, p, X_out, X_in, n)
    integer, intent(in) :: n
    real(kind=rp), intent(in) :: f_min, f_max, p
    real(kind=rp), dimension(n), intent(out) :: x_out
    real(kind=rp), dimension(n), intent(in) :: x_in
 
    x_out = simp_mapping_kernel(f_min, f_max, p, x_in)
 
  end subroutine simp_mapping_apply_cpu
 
  subroutine simp_mapping_apply_backward_cpu(f_min, f_max, p, &
       sens_out, sens_in, X_in, n)
    integer, intent(in) :: n
    real(kind=rp), intent(in) :: f_min, f_max, p
    real(kind=rp), dimension(n), intent(out) :: sens_out
    real(kind=rp), dimension(n), intent(in) :: sens_in
    real(kind=rp), dimension(n), intent(in) :: x_in
 
    sens_out = simp_mapping_backward_kernel(f_min, f_max, p, sens_in, x_in)
 
  end subroutine simp_mapping_apply_backward_cpu
 
  elemental function simp_mapping_kernel(f_min, f_max, p, X_in) result(X_out)
    real(kind=rp), intent(in) :: f_min, f_max, p
    real(kind=rp), intent(in) :: x_in
    real(kind=rp) :: x_out
 
    x_out = f_min + (f_max - f_min) * x_in**p
 
  end function simp_mapping_kernel
 
  elemental function simp_mapping_backward_kernel(f_min, f_max, p, &
       sens_in, X_in) result(sens_out)
    real(kind=rp), intent(in) :: f_min, f_max, p
    real(kind=rp), intent(in) :: sens_in, x_in
    real(kind=rp) :: sens_out
 
    sens_out = sens_in * (f_max - f_min) * p * x_in**(p - 1.0_rp)
 
  end function simp_mapping_backward_kernel
 
end module simp_mapping_cpu