00001 00005 /* 00006 * Author: Steven Ludtke, 04/10/2003 (sludtke@bcm.edu) 00007 * Copyright (c) 2000-2006 Baylor College of Medicine 00008 * 00009 * This software is issued under a joint BSD/GNU license. You may use the 00010 * source code in this file under either license. However, note that the 00011 * complete EMAN2 and SPARX software packages have some GPL dependencies, 00012 * so you are responsible for compliance with the licenses of these packages 00013 * if you opt to use BSD licensing. The warranty disclaimer below holds 00014 * in either instance. 00015 * 00016 * This complete copyright notice must be included in any revised version of the 00017 * source code. Additional authorship citations may be added, but existing 00018 * author citations must be preserved. 00019 * 00020 * This program is free software; you can redistribute it and/or modify 00021 * it under the terms of the GNU General Public License as published by 00022 * the Free Software Foundation; either version 2 of the License, or 00023 * (at your option) any later version. 00024 * 00025 * This program is distributed in the hope that it will be useful, 00026 * but WITHOUT ANY WARRANTY; without even the implied warranty of 00027 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 00028 * GNU General Public License for more details. 00029 * 00030 * You should have received a copy of the GNU General Public License 00031 * along with this program; if not, write to the Free Software 00032 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA 00033 * 00034 * */ 00035 00036 #ifndef eman__interp_h__ 00037 #define eman__interp_h__ 00038 00039 #include <cmath> 00040 00041 namespace EMAN 00042 { 00043 00047 class Interp 00048 { 00049 public: 00050 static float get_hyperg(int i) 00051 { 00052 return HYPERGEOM[i]; 00053 } 00054 00055 static float hyperg(float v) 00056 { 00057 if (v < 0 || v > 4.998f) { 00058 return 0; 00059 } 00060 float r = v / 0.001f; 00061 int a = (int) floor(r); 00062 r -= a; 00063 return (float)(-(HYPERGEOM[a] * (1.0 - r) + HYPERGEOM[a + 1] * r)); 00064 } 00065 00066 static float *get_gimx() 00067 { 00068 if (!gimx) { 00069 init_gimx(); 00070 } 00071 return gimx; 00072 } 00073 00074 private: 00075 static void init_gimx(); 00076 00077 private: 00078 static float HYPERGEOM[]; 00079 static float *gimx; 00080 }; 00081 } 00082 00083 #endif