horner계산을 이용한 테일러급수 계산 및 loop unrolling
void Eval_Poly_Naive(double y[], double x[], int n_x, double a[], int deg) { //p = c[0] + c[1] * x + c[2] * pow(x, 2.0) + c[3] * //pow(x, 3.0) + c[4] * pow(x, 4.0); for (int j = 0; j < deg; j++) { y[j] = a[0]; for (int i = 1; i < deg; i++) { y[j] += a[i] * pow(x[j], i); } printf("naive ret %lf\n", y[j]); } } void Eval_Poly_Horner(double y[], double x[], int n_x, double a[], int deg) { for (int j..
최근접 점의 쌍 찾기
#include #include # define SWAP(x, y, temp) ( (temp)=(x), (x)=(y), (y)=(temp) ) void merging(double* X, unsigned* Xid, unsigned* TMP, unsigned left, unsigned mid, unsigned right) { int i, j, k, l; i = left; j = mid + 1; k = left; /* 분할 정렬된 list의 합병 */ while (i