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
75
|
function writeTestcase(file,usenfsft,usenfft,cutoff,usefpt,threshold,...
bandwidth,theta,phi,f)
%WRITETESTCASE - Write an iterS2 testcase definition to a file
% WRITETESTCASE(FILE, USENFSFT, USENFFT, CUTOFF, USEFPT, THRESHOLD,
% BANDWIDTH, THETA, PHI, F)
%
% Copyright (c) 2002, 2017 Jens Keiner, Stefan Kunis, Daniel Potts
% Copyright (c) 2002, 2017 Jens Keiner, Stefan Kunis, Daniel Potts
%
% This program is free software; you can redistribute it and/or modify it under
% the terms of the GNU General Public License as published by the Free Software
% Foundation; either version 2 of the License, or (at your option) any later
% version.
%
% This program is distributed in the hope that it will be useful, but WITHOUT
% ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
% FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
% details.
%
% You should have received a copy of the GNU General Public License along with
% this program; if not, write to the Free Software Foundation, Inc., 51
% Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
%
% Write NFSFT usage flag.
fprintf(file,'nfsft=%d\n',usenfsft);
if (usenfsft == true)
% Write NFFT usage flag.
fprintf(file,'nfft=%d\n',usenfft);
if (usenfft == true)
% Write NFFT cut-off parameter.
fprintf(file,'cutoff=%d\n',cutoff);
end
% Write FPT usage flag.
fprintf(file,'fpt=%d\n',usefpt);
if (usefpt == true)
% Write FPT threshold.
fprintf(file,'threshold=%e\n',threshold);
end
end
% Write bandwidth
fprintf(file,'bandwidth=%d\n',bandwidth);
% Write number of nodes.
fprintf(file,'nodes=%d\n',length(theta));
% Write nodes and function values.
for j=1:length(theta)
% Write node j and corresponding function value.
fprintf(file,'%f %f %f %f\n',theta(j),phi(j),real(f(j)),imag(f(j)));
end
% Write number of nodes.
fprintf(file,'nodes_eval=%d\n',m*n);
% Write nodes and function values.
for j=1:length(theta)
for k=1:length(phi)
fprintf(file,'%f %f\n',theta(j),phi(k));
end
end
% End of function
return;
|
