enrichment.h

Go to the documentation of this file.

#ifndef PYNE_B3ANNCKDQ5HEJLI33RPZPDNX6A
#define PYNE_B3ANNCKDQ5HEJLI33RPZPDNX6A

#ifndef PYNE_IS_AMALGAMATED
#include "enrichment_symbolic.h"
#endif

/************************************************/
/*** Enrichment Component Class and Functions ***/
/************************************************/

namespace pyne {
namespace enrichment {

  Cascade _fill_default_uranium_cascade();
  extern Cascade default_uranium_cascade;

  double feed_per_prod(double x_feed, double x_prod, double x_tail);
  double feed_per_tail(double x_feed, double x_prod, double x_tail);
  double prod_per_feed(double x_feed, double x_prod, double x_tail);
  double prod_per_tail(double x_feed, double x_prod, double x_tail);
  double tail_per_feed(double x_feed, double x_prod, double x_tail);
  double tail_per_prod(double x_feed, double x_prod, double x_tail);
  double value_func(double x);
  double swu_per_feed(double x_feed, double x_prod, double x_tail);
  double swu_per_prod(double x_feed, double x_prod, double x_tail);
  double swu_per_tail(double x_feed, double x_prod, double x_tail);

  double alphastar_i(double alpha, double Mstar, double M_i);


  Cascade solve_numeric(Cascade & orig_casc, double tolerance=1.0E-7,
                                             int max_iter=100);
  void _recompute_nm(Cascade & casc, double tolerance=1.0E-7);
  void _recompute_prod_tail_mats(Cascade & casc);
  Cascade _norm_comp_secant(Cascade & casc, double tolerance=1.0E-7, int max_iter=100);
  double _deltaU_i_OverG(Cascade & casc, int i);

  Cascade multicomponent(Cascade & orig_casc, char * solver,
                         double tolerance=1.0E-7, int max_iter=100);
  Cascade multicomponent(Cascade & orig_casc, std::string solver="symbolic",
                         double tolerance=1.0E-7, int max_iter=100);

  class EnrichmentInfiniteLoopError: public std::exception
  {
    virtual const char* what() const throw()
    {
      return "Inifinite loop found while calculating enrichment cascade.";
    }
  };

  class EnrichmentIterationLimit: public std::exception
  {
    virtual const char* what() const throw()
    {
      return "Iteration limit hit durring enrichment calculation.";
    }
  };

  class EnrichmentIterationNaN: public std::exception
  {
    virtual const char* what() const throw()
    {
      return "Iteration has hit a point where some values are not-a-number.";
    }
  };

// end enrichment
}
// end pyne
}

#endif