Imaginary base numbers: Difference between revisions

Line 2,951:

=={{header|Phix}}==

<lang Phix>~~include complex.e~~

with javascript_semantics

include complex.e

function base2(atom num, integer radix, precision = -8)

if radix<-36 or radix>-2 then throw("radix out of range (-2..-36)") end if

function base2(atom num, integer radix, precision = -8)

sequence result

if radix<-36 or radix>-2 then throw("radix out of range (-2..-36)") end if

if num=0 then

sequence result

result = {"0",""}

if num=0 then

else

result = {"0",""}

integer place = 0

else

result = ""

integer place = 0

atom v = num

result = ""

atom upper_bound = 1/(1-radix),

atom v = num

lower_bound = radix*upper_bound

atom upper_bound = 1/(1-radix),

while not(lower_bound <= v) or not(v < upper_bound) do

lower_bound = radix*upper_bound

place += 1

while not(lower_bound <= v) or not(v < upper_bound) do

v = num/power(radix,place)

place += 1

end while

v = num/power(radix,place)

end while

while (v or place > 0) and (place > precision) do

integer digit = floor(radix*v - lower_bound)

while (v or place > 0) and (place > precision) do

v = (radix*v - digit)

integer digit = floor(radix*v - lower_bound)

if place=0 and not find('.',result) then result &= '.' end if

v = (radix*v - digit)

result &= digit+iff(digit>9?'a'-10:'0')

if place=0 and not find('.',result) then result &= '.' end if

place -= 1

result &= digit+iff(digit>9?'a'-10:'0')

end while

place -= 1

integer dot = find('.',result)

end while

if dot then

integer dot = find('.',result)

result = trim_tail(result,'0')

if dot then

result = {result[1..dot-1],result[dot+1..$]}

result = trim_tail(result,'0')

else

result = {result[1..dot-1],result[dot+1..$]}

result = {result,""}

else

end if

result = {result,""}

end if

end if

return result

end if

end function

return result

end function

function zip(string a, string b)

integer ld = length(a)-length(b)

function zip(string a, string b)

if ld!=0 then

integer ld = length(a)-length(b)

if ld>0 then

if ld!=0 then

b &= repeat('0',ld)

if ld>0 then

else

b &= repeat('0',ld)

a &= repeat('0',abs(ld))

else

end if

a &= repeat('0',abs(ld))

end if

end if

string res = ""

end if

for i=1 to length(a) do

string res = ""

res &= a[i]&b[i]

for i=1 to length(a) do

end for

res &= a[i]&b[i]

res = trim_tail(res,'0')

end for

if res="" then res = "0" end if

res = trim_tail(res,'0')

return res

if res="" then res = "0" end if

end function

return res

end function

function base(complexn num, integer radix, precision = -8)

function base(complexn num, integer radix, precision = -8)

integer absrad = abs(radix),

radix2 = -power(radix,2)

integer absrad = abs(radix),

if absrad<2 or absrad>6 then throw("base radix out of range") end if

radix2 = -power(radix,2)

if absrad<2 or absrad>6 then throw("base radix out of range") end if

atom {re, im} = {complex_real(num), complex_imag(num)}

string {re_wh, re_fr} = base2(re, radix2, precision),

atom {re, im} = {complex_real(num), complex_imag(num)}

{im_wh, im_fr} = base2(im/radix, radix2, precision)

string {re_wh, re_fr} = base2(re, radix2, precision),

{im_wh, im_fr} = base2(im/radix, radix2, precision)

string whole = reverse(zip(reverse(re_wh), reverse(im_wh))),

fraction = zip(im_fr, re_fr)

string whole = reverse(zip(reverse(re_wh), reverse(im_wh))),

if fraction!="0" then whole &= '.'&fraction end if

fraction = zip(im_fr, re_fr)

return whole

if fraction!="0" then whole &= '.'&fraction end if

end function

return whole

end function

function parse_base(string str, integer radix)

function parse_base(string str, integer radix)

complexn fraction = 0

complexn fraction = 0

integer dot = find('.',str)

if dot then

integer dot = find('.',str)

string fr = str[dot+1..$]

if dot then

for i=1 to length(fr) do

string fr = str[dot+1..$]

integer c = fr[i]

for i=1 to length(fr) do

c -= iff(c>='a'?'a'-10:'0')

