projects
/
unix-history
/ blob
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
tags
|
clone url
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
blame
|
history
|
raw
|
HEAD
date and time created 83/08/05 13:35:20 by sam
[unix-history]
/
usr
/
src
/
old
/
libm
/
libom
/
pow.c
/* @(#)pow.c 4.2 %G% */
/*
computes a^b.
uses log and exp
*/
#include <errno.h>
int
errno
;
double
log
(),
exp
();
double
pow
(
arg1
,
arg2
)
double
arg1
,
arg2
;
{
double
temp
;
long
l
;
#ifdef vax
asm
(
" bispsw $0xe0"
);
#endif
if
(
arg1
<=
0
.) {
if
(
arg1
==
0
.) {
if
(
arg2
<=
0
.)
goto
domain
;
return
(
0
.);
}
l
=
arg2
;
if
(
l
!=
arg2
)
goto
domain
;
temp
=
exp
(
arg2
*
log
(-
arg1
));
if
(
l
&
1
)
temp
= -
temp
;
return
(
temp
);
}
return
(
exp
(
arg2
*
log
(
arg1
)));
domain
:
errno
=
EDOM
;
return
(
0
.);
}