Rev | Line | |
---|
[17802] | 1 | Index: ../trunk-jpl/src/c/shared/Numerics/NewtonSolveDnorm.cpp
|
---|
| 2 | ===================================================================
|
---|
| 3 | --- ../trunk-jpl/src/c/shared/Numerics/NewtonSolveDnorm.cpp (revision 17739)
|
---|
| 4 | +++ ../trunk-jpl/src/c/shared/Numerics/NewtonSolveDnorm.cpp (revision 17740)
|
---|
| 5 | @@ -31,7 +31,7 @@
|
---|
| 6 | while(true){
|
---|
| 7 |
|
---|
| 8 | /*Newton step*/
|
---|
| 9 | - 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));
|
---|
| 10 | + 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));
|
---|
| 11 |
|
---|
| 12 | if( fabs(y2-y1)/(fabs(y2))<threshold ){
|
---|
| 13 | break;
|
---|
Note:
See
TracBrowser
for help on using the repository browser.