Exercise 4-5 (mathops - Calculator)

Chapter_4     Exercise_4-4

Exercise_4-6

---

---

Exercise 4-5     K&R, p. 79

Exercise 4-5. Add access to library functions like sin(), exp(), and pow(). See <math.h> in Appendix B, Section 4.

mathops.c         download

#include <stdio.h> // for printf(), getchar(), EOF
#include <stdlib.h> // for atof()
#include <math.h> // for sin(), exp(), pow()

#define MAXOP 100 // max size of operand or operator
#define NUMBER '0' // signal that a number has been found

int getop(char []); // get the next operand or operator
void push(double); // push operand on stack
double pop(void); // pop operand from stack

int main()
{
int type, i1, i2;
double op2; // operand
char s[MAXOP]; // operand or operator

while((type = getop(s)) != EOF && type != '\0')
{ // nonempty line ending with EOF will not be executed
switch(type)
{
case NUMBER :
push(atof(s));
break;
case '+' :
push(pop() + pop());
break;
case '-' :
op2 = pop();
push(pop() - op2);
break;
case '*' :
push(pop() * pop());
break;
case '/' :
op2 = pop();
if (op2 != 0.0)
{push(pop() / op2);}
else {printf("Error: zero divisor\n");}
break;
case '%' :
i2 = pop(); // automatic conversions
i1 = pop(); // from double to int
if (i2 != 0.0)
{push(i1 % i2);} // auto conversion from int to double, arg of push()
else {printf("Error: zero divisor\n");}
break;
case 's' : // sin()
push(sin(pop()));
break;
case 'e' : // exp(x), exponential function e^x
push(exp(pop()));
break;
case 'p' : // pow(x,y)
op2 = pop();
push(pow(pop(), op2));
break;
case '\n' :
printf("\t%.8g\n", pop()); // top of stack
break;
default :
printf("Unknown command: %s\n", s);
break;
}
}

return 0;
}

#define MAXVAL 100 // stack size
// stack pointer (0 is top of stack):
int sp = 0; // next free stack position
double val[MAXVAL]; // value stack

void push(double f) // push f on (top of) stack
{
if (sp < MAXVAL) {val[sp++] = f;}
else {printf("Error: stack full, can't push %g\n", f);}
}

double pop(void) // pop and return top value from stack
{
if (sp > 0) {return val[--sp];}
else
{
printf("Error: stack empty\n");
return 0.0;
}
}

#include <ctype.h> // for isdigit()

int getch(void);
void ungetch(int);

int getop(char s[]) // get next operator or numeric operand
{
int c, i;

while ((s[0] = c = getch()) == ' ' || c == '\t') // skip beginning ' ', '\t'
{} // last value read in s[0] is not ' ' or '\t', but could be '\n' or EOF
s[1] = '\0'; // end string

if (!isdigit(c) && c != '.' && c != '+' && c != '-')
{return c;} // not a number, probably an operator (or '\n' or '\0' or EOF)
// here c is digit or '.' or '+' or '-'

i = 0;

if (c == '+' || c == '-') // unary or binary sign
{
c = getch();
if (!isdigit(c) && c != '.')
{
ungetch(c);
return s[0]; // binary '+' or '-'
}
else // c is digit or '.'
{
if (s[i] == '-') {i++;} // s[0]
// if (s[i] == '+') {} // skip adding '+' to s[]
s[i] = c; // digit or '.'
}
}

if (isdigit(c)) // collect integer part
{ // s[i] contains c (digit)
while (isdigit(s[++i] = c = getch()))
{}
} // here s[i] is not a digit
if (c == '.') // collect fraction part
{ // s[i] contains c (s[i] == '.')
while (isdigit(s[++i] = c = getch()))
{}
} // here s[i] is not a digit
s[i] = '\0'; // end string containing number (operand)

if (c != EOF) {ungetch(c);} // cannot ungetch(EOF) as EOF is not a char
else {ungetch('\0');} // EOF or '\0' returned by getop() will end the program

return NUMBER;
}

#define BUFSIZE 100

char buf[BUFSIZE]; // buffer for ungetch()
// buf[] pointer (0 is top of buf):
int bufp = 0; // next free position in buf[]

int getch(void) // get a (possibly pushed-back) character
{
return (bufp > 0) ? buf[--bufp] : getchar();
}
// push character back on input (make it available for the next getch()):
void ungetch(int c)
{
if (bufp >= BUFSIZE)
{printf("ungetch(): too many characters\n");}
else {buf[bufp++] = c;}
}
/*
gcc mathops.c -o mathops -lm // link math library
./mathops
0 s // sin(0)
0
3.14 s // sin(pi)
0.0015926529
1.57 s // sin(pi/2)
0.99999968
1 e // e^1
2.7182818
2 e // e^2
7.3890561
3 2 p // 3^2
9
2 3 p // 2^3
8
1 2 +
3
1 2 + 3 * 4 /
2.25
-1 2 +
1
1 2 -
-1 // 1-2
-1 -2 +
-3
-1 2 + 3 *
3
2 0 /
Error: zero divisor
2 // pop() called by '\n' returns first operand, still on stack
1 2 /
0.5
2 3 %
2
4 2 %
0
4 1 %
0
4 3 %
1
4 0 %
Error: zero divisor
Error: stack empty // nothing pushed, '\n' calls pop():
0 // value returned by pop() for an empty stack
a b +
Unknown command: a
Unknown command: b
Error: stack empty // nothing pushed, pop() called by +
Error: stack empty // pop() called by '\n'
0 // value returned by pop() for an empty stack
1 2 9
9 // pop() called by '\n'
+
3 // 1 + 2
1 2 c
Unknown command: c
2 // pop() called by '\n'
3 +
4 // 1 + 3
CTRL^D (EOF) or CTRL^C in Linux, CTRL^Z + Enter in Windows

./mathops < compute.txt
3 // 1 2 +
0 // 1 2 * 3 + 5 -
1 // -1 2 +
1 // -1 -2 -
1 // 4 3 %
Error: zero divisor // 2 0 %
Error: stack empty // nothing pushed, '\n' calls pop():
0 // value returned by pop() for an empty stack
0.33333333 // 1 3 /
Error: zero divisor // 1 0 /
1 // pop() called by '\n' returns first operand, still on stack
*/

Note:  See Exercise_4-3 for the file compute.txt.

---

---

Chapter_4     Exercise_4-4

BACK_TO_TOP

Exercise_4-6