comsol_get_DBC_ind.m

function [com_DBC_ind, com_DBC_value, com_DOF_ind, ud, Null] = comsol_get_DBC_ind(cmodel)
%function [com_DBC_ind, com_DBC_value, com_DOF_ind, ud, Null] = comsol_get_DBC_ind(cmodel)
%
% function returning all the COMSOL_indices of dofs with Dirichlet Boundary
% Conditions and the according values
% Important: com_DBC_ind are the indices COMSOL uses (and not the rbmatlab-ordering)!!!
%
% GERMAN:
% Funktion, die die COMSOL-Indizes der Dofs mit Dirichlet Boundary Conditions
% zurückgibt!
%
% Oliver Zeeb
% 10.09.2012

str = mphmatrix(cmodel,'sol1','Out',{'Null','ud'});
com_DBC_ind = find(sum(str.Null,2)==0);
ud = str.ud; %particular solution --> DBC_values!
Null = str.Null;
com_DBC_value = ud(com_DBC_ind);
com_DOF_ind = find(sum(str.Null,2)~=0);