/*
 * Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
 * See https://llvm.org/LICENSE.txt for license information.
 * SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
 *
 * Part of the f2008 bessel_jn tranformational intrinsic test
 */

#include <stdio.h>
#include <math.h>

void
get_expected_f(float src1, float expct[], int n1, int n2, int n)
{
    int i;
    int order;

        for(order = n1; order <= n2; order++ ) {
#if  defined( _WIN64) || defined(_WIN32) || defined(osx86)
            switch(order) {
            case 0:
                expct[order-n1] = (float)j0(src1);
                break;
            case 1:
                expct[order-n1] = (float)j1(src1);
                break;
            default:
                expct[order-n1] = (float)jn(order, src1);
                break;
            }
#else
            switch(order) {
            case 0:
                expct[order-n1] = (float)j0f(src1);
                break;
            case 1:
                expct[order-n1] = (float)j1f(src1);
                break;
            default:
                expct[order-n1] = (float)jnf(order, src1);
                break;
            }
#endif
            /*printf("%d) jnf(%e) = %e\n", i, src1, expct[i+order-n1]);*/
        }
}

void
get_expected_d(double src1, double expct[], int n1, int n2, int n)
{
    int i;
    int order;

        for(order = n1; order <= n2; order++ ) {
            switch(order) {
            case 0:
                expct[order-n1] = j0(src1);
                break;
            case 1:
                expct[order-n1] = j1(src1);
                break;
            default:
                expct[order-n1] = jn(order, src1);
                break;
            }
            /*printf("%d) jn(%e) = %e\n", i, src1, expct[i+order-n1]);*/
        }
}