projects
/
ctsim.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
r386: no message
[ctsim.git]
/
include
/
filter.h
diff --git
a/include/filter.h
b/include/filter.h
index ca6e67b0314e405d91bbfbf8711bd2353a37f78e..7b696b29e8553fd89e99d2b7dd67fa53a72a1a47 100644
(file)
--- a/
include/filter.h
+++ b/
include/filter.h
@@
-9,7
+9,7
@@
** This is part of the CTSim program
** Copyright (C) 1983-2000 Kevin Rosenberg
**
** This is part of the CTSim program
** Copyright (C) 1983-2000 Kevin Rosenberg
**
-** $Id: filter.h,v 1.2
2 2000/12/16 06:12:4
7 kevin Exp $
+** $Id: filter.h,v 1.2
3 2001/01/12 21:53:2
7 kevin Exp $
**
** This program is free software; you can redistribute it and/or modify
** it under the terms of the GNU General Public License (version 2) as
**
** This program is free software; you can redistribute it and/or modify
** it under the terms of the GNU General Public License (version 2) as
@@
-116,6
+116,12
@@
class SignalFilter {
static const char* convertDomainIDToName (const int idDomain);
static const char* convertDomainIDToTitle (const int idDomain);
static const char* convertDomainIDToName (const int idDomain);
static const char* convertDomainIDToTitle (const int idDomain);
+ static double sinc (double x)
+ { return (fabs(x) > F_EPSILON ? (sin (x) / x) : 1.0); }
+
+ static double sinc (double x, double mult)
+ { return (fabs(x) > F_EPSILON ? (sin (x * mult) / x) : 1.0); }
+
private:
int m_nFilterPoints;
double m_dBandwidth;
private:
int m_nFilterPoints;
double m_dBandwidth;
@@
-152,8
+158,6
@@
class SignalFilter {
double spatialResponseAnalytic (double x) const;
double frequencyResponse (double u) const;
double spatialResponseAnalytic (double x) const;
double frequencyResponse (double u) const;
- static double sinc (double x, double mult)
- { return (fabs(x) > F_EPSILON ? (sin (x * mult) / x) : 1.0); }
static double integral_abscos (double u, double w)
{ return (fabs (u) > F_EPSILON ? (cos (u * w) - 1) / (u * u) + w / u * sin (u * w) : (w * w / 2)); }
};
static double integral_abscos (double u, double w)
{ return (fabs (u) > F_EPSILON ? (cos (u * w) - 1) / (u * u) + w / u * sin (u * w) : (w * w / 2)); }
};