diff --git a/fastlib/u/dongryel/series_expansion/mult_farfield_expansion.h b/fastlib/u/dongryel/series_expansion/mult_farfield_expansion.h
index 708b4a6974..83789be09e 100644
--- a/fastlib/u/dongryel/series_expansion/mult_farfield_expansion.h
+++ b/fastlib/u/dongryel/series_expansion/mult_farfield_expansion.h
@@ -108,7 +108,7 @@ class MultFarFieldExpansion {
    * Evaluates the far-field coefficients at the given point
    */
   double EvaluateField(Matrix* data=NULL, int row_num=-1,
-		       Vector* x_q=NULL) const;
+		       Vector* x_q=NULL, int order=-1) const;
   
   /**
    * Evaluates the two-way convolution mixed with exhaustive computations
@@ -293,13 +293,13 @@ void MultFarFieldExpansion<TKernel, TKernelAux>::RefineCoeffs
 
 template<typename TKernel, typename TKernelAux>
 double MultFarFieldExpansion<TKernel, TKernelAux>::
-  EvaluateField(Matrix* data, int row_num, Vector* x_q) const {
+  EvaluateField(Matrix* data, int row_num, Vector* x_q, int order) const {
   
   // dimension
   int dim = sea_->get_dimension();
 
   // total number of coefficients
-  int total_num_coeffs = sea_->get_total_num_coeffs(order_);
+  int total_num_coeffs = sea_->get_total_num_coeffs(order);
 
   // square root times bandwidth
   double bandwidth_factor = ka_.BandwidthFactor(kernel_.bandwidth_sq());