get_partial_mean_all Function

  function get_partial_mean_all(this, upstream_grad) result(output)
    implicit none
    class(array_type), intent(inout) :: this
    type(array_type), intent(in) :: upstream_grad
    type(array_type) :: output

    real(real32) :: rtmp1
    type(array_type), pointer :: ptr

    ! Calculate the number of elements that were averaged
    rtmp1 = real(product(shape(this%left_operand%val)), real32)

    if(this%is_forward)then
       ptr => sum( upstream_grad ) / rtmp1
    else
       ptr => spread( &
            upstream_grad / rtmp1, &
            dim=1, &
            index=1, &
            ncopies= size(this%left_operand%val, 1) &
       )
       ptr => spread( &
            ptr, &
            dim=2, &
            index=1, &
            ncopies= size(this%left_operand%val, 2) &
       )
    end if
    call output%assign_and_deallocate_source(ptr)
  end function get_partial_mean_all