![]() |
NFFT
3.3.1
|
00001 /* 00002 * Copyright (c) 2002, 2016 Jens Keiner, Stefan Kunis, Daniel Potts 00003 * 00004 * This program is free software; you can redistribute it and/or modify it under 00005 * the terms of the GNU General Public License as published by the Free Software 00006 * Foundation; either version 2 of the License, or (at your option) any later 00007 * version. 00008 * 00009 * This program is distributed in the hope that it will be useful, but WITHOUT 00010 * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS 00011 * FOR A PARTICULAR PURPOSE. See the GNU General Public License for more 00012 * details. 00013 * 00014 * You should have received a copy of the GNU General Public License along with 00015 * this program; if not, write to the Free Software Foundation, Inc., 51 00016 * Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. 00017 */ 00018 00019 #include "infft.h" 00020 00022 void Y(cp_complex)(C *x, C *y, INT n) 00023 { 00024 INT k; 00025 00026 for (k = 0; k < n; k++) 00027 x[k] = y[k]; 00028 } 00029 00031 void Y(cp_double)(R *x, R *y, INT n) 00032 { 00033 INT k; 00034 00035 for (k = 0; k < n; k++) 00036 x[k] = y[k]; 00037 } 00038 00040 void Y(cp_a_complex)(C *x, R a, C *y, INT n) 00041 { 00042 INT k; 00043 00044 for (k = 0; k < n; k++) 00045 x[k] = a * y[k]; 00046 } 00047 00049 void Y(cp_a_double)(R *x, R a, R *y, INT n) 00050 { 00051 INT k; 00052 00053 for (k = 0; k < n; k++) 00054 x[k] = a * y[k]; 00055 } 00056 00057 00059 void Y(cp_w_complex)(C *x, R *w, C *y, INT n) 00060 { 00061 INT k; 00062 00063 for (k = 0; k < n; k++) 00064 x[k] = w[k]*y[k]; 00065 } 00066 00068 void Y(cp_w_double)(R *x, R *w, R *y, INT n) 00069 { 00070 INT k; 00071 00072 for (k = 0; k < n; k++) 00073 x[k] = w[k] * y[k]; 00074 }