get_assembled_matrices.m

function str = get_assembled_matrices(model)
%function str = get_assembled_matrices(model)
%
% This function returns the assembled matrices from comsol.
% Therefore an "assemble-feature" model.sol('sol1').feature('a1')
% must exist in the comsol model!
% 
% Input: 
%    -model = a comsol model
% Output:
%    -str = a struct containing the assemble matrices that were extracted
%    from Comsol
%
% Oliver Zeeb, 2013/11/21

model.sol('sol1').feature('v1').set('scaleval', '1');
model.sol('sol1').feature('v1').set('scalemethod', 'manual');
model.sol('sol1').feature('t1').active(false);
model.sol('sol1').feature('a1').active(true);
model.sol('sol1').feature('a1').set('Nullf', true);
model.sol('sol1').feature('a1').set('Null', true);
model.sol('sol1').feature('a1').set('Lc', true);
model.sol('sol1').feature('a1').set('ud', true);
model.sol('sol1').feature('a1').set('Kc', true);
model.sol('sol1').feature('a1').set('Ec', true);
model.sol('sol1').feature('a1').set('uscale', true);
model.sol('sol1').feature('a1').set('Dc', true);

model.study('std1').run; 

str_Null = mphmatrix(model,'sol1','out',{'Null'});
str.Null = str_Null.Null;

%create Matlab-sparse matrix from the vectors extracted from comsol
str.DcCol = 1+double(model.sol('sol1').feature('a1').getSparseMatrixCol('Dc')); 
          %"1+..." is needed since the indices returned by comsol start at 0, not at 1 as in matlab
str.DcRow = 1+double(model.sol('sol1').feature('a1').getSparseMatrixRow('Dc'));
str.DcVal = model.sol('sol1').feature('a1').getSparseMatrixVal('Dc');
str.Dc = sparse(size(str.Null,2),size(str.Null,2));
for k=1:length(str.DcRow)
    str.Dc(str.DcRow(k), str.DcCol(k)) = str.DcVal(k);
end
str.KcCol = 1+double(model.sol('sol1').feature('a1').getSparseMatrixCol('Kc'));
str.KcRow = 1+double(model.sol('sol1').feature('a1').getSparseMatrixRow('Kc'));
str.KcVal = model.sol('sol1').feature('a1').getSparseMatrixVal('Kc');
str.Kc = sparse(size(str.Null,2),size(str.Null,2));
for k=1:length(str.KcRow)
    str.Kc(str.KcRow(k), str.KcCol(k)) = str.KcVal(k);
end
str.Lc = model.sol('sol1').feature('a1').getVector('Lc');
str.ud = model.sol('sol1').feature('a1').getVector('ud');
str.uscale = model.sol('sol1').feature('a1').getVector('uscale');
end