NFFT  3.4.1
voronoi.c
1 /*
2  * Copyright (c) 2002, 2017 Jens Keiner, Stefan Kunis, Daniel Potts
3  *
4  * This program is free software; you can redistribute it and/or modify it under
5  * the terms of the GNU General Public License as published by the Free Software
6  * Foundation; either version 2 of the License, or (at your option) any later
7  * version.
8  *
9  * This program is distributed in the hope that it will be useful, but WITHOUT
10  * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
11  * FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
12  * details.
13  *
14  * You should have received a copy of the GNU General Public License along with
15  * this program; if not, write to the Free Software Foundation, Inc., 51
16  * Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
17  */
18 
19 #include "api.h"
20 #include "cstripack.h"
21 
23 void Y(voronoi_weights_1d)(R *w, R *x, const INT M)
24 {
25  INT j;
26 
27  w[0] = (x[1]-x[0])/K(2.0);
28 
29  for(j = 1; j < M-1; j++)
30  w[j] = (x[j+1]-x[j-1])/K(2.0);
31 
32  w[M-1] = (x[M-1]-x[M-2])/K(2.0);
33 }