function test:float[](floatarray#[])
  local farr1:int ptr=int ptr varptr floatarray[0]
  local farr2:int ptr=int ptr memalloc(floatarray.length*4)
' do something here
  return float[] varptr farr2 '?????
endfunction

function radixsortfloat#[](floatarray#[])
  local floatarray2#[floatarray.length]
  local farr3%ptr=int ptr varptr floatarray[0]
  local farr4%ptr=int ptr varptr floatarray2[0]
'  ...
  ret farr2

extern "C"
  function memset(m:byte ptr,val%,sz%)
endextern

function radixsortfloat#[](arr#[])
  local pos%
  local arr2#[arr.length]',swp#[]
  local arr3%ptr=int ptr varptr arr[0]
  local arr4%ptr=int ptr varptr arr2[0]
  local cnt%[256]
  repeat
    for local a%=0 until arr.length
      local m%=(arr3[a]~$80000000) shr pos&255 ~$7f shr (~arr3[a] shr 28&$8) 'adjust negative numbers
'      local m%=(arr3[a]~$80000000) shr pos&255 ~$7f*(arr3[a]<0) 'adjust negative numbers
      cnt[m]:+1
    next

    for local a%=1 until 256
      cnt[a]:+cnt[a-1]
    next

    for local a%=arr.length-1 to 0 step-1
      local m%=(arr3[a]~$80000000) shr pos&255 ~$7f shr (~arr3[a] shr 28&$8) 'adjust negative numbers
'      local m%=(arr3[a]~$80000000) shr pos&255 ~$7f*(arr3[a]<0) 'adjust negative numbers
      cnt[m]:-1
      arr4[cnt[m]]=arr3[a]
    next
    pos:+8
    if (pos>=32) then return arr2
    local swp2#[]=arr
    arr=arr2
    arr2=swp2
    arr3=int ptr varptr arr[0]
    arr4=int ptr varptr arr2[0]
'    local swp%ptr=arr3
'    arr3=arr4
'    arr4=swp

    setmem byte ptr cnt,0,256*4
  forever
endfunction

For local a%=0 until 256
  cnt[a]=0
next

;Radix.s

format MS COFF
section "code" code align 4
;in esp+4 int ptr
;in esp+8 val
;in esp+12 size
public _fillint
_fillint:
  mov ecx,[esp+12] ;size
  or ecx,ecx
 jz .out
  mov eax,[esp+8] ;val
   push edi
  mov edi,[esp+4+4] ;data
  rep stosd
   pop edi
.out:
ret

import "radix.s"
extern
  function fillint(iarr%ptr,val%,sz%)
endextern

function radixsortint%[](arr%[])
  local pos%
  local arr2%[arr.length]
  local cnt%[256]
  repeat
    for local a%=0 until arr.length
      local m%=(arr[a]) shr pos&255 ~$80 'adjust negative numbers
      cnt[m]:+1
    next

    for local a%=1 until 256
      cnt[a]:+cnt[a-1]
    next

    for local a%=arr.length-1 to 0 step-1
      local m%=(arr[a]) shr pos&255 ~$80 'adjust negative numbers
      cnt[m]:-1
      arr2[cnt[m]]=arr[a]
    next
    pos:+8
    if (pos>=32) then return arr2
    local swp%[]=arr
    arr=arr2
    arr2=swp

    fillint int ptr cnt,0,256
  forever
endfunction

function radixsortfloat#[](arr#[])
  local pos%
  local arr2#[arr.length]',swp#[]
  local arr3%ptr=int ptr varptr arr[0]
  local arr4%ptr=int ptr varptr arr2[0]
  local cnt%[256]
  repeat
    for local a%=0 until arr.length
      local m%=-arr3[a] shr pos&255 ~$7f shr (arr3[a] shr 28&$8) 'adjust negative numbers
      cnt[m]:+1
    next

    for local a%=1 until 256
      cnt[a]:+cnt[a-1]
    next

    for local a%=0 until arr.length
      local m%=-arr3[a] shr pos&255 ~$7f shr (arr3[a] shr 28&$8) 'adjust negative numbers
      cnt[m]:-1
      arr4[cnt[m]]=arr3[a]
    next
    pos:+8
    if (pos>=32) then return arr2
    local swp2#[]=arr
    arr=arr2
    arr2=swp2
    arr3=int ptr varptr arr[0]
    arr4=int ptr varptr arr2[0]

    fillint int ptr cnt,0,256
  forever
endfunction

function radixsortstr$[](arr$[],maxlen%=-1)
  if (maxlen<0) then for local a%=0 until arr.length
    if (maxlen<arr[a].length) then maxlen=arr[a].length
  next
  if not maxlen then return arr
  maxlen:-1
  local arr2$[arr.length]
  local cnt%[256]
  repeat
    for local a%=0 until arr.length
      local m%
      if arr[a] then
        m=arr[a][min(maxlen,arr[a].length-1)]&255
        if (m>=asc"A") and (m<=asc"Z") then m:+asc"a"-asc"A"
      endif
      cnt[m]:+1
    next

    for local a%=1 until 256
      cnt[a]:+cnt[a-1]
    next

    for local a%=0 until arr.length
      local m%
      if arr[a] then
        m=arr[a][min(maxlen,arr[a].length-1)]&255
        if (m>=asc"A") and (m<=asc"Z") then m:+asc"a"-asc"A"
      endif
      cnt[m]:-1
      arr2[cnt[m]]=arr[a]
    next
    maxlen:-1
    if (maxlen<0) then return arr2
    local swp$[]=arr
    arr=arr2
    arr2=swp

    fillint int ptr cnt,0,256
  forever
endfunction

const maxint%=$7fffffff
global n%[10]
seedrnd 10
for local a%=0 until n.length
  n[a]=int(rnd(maxint))+int(rnd(maxint)) 'trick 17 :D times 2 would mean you only have equal numbers, no odds
next
for local a%=0 until n.length
  print n[a]
next
print
n=radixsortint(n)
for local a%=0 until n.length
  print n[a]
next
input

global n4#[10]
'n4[0]=5.5
'print float2fpu(n4[0]) 'varptr n4[0]
'n4[0]=-5.5
'print float2fpu(n4[0])
seedrnd 15
for local a%=0 until n4.length
  n4[a]=rnd(10)-5
next
for local a%=0 until n4.length
  print n4[a]
next
print
n4=radixsortfloat(n4)
'qsortfloat n4,0,n4.length-1
for local a%=0 until n4.length
  print n4[a]
next
input

global n2$[10],n3$[10]
seedrnd 118
for local a%=0 until n2.length
for local b%=0 until rnd(7)+1
  n2[a]:+chr(asc("A")+32*int(rnd(2))+rnd(24))
next;next
n2[1]="va"
for local a%=0 until n2.length
  n3[a]=n2[a]
  print n2[a]
next
print
'qsortstr n3,0,n3.length-1
'n3.sort
'n2=radixsortstrlen(n2)
n2=radixsortstr(n2)
for local a%=0 until n2.length
  print n2[a]+":"+n3[a]
next
input
end

function radixsortint%[](arr%[])
  local cnt%[256]
  local pos%
  local arr2%[arr.length]
  repeat
    for local a%=0 until arr.length
      local m%=(arr[a]~$80000000) shr pos&255 'adjust negative numbers
      cnt[m]:+1
    next

    for local a%=1 until 256
      cnt[a]:+cnt[a-1]
    next

    for local a%=arr.length-1 to 0 step-1
      local m%=(arr[a]~$80000000) shr pos&255 'adjust negative numbers
      cnt[m]:-1
      arr2[cnt[m]]=arr[a]
    next

    pos:+8
    if (pos=32) then return arr2
    local swp%[]=arr
    arr=arr2
    arr2=swp

    for local a%=0 until 256
      cnt[a]=0
    next
  forever
endfunction

function radixsortfloat#[](source#[])
  local cnt%[256]
  local temp#[]
  local negativecount%
  local pos%,m%
  local result:float[source.length]
'  local srclen%=source.length-1
  local sourceptr%ptr=int ptr varptr source[0]
  local resultptr%ptr=int ptr varptr result[0]
  repeat
    for local a%=0 until source.length
      m=sourceptr[a] shr pos&$ff
      cnt[m]:+1
      if (pos=24) then negativecount:+m shr 7
    next

    for local a%=1 until 256
      cnt[a]:+cnt[a-1]
    next

    for local a%=source.length-1 to 0 step-1
      m=sourceptr[a] shr pos&$ff
      cnt[m]:-1

      if (pos=24) then
        if (m shr 7) then resultptr[srclen-cnt[m]]=sourceptr[a]..
                     else resultptr[cnt[m]+negativecount]=sourceptr[a]
      else
        resultptr[cnt[m]]=sourceptr[a]
      endif
    next

    pos:+8
    if (pos=32) then return result
    temp=source
    source=result
    result=temp

    sourceptr=int ptr varptr source[0]
    resultptr=int ptr varptr result[0]

    for local a%=0 until 256
      cnt[a]=0
    next
  forever
endfunction

function radixsortstr$[](arr$[],maxlen%=-1)
  if (maxlen<0) then for local a%=0 until arr.length
    if (maxlen<arr[a].length) then maxlen=arr[a].length
  next
  if not maxlen then return arr
  maxlen:-1
  local arr2$[arr.length]
  local cnt%[256]
  repeat
    for local a%=0 until arr.length
      local m%
      if arr[a] then
        m=arr[a][min(maxlen,arr[a].length-1)]&255
        if (m>=asc"A") and (m<=asc"Z") then m:+asc"a"-asc"A"
      endif
      cnt[m]:+1
    next

    for local a%=1 until 256
      cnt[a]:+cnt[a-1]
    next

    for local a%=0 until arr.length
      local m%
      if arr[a] then
        m=arr[a][min(maxlen,arr[a].length-1)]&255
        if (m>=asc"A") and (m<=asc"Z") then m:+asc"a"-asc"A"
      endif
      cnt[m]:-1
      arr2[cnt[m]]=arr[a]
    next

    if not maxlen then return arr2
    maxlen:-1
    local swp$[]=arr
    arr=arr2
    arr2=swp

    for local a%=0 until 256
      cnt[a]=0
    next
  forever
endfunction