non-overflowing-multiplaction.thump

// 21*23/31=15 (.5806451613)

pi $0, si 0
pi $1, si 1

pi $fac1, si 21
pi $fac2, si 23
pi $div, si 31

@muldiv
xi 0, yi 0, zi 4, s $fac1, sf $fac2

pi $result, si 0
pi $_x, si 0
pi $_r, si 0
@calc
x $_x, yi 0, zi 4, g $v
pi $_f, s $v, div $div
pi $result, add $_f
pi $_r2, s $v, mod $div
pi $_r, add $_r2
pi $_f, s $_r, div $div
pi $_r, mod $div
pi $result, add $_f
pi $_x, add $1
pi $v, seq $0, j @calc

l $result