source: issm/oecreview/Archive/16554-17801/ISSM-17739-17740.diff@ 17802

Last change on this file since 17802 was 17802, checked in by Mathieu Morlighem, 11 years ago

Added archives
File size: 656 bytes

  • ../trunk-jpl/src/c/shared/Numerics/NewtonSolveDnorm.cpp

     
    3131        while(true){
    3232
    3333                /*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));
    3535
    3636                if( fabs(y2-y1)/(fabs(y2))<threshold ){
    3737                        break;
Note: See TracBrowser for help on using the repository browser.