使用Brobdingnag包处理矩阵

Question

我需要构建一个包含极小条目的矩阵。 到目前为止，我意识到定义我需要的矩阵类型的最快方法是：

1. 定义坐标的矢量化函数：

   func = function(m,n){...}
   

2. 使用outer：

   组合每个可能的坐标

   matrix = outer(1:100,1:100,FUN=func)
   

我必须处理使用brob数字在func环境中工作的非常小的数字，因此它的输出将与brob的类型相同：

typeof(func(0:100,0:100) )
[1] "S4"

如果我在我的函数func中直接插入两个向量0：100，它会返回一个brobs向量，但如果我尝试将它与外部一起使用，我会得到错误：

Error in outer(1:100, 1:100, FUN = func) : invalid first argument

我想这是因为包Brobdingnag可以某种方式处理向量而不是矩阵。这样对吗？有没有办法让它发挥作用？