torch_tensor_to_array_real64_4d Subroutine

  subroutine torch_tensor_to_array_real64_4d(tensor, data_out, sizes)
    use, intrinsic :: iso_c_binding, only : c_int, c_int64_t, c_loc
    use, intrinsic :: iso_fortran_env, only : real64, int64
    type(torch_tensor), intent(in) :: tensor  !! Returned tensor
    real(kind=real64), pointer, intent(out) :: data_out(:,:,:,:)  !! Pointer to tensor data
    integer, optional, intent(in) :: sizes(4)  !! Number of entries for each rank
    integer(kind=int64), allocatable :: my_shape(:)  !! Number of entries for each rank

    ! Local data
    integer(c_int), parameter :: dtype = torch_kFloat64  !! Data type
    type(c_ptr) :: cptr

    my_shape = tensor%get_shape()

    if (present(sizes)) then
      if (.not. all(my_shape == sizes)) then
        write(*,*) 'Error :: sizes argument does not match shape of tensor'
        write(*,'(A, 4(I0, " "), A)') 'sizes        :: [ ', sizes(:), ']'
        write(*,'(A, 4(I0, " "), A)') 'tensor shape :: [ ', my_shape(:), ']'
        stop 1
      end if
    end if

    ! Have the data_out array point to the Tensor data
    cptr = torch_to_blob_c(tensor%p, dtype)
    call c_f_pointer(cptr, data_out, my_shape)

  end subroutine torch_tensor_to_array_real64_4d