楼主: B.N.F
1902 1

Dow loop快还是temp array快? [推广有奖]

  • 1关注
  • 0粉丝

博士生

54%

还不是VIP/贵宾

-

威望
0
论坛币
2124 个
通用积分
247.2645
学术水平
15 点
热心指数
15 点
信用等级
15 点
经验
4683 点
帖子
208
精华
0
在线时间
203 小时
注册时间
2012-8-19
最后登录
2024-3-3

10论坛币
  1. data work.a;
  2. array a(10) a1-a10;
  3. do i=1 to 1000000;
  4. do j=1 to 10;
  5. a(j)=normal(111);
  6. end;
  7. output;
  8. end;
  9. run;
  10. data _null_;
  11. sasfile work.a load;
  12. /*读入work.a的数据用于跟array部分一致的基准基础*/
  13. array x[1000000,10] _temporary_;
  14. array a[10] a1-a10;
  15. do i=1 to 1000000;
  16. set work.a;
  17. do j=1 to 10;
  18. x[i,j]=a[j];
  19. end;
  20. end;
  21. /* 代码正文 */
  22. do m=1 to 100;
  23. do until (eof);
  24. set work.a (firstobs=11) end=eof;
  25. y=20;
  26. a11=lag(a1);
  27. a22=lag2(a2);
  28. a55=lag5(a5);
  29. a77=lag7(a7);
  30. a1010=lag10(a10);
  31. pi = constant('PI');
  32. y1 = (a11-(5.1/(4*pi**2))*a22*a1010+5*a55/pi-6)**2;
  33. y2 = 10*(1-1/(8*pi))*cos(a77);
  34. rec= (y1+y2+10);
  35. end;
  36. put _all_;
  37. end;
  38. run;
  39. /*加array temp查看temp array在pdv的作用和效率*/
  40. data _null_;
  41. /*数据准备成array*/
  42. array x[1000000,10] _temporary_;
  43. array a[10] a1-a10;
  44. do i=1 to 1000000;
  45. set work.a;
  46. do j=1 to 10;
  47. x[i,j]=a[j];
  48. end;
  49. end;
  50. /*代码正文*/
  51. do m=1 to 100;
  52. do j=11 to 1000000;/*遍历全数据array*/
  53. y=20;
  54. pi = constant('PI');
  55. y1 = (x(j-1,1)-(5.1/(4*pi**2))*x(j-2,2)*x(j-10,10)+5*x(j-5,5)/pi-6)**2;
  56. y2 = 10*(1-1/(8*pi))*cos(x(j-7,7));
  57. rec= (y1+y2+10);
  58. end;
  59. end;
  60. run;
  61. sasfile work.a close;
复制代码

几个问题:
1、我用sasfile的dow loop的快,还是temp array快?因为考虑到一个是如果是temp array的话,数据全部进内存,但是如果读数据集的话会不会慢?另一个是如果temp array进内存的话,那么PDV是不是上面的情况dow loop dataset的PDV就会比temp array的短很多,那么会快不会快?
2、为什么dow loop的代码我想循环100遍set 不行?
3、有没有更好的方法,因为得重复计算data步,也就是说每计算一次都得set一遍,所以考虑array,但是array每次都得赋值。

关键词:array Temp loop ARR TEM normal null
沙发
B.N.F 发表于 2014-8-19 16:15:37 |只看作者 |坛友微信交流群
顶一下

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-20 09:16