d_diagonal Function

pure function d_diagonal(A) result(dd)

Arguments

Type IntentOptional Attributes Name
real(kind=8), intent(in) :: A(:,:)

Return Value real(kind=8), dimension(size(A,1))


Source Code

pure function d_diagonal(A) result(dd)
  real(8),intent(in)           :: A(:,:)
  real(8),dimension(size(A,1)) :: dd
  integer                      :: i
  do i = 1,size(A,1)
     dd(i) = A(i,i)
  end do
end function d_diagonal