stack32 out_stack
stack32 u_stack
int32 eax
int32 ebx
int32 esi
int32 edi
int32 ebp
int32 out
int32 u

int32 word0
int32 word1
int32 word2
int32 word3
int32 word4
float80 out0
float80 out1
float80 out2
float80 out3
float80 out4
float80 carry0
float80 carry1
float80 carry2
float80 carry3
float80 carry4
stack64 d0
stack64 d1
stack64 d2
stack64 d3
stack64 d4

enter gaudry_pm_todouble
input out_stack
input u_stack
caller eax
caller ebx
caller esi
caller edi
caller ebp

  u = u_stack

  word0 = *(uint32 *) (u + 0)
  d0 top = 0x43300000
  inplace d0 bottom = word0

  word1 = *(uint32 *) (u + 4)
  word1 &= 0xffffff
  d1 top = 0x45300000
  inplace d1 bottom = word1

  word2 = *(uint32 *) (u + 7)
  word2 &= 0xffffff
  d2 top = 0x46b00000
  inplace d2 bottom = word2

  word3 = *(uint32 *) (u + 10)
  word3 &= 0xffffff
  d3 top = 0x48300000
  inplace d3 bottom = word3

  word4 = *(uint32 *) (u + 12)
  (uint32) word4 >>= 8
  d4 top = 0x49b00000
  inplace d4 bottom = word4

  out = out_stack

  out4 = *(float64 *) &d4
  out4 -= *(float64 *) &gaudry_pm_in4offset
  carry4 = *(float64 *) &gaudry_pm_alpha127
  carry4 += out4
  carry4 -= *(float64 *) &gaudry_pm_alpha127
  out4 -= carry4

  out0 = *(float64 *) &d0
  out0 -= *(float64 *) &gaudry_pm_in0offset
  carry4 *= *(float64 *) &gaudry_pm_scale
  out0 += carry4
  carry0 = *(float64 *) &gaudry_pm_alpha26
  carry0 += out0
  carry0 -= *(float64 *) &gaudry_pm_alpha26
  out0 -= carry0
  *(float64 *) (out + 0) = out0

  out1 = *(float64 *) &d1
  out1 -= *(float64 *) &gaudry_pm_in1offset
  out1 += carry0
  carry1 = *(float64 *) &gaudry_pm_alpha51
  carry1 += out1
  carry1 -= *(float64 *) &gaudry_pm_alpha51
  out1 -= carry1
  *(float64 *) (out + 8) = out1

  out2 = *(float64 *) &d2
  out2 -= *(float64 *) &gaudry_pm_in2offset
  out2 += carry1
  carry2 = *(float64 *) &gaudry_pm_alpha77
  carry2 += out2
  carry2 -= *(float64 *) &gaudry_pm_alpha77
  out2 -= carry2
  *(float64 *) (out + 16) = out2

  out3 = *(float64 *) &d3
  out3 -= *(float64 *) &gaudry_pm_in3offset
  out3 += carry2
  carry3 = *(float64 *) &gaudry_pm_alpha102
  carry3 += out3
  carry3 -= *(float64 *) &gaudry_pm_alpha102
  out3 -= carry3
  *(float64 *) (out + 24) = out3

  out4 += carry3
  *(float64 *) (out + 32) = out4

leave