get_partial_gaussian_val Subroutine

  pure subroutine get_partial_gaussian_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

    output = -upstream_grad * this%val * &
         this%left_operand%val / ( this%right_operand%val(2,1)**2._real32 )

    !(this%val - this%right_operand%val(1,1)) * &
    !     exp( -0.5_real32 * ( &
    !          (this%val - this%right_operand%val(1,1)) / &
    !          this%right_operand%val(2,1) &
    !     )**2 )
  end subroutine get_partial_gaussian_val