1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
|
classdef nfsftTestcaseInitDelegate
%NFSFTTESTCASEDELEGATEINIT Summary of this class goes here
% Detailed explanation goes here
properties
name = '';
txt = '';
flags = 0;
kappa = 1000;
fpt_flags = 0;
nfft_flags = [];
nfft_cutoff_m = [];
end
methods
function h = nfsftTestcaseInitDelegate(name, flags, fpt_flags, nfft_flags, nfft_cutoff_m)
h.name = name;
h.txt = h.name;
if exist('flags','var') && ~isempty(flags)
h.flags = flags;
if bitand(flags,NFSFT_USE_DPT) ~= 0
h.txt = sprintf('%s DPT', h.name);
end
end
if exist('fpt_flags','var') && ~isempty(fpt_flags)
h.fpt_flags = fpt_flags;
end
if exist('nfft_flags','var') && ~isempty(nfft_flags)
h.nfft_flags = nfft_flags;
end
if exist('nfft_cutoff_m','var') && ~isempty(nfft_cutoff_m)
h.nfft_cutoff_m = nfft_cutoff_m;
end
end
function nfsft_precompute(h, N)
if ~isempty(h.fpt_flags)
nfsft_precompute(N,h.kappa,h.flags,h.fpt_flags);
elseif ~isempty(h.flags)
nfsft_precompute(N,h.kappa,h.flags);
elseif ~isempty(h.kappa)
nfsft_precompute(N,h.kappa);
else
nfsft_precompute(N);
end
end
function [h, plan] = init(h, N, M)
switch h.name
case 'init'
h.nfsft_precompute(N);
plan = nfsft_init(N,M);
case 'init_advanced'
h.nfsft_precompute(N);
plan = nfsft_init_advanced(N,M,h.flags);
case 'init_guru'
h.nfsft_precompute(N);
plan = nfsft_init_guru(N,M,h.flags,h.nfft_flags,h.nfft_cutoff_m);
case 'init_class'
if isempty(h.flags) && isempty(h.nfft_cutoff_m) && isempty(h.fpt_flags) && isempty(h.nfft_flags)
plan = nfsft(N,M);
elseif ~isempty(h.flags)
plan = nfsft(N,M,h.flags);
else
plan = nfsft(N,M,h.flags,h.kappa,h.nfft_cutoff_m,h.fpt_flags,h.nfft_flags);
end
otherwise
error('Unknown init name: %s', h.name);
end
end
end
end
|
