QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 9000|回复: 0
打印 上一主题 下一主题

Lu数值计算扩展动态库LuMath & 小矩阵大运算量效率测试

[复制链接]
字体大小: 正常 放大
forcal 实名认证       

45

主题

3

听众

282

积分

升级  91%

  • TA的每日心情
    难过
    2012-8-27 18:22
  • 签到天数: 1 天

    [LV.1]初来乍到

    跳转到指定楼层
    1#
    发表于 2011-11-10 21:43 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    Lu数值计算扩展动态库LuMath V1.0:http://www.forcal.net/sysm/lu1/luhtm/lumath.htm
    6 r1 d& \" P- a: _" e) |# O: ^$ C5 }' t1 f
    LuMath 库的数值算法是以矩阵为基础,并进行了运算符重载,具有内存消耗低、执行效率高、代码简洁、实用性强的特点。( m1 L' N+ U4 b! S' k2 g3 i

    - _& K: |8 |. e7 ?& _例子:
    2 y( {6 C( l+ v! u9 X+ D2 B* \8 `0 q& r; P$ ?5 l9 H6 J+ [  J& c
    matlab2009a代码:
    1. clear all
      : f1 g  `1 g. C; z$ m
    2. clc( E! L$ R0 @9 k9 f
    3. tic
      % d3 T- g* D9 I+ l8 x
    4. k = zeros(5,5); % //生成5×5全0矩阵
      / P5 W$ Z/ c, a4 n) A\" F
    5. % 循环计算以下程序段100000次:) O) x# D+ ~% H9 P% g- H8 p
    6. for m = 1:100000  ]9 Q2 q' H& x) _3 d
    7.     a = rand(5,7);
      - t% O6 D3 `6 {: s7 j
    8.     b = rand(7,5);%//生成5×7矩阵a,7×5矩阵b,用0~1之间的随机数初始化
      , `/ R. W( J. }3 W7 A% f
    9.     k = k + a * b + a(1:5, 2:6) * b(2:6, 1:5) - a(:, 7) * b(3, :);
      8 \+ L! M; V  L/ r3 D
    10. end2 l: T/ x+ v/ y* k8 l
    11. k
      3 K9 J1 W* Q\" x- c
    12. toc
    复制代码
    多次运行,结果约为:
    1. k =1 F6 `9 T8 e* }

    2. . B- K4 u& Y' x2 F\" c( V, F6 W\" N
    3.   1.0e+005 *8 x) _1 s' g+ h) a- X; f  E7 Y& |

    4. ) L# w5 g3 @  f  @
    5.     2.7530    2.7525    2.7509    2.7539    2.75291 K; {$ l# g: {2 J
    6.     2.7521    2.7494    2.7493    2.7542    2.7476
      ( S1 s. W( b5 s) |0 K5 r7 b4 N$ \% c
    7.     2.7549    2.7524    2.7527    2.7564    2.7517
      ( [\" g% u/ ?- q! D2 ]
    8.     2.7530    2.7523    2.7533    2.7525    2.7502
      1 \2 B0 X8 ~+ Y8 `% }- [1 G
    9.     2.7527    2.7521    2.7505    2.7519    2.7523
      0 [2 b7 n- F+ L% ~, w
    10. . h6 W2 |1 W. M
    11. Elapsed time is 2.126704 seconds.
    复制代码
    Lu代码:
    1. !!!using["math"];
    2. ( L9 [9 e/ d+ s$ a% D
    3. main(:t0,k,i,a,b)=3 W- t2 v' C3 N3 O  L* |0 a* ]
    4. {$ ^) p1 e; d\\" ]! j+ s/ [9 t6 s2 v* X
    5.     t0=clock(),2 {( G! b: H' p, ~9 a! I
    6.     k=zeros(5,5),         //生成5×5矩阵k,初始化为0
    7. % q3 \5 X. K7 w( o4 g
    8.     i=0,(++i<=100000).while{//循环计算100000次- L8 m- ]8 @9 B
    9.         a=rand(5,7), b=rand(7,5), //生成5×7矩阵a,7×5矩阵b,用0~1之间的随机数初始化. ~- |7 }, v6 `1 z6 i: O
    10.         k.=k+a*b+a(0,4:1,5)*b(1,5:0,4)-a(all:6)*b(2:all) //计算k=k+a*b+a(0,4:1,5)*b(1,5:0,4)+a(all:6)*b(3:all)
    11. ; o' C5 |% {9 ]4 Y, t
    12.     },
    13. ( e$ m& Q5 d8 |+ b; R0 \9 v2 @7 c6 _
    14.     k.outa(),             //输出矩阵k
    15. % f* t/ `' u7 ~
    16.     [clock()-t0]/1000.    //得到计算时间,秒; `/ _! b4 _: B' T9 U6 K' ^/ v8 S
    17. };
    多次运行,结果约为:
    1.          275024         274895         275054         275083         2750242 d$ i4 e\" V6 w\" d+ l8 X+ d
    2.          274965         275036         275091         275112         2749319 Q/ R1 Y& f6 w) I5 O1 K. L
    3.          274975         275012         275020         274963         2749395 ~. p- f% j( i
    4.          274981         274976         275117         275065         275044+ y% l\" z! N5 u' `# G( a. f9 O+ @
    5.          274920         274906         275005         274972         274920
      # }1 }5 P$ h$ P& k6 y
    6. 1.25 秒
    复制代码
    Lu快于matlab的原因应该仍在于内存管理效率较高。
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-9-1 04:02 , Processed in 0.581499 second(s), 49 queries .

    回顶部