1.1 --- a/tools/decompose_2.cpp	Fri Jul 01 11:04:42 2011 +0200
     1.2 +++ b/tools/decompose_2.cpp	Fri Jul 01 11:25:13 2011 +0200
     1.3 @@ -54,7 +54,17 @@
     1.4  		 std::back_inserter(masses),
     1.5  		 OT::null_output_iterator(),
     1.6  		 0.001, 30);
     1.7 -  decompose_2(positions, masses, decomposition, k, L, threshold);
     1.8 +#if 1
     1.9 +  std::ofstream os("/tmp/test.diracs");
    1.10 +  for (size_t i = 0; i < positions.size(); ++i)
    1.11 +    {
    1.12 +      os << positions[i].x() << " "
    1.13 +	 << positions[i].y() << " "
    1.14 +	 << masses[i] << "\n";
    1.15 +    }
    1.16 +#endif
    1.17 +
    1.18 +  decompose_2(positions, masses, decomposition, L, k, threshold);
    1.19  }
    1.20  
    1.21  bool
    1.22 @@ -131,7 +141,7 @@
    1.23  	  std::vector<Point> positions;
    1.24  	  std::vector<double> masses;
    1.25  	  load_discrete_measure_2 (positions, masses, source);
    1.26 -	  decompose_2(positions, masses, decomposition, N, L, k);
    1.27 +	  decompose_2(positions, masses, decomposition, L, k);
    1.28  	}
    1.29        else
    1.30  	{