1 /*

     2  * LIBOIL - Library of Optimized Inner Loops

     3  * Copyright (c) 2003,2004 David A. Schleef <ds@schleef.org>

     4  * All rights reserved.

     5  *

     6  * Redistribution and use in source and binary forms, with or without

     7  * modification, are permitted provided that the following conditions

     8  * are met:

     9  * 1. Redistributions of source code must retain the above copyright

    10  *    notice, this list of conditions and the following disclaimer.

    11  * 2. Redistributions in binary form must reproduce the above copyright

    12  *    notice, this list of conditions and the following disclaimer in the

    13  *    documentation and/or other materials provided with the distribution.

    14  * 

    15  * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR

    16  * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED

    17  * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE

    18  * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT,

    19  * INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES

    20  * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR

    21  * SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)

    22  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT,

    23  * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING

    24  * IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE

    25  * POSSIBILITY OF SUCH DAMAGE.

    26  */

    27 //Portions Copyright (c)  2008-2009 Nokia Corporation and/or its subsidiary(-ies). All rights reserved. 

    28 

    29 #ifdef HAVE_CONFIG_H

    30 #include "config.h"

    31 #endif

    32 

    33 #include <liboil/liboilfunction.h>

    34 #include "liboil/simdpack/simdpack.h"

    35 

    36 #define ABS(x) ((x)>0 ? (x) : -(x))

    37 

    38 static void

    39 abs_u16_s16_unroll2 (uint16_t * dest, int dstr, int16_t * src, int sstr, int n)

    40 {

    41   while (n & 1) {

    42     *dest = ABS (*src);

    43     OIL_INCREMENT (dest, dstr);

    44     OIL_INCREMENT (src, sstr);

    45     n--;

    46   }

    47   while (n > 0) {

    48     *dest = ABS (*src);

    49     OIL_INCREMENT (dest, dstr);

    50     OIL_INCREMENT (src, sstr);

    51     *dest = ABS (*src);

    52     OIL_INCREMENT (dest, dstr);

    53     OIL_INCREMENT (src, sstr);

    54     n -= 2;

    55   }

    56 }

    57 OIL_DEFINE_IMPL (abs_u16_s16_unroll2, abs_u16_s16);

    58 

    59 static void

    60 abs_u16_s16_unroll4 (uint16_t * dest, int dstr, int16_t * src, int sstr, int n)

    61 {

    62   while (n & 3) {

    63     *dest = ABS (*src);

    64     OIL_INCREMENT (dest, dstr);

    65     OIL_INCREMENT (src, sstr);

    66     n--;

    67   }

    68   while (n > 0) {

    69     *dest = ABS (*src);

    70     OIL_INCREMENT (dest, dstr);

    71     OIL_INCREMENT (src, sstr);

    72     *dest = ABS (*src);

    73     OIL_INCREMENT (dest, dstr);

    74     OIL_INCREMENT (src, sstr);

    75     *dest = ABS (*src);

    76     OIL_INCREMENT (dest, dstr);

    77     OIL_INCREMENT (src, sstr);

    78     *dest = ABS (*src);

    79     OIL_INCREMENT (dest, dstr);

    80     OIL_INCREMENT (src, sstr);

    81     n -= 4;

    82   }

    83 }

    84 OIL_DEFINE_IMPL (abs_u16_s16_unroll4, abs_u16_s16);

    85 

    86 static void

    87 abs_u16_s16_fast (uint16_t * dest, int dstr, int16_t * src, int sstr, int n)

    88 {

    89   int i;

    90   int16_t x, y;

    91 

    92   for (i = 0; i < n; i++) {

    93     x = *src;

    94     y = ((x >> 15) & x);

    95     *dest = x - y - y;

    96     OIL_INCREMENT (dest, dstr);

    97     OIL_INCREMENT (src, sstr);

    98   }

    99 }

   100 

   101 OIL_DEFINE_IMPL (abs_u16_s16_fast, abs_u16_s16);

   102 

   103 

   104 /* Similar algorithm as clipping, but for absolute value.  Note the stunt

   105  * pulled to get around 2*y rolling over.  */

   106 

   107 static void

   108 abs_u32_s32_fast (uint32_t * dest, int dstr, int32_t * src, int sstr, int n)

   109 {

   110   int i;

   111   int32_t x, y;

   112 

   113   for (i = 0; i < n; i++) {

   114     x = OIL_GET(src, i * sstr, int32_t);

   115     y = ((x >> 31) & x);

   116     OIL_GET(dest, i * dstr, uint32_t) = x - y - y;

   117   }

   118 }

   119 

   120 OIL_DEFINE_IMPL (abs_u32_s32_fast, abs_u32_s32);

   121 

   122 

   123 

   124 #ifdef	__SYMBIAN32__

   125  

   126 OilFunctionImpl* __oil_function_impl_abs_u16_s16_unroll2() {

   127 		return &_oil_function_impl_abs_u16_s16_unroll2;

   128 }

   129 #endif

   130 

   131 #ifdef	__SYMBIAN32__

   132  

   133 OilFunctionImpl* __oil_function_impl_abs_u16_s16_unroll4() {

   134 		return &_oil_function_impl_abs_u16_s16_unroll4;

   135 }

   136 #endif

   137 

   138 #ifdef	__SYMBIAN32__

   139  

   140 OilFunctionImpl* __oil_function_impl_abs_u16_s16_fast() {

   141 		return &_oil_function_impl_abs_u16_s16_fast;

   142 }

   143 #endif

   144 

   145 #ifdef	__SYMBIAN32__

   146  

   147 OilFunctionImpl* __oil_function_impl_abs_u32_s32_fast() {

   148 		return &_oil_function_impl_abs_u32_s32_fast;

   149 }

   150 #endif

   151