#include <util.h>
Inheritance diagram for EMAN::Util::FakeKaiserBessel:
Public Member Functions | |
FakeKaiserBessel (float alpha, int K, float r_, float v_, int N_, float vtable_=0.f, int ntable_=5999) | |
float | sinhwin (float x) const |
Kaiser-Bessel Sinh window function. | |
float | i0win (float x) const |
Kaiser-Bessel I0 window function. | |
void | build_I0table () |
2*pi*alpha*r*vadjust |
Definition at line 386 of file util.h.
EMAN::Util::FakeKaiserBessel::FakeKaiserBessel | ( | float | alpha, | |
int | K, | |||
float | r_, | |||
float | v_, | |||
int | N_, | |||
float | vtable_ = 0.f , |
|||
int | ntable_ = 5999 | |||
) | [inline] |
void Util::FakeKaiserBessel::build_I0table | ( | ) | [virtual] |
2*pi*alpha*r*vadjust
Reimplemented from EMAN::Util::KaiserBessel.
Definition at line 2068 of file util_sparx.cpp.
References EMAN::Util::KaiserBessel::facadj, EMAN::Util::KaiserBessel::fltb, EMAN::Util::KaiserBessel::i0table, EMAN::Util::KaiserBessel::K, EMAN::Util::KaiserBessel::N, EMAN::Util::KaiserBessel::ntable, EMAN::Util::round(), sqrt(), and EMAN::Util::KaiserBessel::vadjust.
02068 { 02069 i0table.resize(ntable+1); // i0table[0:ntable] 02070 int ltab = int(round(float(ntable)/1.1f)); 02071 fltb = float(ltab)/(K/2); 02072 float val0 = sqrt(facadj)*static_cast<float>(gsl_sf_bessel_I1(facadj)); 02073 for (int i=ltab+1; i <= ntable; i++) i0table[i] = 0.f; 02074 for (int i=0; i <= ltab; i++) { 02075 float s = float(i)/fltb/N; 02076 if (s < vadjust) { 02077 float rt = sqrt(1.f - pow(s/vadjust, 2)); 02078 i0table[i] = sqrt(facadj*rt)*static_cast<float>(gsl_sf_bessel_I1(facadj*rt))/val0; 02079 } else { 02080 i0table[i] = 0.f; 02081 } 02082 } 02083 }
float Util::FakeKaiserBessel::i0win | ( | float | x | ) | const [virtual] |
Kaiser-Bessel I0 window function.
Reimplemented from EMAN::Util::KaiserBessel.
Definition at line 2059 of file util_sparx.cpp.
References EMAN::Util::KaiserBessel::facadj, sqrt(), and EMAN::Util::KaiserBessel::vadjust.
02059 { 02060 float val0 = sqrt(facadj)*float(gsl_sf_bessel_I1(facadj)); 02061 float absx = fabs(x); 02062 if (absx > vadjust) return 0.f; 02063 float rt = sqrt(1.f - pow(absx/vadjust, 2)); 02064 float res = sqrt(facadj*rt)*float(gsl_sf_bessel_I1(facadj*rt))/val0; 02065 return res; 02066 }
float Util::FakeKaiserBessel::sinhwin | ( | float | x | ) | const [virtual] |
Kaiser-Bessel Sinh window function.
Reimplemented from EMAN::Util::KaiserBessel.
Definition at line 2085 of file util_sparx.cpp.
References EMAN::Util::KaiserBessel::alphar, EMAN::Util::KaiserBessel::fac, and sqrt().
02085 { 02086 float val0 = sinh(fac)/fac; 02087 float absx = fabs(x); 02088 if (0.0 == x) { 02089 float res = 1.0f; 02090 return res; 02091 } else if (absx == alphar) { 02092 return 1.0f/val0; 02093 } else if (absx < alphar) { 02094 float rt = sqrt(1.0f - pow((x/alphar), 2)); 02095 float facrt = fac*rt; 02096 float res = (sinh(facrt)/facrt)/val0; 02097 return res; 02098 } else { 02099 float rt = sqrt(pow((x/alphar),2) - 1.f); 02100 float facrt = fac*rt; 02101 float res = (sin(facrt)/facrt)/val0; 02102 return res; 02103 } 02104 }