stack32 out_stack
int32 out
stack32 u_stack
int32 u
stack32 v_stack
int32 v
int32 eax
int32 ebx
int32 esi
int32 edi
int32 ebp
float80 out0
float80 out1
float80 out2
float80 out3
float80 out4
float80 carry0
float80 carry1
float80 carry2
float80 carry3
float80 carry4
float80 u0v0
float80 u0v1
float80 u0v2
float80 u0v3
float80 u0v4
float80 u1v0
float80 u1v1
float80 u1v2
float80 u1v3
float80 u1v4
float80 u2v0
float80 u2v1
float80 u2v2
float80 u2v3
float80 u2v4
float80 u3v0
float80 u3v1
float80 u3v2
float80 u3v3
float80 u3v4
float80 u4v0
float80 u4v1
float80 u4v2
float80 u4v3
float80 u4v4

enter gaudry_pm_mult
input out_stack
input u_stack
input v_stack

caller eax
caller ebx
caller esi
caller edi
caller ebp

  u = u_stack
  v = v_stack


  out3 = *(float64 *) (u + 32)
  out3 *= *(float64 *) (v + 32)

  out3 *= *(float64 *) &gaudry_pm_scale

  u0v3 = *(float64 *) (u + 0)
  u0v3 *= *(float64 *) (v + 24)
  out3 += u0v3

  u1v2 = *(float64 *) (u + 8)
  u1v2 *= *(float64 *) (v + 16)
  out3 += u1v2

  u2v1 = *(float64 *) (u + 16)
  u2v1 *= *(float64 *) (v + 8)
  out3 += u2v1

  u3v0 = *(float64 *) (u + 24)
  u3v0 *= *(float64 *) (v + 0)
  out3 += u3v0



    out4 = *(float64 *) (u + 0)
    out4 *= *(float64 *) (v + 32)
  
    u1v3 = *(float64 *) (u + 8)
    u1v3 *= *(float64 *) (v + 24)
    out4 += u1v3
  
    u2v2 = *(float64 *) (u + 16)
    u2v2 *= *(float64 *) (v + 16)
    out4 += u2v2
  
    u3v1 = *(float64 *) (u + 24)
    u3v1 *= *(float64 *) (v + 8)
    out4 += u3v1
  
    u4v0 = *(float64 *) (u + 32)
    u4v0 *= *(float64 *) (v + 0)
    out4 += u4v0


            carry3 = *(float64 *) &gaudry_pm_alpha102
            carry3 += out3
            carry3 -= *(float64 *) &gaudry_pm_alpha102
          
            out4 += carry3
            out3 -= carry3


      out0 = *(float64 *) (u + 8)
      out0 *= *(float64 *) (v + 32)
    
      u2v3 = *(float64 *) (u + 16)
      u2v3 *= *(float64 *) (v + 24)
      out0 += u2v3
    
      u3v2 = *(float64 *) (u + 24)
      u3v2 *= *(float64 *) (v + 16)
      out0 += u3v2
    
      u4v1 = *(float64 *) (u + 32)
      u4v1 *= *(float64 *) (v + 8)
      out0 += u4v1


              carry4 = *(float64 *) &gaudry_pm_alpha127
              carry4 += out4
              carry4 -= *(float64 *) &gaudry_pm_alpha127
            
              out0 += carry4
              out4 -= carry4
            
            
              out0 *= *(float64 *) &gaudry_pm_scale
            
              u0v0 = *(float64 *) (u + 0)
              u0v0 *= *(float64 *) (v + 0)
              out0 += u0v0



        out1 = *(float64 *) (u + 16)
        out1 *= *(float64 *) (v + 32)
      
        u3v3 = *(float64 *) (u + 24)
        u3v3 *= *(float64 *) (v + 24)
        out1 += u3v3
      
        u4v2 = *(float64 *) (u + 32)
        u4v2 *= *(float64 *) (v + 16)
        out1 += u4v2
      
        out1 *= *(float64 *) &gaudry_pm_scale
      
        u0v1 = *(float64 *) (u + 0)
        u0v1 *= *(float64 *) (v + 8)
        out1 += u0v1
      
        u1v0 = *(float64 *) (u + 8)
        u1v0 *= *(float64 *) (v + 0)
        out1 += u1v0


                carry0 = *(float64 *) &gaudry_pm_alpha26
                carry0 += out0
                carry0 -= *(float64 *) &gaudry_pm_alpha26
              
                out1 += carry0
                out0 -= carry0



          out2 = *(float64 *) (u + 24)
          out2 *= *(float64 *) (v + 32)
        
          u4v3 = *(float64 *) (u + 32)
          u4v3 *= *(float64 *) (v + 24)
          out2 += u4v3
        
          out2 *= *(float64 *) &gaudry_pm_scale
        
          u0v2 = *(float64 *) (u + 0)
          u0v2 *= *(float64 *) (v + 16)
          out2 += u0v2
        
          u1v1 = *(float64 *) (u + 8)
          u1v1 *= *(float64 *) (v + 8)
          out2 += u1v1
        
          u2v0 = *(float64 *) (u + 16)
          u2v0 *= *(float64 *) (v + 0)
          out2 += u2v0


  out = out_stack


                  carry1 = *(float64 *) &gaudry_pm_alpha51
                  carry1 += out1
                  carry1 -= *(float64 *) &gaudry_pm_alpha51
                
                  out2 += carry1
                  out1 -= carry1


                    carry2 = *(float64 *) &gaudry_pm_alpha77
                    carry2 += out2
                    carry2 -= *(float64 *) &gaudry_pm_alpha77
                  
                    out3 += carry2
                    out2 -= carry2


                      carry3 = *(float64 *) &gaudry_pm_alpha102
                      carry3 += out3
                      carry3 -= *(float64 *) &gaudry_pm_alpha102
                    
                      out4 += carry3
                      out3 -= carry3

  *(float64 *) (out + 0) = out0
  *(float64 *) (out + 8) = out1
  *(float64 *) (out + 16) = out2
  *(float64 *) (out + 24) = out3
  *(float64 *) (out + 32) = out4



leave