Continuation (延续)

首选是说下尾递归.如果一个函数呈现下列情况,且中间没有再次递归使用自己,可以认为是一个尾递归.

R function(a,b)

{

　　......

　　return function(c,d);
}

尾递归可以很轻松的改为一个循环的,如果编译器(解释器)支持的话,可以直接优化的.一般的FP都支持尾递归的.

如果像树,Fibonacc这种结构来使用递归的,好像很难转化为尾递归的,这时候使用到的另一种方法是Continuation ,Continuation 可以是看做是一个尾递归的,不过得把 b,d 看做函数它们有着相同的函数签名.

看个例子把

let rec fib n = 
    match n with
    | 0|1 -> n    
    |_ -> fib (n-1)+fib (n-2)

简单的fibonacc程序,如果这个改成continuation结构的话如下:

let fibonacci_cps n = 
  let rec fibonacci_cont a cont = 
    if a <= 2 then cont 1 
    else 
      fibonacci_cont (a - 2) (fun x -> 
        fibonacci_cont (a - 1) (fun y ->  
          cont(x + y)))            
  fibonacci_cont n (fun x -> x)

fibonacci_cont 函数中展开了结果函数调用.原来的程序是将函数放在了栈上,而这个函数相当于在堆上生成了很多的函数对象.其实程序的复杂度没有变化的.

再看下C++中的continuation,本来想使用c++0x来重写上面的例子的,但是发现不行主要是c++lambda支持"对lambda外部局部变量引用无效",二即使写出来最终的结果也是放在堆上的.

当然这个时候有个解决方案,使用传说中的tbb.

tbb中task是带continuation应用的.在讨论tbb前,我们先看下 return fib(n-1)+fib(n-2)的含义.在栈上.调用函数将n-1 和 n-2分别给了函数 fib fib,在这两个函数结束之前是不能够运行本函数的,好了这个时候我们就发现了我们需要在堆上模拟这种情况.一个函数对象,它能够调用2个函数,并且在这2个函数调用之前它自己不能运行.

好了直接上代码了,代码中有解释的

struct FibContinuation: public task {
   long* const sum;
   long x, y;
   FibContinuation( long* sum_ ) : sum(sum_) {}
   task* execute() {
   *sum = x+y;        //fib(n-1)+fib(n-2)
   return NULL;
   }
};
struct FibTask: public task {
   const long n;
   long* const sum;
   FibTask( long n_, long* sum_ ) : n(n_), sum(sum_) {}

   task* execute() {
      if( n<CutOff ) {
         *sum = SerialFib(n);
         return NULL;
      } else {
         FibContinuation& c =
         *new( allocate_continuation() ) FibContinuation(sum);       //生成函数对象
         FibTask& a = *new( c.allocate_child() ) FibTask(n-2,&c.x);  //fib(n-2)
         FibTask& b = *new( c.allocate_child() ) FibTask(n-1,&c.y);  //fib(n-1)
         // Set ref_count to "two children plus one for the wait".
         c.set_ref_count(2);                                         //设置要等两个函数
         c.spawn( b );                                               //不阻塞地运行fib(n-1)
         c.spawn( a );                                               //不阻塞的运行fib(n-2)
         return NULL;
      }
   }
};