00001 #ifndef UTIL_H 00002 #define UTIL_H 00003 00004 #include <complex.h> 00005 #include "../../../include/util.h" 00006 00014 inline int ngpt(int n); 00015 00016 #define SWAPCT(x,y,temp) {temp=(x); (x)=(y); (y)=temp;} 00017 00018 void vpr_c_hat (complex **x, int n, char * text); 00019 00020 double error_complex_inf_old(complex *x0, complex *x, int n); 00021 double error_complex_inf_r_old(complex *x0, complex *x, int n); 00022 double error_complex_1_old(complex *a, complex *b, int n); 00023 double error_complex_2_old(complex *a, complex *b, int n); 00024 double error_complex3_old(complex *a, complex *b, int n); 00025 double err_f_hat_infty_old(complex **f_hat, complex **f_hat2, int M); 00026 double err_f_hat_1_old(complex **f_hat, complex **f_hat2, int M); 00027 double err_f_hat_2_old(complex **f_hat, complex **f_hat2, int M); 00028 double norm_complex_inf_old(complex *x, int n); 00029 double norm_complex_1_old(complex *x, int n); 00030 double norm_complex_2_old(complex *x, int n); 00031 double norm_f_hat_1_old(complex **f_hat, int M); 00032 double norm_nfft_1_old(complex *f_hat, int M); 00033 double norm_f_hat_1s_old(complex **f_hat, int M); 00034 int min_old(const int a, const int b); 00035 int max_old(const int a, const int b); 00036 long int max_oldl(const long int a, const long int b); 00037 00038 void copyc_hat_old(complex **x,complex **y, int M); 00039 void copyc_w_hat_old(complex **x, double **w, complex **y, int M); 00040 double dotproductc_w_hat_old(complex **x, double **w, int M); 00041 double dotproductc_hat_old(complex** x, int M); 00042 void updatec_xpay_hat_old(complex **x,double a, complex **y, int M); 00043 void updatec_xpawy_hat_old(complex **x, double a, double **w, complex **y, int M); 00044 void updatec_axpy_hat_old(complex **x,double a, complex **y, int M); 00045 void updatec_axpby_hat_old(complex **x, double a, complex **y, double b, int M); 00046 void updatec_axpy_2_old(complex* x,double a, complex* y, int n); 00047 void updatec_xpay_2_old(complex* x,double a, complex* y, int n); 00048 void updatec_axpby_2(complex* x, double a, complex* y, double b, int n); 00049 void updatec_xpawy_2_old(complex* x,double a, double* w, complex* y, int n); 00050 void uvxpwy_old(complex* u, complex* x, double* v, complex* y, double* w, int n); 00051 void auvxpwy_old(double a, complex* u, complex* x, double* v, complex* y, double* w, int n); 00052 void abuvxpwy_old(double a, double b, complex* u, complex* x, double* v, complex* y, double* w, int n); 00053 void normalize_f_hat_old(complex **f_hat, int M); 00054 int abs(const int a); 00055 #endif