Changeset 17740
- Timestamp:
- 04/15/14 16:20:35 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/shared/Numerics/NewtonSolveDnorm.cpp
r17650 r17740 32 32 33 33 /*Newton step*/ 34 y2 = y1 - (c1*pow(pow(10.,y1),s-1.) + c2*pow(10.,y1) - c3)/((s-1)*c1*log(10 )*pow(pow(10.,y1),s-1.) + c2*log(10.)*pow(10.,y1));34 y2 = y1 - (c1*pow(pow(10.,y1),s-1.) + c2*pow(10.,y1) - c3)/((s-1)*c1*log(10.)*pow(pow(10.,y1),s-1.) + c2*log(10.)*pow(10.,y1)); 35 35 36 36 if( fabs(y2-y1)/(fabs(y2))<threshold ){
Note:
See TracChangeset
for help on using the changeset viewer.