Compiler Explorer

Source code

module cobyla_c_mod
implicit none
private
public :: cobyla_c
integer, parameter :: RP = 4

contains

subroutine cobyla_c(array_of_double_ptr, num_elements, array_of_double_value) bind(C)
use, intrinsic :: iso_c_binding
implicit none
type(C_PTR), value :: array_of_double_ptr
integer(C_INT), value :: num_elements
real(C_DOUBLE) :: array_of_double_value(num_elements)
real(C_DOUBLE), pointer :: array_of_double_ptr_loc_intermediate(:)
real(RP), target :: array_of_double_ptr_loc_intermediate2(num_elements)
real(RP), pointer :: array_of_double_ptr_loc(:)
real(RP) :: array_of_double_value_loc(num_elements)

! Deal with the pass by value one
array_of_double_value_loc = real(array_of_double_value, kind(array_of_double_value_loc))
print *, "from value: ", array_of_double_value_loc

! Deal with the ptr
if (C_ASSOCIATED(array_of_double_ptr)) then
  print *, "associated array_of_double_ptr_loc_intermediate: ", associated(array_of_double_ptr_loc_intermediate)
  call C_F_POINTER(array_of_double_ptr, array_of_double_ptr_loc_intermediate, shape=[num_elements])
  print *, "ptr is cast", array_of_double_ptr_loc_intermediate
  array_of_double_ptr_loc_intermediate2 = real(array_of_double_ptr_loc_intermediate, kind(array_of_double_ptr_loc_intermediate2))
  print *, "test", array_of_double_ptr_loc_intermediate2
  array_of_double_ptr_loc => array_of_double_ptr_loc_intermediate2
  print *, "from ptr: ", array_of_double_ptr_loc
else
  array_of_double_ptr_loc => null()
end if

call check_if_present(array_of_double_ptr_loc)

end subroutine

subroutine check_if_present(array_of_values)
implicit none
real(RP), optional :: array_of_values(:)
print *, "value is present? ", present(array_of_values)
end subroutine
end module