Homec4science
Diffusion SPClibs (feature/fft-split) Browse

spclibs/matlab/pde1dp_cmpl_dft.m
93a1e9712d04feature/fft-split

Last Change
Raw File

pde1dp_cmpl_dft.m
View Options

%
% @file pde1dp_cmpl_dft.m
%
% @brief
%
% @copyright
% Copyright (©) 2021 EPFL (Ecole Polytechnique Fédérale de Lausanne)
% SPC (Swiss Plasma Center)
%
% SPClibs is free software: you can redistribute it and/or modify it under
% the terms of the GNU Lesser General Public License as published by the Free
% Software Foundation, either version 3 of the License, or (at your option)
% any later version.
%
% SPClibs 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 Lesser General Public License
% along with this program. If not, see <https://www.gnu.org/licenses/>.
%
% @authors
% (in alphabetical order)
% @author Trach-Minh Tran <trach-minh.tran@epfl.ch>
%
file='pde1dp_cmpl_dft.h5';
xgrid=hdf5read(file,'/xgrid');
nx = size(xgrid)-1;
mode=0:nx-1;
rhs=h5Complex(file, '/rhs');
sol=h5Complex(file, '/sol');
rhs_fft=h5Complex(file, '/rhs_fft');
sol_fft=h5Complex(file, '/sol_fft');
mat=h5Complex(file, '/mat');
x=hdf5read(file,'/x');
solana=h5Complex(file,'/solana');
solcal=h5Complex(file,'/solcal');
err=hdf5read(file, '/err');
figure
subplot(211)
plot(x, real(solana),x,imag(solana),x,real(solcal),'o',x, ...
imag(solcal), '*')
legend('Exact Real', 'Exact Imag', 'Calc. Real', 'Calc. Imag')
xlabel('X'); ylabel('SOL');
subplot(212)
plot(x, err, 'o-');
xlabel('X'); ylabel('|Error|')
grid on
figure
subplot(311)
plot(mode, real(mat), 'o', mode, imag(mat), '*')
xlabel('mode'); ylabel('DFT of MAT')
legend('Real', 'Imag')
grid on
subplot(312)
plot(mode, real(rhs_fft), 'o', mode, imag(rhs_fft), '*')
xlabel('mode'); ylabel('DFT of RHS')
legend('Real', 'Imag')
grid on
subplot(313)
plot(mode, real(sol_fft), 'o', mode, imag(sol_fft), '*')
xlabel('mode'); ylabel('DFT of SOL')
legend('Real', 'Imag')
grid on
c4science · Help