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
|
function tsgLoadHCoefficients(lGrid, mValues)
%
% tsgLoadHCoefficients(lGrid, mValues)
%
% loads the values of the target function at the needed points
%
% INPUT:
%
% lGrid: a grid list created by tsgMake***(...) command
%
% mValues: a matrix with dimension [num_needed_points, iOut]
% each row corresponds to the values of the hierarchical
% coefficients at the corresponding needed point.
% The order and leading dimension must match the points
% obtained form tsgGetNeededPoints(...) command
%
% OUTPUT:
%
% The grid file associated with lGrid is modified
%
[sFiles, sTasGrid] = tsgGetPaths();
[sFileG, sFileX, sFileV, sFileO, sFileW, sFileC] = tsgMakeFilenames(lGrid);
sCommand = [sTasGrid,' -setcoefficients'];
sCommand = [sCommand, ' -gridfile ', sFileG];
if (strcmp(lGrid.sType, 'fourier'))
mRealMat = zeros(size(mValues, 1), 2 * size(mValues, 2));
mRealMat(:, 1:2:end) = real(mValues);
mRealMat(:, 2:2:end) = imag(mValues);
mValues = mRealMat;
end
tsgWriteMatrix(sFileV, mValues);
lClean.sFileV = 1;
sCommand = [sCommand, ' -vf ', sFileV];
[status, cmdout] = system(sCommand);
if (max(size(strfind(cmdout, 'ERROR'))) ~= 0)
disp(cmdout);
error('The tasgrid execurable returned an error, see above');
return;
else
if (~ isempty(cmdout))
fprintf(1,['Warning: Command had non-empty output:\n']);
disp(cmdout);
end
end
tsgCleanTempFiles(lGrid, lClean);
end
|
