避免在嵌套循环中声明过多相同的Rcpp列表元素。
创始人
2024-12-17 11:01:11
0

可以将列表中相同的元素声明放在外层循环中,减少重复声明。例如:

#include 
using namespace Rcpp;

// 嵌套循环
// 在内层循环中重复声明相同的Rcpp列表元素
// 需要优化
// [[Rcpp::export]]
NumericVector nestedLoopWithRepeatedListElements(NumericMatrix x, NumericMatrix y) {
  int n = x.nrow(), m = y.ncol();
  NumericVector ans(n);
  for(int i = 0; i < n; i++) {
    Rcpp::List tmp(m); // 重复声明
    for(int j = 0; j < m; j++) {
      tmp[j] = NumericVector(m); // 重复声明
    }
    for(int j = 0; j < m; j++) {
      for(int k = 0; k < m; k++) {
        tmp[j][k] = x(i, j) * y(k, j); // 计算
      }
    }
    for(int j = 0; j < m; j++) {
      ans[i] += sum(tmp[j]); // 累加
    }
  }
  return ans;
}

// 优化后的代码
// 在外层循环中声明列表元素,避免重复声明
// [[Rcpp::export]]
NumericVector nestedLoopWithOptimizedListElements(NumericMatrix x, NumericMatrix y) {
  int n = x.nrow(), m = y.ncol();
  NumericVector ans(n);
  Rcpp::List tmp(m); // 声明列表元素
  for(int j = 0; j < m; j++) {
    tmp[j] = NumericVector(m); // 声明列表元素
  }
  for(int i = 0; i < n; i++) {
    for(int j = 0; j < m; j++) {
      for(int k = 0; k < m; k++) {
        tmp[j][k] = x(i, j) * y(k, j); // 计算
      }
    }
    for(int j = 0; j < m; j++) {
      ans[i] += sum(tmp[j]); // 累加
    }
  }
  return ans;
}

相关内容

热门资讯

iwatch怎么连接安卓系统,... 你有没有想过,那款时尚又实用的iWatch,竟然只能和iPhone好上好?别急,今天就来给你揭秘,怎...
安卓系统怎么连不上carlif... 安卓系统无法连接CarLife的原因及解决方法随着智能手机的普及,CarLife这一车载互联功能为驾...
oppo手机安卓系统换成苹果系... OPPO手机安卓系统换成苹果系统:现实吗?如何操作?随着智能手机市场的不断发展,用户对于手机系统的需...
安卓平板改windows 系统... 你有没有想过,你的安卓平板电脑是不是也能变身成Windows系统的超级英雄呢?想象在同一个设备上,你...
iphone系统与安卓系统更新... 最近是不是你也遇到了这样的烦恼?手机更新系统总是失败,急得你团团转。别急,今天就来给你揭秘为什么iP...
安卓系统上滑按键,便捷生活与高... 你有没有发现,现在手机屏幕越来越大,操作起来却越来越方便了呢?这都得归功于安卓系统上的那些神奇的上滑...
安卓系统连接耳机模式,蓝牙、有... 亲爱的手机控们,你们有没有遇到过这种情况:手机突然变成了“耳机模式”,明明耳机没插,声音却只从耳机孔...
希沃系统怎么装安卓系统,解锁更... 亲爱的读者们,你是否也像我一样,对希沃一体机上的安卓系统充满了好奇呢?想象在教室里,你的希沃一体机不...
安装了Anaconda之后找不... 在安装Anaconda后,如果找不到Jupyter Notebook,可以尝试以下解决方法:检查环境...
安卓换鸿蒙系统会卡吗,体验流畅... 最近手机圈可是热闹非凡呢!不少安卓用户都在议论纷纷,说鸿蒙系统要来啦!那么,安卓手机换上鸿蒙系统后,...