get_partial_max_right_val Subroutine

  pure subroutine get_partial_max_right_val(this, upstream_grad, output)
    implicit none
    class(array_type), intent(in) :: this
    real(real32), dimension(:,:), intent(in) :: upstream_grad
    real(real32), dimension(:,:), intent(out) :: output

    integer :: i, j

    output = 0._real32
    do concurrent(i=1:size(upstream_grad,1), j=1:size(upstream_grad,2))
       if( abs( this%val(i,j) - this%right_operand%val(i,j) ) .lt. 1.E-6_real32 ) then
          output(i,j) = upstream_grad(i,j)
       end if
    end do

  end subroutine get_partial_max_right_val