Commit 9bacca3d authored by thomas.forbriger's avatar thomas.forbriger Committed by thomas.forbriger
Browse files

waterlevel in omega

This is a legacy commit from before 2015-03-01.
It may be incomplete as well as inconsistent.
See COPYING.legacy and README.history for details.


SVN Path:     http://gpitrsvn.gpi.uni-karlsruhe.de/repos/TFSoftware/trunk
SVN Revision: 1036
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent becb5174
# this is <Makefile>
# ----------------------------------------------------------------------------
# $Id: Makefile,v 1.2 2002-11-11 15:27:00 forbrig Exp $
# $Id: Makefile,v 1.3 2002-11-13 18:51:51 forbrig Exp $
#
# Copyright (c) 2002 by Thomas Forbriger (IMG Frankfurt)
#
......@@ -14,7 +14,7 @@
# ============================================================================
#
all:
all: libfourier.a libfourier.doc
flist: Makefile $(wildcard *.f *.inc)
echo $^ | tr ' ' '\n' | sort > $@
......
c this is <polesnzeroes.f>
c ----------------------------------------------------------------------------
c ($Id: polesnzeroes.f,v 1.3 2002-11-12 08:24:35 forbrig Exp $)
c ($Id: polesnzeroes.f,v 1.4 2002-11-13 18:51:51 forbrig Exp $)
c
c Copyright (c) 2002 by Thomas Forbriger (IMG Frankfurt)
c
......@@ -8,6 +8,7 @@ c pole and zero access subroutines
c
c REVISIONS and CHANGES
c 11/11/2002 V1.0 Thomas Forbriger
c 13/11/2002 V1.1 introduce waterlevel to omega
c
cS
c ============================================================================
......@@ -17,7 +18,8 @@ c
c print library code revision
c
cE
print *,'$Id: polesnzeroes.f,v 1.3 2002-11-12 08:24:35 forbrig Exp $'
print *,
& '$Id: polesnzeroes.f,v 1.4 2002-11-13 18:51:51 forbrig Exp $'
return
end
c
......@@ -229,17 +231,20 @@ c
cS
c----------------------------------------------------------------------
c
double complex function fou_eval(omega)
double complex function fou_eval(om)
c
c evaluate filter response coefficient at real frequency omega
c evaluate filter response coefficient at real frequency om
c the function applies a waterlevel of 1.d-30 to om
c
double precision omega
double precision om
c
include 'polesnzeros.inc'
c
cE
double complex value, thenumer, thedenom
double complex value, thenumer, thedenom, omega
integer i
c
omega=max(1.d-30,om)
c
thenumer=(1.d0,0.d0)
do i=1,fourier_nzeros
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment