Homec4science
Diffusion SPClibs (feature/fft-split) Browse

spclibs/matlab/pde2d_mumps.m
03540012fedefeature/fft-split

Last Change
Raw File

pde2d_mumps.m
View Options

%
% @file pde2d_mumps.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='pde2d_mumps.h5';
nr=hdf5read(file,'/', 'NX');
nth=hdf5read(file,'/', 'NY');
NIDBAS1=hdf5read(file,'/','NIDBAS1');
NIDBAS2=hdf5read(file,'/','NIDBAS2');
MBESS=hdf5read(file,'/','MBESS');
r=hdf5read(file,'/xgrid');
t=hdf5read(file,'/ygrid');
sol=hdf5read(file,'/sol')';
solexact=hdf5read(file,'/solana')';
err=hdf5read(file,'/errors')';
solr=hdf5read(file,'/derivx')';
solt=hdf5read(file,'/derivy')';
LABEL=sprintf('nr = %d, ntheta = %d, nidbas = (%d,%d), mbess = %d', nr, nth, ...
NIDBAS1, NIDBAS2, MBESS);
figure
subplot(211)
plot(r, sol(nth/2,:), 'o', r, solexact(nth/2,:))
xlabel('r')
grid on
title(LABEL)
subplot(212)
if MBESS == 0
plot(t, sol(:,1), 'o', t, solexact(:,1))
else
plot(t, sol(:,nr/2), 'o', t, solexact(:,nr/2))
end
xlabel('\theta')
grid on
% $$$ if verLessThan('matlab', '7.9');
% $$$ n = hdf5read(file,'/MAT/RANK');
% $$$ nnz = hdf5read(file,'/MAT/NNZ');
% $$$ nlsym = hdf5read(file,'/MAT/NLSYM');
% $$$ else
% $$$ n = hdf5read(file,'/MAT/', 'RANK');
% $$$ nnz = hdf5read(file,'/MAT/', 'NNZ');
% $$$ nlsym = hdf5read(file,'/MAT/', 'NLSYM');
% $$$ end
% $$$
% $$$ cols=hdf5read(file, '/MAT/cols');
% $$$ irow=hdf5read(file, '/MAT/irow');
% $$$ val=hdf5read(file, '/MAT/val');
% $$$ perm=hdf5read(file, '/MAT/perm');
% $$$
% $$$ rows = zeros(nnz,1);
% $$$ cols=double(cols);
% $$$ irow=double(irow);
% $$$ perm=double(perm);
% $$$
% $$$ for i=1:n
% $$$ s = irow(i);
% $$$ e = irow(i+1)-1;
% $$$ rows(s:e) = i;
% $$$ end
% $$$
% $$$ mat = sparse(rows,cols,val);
% $$$ figure
% $$$ subplot(121);
% $$$ spy(mat(perm,perm));
% $$$ title('Matrix structure')
% $$$ subplot(122);
% $$$ spy(chol(mat(perm,perm)));
% $$$ title('Factor L^T')
c4science · Help