A <- matrix( c(1:48),ncol = 16 ); 
B <- matrix( c(1:80),ncol = 16 ); 
 
mutual <- function(delta, a, FUN = "*"){
  delta  <- as.matrix( delta );
  a      <- as.matrix( a     );
  result <- matrix( outer(delta, a, FUN), nrow = nrow(delta) )
  return( result )
}
# > A
#      [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16]
# [1,]    1    4    7   10   13   16   19   22   25    28    31    34    37    40    43    46
# [2,]    2    5    8   11   14   17   20   23   26    29    32    35    38    41    44    47
# [3,]    3    6    9   12   15   18   21   24   27    30    33    36    39    42    45    48
# > B
#      [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16]
# [1,]    1    6   11   16   21   26   31   36   41    46    51    56    61    66    71    76
# [2,]    2    7   12   17   22   27   32   37   42    47    52    57    62    67    72    77
# [3,]    3    8   13   18   23   28   33   38   43    48    53    58    63    68    73    78
# [4,]    4    9   14   19   24   29   34   39   44    49    54    59    64    69    74    79
# [5,]    5   10   15   20   25   30   35   40   45    50    55    60    65    70    75    80
lapply( 1:ncol(A) , function(x) mutual( A[,x],B[,x] ) )
# > lapply( 1:ncol(A) , function(x) mutual( A[,x],B[,x] ) )
# [[1]]
#      [,1] [,2] [,3] [,4] [,5]
# [1,]    1    2    3    4    5
# [2,]    2    4    6    8   10
# [3,]    3    6    9   12   15
# 
# [[2]]
#      [,1] [,2] [,3] [,4] [,5]
# [1,]   24   28   32   36   40
# [2,]   30   35   40   45   50
# [3,]   36   42   48   54   60
# 
# [[3]]
#      [,1] [,2] [,3] [,4] [,5]
# [1,]   77   84   91   98  105
# [2,]   88   96  104  112  120
# [3,]   99  108  117  126  135   ......
 
 

QSA8LSBtYXRyaXgoIGMoMTo0OCksbmNvbCA9IDE2ICk7IApCIDwtIG1hdHJpeCggYygxOjgwKSxuY29sID0gMTYgKTsgCgptdXR1YWwgPC0gZnVuY3Rpb24oZGVsdGEsIGEsIEZVTiA9ICIqIil7CiAgZGVsdGEgIDwtIGFzLm1hdHJpeCggZGVsdGEgKTsKICBhICAgICAgPC0gYXMubWF0cml4KCBhICAgICApOwogIHJlc3VsdCA8LSBtYXRyaXgoIG91dGVyKGRlbHRhLCBhLCBGVU4pLCBucm93ID0gbnJvdyhkZWx0YSkgKQogIHJldHVybiggcmVzdWx0ICkKfQojID4gQQojICAgICAgWywxXSBbLDJdIFssM10gWyw0XSBbLDVdIFssNl0gWyw3XSBbLDhdIFssOV0gWywxMF0gWywxMV0gWywxMl0gWywxM10gWywxNF0gWywxNV0gWywxNl0KIyBbMSxdICAgIDEgICAgNCAgICA3ICAgMTAgICAxMyAgIDE2ICAgMTkgICAyMiAgIDI1ICAgIDI4ICAgIDMxICAgIDM0ICAgIDM3ICAgIDQwICAgIDQzICAgIDQ2CiMgWzIsXSAgICAyICAgIDUgICAgOCAgIDExICAgMTQgICAxNyAgIDIwICAgMjMgICAyNiAgICAyOSAgICAzMiAgICAzNSAgICAzOCAgICA0MSAgICA0NCAgICA0NwojIFszLF0gICAgMyAgICA2ICAgIDkgICAxMiAgIDE1ICAgMTggICAyMSAgIDI0ICAgMjcgICAgMzAgICAgMzMgICAgMzYgICAgMzkgICAgNDIgICAgNDUgICAgNDgKIyA+IEIKIyAgICAgIFssMV0gWywyXSBbLDNdIFssNF0gWyw1XSBbLDZdIFssN10gWyw4XSBbLDldIFssMTBdIFssMTFdIFssMTJdIFssMTNdIFssMTRdIFssMTVdIFssMTZdCiMgWzEsXSAgICAxICAgIDYgICAxMSAgIDE2ICAgMjEgICAyNiAgIDMxICAgMzYgICA0MSAgICA0NiAgICA1MSAgICA1NiAgICA2MSAgICA2NiAgICA3MSAgICA3NgojIFsyLF0gICAgMiAgICA3ICAgMTIgICAxNyAgIDIyICAgMjcgICAzMiAgIDM3ICAgNDIgICAgNDcgICAgNTIgICAgNTcgICAgNjIgICAgNjcgICAgNzIgICAgNzcKIyBbMyxdICAgIDMgICAgOCAgIDEzICAgMTggICAyMyAgIDI4ICAgMzMgICAzOCAgIDQzICAgIDQ4ICAgIDUzICAgIDU4ICAgIDYzICAgIDY4ICAgIDczICAgIDc4CiMgWzQsXSAgICA0ICAgIDkgICAxNCAgIDE5ICAgMjQgICAyOSAgIDM0ICAgMzkgICA0NCAgICA0OSAgICA1NCAgICA1OSAgICA2NCAgICA2OSAgICA3NCAgICA3OQojIFs1LF0gICAgNSAgIDEwICAgMTUgICAyMCAgIDI1ICAgMzAgICAzNSAgIDQwICAgNDUgICAgNTAgICAgNTUgICAgNjAgICAgNjUgICAgNzAgICAgNzUgICAgODAKbGFwcGx5KCAxOm5jb2woQSkgLCBmdW5jdGlvbih4KSBtdXR1YWwoIEFbLHhdLEJbLHhdICkgKQojID4gbGFwcGx5KCAxOm5jb2woQSkgLCBmdW5jdGlvbih4KSBtdXR1YWwoIEFbLHhdLEJbLHhdICkgKQojIFtbMV1dCiMgICAgICBbLDFdIFssMl0gWywzXSBbLDRdIFssNV0KIyBbMSxdICAgIDEgICAgMiAgICAzICAgIDQgICAgNQojIFsyLF0gICAgMiAgICA0ICAgIDYgICAgOCAgIDEwCiMgWzMsXSAgICAzICAgIDYgICAgOSAgIDEyICAgMTUKIyAKIyBbWzJdXQojICAgICAgWywxXSBbLDJdIFssM10gWyw0XSBbLDVdCiMgWzEsXSAgIDI0ICAgMjggICAzMiAgIDM2ICAgNDAKIyBbMixdICAgMzAgICAzNSAgIDQwICAgNDUgICA1MAojIFszLF0gICAzNiAgIDQyICAgNDggICA1NCAgIDYwCiMgCiMgW1szXV0KIyAgICAgIFssMV0gWywyXSBbLDNdIFssNF0gWyw1XQojIFsxLF0gICA3NyAgIDg0ICAgOTEgICA5OCAgMTA1CiMgWzIsXSAgIDg4ICAgOTYgIDEwNCAgMTEyICAxMjAKIyBbMyxdICAgOTkgIDEwOCAgMTE3ICAxMjYgIDEzNSAgIC4uLi4uLgoK