Merge pull request #17 from HaibingShao/react_trans_GIA

change the minimization solve method applied in OPS module. This will sl...
including activity models( Debye-Hückel and  Davies equations) for GIA
including activity models for GIA, part2
implemented starting value search to avoid negative concentrations. part1
using tilded xis to calculate concentrations. starting value search part2.
call update_lnK function
get the right set of lnK for each node and pass it to local problem
use the updated lnK in local problem
calculate AI vector only once at the beginning of local problem
clean up local problem
improve efficiency by moving AI vector calculation to the beginning of local problem.
fix equation 3.62.
passing vec_AI
changing quitting criteria in local problem to get a better convergence
fixing the line search algorithm
use solve_minimization function in vprime
use constant molecular diffusion coefficient (1.0E-9)
correcting RHS in global residual
use JH back transformation to concentration in global jacobian and  use the updated lnK values
change the quitting criteria for line search