小编yan*_*ick的帖子

我有两个平方矩阵A和B.A是对称的,B是对称的正定.我想计算$ trace(AB ^ { - 1})$.现在,我计算B的Cholesky分解,求解方程$ A = CB $中的C并总结对角元素.

有更有效的方法吗？

我计划使用Eigen.如果矩阵稀疏(A通常是对角线,B通常是带对角线),你能提供一个实现吗？

假设我们有两种类

• 输入类 Input
  • 定义一种类型 result_type
  • 定义 set(result_type)
• 输出类 Output
  • 定义一种类型 result_type
  • 定义 result_type get() const
  • 有许多Input类作为成员变量,其输出依赖于它

给定输出类和几个输入类(任意数),请考虑以下过程:

• 遍历每个输入类并set()使用适当的值调用(事先定义)
• 调用get()输出类并收集结果.

此过程可视为对函数的调用,该函数将输入的值作为参数返回输出值.

编写在一般情况下构造这种可变参数函数的仿函数.

约束是:C++(很可能是C++ 11),可能不同Input::result_type的任意数量的输入类.注意,Input::result_type不一定与Output::result_type.目标首先应该是效率,但如果代码优雅且可读,则会有很大的好处.

细节:对于那些想知道如何Output相关的人Input,可以想象Input有一个result_type get() const方法,它返回你提供的任何东西set().Output然后有一个构造函数,它接受各种Inputs,并将它们(或它们的引用)存储为成员变量.Output::get()然后通过使用其输入get()方法的返回值进行一些数学运算,并返回一些类型的结果Output::result_type.

我适合使用广义加模式gam从mgcv包.我有一个数据表,包含我的因变量Y,一个自变量X,其他自变量Oth和一个两级因子Fac.我想适合以下模型

Y ~ s(X) + Oth

但是,附加约束条件是该s(X)术语仅适用于因子的两个等级之一Fac==1.其他术语Oth应与整个数据相符.

我尝试过探索,s(X,by=Fac)但这会偏向于适合Oth.换句话说,我想表达的信念X涉及到Y只有Fac==1,否则就没有意义模型X.

是否有免费提供的GPL库或一段代码实现了虚构的错误函数:

erfi(x)=-i*erf(i*x)

哪里x有任何复数(或至少是实数)并且i是虚数？