deye_indices Function

pure function deye_indices(i, j) result(a)

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: i
integer, intent(in) :: j

Return Value real(kind=8)


Source Code

pure function deye_indices(i,j) result(a)
  integer, intent(in) :: i,j
  real(8)             :: a
  a = 0d0
  if(i==j)a=1d0
end function deye_indices