1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
%% Fixes for c-code generation
%% Set PI explicitly here to avoid later conflict with cc headers
%% if PI is not already a number (i.e. on rounded has not been set)
IF NOT NUMBERP (pi) THEN LET PI = 3.14159$
ON ROUNDED$ % No integer output
%% Changes x^y to pow(x,y)
OPERATOR pow$
FOR ALL x,y LET x^y = pow(x,y)$ % Use the pow function
|
|
>
>
|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
|
%% Fixes for c-code generation
%% Set PI explicitly here to avoid later conflict with cc headers
%% if PI is not already a number (i.e. on rounded has not been set)
IF NOT NUMBERP (pi) THEN LET PI = 3.14159$
ON ROUNDED$ % No integer output
%% Changes x^y to pow(x,y)
OPERATOR pow$
FOR ALL x,y LET x^y = pow(x,y)$ % Use the pow function
%% Derivatives
FOR ALL f,g,x LET df(pow(f,g),x)=
pow(f,g-1) * (df(f,x)*g + df(g,x)*f*log(f))$
%% Special cases
FOR ALL x LET pow(x,0) = 1$
FOR ALL x LET pow(x,1) = x$
FOR ALL x,y,z LET pow(x,y)*pow(x,z) = pow(x,y+z)$
FOR ALL x,y,z LET pow(pow(x,y),z) = pow(x,y*z)$
OPERATOR fabs$
FOR ALL x let abs(x) = fabs(x)$
END$
|