Algebrite/sources/qdiv.ts at master · davidedc/Algebrite

History

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

import { Constants, MZERO, Num } from '../runtime/defs';

import { stop } from '../runtime/run';

import { makeSignSameAs } from './bignum';

import { mgcd } from './mgcd';

import { mdiv, mmul } from './mmul';

// Divide rational numbers

//

// Input: p1 dividend

// p2 divisor

//

// Output: quotient

export function qdiv(p1: Num, p2: Num): Num {

// zero?

if (MZERO(p2.q.a)) {

stop('divide by zero');

}

if (MZERO(p1.q.a)) {

return Constants.zero;

}

const aa = mmul(p1.q.a, p2.q.b);

const bb = mmul(p1.q.b, p2.q.a);

let c = mgcd(aa, bb);

c = makeSignSameAs(c, bb);

return new Num(mdiv(aa, c), mdiv(bb, c));

}

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

qdiv.ts

qdiv.ts

Files

qdiv.ts

Latest commit

History

qdiv.ts

File metadata and controls