UNPKG

algebrite

Version:

Computer Algebra System in Coffeescript

github.com/davidedc/Algebrite

davidedc/Algebrite

62 lines (41 loc) • 1.6 kB
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
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
# Add rational numbers # # Input: tos-2 addend # # tos-1 addend # # Output: sum on stack qadd = -> # a, qadd_ab, b, qadd_ba, c are all bigNum # we are adding the fractions qadd_frac1 + qadd_frac2 i.e. # qadd_frac1.q.a/qadd_frac1.q.b + qadd_frac2.q.a/qadd_frac2.q.b qadd_frac2 = pop() qadd_frac1 = pop() qadd_ab = mmul(qadd_frac1.q.a, qadd_frac2.q.b) qadd_ba = mmul(qadd_frac1.q.b, qadd_frac2.q.a) qadd_numerator = madd(qadd_ab, qadd_ba) #mfree(qadd_ab) #mfree(qadd_ba) # zero? if (MZERO(qadd_numerator)) #console.log "qadd IS ZERO" #mfree(qadd_numerator) push(zero) return qadd_denominator = mmul(qadd_frac1.q.b, qadd_frac2.q.b) gcdBetweenNumeratorAndDenominator = mgcd(qadd_numerator, qadd_denominator) #console.log "gcd("+qadd_numerator+","+qadd_denominator+"): " + gcdBetweenNumeratorAndDenominator gcdBetweenNumeratorAndDenominator = makeSignSameAs(gcdBetweenNumeratorAndDenominator,qadd_denominator) #console.log "qadd qadd_denominator: " + qadd_denominator #console.log "qadd gcdBetweenNumeratorAndDenominator: " + gcdBetweenNumeratorAndDenominator resultSum = new U() resultSum.k = NUM resultSum.q.a = mdiv(qadd_numerator, gcdBetweenNumeratorAndDenominator) resultSum.q.b = mdiv(qadd_denominator, gcdBetweenNumeratorAndDenominator) #console.log "qadd resultSum.q.a: " + resultSum.q.a #console.log "qadd resultSum.q.b: " + resultSum.q.b #mfree(qadd_numerator) #mfree(qadd_denominator) #mfree(gcdBetweenNumeratorAndDenominator) push(resultSum) #console.log "qadd result: " + resultSum