git.gag.com Git - debian/gnuradio/blob - gnuradio-core/src/lib/filter/3dnow_float_dotprod_really_simple.S

   1 #
   2 # Copyright 2002 Free Software Foundation, Inc.
   3 #
   4 # This file is part of GNU Radio
   5 #
   6 # GNU Radio is free software; you can redistribute it and/or modify
   7 # it under the terms of the GNU General Public License as published by
   8 # the Free Software Foundation; either version 3, or (at your option)
   9 # any later version.
  10 #
  11 # GNU Radio is distributed in the hope that it will be useful,
  12 # but WITHOUT ANY WARRANTY; without even the implied warranty of
  13 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  14 # GNU General Public License for more details.
  15 #
  16 # You should have received a copy of the GNU General Public License
  17 # along with GNU Radio; see the file COPYING.  If not, write to
  18 # the Free Software Foundation, Inc., 51 Franklin Street,
  19 # Boston, MA 02110-1301, USA.
  20 #
  21
  22
  23 # input and taps are guarenteed to be 16 byte aligned.
  24 # n_4_float_blocks is != 0
  25 #
  26 #
  27 #  float
  28 #  sse_float_dotprod (const float *input,
  29 #                  const float *taps, unsigned n_4_float_blocks)
  30 #  {
  31 #    float sum0 = 0;
  32 #    float sum1 = 0;
  33 #    float sum2 = 0;
  34 #    float sum3 = 0;
  35 #
  36 #    do {
  37 #
  38 #      sum0 += input[0] * taps[0];
  39 #      sum1 += input[1] * taps[1];
  40 #      sum2 += input[2] * taps[2];
  41 #      sum3 += input[3] * taps[3];
  42 #
  43 #      input += 4;
  44 #      taps += 4;
  45 #
  46 #    } while (--n_4_float_blocks != 0);
  47 #
  48 #
  49 #    return sum0 + sum1 + sum2 + sum3;
  50 #  }
  51 #
  52
  53
  54         .file   "3dnow_float_dotprod_really_simple.s"
  55         .version        "01.01"
  56 .text
  57         .p2align 4
  58 .globl sse_float_dotprod
  59         .type    sse_float_dotprod,@function
  60 sse_float_dotprod:
  61         pushl   %ebp
  62         movl    %esp, %ebp
  63         movl    8(%ebp), %edx
  64         movl    12(%ebp), %eax
  65         movl    16(%ebp), %ecx
  66
  67
  68         # The plan is to get it computing the correct answer, and
  69         # then to unroll and schedule the inner loop.
  70
  71         pxor    %mm4, %mm4              # mm4 = 0 0
  72         shll    $1, %ecx                # count * 2
  73
  74         .p2align 4
  75 .Loop1:
  76         movq    (%eax), %mm0
  77         pfmul   (%edx), %mm0
  78         pfadd   %mm0, %mm4
  79         addl    $8, %edx
  80         addl    $8, %eax
  81         decl    %ecx
  82         jne     .Loop1
  83
  84         # at this point mm4 contains partial sums
  85
  86         pfacc   %mm4, %mm4
  87         movd    %mm4, 16(%ebp)
  88         femms
  89         flds    16(%ebp)
  90
  91         popl    %ebp
  92         ret
  93 .Lfe1:
  94         .size    sse_float_dotprod,.Lfe1-sse_float_dotprod
  95         .ident  "Hand coded x86 3DNow! assembly"
  96
  97 #if defined(__linux__) && defined(__ELF__)
  98 .section .note.GNU-stack,"",%progbits
  99 #endif