Я имею матрицу n x p и хотел бы вычислить матрицу B n x n, определенную как

B[i, j] = f(A[i,], A[j,])

где f - это функция, которая принимает аргументы соответствующей размерности. Есть ли хитрый трюк, чтобы вычислить это в R? f является симметричным и положительно определенным (если это может помочь в вычислении).

РЕДАКТИРОВАТЬ: Praneet попросил указать f. Это хороший момент. Хотя я думаю, что было бы интересно иметь эффективное решение для любой функции, я получил бы много пользы от эффективных вычислений в важном случае, когда f (x, y) равно base :: norm (xy, type = 'F'). ).