60 lines
664 B
C++
60 lines
664 B
C++
// Subtract rational numbers
|
|
//
|
|
// Input: tos-2 minuend
|
|
//
|
|
// tos-1 subtrahend
|
|
//
|
|
// Output: difference on stack
|
|
|
|
#include "stdafx.h"
|
|
#include "defs.h"
|
|
|
|
void
|
|
qsub(void)
|
|
{
|
|
unsigned int *a, *ab, *b, *ba, *c;
|
|
|
|
save();
|
|
|
|
p2 = pop();
|
|
p1 = pop();
|
|
|
|
ab = mmul(p1->u.q.a, p2->u.q.b);
|
|
ba = mmul(p1->u.q.b, p2->u.q.a);
|
|
|
|
a = msub(ab, ba);
|
|
|
|
mfree(ab);
|
|
mfree(ba);
|
|
|
|
// zero?
|
|
|
|
if (MZERO(a)) {
|
|
mfree(a);
|
|
push(_zero);
|
|
restore();
|
|
return;
|
|
}
|
|
|
|
b = mmul(p1->u.q.b, p2->u.q.b);
|
|
|
|
c = mgcd(a, b);
|
|
|
|
MSIGN(c) = MSIGN(b);
|
|
|
|
p1 = alloc();
|
|
|
|
p1->k = NUM;
|
|
|
|
p1->u.q.a = mdiv(a, c);
|
|
p1->u.q.b = mdiv(b, c);
|
|
|
|
mfree(a);
|
|
mfree(b);
|
|
mfree(c);
|
|
|
|
push(p1);
|
|
|
|
restore();
|
|
}
|