Commit graph

3043 commits

Author SHA1 Message Date
Harry Moffat
14bf2b1ee2 Started putting these files into the standard form.
Started eliminating doxygen errors.
2011-01-10 23:22:19 +00:00
Harry Moffat
f727cd44af Added in a more robust treatment of handling limits 2011-01-05 03:55:13 +00:00
Harry Moffat
a36864543a Got rid of a print statement. 2011-01-04 23:02:07 +00:00
Harry Moffat
b1750d579a Added return variables to a lot of Residual evaluation functions.
These are now used by the nonlinear solver to bail if the residual
says to bail.
2011-01-04 22:04:59 +00:00
Harry Moffat
5395e214b8 Fixed a doxygen error. 2010-12-31 23:54:29 +00:00
Harry Moffat
eb35f9e253 Some changes to reduce extra printing
Some changes due to a bad algorithm - Known - need to fix this.
2010-12-31 23:20:55 +00:00
Harry Moffat
8c4d47f59c Fixed an error in the phase stability calculation 2010-12-18 18:26:16 +00:00
Harry Moffat
fc7b7e8dc6 Initial attempt at exposing the phase stability routine to the
outside without actually doing an equilibrium solution.
2010-12-16 22:44:29 +00:00
Harry Moffat
32816b1fa7 Incremental commit of new ability to determine phase stability
without actually solving the equilibrium equations.
2010-12-15 20:20:13 +00:00
Harry Moffat
e10fa97af0 Added the phaseIsStable idea 2010-12-15 02:54:08 +00:00
Harry Moffat
73d8df7f4c Tightend up the tolerances on the nonlinear solvers. 2010-12-15 02:15:27 +00:00
Harry Moffat
ec8bbc5be2 Fixed the help section to be in line with what the code actually does. 2010-12-14 23:56:12 +00:00
Harry Moffat
d5f44a5e6b ResidJacEval:
Added in support for returning a failure condition from calcDeltaSolnVariables()


NonLinearSolver:
   Added in support for returning a failure condition from bueler_jac().

   Added in a call to return an editable bounds constraint vector.

DenseMatrix:
   Added in comments about error handling from LAPACK routines.


mdp_allo:
   Added in a new routine, mdp_zero_int_1

thermo.h
electrolyteThermo.h
   Started to add in ifdef blocks.  These are not handled correctly currently.
2010-12-07 18:05:25 +00:00
Harry Moffat
44e03fcf79 Basic upgrade to the DenseMatrix Class.
We've added a capability to change the handling of singular matrix
    conditions.
     Options exist to either throw an exception or return an error code.
     Options exist to print out a warning or to not print out anything.
