mod.coffee

Eval_mod = ->
	push(cadr(p1))
	Eval()
	push(caddr(p1))
	Eval()
	mod()

mod = ->
	n = 0

	save()

	p2 = pop()
	p1 = pop()

	if (isZeroAtomOrTensor(p2))
		stop("mod function: divide by zero")

	if (!isNumericAtom(p1) || !isNumericAtom(p2))
		push_symbol(MOD)
		push(p1)
		push(p2)
		list(3)
		restore()
		return

	if (isdouble(p1))
		push(p1)
		n = pop_integer()
		if (isNaN(n))
			stop("mod function: cannot convert float value to integer")
		push_integer(n)
		p1 = pop()

	if (isdouble(p2))
		push(p2)
		n = pop_integer()
		if (isNaN(n))
			stop("mod function: cannot convert float value to integer")
		push_integer(n)
		p2 = pop()

	if (!isinteger(p1) || !isinteger(p2))
		stop("mod function: integer arguments expected")

	p3 = new U()
	p3.k = NUM
	p3.q.a = mmod(p1.q.a, p2.q.a)
	p3.q.b = mint(1)
	push(p3)

	restore()