EMAN2
interp.h
Go to the documentation of this file.
1/*
2 * Author: Steven Ludtke, 04/10/2003 (sludtke@bcm.edu)
3 * Copyright (c) 2000-2006 Baylor College of Medicine
4 *
5 * This software is issued under a joint BSD/GNU license. You may use the
6 * source code in this file under either license. However, note that the
7 * complete EMAN2 and SPARX software packages have some GPL dependencies,
8 * so you are responsible for compliance with the licenses of these packages
9 * if you opt to use BSD licensing. The warranty disclaimer below holds
10 * in either instance.
11 *
12 * This complete copyright notice must be included in any revised version of the
13 * source code. Additional authorship citations may be added, but existing
14 * author citations must be preserved.
15 *
16 * This program is free software; you can redistribute it and/or modify
17 * it under the terms of the GNU General Public License as published by
18 * the Free Software Foundation; either version 2 of the License, or
19 * (at your option) any later version.
20 *
21 * This program is distributed in the hope that it will be useful,
22 * but WITHOUT ANY WARRANTY; without even the implied warranty of
23 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
24 * GNU General Public License for more details.
25 *
26 * You should have received a copy of the GNU General Public License
27 * along with this program; if not, write to the Free Software
28 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
29 *
30 * */
31
32#ifndef eman__interp_h__
33#define eman__interp_h__
34
35#include <cmath>
36
37namespace EMAN
38{
39
43 class Interp
44 {
45 public:
46 static float get_hyperg(int i)
47 {
48 return HYPERGEOM[i];
49 }
50
51 static float hyperg(float v)
52 {
53 if (v < 0 || v > 4.998f) {
54 return 0;
55 }
56 float r = v / 0.001f;
57 int a = (int) floor(r);
58 r -= a;
59 return (float)(-(HYPERGEOM[a] * (1.0 - r) + HYPERGEOM[a + 1] * r));
60 }
61
62 static float *get_gimx()
63 {
64 if (!gimx) {
65 init_gimx();
66 }
67 return gimx;
68 }
69
70 private:
71 static void init_gimx();
72
73 private:
74 static float HYPERGEOM[];
75 static float *gimx;
76 };
77}
78
79#endif
Interp defines the interpolation function used to generate a e^-x^4 function in real space.
Definition: interp.h:44
static float hyperg(float v)
Definition: interp.h:51
static float * gimx
Definition: interp.h:75
static float HYPERGEOM[]
Definition: interp.h:74
static void init_gimx()
Definition: interp.cpp:1711
static float * get_gimx()
Definition: interp.h:62
static float get_hyperg(int i)
Definition: interp.h:46
E2Exception class.
Definition: aligner.h:40