Go to the source code of this file.
|
typedef void(* | error_print_t )(long, long, const valarray< double > &, const valarray< double > &) |
|
|
bool | newton_step (GroupMap &mole_map, const valarray< double > &oldmols, valarray< double > &newmols, realnum *eqerror, realnum *error, const long n, double *rlimit, double *rmax, valarray< double > &escale, void(*jacobn)(GroupMap &mole_map, const valarray< double > &b2vec, double *const ervals, double *const amat, const bool lgJac, bool *lgConserved)) |
|
int32 | solve_system (const valarray< double > &a, valarray< double > &b, long int n, error_print_t error_print) |
|
typedef void(* error_print_t)(long, long, const valarray< double > &, const valarray< double > &) |
bool newton_step |
( |
GroupMap & |
mole_map, |
|
|
const valarray< double > & |
oldmols, |
|
|
valarray< double > & |
newmols, |
|
|
realnum * |
eqerror, |
|
|
realnum * |
error, |
|
|
const long |
n, |
|
|
double * |
rlimit, |
|
|
double * |
rmax, |
|
|
valarray< double > & |
escale, |
|
|
void(*)(GroupMap &mole_map, const valarray< double > &b2vec, double *const ervals, double *const amat, const bool lgJac, bool *lgConserved) |
jacobn |
|
) |
| |
Take one Newton step
- Parameters
-
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()().