Homec4science
Diffusion SPClibs (feature/refactoring) Browse

spclibs/matlab/gbmat.m
93a1e9712d04feature/refactoring

Last Change
Raw File

gbmat.m
View Options

%
% @file gbmat.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.h5';
%
% Attributes of GB matrix
%
NR=hdf5read(file,'/NX'); NTH=hdf5read(file,'/NY');
NIDBAS1=hdf5read(file,'/NIDBAS1');
NIDBAS2=hdf5read(file,'/NIDBAS2');
%mat='/MAT1';
rank=hdf5read(file,strcat(mat,'/RANK'));
ku=hdf5read(file,strcat(mat,'/KU'));
kl=hdf5read(file,strcat(mat,'/KL'));
gb_mat=hdf5read(file,mat);
rhs0=hdf5read(file,'/RHS');
sol0=hdf5read(file,'/SOL');
%
m=rank; n=rank;
a = zeros(m,n);
for i=1:m
jmin = max(1,i-kl);
jmax = min(n,i+ku);
for j=jmin:jmax
ib = kl+ku+i-j+1;
a(i,j)=gb_mat(ib,j);
end
end
%clear gb_mat;
S = sparse(a);
%clear a;
figure
spy(S);
LABEL=sprintf('nr = %d, ntheta = %d, nidbas = (%d,%d), rank = %d', NR, NTH, ...
NIDBAS1, NIDBAS2, rank);
title(LABEL)
c4science · Help