indentation
This commit is contained in:
parent
cf734f60ac
commit
a2d0d8031f
File diff suppressed because it is too large
Load Diff
|
|
@ -24,7 +24,7 @@ extern double asinh(), acosh(), atanh();
|
||||||
|
|
||||||
double PTfudge_factor;
|
double PTfudge_factor;
|
||||||
|
|
||||||
#define MODULUS(NUM,LIMIT) ((NUM) - ((int) ((NUM) / (LIMIT))) * (LIMIT))
|
#define MODULUS(NUM,LIMIT) ((NUM) - ((int) ((NUM) / (LIMIT))) * (LIMIT))
|
||||||
|
|
||||||
double
|
double
|
||||||
PTabs(double arg)
|
PTabs(double arg)
|
||||||
|
|
@ -60,9 +60,9 @@ double
|
||||||
PTdivide(double arg1, double arg2)
|
PTdivide(double arg1, double arg2)
|
||||||
{
|
{
|
||||||
if (arg2 >= 0.0)
|
if (arg2 >= 0.0)
|
||||||
arg2 += PTfudge_factor;
|
arg2 += PTfudge_factor;
|
||||||
else
|
else
|
||||||
arg2 -= PTfudge_factor;
|
arg2 -= PTfudge_factor;
|
||||||
|
|
||||||
if (arg2 == 0.0)
|
if (arg2 == 0.0)
|
||||||
return (HUGE);
|
return (HUGE);
|
||||||
|
|
@ -74,11 +74,11 @@ double
|
||||||
PTpower(double arg1, double arg2)
|
PTpower(double arg1, double arg2)
|
||||||
{
|
{
|
||||||
if (arg1 < 0.0) {
|
if (arg1 < 0.0) {
|
||||||
if (fabs(arg2 - ((int) arg2)) / (arg2 + 0.001) < 0.000001) {
|
if (fabs(arg2 - ((int) arg2)) / (arg2 + 0.001) < 0.000001) {
|
||||||
arg2 = (int) arg2;
|
arg2 = (int) arg2;
|
||||||
} else {
|
} else {
|
||||||
arg1 = -arg1;
|
arg1 = -arg1;
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
return (pow(arg1, arg2));
|
return (pow(arg1, arg2));
|
||||||
}
|
}
|
||||||
|
|
@ -108,7 +108,7 @@ PTacosh(double arg)
|
||||||
return (acosh(arg));
|
return (acosh(arg));
|
||||||
#else
|
#else
|
||||||
if (arg < 1.0)
|
if (arg < 1.0)
|
||||||
arg = 1.0;
|
arg = 1.0;
|
||||||
return (log(arg + sqrt(arg*arg-1.0)));
|
return (log(arg + sqrt(arg*arg-1.0)));
|
||||||
#endif
|
#endif
|
||||||
}
|
}
|
||||||
|
|
@ -142,9 +142,9 @@ PTatanh(double arg)
|
||||||
return (atanh(arg));
|
return (atanh(arg));
|
||||||
#else
|
#else
|
||||||
if (arg < -1.0)
|
if (arg < -1.0)
|
||||||
arg = -1.0 + PTfudge_factor + 1e-10;
|
arg = -1.0 + PTfudge_factor + 1e-10;
|
||||||
else if (arg > 1.0)
|
else if (arg > 1.0)
|
||||||
arg = 1.0 - PTfudge_factor - 1e-10;
|
arg = 1.0 - PTfudge_factor - 1e-10;
|
||||||
return (log((1.0 + arg) / (1.0 - arg)) / 2.0);
|
return (log((1.0 + arg) / (1.0 - arg)) / 2.0);
|
||||||
#endif
|
#endif
|
||||||
}
|
}
|
||||||
|
|
@ -153,11 +153,11 @@ double
|
||||||
PTustep(double arg)
|
PTustep(double arg)
|
||||||
{
|
{
|
||||||
if (arg < 0.0)
|
if (arg < 0.0)
|
||||||
return 0.0;
|
return 0.0;
|
||||||
else if (arg > 0.0)
|
else if (arg > 0.0)
|
||||||
return 1.0;
|
return 1.0;
|
||||||
else
|
else
|
||||||
return 0.5; /* Ick! */
|
return 0.5; /* Ick! */
|
||||||
}
|
}
|
||||||
|
|
||||||
/* MW. PTcif is like "C" if - 0 for (arg<=0), 1 elsewhere */
|
/* MW. PTcif is like "C" if - 0 for (arg<=0), 1 elsewhere */
|
||||||
|
|
@ -166,11 +166,11 @@ double
|
||||||
PTustep2(double arg)
|
PTustep2(double arg)
|
||||||
{
|
{
|
||||||
if (arg <= 0.0)
|
if (arg <= 0.0)
|
||||||
return 0.0;
|
return 0.0;
|
||||||
else if (arg <= 1.0)
|
else if (arg <= 1.0)
|
||||||
return arg;
|
return arg;
|
||||||
else /* if (arg > 1.0) */
|
else /* if (arg > 1.0) */
|
||||||
return 1.0;
|
return 1.0;
|
||||||
}
|
}
|
||||||
|
|
||||||
double
|
double
|
||||||
|
|
@ -213,9 +213,9 @@ double
|
||||||
PTuramp(double arg)
|
PTuramp(double arg)
|
||||||
{
|
{
|
||||||
if (arg < 0.0)
|
if (arg < 0.0)
|
||||||
return 0.0;
|
return 0.0;
|
||||||
else
|
else
|
||||||
return arg;
|
return arg;
|
||||||
}
|
}
|
||||||
|
|
||||||
double
|
double
|
||||||
|
|
|
||||||
Loading…
Reference in New Issue