在R中的矩阵行上的两级循环

Question

我需要将下面的迭代转换为矩阵运算。

a = certain_matrix(3,5)
mat = matrix(0, 3,3)

for ( i in 1:3 )
  for ( j in 1:3 )
    mat[i,j] = certain_distance(a[i,], a[j,])

我大概需要使用apply（）。但是，我无法掌握

中的行索引

apply(a, 1, function(x) XXX)

那么如何快速执行该操作？