EMAN2
interp.h
Go to the documentation of this file.
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