サボタージュ禁止のおさぼり日記

これまでと同様です. 説明しないので, 気になったらコピペしてみてください. おまじない library (Rcpp) library (inline) library (RcppArmadillo) library (Matrix) # for sp_mat コレスキー分解 # コレスキー分解 src <- 'arma::mat X = arma::randu(5,5); arma::mat Y = X.t() * X; arma::mat R = arma::chol(Y); arma…