/*
 * Copyright (c) 2000, NVIDIA CORPORATION.  All rights reserved.
 *
 * 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
 *
 *     http://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.
 *
 * Part of the f2008 bessel_yn 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)y0(src1);
                break;
            case 1:
                expct[order-n1] = (float)y1(src1);
                break;
            default:
                expct[order-n1] = (float)yn(order, src1);
                break;
            }
#else
            switch(order) {
            case 0:
                expct[order-n1] = (float)y0f(src1);
                break;
            case 1:
                expct[order-n1] = (float)y1f(src1);
                break;
            default:
                expct[order-n1] = (float)ynf(order, src1);
                break;
            }
#endif
            /*printf("%d) ynf(%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] = y0(src1);
                break;
            case 1:
                expct[order-n1] = y1(src1);
                break;
            default:
                expct[order-n1] = yn(order, src1);
                break;
            }
            /*printf("%d) yn(%e) = %e\n", i, src1, expct[i+order-n1]);*/
        }
}