subtract_real2d Module Function

  module function subtract_real2d(a, b) result(c)
    !! Subtract a real array from an autodiff array
    implicit none
    class(array_type), intent(in), target :: a
    real(real32), dimension(:,:), intent(in) :: b
    type(array_type), pointer :: c

    integer :: i, s

    c => a%create_result()
    do concurrent(s = 1:size(a%val,2), i = 1:size(a%val,1))
       c%val(i,s) = a%val(i,s) - b(i,s)
    end do

    c%get_partial_left => get_partial_add
    c%get_partial_left_val => get_partial_add_val
    if(a%requires_grad)then
       c%requires_grad = .true.
       c%is_forward = a%is_forward
       c%operation = 'subtract_scalar'
       c%left_operand => a
       c%owns_left_operand = a%is_temporary
    end if
  end function subtract_real2d