mole_solve.cpp File Reference

#include "cddefines.h"
#include "deuterium.h"
#include "ionbal.h"
#include "phycon.h"
#include "hmi.h"
#include "dynamics.h"
#include "conv.h"
#include "trace.h"
#include "grainvar.h"
#include "newton_step.h"
#include "h2.h"
#include "mole_priv.h"
#include "mole.h"
#include "dense.h"

Include dependency graph for mole_solve.cpp:

Macros
#define	SMALLABUND   1e-24
#define	ABSLIM   1e-12
#define	ERRLIM   1e-12
#define	MAT(a, I_, J_)

Enumerations
enum	{ PRINTSOL = false }

Functions
void	check_co_ion_converge (void)
STATIC void	funjac (GroupMap &MoleMap, const valarray< double > &b2vec, double *const ervals, double *const amat, const bool lgJac, bool *lgConserve)
STATIC void	mole_h_fixup (void)
STATIC void	grouped_elems (const double bvec[], double mole_elems[])
double	mole_solve ()
STATIC void	mole_eval_dynamic_balance (long int num_total, double *b, bool lgJac, multi_arr< double, 2 > &c)

Macro Definition Documentation

ABSLIM

#define ABSLIM   1e-12

ERRLIM

#define ERRLIM   1e-12

MAT

#define MAT	(		a,
			I_,
			J_ )

Value:

((a)[(I_)*(mole_global.num_compacted)+(J_)])

Referenced by funjac().

SMALLABUND

#define SMALLABUND   1e-24

Referenced by newton_step(), and GroupMap::setup().

Enumeration Type Documentation

anonymous enum

anonymous enum

Enumerator
PRINTSOL

Function Documentation

check_co_ion_converge()

void check_co_ion_converge

(

void

)

References conv, DEBUG_ENTRY, dense, findspecieslocal(), ipCARBON, ipOXYGEN, and SDIV().

Referenced by mole_solve().

Here is the call graph for this function:

funjac()

STATIC void funjac	(	GroupMap &	MoleMap,
		const valarray< double > &	b2vec,
		double *const	ervals,
		double *const	amat,
		const bool	lgJac,
		bool *	lgConserve )

References ASSERT, DEBUG_ENTRY, dynamics, GroupMap::fion, fnzone, fprintf(), get_ptr(), grouped_elems(), groupspecies, ioQQQ, iteration, MAT, mole, mole_eval_dynamic_balance(), mole_global, GroupMap::molElems, nuclide_list, PRINTSOL, SDIV(), STATIC, trace, and GroupMap::updateMolecules().

Referenced by mole_solve().

Here is the call graph for this function:

grouped_elems()

STATIC void grouped_elems	(	const double	bvec[],
		double	mole_elems[] )

References groupspecies, mole_global, nuclide_list, and STATIC.

Referenced by funjac(), and GroupMap::setup().

mole_eval_dynamic_balance()

STATIC void mole_eval_dynamic_balance	(	long int	num_total,
		double *	b,
		bool	lgJac,
		multi_arr< double, 2 > &	c )

References ASSERT, DEBUG_ENTRY, dense, dynamics, iteration, LIMELM, mole, mole_eval_balance(), mole_global, and STATIC.

Referenced by funjac().

Here is the call graph for this function:

mole_h_fixup()

STATIC void mole_h_fixup

(

void

)

References DEBUG_ENTRY, molezone::den, dense, findspecieslocal(), gv, h2, hmi, ipHYDROGEN, mole, mole_global, SMALLFLOAT, and STATIC.

Referenced by mole_solve().

Here is the call graph for this function:

mole_solve()

double mole_solve

(

void

)

mole_solve fills in matrix for heavy elements molecular routines

References ASSERT, BIGFLOAT, check_co_ion_converge(), conv, DEBUG_ENTRY, dense, fixit, fnzone, fprintf(), funjac(), get_ptr(), groupspecies, hmi, ionbal, ioQQQ, ipHYDROGEN, lgElemsConserved(), mole, mole_global, mole_h_fixup(), mole_return_cached_species(), newton_step(), nuclide_list, nzone, phycon, GroupMap::setup(), trace, and GroupMap::updateMolecules().

Referenced by mole_drive().

Here is the call graph for this function: