NFFT  3.3.1
vector2.c
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 }