7楼: Originally posted by chyanog at 2013-08-03 20:11:32
Compile一点也不鸡肋，并不是只能编译一些简单函数，一些数值计算程序可以提速30~100倍以上，但使用一些高级的函数时（NSolve，NIntegrate等）就没必要Compile了，不会有加速效果，
使用CompilationTarget -> & ...

12楼: Originally posted by walk1997 at 2013-08-07 11:32:38
test3 :=
  Position;

可以进一步提高下

13楼: Originally posted by zdzqy at 2013-08-07 11:47:40
这个代码复制过来貌似乱了，能否给我发个nb的文件？...

12楼: Originally posted by walk1997 at 2013-08-07 11:32:38
test3 :=
  Position;

可以进一步提高下

SparseArray[ 1 - Unitize@  FractionalPart@Sqrt@N@Array[#1^2 + #2^2 &, {1000, 1000}]][ "NonzeroPositions"] // Length // AbsoluteTiming

(*---------------------------*)
SparseArray[ 1 - Unitize@FractionalPart@Sqrt@N@(#1^2 + #2^2) & @@ {#, Transpose@#} &[ConstantArray[Range[1000], 1000]]][
   "NonzeroPositions"] // Length // AbsoluteTiming

15楼: Originally posted by chyanog at 2013-08-07 11:56:26
可以提速的，向量化的方法也很快，和Compile后的速度可以一比，但是一个缺点是比较耗内存


SparseArray // Length // AbsoluteTiming

(*---------------------------*)
SparseArray[ 1 - Unitize@Fraction ...

15楼: Originally posted by chyanog at 2013-08-07 11:56:26
可以提速的，向量化的方法也很快，和Compile后的速度可以一比，但是一个缺点是比较耗内存


SparseArray // Length // AbsoluteTiming

(*---------------------------*)
SparseArray[ 1 - Unitize@Fraction ...

Clear["`*"]
mandel[a_, it_: 100] :=
  Module[{z = 0. I, c = a[[2]] + a[[1]] I},
   Catch[Do[If[Abs[z] > 2, Throw@i];
     z = z^2 + c;, {i, it}]; it]
   ];

mandelc =  Compile[{{a, _Real, 1}, {it, _Integer}},
   Module[{z = 0. I, c = a[[2]] + a[[1]] I},
    Catch[Do[If[Abs[z] > 2, Throw@i];
      z = z^2 + c;, {i, it}]; it]
    ], RuntimeAttributes -> {Listable}, RuntimeOptions -> "Speed"];

d1 = mandelc[Outer[List, Range[-1., 1, 2./199], Range[-2., 0.5, 2.5/199]], 100]; // AbsoluteTiming
d2 = Table[ mandel[{i, j}], {i, -1., 1, 2./199}, {j, -2., 0.5,  2.5/199}]; // AbsoluteTiming

d1 == d2
Colorize[#, ColorFunction -> Hue] & /@ {d1, d2}

17楼: Originally posted by walk1997 at 2013-08-07 12:12:55
btw:
  想到一个问题 严格来说 这样求解出来可能会有错误
因为有可能a^2+b^2//N//Sqrt 可能非常非常接近整数
不知道在机器精度下 第一次出现错误的解会在哪里？
不知道存在不？...

10楼: Originally posted by walk1997 at 2013-08-06 23:58:27
学习了  非常好的例子
函数内部的运算量/次数很大的时候 加速比较明显
也和运算的种类有关 编译成C 和编译成 WMV 也差个几倍...