- 在线时间
- 13 小时
- 最后登录
- 2013-12-8
- 注册时间
- 2010-5-13
- 听众数
- 3
- 收听数
- 0
- 能力
- 0 分
- 体力
- 399 点
- 威望
- 11 点
- 阅读权限
- 30
- 积分
- 282
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 97
- 主题
- 45
- 精华
- 0
- 分享
- 0
- 好友
- 1
升级   91% TA的每日心情 | 难过 2012-8-27 18:22 |
---|
签到天数: 1 天 [LV.1]初来乍到
 |
Lu数值计算扩展动态库LuMath V1.0:http://www.forcal.net/sysm/lu1/luhtm/lumath.htm @! Q0 |8 a0 f1 ]+ t" a3 _
+ M' I+ R# n* y: p4 V
LuMath 库的数值算法是以矩阵为基础,并进行了运算符重载,具有内存消耗低、执行效率高、代码简洁、实用性强的特点。
+ M8 |/ i a8 K9 n
# c& v S# t/ w$ a例子:
4 n2 O5 U% m: A5 Z" v/ {; w9 F( y7 u- [! r+ u
matlab2009a代码:- clear all
3 r6 G, n/ T7 `0 u7 d( q' z7 ] - clc# _+ C) A# d5 K
- tic6 w; i1 [# m, f: b
- k = zeros(5,5); % //生成5×5全0矩阵\" w) u& \: H/ w\" d. ^- \. _
- % 循环计算以下程序段100000次:
: a4 ]. j& v4 o6 b7 N& e - for m = 1:100000; Q' d5 K\" b. |\" r+ f) f. ?
- a = rand(5,7);2 c/ _\" ~- S* m; X
- b = rand(7,5);%//生成5×7矩阵a,7×5矩阵b,用0~1之间的随机数初始化: ?1 J, z O. A. S+ Q! z7 Q* C- w, S
- k = k + a * b + a(1:5, 2:6) * b(2:6, 1:5) - a(:, 7) * b(3, :);
( W. e) S+ y8 t3 ?. q: y) q/ z - end' G s+ _% u! f: v2 C6 K
- k
% `0 V7 Q$ M; x2 J8 o' K. }% P - toc
复制代码 多次运行,结果约为:- k =
* f6 j8 x1 `) {1 U. [
. @6 J. }4 x* K1 R* ?5 a& b- 1.0e+005 *, Q2 r7 m) D5 D5 [1 v5 D9 K
( N& O3 b0 h( y- 2.7530 2.7525 2.7509 2.7539 2.7529
0 {2 Y- F5 x( b/ u/ k1 C5 [ - 2.7521 2.7494 2.7493 2.7542 2.7476
( u+ c$ m2 V9 T\" B( D, o& x) |8 q - 2.7549 2.7524 2.7527 2.7564 2.7517
( W$ u+ q( d4 j+ w3 ] - 2.7530 2.7523 2.7533 2.7525 2.7502
3 Z9 g$ d' ?8 q: _ - 2.7527 2.7521 2.7505 2.7519 2.75238 c, j2 d3 i1 f- R* }
- 5 y5 X: g& o* Q+ q P
- Elapsed time is 2.126704 seconds.
复制代码 Lu代码: - !!!using["math"];
- 0 q3 w0 x9 z) P( L v O
- main(:t0,k,i,a,b)=
- / T) r, A; K5 Q! N
- {9 e; I, Q$ J1 h8 j* `4 G% P ]
- t0=clock(),( [, m+ f6 B2 z1 Q4 Z
- k=zeros(5,5), //生成5×5矩阵k,初始化为0
- 4 R2 |! c0 W1 `6 b\\" A1 Q! P
- i=0,(++i<=100000).while{//循环计算100000次8 N X0 s, w; s8 a# g\\" k7 t0 f
- a=rand(5,7), b=rand(7,5), //生成5×7矩阵a,7×5矩阵b,用0~1之间的随机数初始化
- 7 E- J- o3 D3 e! ~/ _
- 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)
- 1 O' W5 e\\" D8 Y* j* h
- },
- 2 P9 e+ ?: Y. P2 z
- k.outa(), //输出矩阵k
- . x\\" W; [4 I: P5 `6 q
- [clock()-t0]/1000. //得到计算时间,秒
- ' R7 H! h+ V$ k& R+ U1 x4 C; R
- };
多次运行,结果约为:- 275024 274895 275054 275083 275024
# c6 A) C3 t\" G& |8 W9 D S\" K - 274965 275036 275091 275112 274931& K- A5 P |8 F5 ~$ m: x4 R
- 274975 275012 275020 274963 274939
! ~, I2 o9 q, t\" Y# v/ x - 274981 274976 275117 275065 275044
3 T/ r& }$ R% R# J; k. S, K; c - 274920 274906 275005 274972 274920
( U* V _& D\" L1 ~/ U7 X: J - 1.25 秒
复制代码 Lu快于matlab的原因应该仍在于内存管理效率较高。 |
zan
|