newton_step.cpp File Reference

#include "cddefines.h"
#include "newton_step.h"
#include "conv.h"
#include "thirdparty.h"
#include "mole.h"
#include "mole_priv.h"
#include "trace.h"
#include "save.h"

Include dependency graph for newton_step.cpp:

Macros
#define	MAT(a, I_, J_)   ((a)[(I_)*(n)+(J_)])

Enumerations
enum	{ PRINTSOL = false }

Functions
STATIC void	mole_system_error (long n, long merror, const valarray< double > &a, const valarray< double > &b)
bool	newton_step (GroupMap &MoleMap, const valarray< double > &b0vec, valarray< double > &b2vec, realnum *eqerror, realnum *error, const long n, double *rlimit, double *rmax, valarray< double > &escale, void(*jacobn)(GroupMap &MoleMap, const valarray< double > &b2vec, double *const ervals, double *const amat, const bool lgJac, bool *lgConserve))
int32	solve_system (const valarray< double > &a, valarray< double > &b, long int n, error_print_t error_print)

Macro Definition Documentation

#define MAT	(		a,
			I_,
			J_
	)		((a)[(I_)*(n)+(J_)])

Referenced by newton_step().

Enumeration Type Documentation

anonymous enum

Enumerator
PRINTSOL

Function Documentation

STATIC void mole_system_error	(	long	n,
		long	merror,
		const valarray< double > &	a,
		const valarray< double > &	b
	)

References fprintf(), groupspecies, ioQQQ, and mole_print_species_reactions().

Referenced by newton_step().

Here is the call graph for this function:

bool newton_step	(	GroupMap &	MoleMap,
		const valarray< double > &	b0vec,
		valarray< double > &	b2vec,
		realnum *	eqerror,
		realnum *	error,
		const long	n,
		double *	rlimit,
		double *	rmax,
		valarray< double > &	escale,
		void(*)(GroupMap &MoleMap, const valarray< double > &b2vec, double *const ervals, double *const amat, const bool lgJac, bool *lgConserve)	jacobn
	)

References conv, DEBUG_ENTRY, e1(), e2(), fprintf(), get_ptr(), groupspecies, ioQQQ, t_trace::lgTrace, MAT, MIN2, mole_system_error(), t_conv::register_(), SMALLABUND, solve_system(), and trace.

Referenced by mole_solve().

Here is the call graph for this function:

int32 solve_system	(	const valarray< double > &	a,
		valarray< double > &	b,
		long int	n,
		error_print_t	error_print
	)

References ASSERT, DEBUG_ENTRY, fprintf(), get_ptr(), getrf_wrapper(), getrs_wrapper(), ioQQQ, maxe, and NULL.

Referenced by find_solution(), newton_step(), and Atom_LevelN::operator()().

Here is the call graph for this function: