always check quantity against eps_rsrc() instead of zero

src/material.cc

@@ -141,7 +141,7 @@ void Material::Transmute(Composition::Ptr c) {
 }
 
 Resource::Ptr Material::PackageExtract(double qty, std::string new_package_name) {
-  if (qty > qty_) {
+  if ((qty - qty_) > eps_rsrc()) {
     throw ValueError("Attempted to extract more quantity than exists.");
   }
 
@@ -155,7 +155,7 @@ Resource::Ptr Material::PackageExtract(double qty, std::string new_package_name)
   // this call to res_tracker must come first before the parent resource 
   // state id gets modified
   other->tracker_.Package(&tracker_);
-  if (qty_ > cyclus::eps_rsrc()) {
+  if (qty_ > eps_rsrc()) {
     tracker_.Modify();
   }
   return boost::static_pointer_cast<Resource>(other);

src/package.cc

@@ -1,5 +1,6 @@
 #include "package.h"
 #include "error.h"
+#include "cyc_limits.h"
 
 namespace cyclus {
 
@@ -30,7 +31,7 @@ Package::Ptr& Package::unpackaged() {
 }
 
 std::pair<double, int> Package::GetFillMass(double qty) {
-  if (qty < fill_min_) {
+  if ((qty - fill_min_) < -eps_rsrc()) {
     // less than one pkg of material available
     return std::pair<double, int> (0, 0);
   }

src/resource.h

@@ -6,6 +6,7 @@
 #include <boost/shared_ptr.hpp>
 
 #include "package.h"
+#include "cyc_limits.h"
 
 class SimInitTest;
 
@@ -151,7 +152,7 @@ std::vector<typename T::Ptr> Resource::Package(Package::Ptr pkg) {
   int approx_num_pkgs = fill.second;
   Package::ExceedsSplitLimits(approx_num_pkgs);
 
-  while (quantity() > 0 && quantity() >= pkg->fill_min()) {
+  while (quantity() > 0 && (quantity() - pkg->fill_min()) >= -eps_rsrc()) {
     double pkg_fill = std::min(quantity(), fill_mass);
     t_pkgd = boost::dynamic_pointer_cast<T>(PackageExtract(pkg_fill, pkg->name()));
     ts_pkgd.push_back(t_pkgd); 

src/toolkit/matl_buy_policy.cc

@@ -318,7 +318,7 @@ void MatlBuyPolicy::AcceptMatlTrades(
   // check if cumulative cap has been reached. If yes, then sample for dormant
   // length and reset cycle_total_inv
   if (use_cumulative_capacity() && (
-    (cumulative_cap_ - cycle_total_inv_) < eps())) {
+    (cumulative_cap_ - cycle_total_inv_) < eps_rsrc())) {
       SetNextDormantTime();
       LGH(INFO3) << "cycle cumulative inventory has been reached. Dormant period will end at " << next_dormant_end_ << std::endl;
       cycle_total_inv_ = 0;

src/toolkit/matl_sell_policy.cc

@@ -265,7 +265,7 @@ void MatlSellPolicy::GetMatlTrades(
     if (shippable_pkgs > 0){
       qty = it->amt;
       LGH(INFO3) << " sending " << qty << " kg of " << it->request->commodity();
-      Material::Ptr mat = buf_->Pop(qty, cyclus::eps_rsrc());
+      Material::Ptr mat = buf_->Pop(qty, eps_rsrc());
       Material::Ptr trade_mat;
 
       // don't go through packaging if you don't need to. packaging always bumps
@@ -274,7 +274,7 @@ void MatlSellPolicy::GetMatlTrades(
       if (package_->name() != mat->package_name()) { // packaging needed
         std::vector<Material::Ptr> mat_pkgd = mat->Package<Material>(package_);
 
-        if (mat->quantity() > eps()) {
+        if (mat->quantity() > eps_rsrc()) {
           // push any extra material that couldn't be packaged back onto buffer
           // don't push unless there's leftover material
           buf_->Push(mat);