2010-11-30 15:33:50 +00:00
Harry Moffat
2e50b8a929 Fixed an error. The m_speciesData XML tree wasn't being created.
I created one in the constructor.
2010-11-18 03:09:45 +00:00
Harry Moffat
1f53f7db72 Fixed minimal python build. 2010-11-15 22:34:21 +00:00
Harry Moffat
3e0d9cef21 Fixed an error that crept in from the recent merge 2010-11-12 22:06:23 +00:00
Harry Moffat
ff3ebb25c2 Merge of the main branch into the LiquidTransportDevelop branch. 2010-11-12 21:37:41 +00:00
Harry Moffat
c3405dad33 Added a file 2010-11-12 20:06:53 +00:00
Harry Moffat
b6e1b14216 Deleted a file 2010-11-12 20:00:56 +00:00
Harry Moffat
fa56309462 Modifications to the Makefiles to make the THREAD_SAFE_CANTERA ifdef
block work again.
2010-11-12 19:18:49 +00:00
Harry Moffat
4eb6e622be Fixed a compiler error on the THREAD_SAFE ifdef 2010-11-12 18:08:51 +00:00
Harry Moffat
56ad9ca9e8 Fixed a compile error in the THREAD_SAFE capability 2010-11-12 18:02:21 +00:00
Harry Moffat
58baed0f2e Fixed issue 17
There was an obvious error in the THREAD_SAFE_CANTERA compile.
2010-11-11 00:19:59 +00:00
Harry Moffat
4942cf30aa Fixed issue #17.
There was an obvious error in the THREAD_SAFE_CANTERA ifdef block compile.
2010-11-11 00:17:45 +00:00
Harry Moffat
54e1280e08 Added hooks in for filters for New step directions and for
new solutions.
2010-11-10 05:01:22 +00:00
Harry Moffat
813a98ddf1 Fixed the residual norm calculation. What was in there before
was necessary but not sufficient. Since there is no good 
theory in tying the residual calculation to convergence short
of calculating the condition number, I've made its magnitude 
a slave to the solution update magnitude.
2010-11-10 01:39:48 +00:00
Harry Moffat
12b6f572e6 Added in a new Residual evaluation type, Base_LaggedSolutionComponents,
to properly handle Picard iterations within the damping algorithm.
2010-11-08 22:22:13 +00:00
Harry Moffat
154b234da9 Added a printout of the jacobian, when an ifdef block is turned on. 2010-11-05 19:44:05 +00:00
Harry Moffat
c44b477285 Fixed an error in the specification of m_deltaBoundsMagnitudes
Made the residual check in dampStep a lot harder to satisfy.
2010-11-02 15:44:09 +00:00
Harry Moffat
7fc2841f85 Fixed an error in vcs_solve_TP. m_tPhaseMoles_new wasn't set
before being used.
2010-10-28 01:33:54 +00:00
Harry Moffat
c7ce446b32 Added a special constructor for setting the element potential directly. 2010-10-27 21:07:51 +00:00
Harry Moffat
9373fe02f7 These routines needed std:: added to the arguments and function bodies
of std c++ library calls.
2010-10-27 21:07:15 +00:00
Harry Moffat
a96f865403 Added FixedChemPotSSTP object.
This is a utility object that will provide a bath species at a certain
chemical potential to the equilibrium solver.
2010-10-27 17:00:46 +00:00
Harry Moffat
13f151a358 Modification to the structure of solve_TP()
Took the special case of SS species deletion in RxnStep and moved
it into the main routine. It now counts as an interation step.
I'm trying to simplify the structure of the program.
2010-10-25 22:12:55 +00:00
Harry Moffat
1bdee90289 Fixed an error that occured when components of multispecies
phases were  zeroed because other phases were being zeroed.
Then the noncomponents of that same multispecies phase could 
comprise the entire phase and the Hessian diagonal would then be
zero. The fix to zero the species seems to work.
2010-10-23 01:35:10 +00:00
Harry Moffat
f26ff94525 Changed NaCl_equil printouts due to rebaseline.
Fixed Makefile error
2010-10-22 22:30:24 +00:00
Harry Moffat
7632281197 Changed table at the start of debug printing 2010-10-22 20:18:10 +00:00
Harry Moffat
e23a780629 Added missing functions from PureFluidPhase 2010-10-20 18:18:05 +00:00
Harry Moffat
550c33c628 Added a getDeltaElectrochemPotentials() member function 2010-10-18 19:27:06 +00:00
Harry Moffat
78b2d5689c Took out an unneeded m_atol call that wasn't necessary and caused
problems with numerical deltas.
2010-10-15 21:02:19 +00:00
Harry Moffat
43c4d1215c Fixed a printLvl error 2010-10-15 18:33:57 +00:00
Harry Moffat
fad1f9f7ff Added hooks for rtol and max newton iterations 2010-10-14 22:34:59 +00:00
Harry Moffat
97f326a628 Tightened the convergence tolerances 2010-10-14 02:44:23 +00:00
Christopher Lueth
9630216f3e Fixed nan problem if m_normResid0 = 0.0 at line 965 2010-10-13 22:06:03 +00:00
Harry Moffat
ae989985eb Added a hook for specifying atol 2010-10-13 16:58:34 +00:00
Harry Moffat
1d0ba2d3b8 Changed the interface for ResidJacEval.
Updated NonlinearSolver with a lot of heuristic algorithm changes.
2010-10-08 17:16:25 +00:00
Harry Moffat
7905ea977f Rebaselined some test problems. 2010-10-04 23:00:26 +00:00
Harry Moffat
452e8a4957 Added DenseMatrix 2010-10-04 22:57:34 +00:00
Harry Moffat
367d1f90b0 Small changes that I found improved things 2010-10-02 04:51:00 +00:00