unit_matrix Subroutine

public subroutine unit_matrix(mat)

Arguments

Type IntentOptional Attributes Name
real(kind=8), intent(out), dimension(:,:) :: mat

Called by

proc~~unit_matrix~~CalledByGraph proc~unit_matrix unit_matrix proc~broyden1~2 broyden1 proc~broyden1~2->proc~unit_matrix

Source Code

  subroutine unit_matrix(mat)
    real(8), dimension(:,:), intent(out) :: mat
    integer :: i,n
    n=min(size(mat,1),size(mat,2))
    mat(:,:)=0.0
    do i=1,n
       mat(i,i)=1.0
    end do
  end subroutine unit_matrix