polynomial.h

Go to the documentation of this file.

/*
 * Copyright (c) 2014 The University of Manchester
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 *     https://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 */
 
#ifndef __POLYNOMIAL_H__
#define __POLYNOMIAL_H__
 
#include "arm_acle.h"
#include "assert.h"
 
static inline int __horner_int_b(
    int *a,
    int x,
    int n)
{
#ifdef __ARM_FEATURE_DSP
    register int r = *a++;
 
    assert(!__reset_and_saturation_occurred());
    for ( ; n > 0; n--) {
        r = __smlawb(r, x, *a++);
    }
    assert(!__saturation_occurred());
    return r;
#else  
    register int64_t t = *a++;
    register int64_t dx = (int64_t) ((int16_t) (x & 0xFFFF));
 
    for ( ; n > 0; n--) {
        t = (t * dx >> 16) + *a++;
    }
    return (int) (t & 0xFFFFFFFF);
#endif /*__ARM_FEATURE_DSP*/
}
 
static inline int __horner_int_t(
    int *a,
    int x,
    int n)
{
#ifdef __ARM_FEATURE_DSP
    register int r = *a++;
 
    assert(!__reset_and_saturation_occurred());
    for ( ; n > 0; n--) {
        r = __smlawt(r, x, *a++);
    }
    assert(!__saturation_occurred());
    return r;
#else  
    register int64_t t = *a++;
    register int64_t dx = (int64_t) (x >> 16);
 
    for ( ; n > 0; n--) {
        t = (t * dx >> 16) + *a++;
    }
    return (int) (t & 0xFFFFFFFF);
#endif /*__ARM_FEATURE_DSP*/
}
 
#endif /*__POLYNOMIAL_H__*/