|
|
|
@ -298,10 +298,10 @@ static double CeilingPow |
|
|
|
double xmax2 = (xmax) * (xmax) ; \ |
|
|
|
double xp = 1.0 , xmp = 1.0 ; \ |
|
|
|
int m , mm ; \ |
|
|
|
double arg , dnm ; \ |
|
|
|
double arg_cp , dnm ; \ |
|
|
|
for ( m = 0 ; m < pw ; m ++ ) { xp *= x2 ; xmp *= xmax2 ; } \ |
|
|
|
arg = xp + xmp ; \ |
|
|
|
dnm = arg ; \ |
|
|
|
arg_cp = xp + xmp ; \ |
|
|
|
dnm = arg_cp ; \ |
|
|
|
if ( pw == 1 || pw == 2 || pw == 4 || pw == 8 ) { \ |
|
|
|
if ( pw == 1 ) { mm = 1 ; \ |
|
|
|
} else if ( pw == 2 ) { mm = 2 ; \ |
|
|
|
@ -311,7 +311,7 @@ static double CeilingPow |
|
|
|
} else { dnm = Fn_Pow( dnm , 1.0 / ( 2.0 * pw ) ) ; } \ |
|
|
|
dnm = 1.0 / dnm ; \ |
|
|
|
y = (x) * (xmax) * dnm ; \ |
|
|
|
dx = (xmax) * xmp * dnm / arg ; \ |
|
|
|
dx = (xmax) * xmp * dnm / arg_cp ; \ |
|
|
|
} |
|
|
|
|
|
|
|
#define Fn_SU_CP2( y , x , xmax , delta , pw , dx , dxmax ) { \ |
|
|
|
@ -334,10 +334,10 @@ static double CeilingPow |
|
|
|
double xmax2 = (xmax) * (xmax) ; \ |
|
|
|
double xp = 1.0 , xmp = 1.0 ; \ |
|
|
|
int m =0, mm =0; \ |
|
|
|
double arg =0.0, dnm =0.0; \ |
|
|
|
double arg_cp2 =0.0, dnm =0.0; \ |
|
|
|
for ( m = 0 ; m < pw ; m ++ ) { xp *= x2 ; xmp *= xmax2 ; } \ |
|
|
|
arg = xp + xmp ; \ |
|
|
|
dnm = arg ; \ |
|
|
|
arg_cp2 = xp + xmp ; \ |
|
|
|
dnm = arg_cp2 ; \ |
|
|
|
if ( pw == 1 || pw == 2 || pw == 4 || pw == 8 ) { \ |
|
|
|
if ( pw == 1 ) { mm = 1 ; \ |
|
|
|
} else if ( pw == 2 ) { mm = 2 ; \ |
|
|
|
@ -347,8 +347,8 @@ static double CeilingPow |
|
|
|
} else { dnm = Fn_Pow( dnm , 1.0 / ( 2.0 * pw ) ) ; } \ |
|
|
|
dnm = 1.0 / dnm ; \ |
|
|
|
y = (x) * (xmax) * dnm ; \ |
|
|
|
dx = (xmax) * xmp * dnm / arg ; \ |
|
|
|
dxmax = (x) * xp * dnm / arg ; \ |
|
|
|
dx = (xmax) * xmp * dnm / arg_cp2 ; \ |
|
|
|
dxmax = (x) * xp * dnm / arg_cp2 ; \ |
|
|
|
} |
|
|
|
|
|
|
|
#define Fn_SU_CP3( y , x , xmax , delta , pw , dx , dxmax , ddelta ) { \ |
|
|
|
|