/[escript]/trunk/paso/src/ReactiveSolver.cpp
ViewVC logotype

Diff of /trunk/paso/src/ReactiveSolver.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 5774 by jfenwick, Fri Apr 24 01:36:26 2015 UTC revision 5775 by sshaw, Thu Jul 30 08:01:06 2015 UTC
# Line 57  err_t ReactiveSolver::solve(double* u, d Line 57  err_t ReactiveSolver::solve(double* u, d
57                  const double F_i = source[i];                  const double F_i = source[i];
58                  const double e_i = exp(x_i);                  const double e_i = exp(x_i);
59                  double u_i = e_i*u_old[i];                  double u_i = e_i*u_old[i];
60                  if (abs(x_i) > EXP_LIM_MIN) {                  if (std::abs(x_i) > EXP_LIM_MIN) {
61                      u_i += F_i/d_ii*(e_i-1.);                      u_i += F_i/d_ii*(e_i-1.);
62                  } else {                  } else {
63                      // second order approximation of (exp(x_i)-1)/x_i                      // second order approximation of (exp(x_i)-1)/x_i

Legend:
Removed from v.5774  
changed lines
  Added in v.5775

  ViewVC Help
Powered by ViewVC 1.1.26