diff --git a/fastlib/trunk/contrib/dongryel/thesis_research/core/gnp/dualtree_dfs.h b/fastlib/trunk/contrib/dongryel/thesis_research/core/gnp/dualtree_dfs.h
index 0f3049e1b2..4a884ccf22 100644
--- a/fastlib/trunk/contrib/dongryel/thesis_research/core/gnp/dualtree_dfs.h
+++ b/fastlib/trunk/contrib/dongryel/thesis_research/core/gnp/dualtree_dfs.h
@@ -252,6 +252,7 @@ class DualtreeDfs {
       TreeType *qnode,
       TreeType *rnode,
       const typename ProblemType::DeltaType &delta,
+      const core::math::Range &squared_distance_range,
       typename ProblemType::ResultType *query_results);
 
     /** @brief Summarize a given pair of query/reference using a
diff --git a/fastlib/trunk/contrib/dongryel/thesis_research/core/gnp/dualtree_dfs_dev.h b/fastlib/trunk/contrib/dongryel/thesis_research/core/gnp/dualtree_dfs_dev.h
index 605085a6b5..bb6793d833 100644
--- a/fastlib/trunk/contrib/dongryel/thesis_research/core/gnp/dualtree_dfs_dev.h
+++ b/fastlib/trunk/contrib/dongryel/thesis_research/core/gnp/dualtree_dfs_dev.h
@@ -276,6 +276,7 @@ bool DualtreeDfs<ProblemType>::CanSummarize_(
   typename ProblemType::TableType::TreeType *qnode,
   typename ProblemType::TableType::TreeType *rnode,
   const typename ProblemType::DeltaType &delta,
+  const core::math::Range &squared_distance_range,
   typename ProblemType::ResultType *query_results) {
 
   typename ProblemType::StatisticType &qnode_stat = qnode->stat();
@@ -284,7 +285,8 @@ bool DualtreeDfs<ProblemType>::CanSummarize_(
   new_summary.ApplyDelta(delta);
 
   return new_summary.CanSummarize(
-           problem_->global(), delta, qnode, rnode, query_results);
+           problem_->global(), delta, squared_distance_range,
+           qnode, rnode, query_results);
 }
 
 template<typename ProblemType>
@@ -348,7 +350,8 @@ bool DualtreeDfs<ProblemType>::DualtreeCanonical_(
                              squared_distance_range);
 
   // If it is prunable, then summarize and return.
-  if(CanSummarize_(qnode, rnode, delta, query_results)) {
+  if(CanSummarize_(
+        qnode, rnode, delta, squared_distance_range, query_results)) {
     Summarize_(qnode, delta, query_results);
     num_deterministic_prunes_++;
     return true;
diff --git a/fastlib/trunk/contrib/dongryel/thesis_research/core/gnp/dualtree_dfs_iterator_dev.h b/fastlib/trunk/contrib/dongryel/thesis_research/core/gnp/dualtree_dfs_iterator_dev.h
index 464a549388..bda3f0c6fd 100644
--- a/fastlib/trunk/contrib/dongryel/thesis_research/core/gnp/dualtree_dfs_iterator_dev.h
+++ b/fastlib/trunk/contrib/dongryel/thesis_research/core/gnp/dualtree_dfs_iterator_dev.h
@@ -153,7 +153,8 @@ void DualtreeDfs<ProblemType>::iterator<IteratorMetricType>::operator++() {
       metric_, engine_->problem_->global(),
       qnode, rnode, squared_distance_range);
     bool prunable = engine_->CanSummarize_(
-                      qnode, rnode, delta, query_results_);
+                      qnode, rnode, delta,
+                      squared_distance_range, query_results_);
 
     if(prunable) {
       engine_->Summarize_(qnode, delta, query_results_);
diff --git a/fastlib/trunk/contrib/dongryel/thesis_research/mlpack/kde/kde_dualtree.h b/fastlib/trunk/contrib/dongryel/thesis_research/mlpack/kde/kde_dualtree.h
index f6a9584465..94c9a0ed0b 100644
--- a/fastlib/trunk/contrib/dongryel/thesis_research/mlpack/kde/kde_dualtree.h
+++ b/fastlib/trunk/contrib/dongryel/thesis_research/mlpack/kde/kde_dualtree.h
@@ -623,6 +623,60 @@ class KdeSummary {
     // For Boost serialization.
     friend class boost::serialization::access;
 
+    template <
+    typename GlobalType, typename DeltaType, typename TreeType,
+             typename ResultType >
+    bool CanSummarizeSeriesExpansion_(
+      const GlobalType &global, const DeltaType &delta,
+      const core::math::Range &squared_distance_range,
+      TreeType *qnode, TreeType *rnode, double right_hand_side,
+      ResultType *query_results, int *order_farfield_to_local,
+      int *order_farfield, int *order_local) const {
+
+      // The far-field expansion of the reference node.
+      const typename GlobalType::KernelAuxType::FarFieldType &
+      farfield_expansion = rnode->stat().farfield_expansion_;
+
+      // The local expansion of the query node.
+      const typename GlobalType::KernelAuxType::LocalType &
+      local_expansion = qnode->stat().local_expansion_;
+
+      // Actual amount of error incurred per each query/ref pair.
+      double actual_err_farfield_to_local = 0;
+      double actual_err_farfield = 0;
+      double actual_err_local = 0;
+
+      // The allowed error per each query/ref pair.
+      double allowed_err =
+        right_hand_side / static_cast<double>(rnode->count());
+
+      // Estimated computational cost.
+      int cost_farfield_to_local = std::numeric_limits<int>::max();
+      int cost_farfield = std::numeric_limits<int>::max();
+      int cost_local = std::numeric_limits<int>::max();
+      int cost_exhaustive = (qnode->count()) * (rnode->count()) *
+                            (qnode->bound().dim());
+      int min_cost = 0;
+
+      // Get the order of approximations.
+      *order_farfield_to_local =
+        global.kernel_aux().OrderForConvertingFromFarFieldToLocal(
+          rnode->bound(), qnode->bound(),
+          squared_distance_range.lo, squared_distance_range.hi,
+          allowed_err, &actual_err_farfield_to_local);
+      *order_farfield =
+        global.kernel_aux().OrderForEvaluatingFarField(
+          rnode->bound(), qnode->bound(),
+          squared_distance_range.lo, squared_distance_range.hi,
+          allowed_err, &actual_err_farfield);
+      *order_local =
+        global.kernel_aux().OrderForEvaluatingLocal(
+          rnode->bound(), qnode->bound(),
+          squared_distance_range.lo, squared_distance_range.hi,
+          allowed_err, &actual_err_local);
+      return false;
+    }
+
   public:
 
     double densities_l_;
@@ -776,6 +830,7 @@ class KdeSummary {
              typename ResultType >
     bool CanSummarize(
       const GlobalType &global, const DeltaType &delta,
+      const core::math::Range &squared_distance_range,
       TreeType *qnode, TreeType *rnode, ResultType *query_results) const {
 
       double left_hand_side = delta.used_error_;
@@ -795,15 +850,14 @@ class KdeSummary {
       // Otherwise, try series expansion.
       else {
 
-        // The far-field expansion of the reference node.
-        const typename GlobalType::KernelAuxType::FarFieldType &
-        farfield_expansion = rnode->stat().farfield_expansion_;
+        int order_farfield_to_local;
+        int order_farfield;
+        int order_local;
 
-        // The local expansion of the query node.
-        const typename GlobalType::KernelAuxType::LocalType &
-        local_expansion = qnode->stat().local_expansion_;
-
-        return false;
+        return CanSummarizeSeriesExpansion_(
+                 global, delta, squared_distance_range,
+                 qnode, rnode, right_hand_side, query_results,
+                 &order_farfield_to_local, &order_farfield, &order_local);
       }
       return false;
     }