integer c = fr[i]

fraction = complex_add(fraction,complex_mul(c,complex_power({0,radix},-i)))

c -= iff(c>='a'?'a'-10:'0')

end for

fraction = complex_add(fraction,complex_mul(c,complex_power({0,radix},-i)))

str = str[1..dot-1]

end for

end if

str = str[1..dot-1]

end if

str = reverse(str)

for i=1 to length(str) do

str = reverse(str)

integer c = str[i]

for i=1 to length(str) do

c -= iff(c>='a'?'a'-10:'0')

integer c = str[i]

fraction = complex_add(fraction,complex_mul(c,complex_power({0,radix},(i-1))))

c -= iff(c>='a'?'a'-10:'0')

end for

fraction = complex_add(fraction,complex_mul(c,complex_power({0,radix},(i-1))))

end for

return fraction

end function

return fraction

end function

constant tests = {{0,2},{1,2},{5,2},{-13,2},{{0,9},2},{{0,-3},2},{{7.75,-7.5}, 2},{.25, 2}, -- base 2i tests

{{5,5}, 2},{{5,5}, 3},{{5,5}, 4},{{5,5}, 5},{{5,5}, 6}, -- same value, positive imaginary bases

constant tests = {{0,2},{1,2},{5,2},{-13,2},{{0,9},2},{{0,-3},2},{{7.75,-7.5}, 2},{.25, 2}, -- base 2i tests

{{5,5},-2},{{5,5},-3},{{5,5},-4},{{5,5},-5},{{5,5},-6}, -- same value, negative imaginary bases

{{5,5}, 2},{{5,5}, 3},{{5,5}, 4},{{5,5}, 5},{{5,5}, 6}, -- same value, positive imaginary bases

{{227.65625,10.859375},4}, -- larger test value

{{5,5},-2},{{5,5},-3},{{5,5},-4},{{5,5},-5},{{5,5},-6}, -- same value, negative imaginary bases

{{-579.8225308641975744,-5296.406378600824},6}} -- phix.rules

{{227.65625,10.859375},4}, -- larger test value

{{-579.8225308641975744,-5296.406378600824},6}} -- phix.rules

-- matches output of Sidef and Raku:

for t=1 to length(tests) do

-- matches output of Sidef and Raku:

{complexn v, integer r} = tests[t]

for t=1 to length(tests) do

string ibase = base(v,r),

{complexn v, integer r} = tests[t]

strv = complex_sprint(v),

string ibase = base(v,r),

strb = complex_sprint(parse_base(ibase, r))

strv = complex_sprint(v),

printf(1,"base(%20s, %2di) = %-10s : parse_base(%12s, %2di) = %s\n",

strb = complex_sprint(parse_base(ibase, r))

{strv, r, ibase, '"'&ibase&'"', r, strb})

printf(1,"base(%20s, %2di) = %-10s : parse_base(%12s, %2di) = %s\n",

end for

{strv, r, ibase, '"'&ibase&'"', r, strb})

end for

-- matches output of Kotlin, Java, Go, D, and C#:

for ri=1 to 2 do -- real then imag

-- matches output of Kotlin, Java, Go, D, and C#:

for i=1 to 16 do

for ri=1 to 2 do -- real then imag

complexn c = iff(ri=1?i:{0,i}),

for i=1 to 16 do

nc = complex_neg(c)

complexn c = iff(ri=1?i:{0,i}),

string sc = complex_sprint(c),

nc = complex_neg(c)

snc = complex_sprint(nc),

string sc = complex_sprint(c),

ib = base(c,2),

snc = complex_sprint(nc),

inb = base(nc,2),

ib = base(c,2),

rc = complex_sprint(parse_base(ib,2)),

inb = base(nc,2),

rnc = complex_sprint(parse_base(inb,2))

rc = complex_sprint(parse_base(ib,2)),

printf(1,"%4s -> %8s -> %4s %4s -> %8s -> %4s\n",

rnc = complex_sprint(parse_base(inb,2))

{sc, ib, rc, snc, inb, rnc })

printf(1,"%4s -> %8s -> %4s %4s -> %8s -> %4s\n",

end for

{sc, ib, rc, snc, inb, rnc })

puts(1,"\n")

end for

end for</lang>

puts(1,"\n")

end for

Matches the output of Sidef and Raku, except for the final line: