Generalised floating point multiplication: Difference between revisions

{{works with|ALGOL 68G|Any - tested with release [http://sourceforge.net/projects/algol68/files/algol68g/algol68g-2.3.3 algol68g-2.3.3].}}

{{wont work with|ELLA ALGOL 68|Any (with appropriate job cards) - tested with release [http://sourceforge.net/projects/algol68/files/algol68toc/algol68toc-1.8.8d/algol68toc-1.8-8d.fc9.i386.rpm/download 1.8-8d] - due to extensive use of '''format'''[ted] ''transput''.}}

# TASK CODE #

# Actual generic mulitplication operator #

 

OP *:= = (REF DIGITS lhs, DIGIT arg)DIGITS: lhs := lhs * INITDIGITS arg;

 

################################################################

OD;

out SHR (UPB s-point)

####################################################################

# A program to test arbitrary length floating point multiplication #

printf($l$)

OD

<pre>

a = +523.23914037494284407864655 +-0++0+.+-0++0+

{{trans|Phix}}

In the interests of brevity many of the comments and all of the commented-out code has been omitted.

 

import (

test("septemvigesimal", septemVigesimal)

multTable()

 

{{out}}

10| +000 | +-000 | +0000 | ++000 | +--000 | +-0000 | +-+000 | +0-000 | +00000 | +0+000 | ++-000 | ++0000 |

</pre>

 

=={{header|Julia}}==

import Base.BigInt, Base.BigFloat, Base.print, Base.+, Base.-, Base.*

 

end

arr = reverse(digits(BigInt(round(x * big"3.0"^MAX_PRECISION)), base=3))

end

 

elseif i < 0

arr = [arr; b1.dig[ldigits1+1:end]; zeros(Int8, -i)]

arr2 = [arr2; b2.dig[ldigits2+1:end]]

end

end

 

a = T("+-0++0+.+-0++0+")

b = T(-436.436)

println("\na * (b - c) = ", String(a * (b - c)), "\n = ", format(BigFloat(a * (b - c))))

 

println(" x|+ (1) |+- (2) |+0 (3) |++ (4) |+-- (5)|+-0 (6)|+-+ (7)|+0- (8)|+e+-(9)|+0+(10)|++-(11)|++0(12)|")

for i in 1:27

 

code_reuse_task(BalancedTernary)

<pre>

a = +-0++0+.+-0++0+ = 523.23914

= -262510.90268

 

x|+ (1) |+- (2) |+0 (3) |++ (4) |+-- (5)|+-0 (6)|+-+ (7)|+0- (8)|+e+-(9)|+0+(10)|++-(11)|++0(12)|

1| +| +-| +0| ++| +--| +-0| +-+| +0-| +00| +0+| ++-| ++0|

27| +000| +-000| +0000| ++000| +--000| +-0000| +-+000| +0-000| +00000| +0+000| ++-000| ++0000|

</pre>

 

=={{header|Phix}}==

0.1 accurate to several million decimal places, but just never quite exact.

 

{{out}}

<pre>

=== multiplication table ===

Without e notation, with hexadecimal across, septemvigesimal down, and balanced ternary contents!

{{out}}

<pre>

* | #1 + | #2 +- | #3 +0 | #4 ++ | #5 +-- | #6 +-0 | #7 +-+ | #8 +0- | #9 +00 | #A +0+ | #B ++- | #C ++0 |

1 | + | | | | | | | | | | | |