zeye_indices Function

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

Arguments

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

Return Value complex(kind=8)

Source Code

pure function zeye_indices(i,j) result(a)
  integer, intent(in) :: i,j
  complex(8)          :: a
  a = zero
  if(i==j)a=one
end function zeye_indices