rbmatlab/1.16.09/ellipt__compliant__gen__model__data_8m_source.html

function model_data = ellipt_compliant_gen_model_data(model);

%function model_data = ellipt_compliant_gen_model_data(model)

%

% computation of mu independet data

% input: model

%

% output:

%      model_data.grid: contains p,e,t

%                       call model_data.grid.p etc

%      model_data.dl: decomposed domain

%


% B. Haasdonk 8.3.2010

% S. Langhof


% alles was geometrie generiert


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%GEOMETRY

X=[];

Y=[];

def=repmat([3;4], 1, model.B1*model.B2);

GEO=[];


%%%%%%%%%%%%%%%%%%%

%%%%% X direction

x1=0:model.B1-1;

x2=1:model.B1;


X=[x1;x2;x2;x1];

X=X*(1/model.B1);


X=repmat(X,1,model.B2);


%%%%%%%%%%%%%%%%%%%%

%%%%% Y direction

y1=1:(model.B2);

y1=repmat(y1,1,model.B1);

y1=sort(y1);


y2=0:(model.B2-1);

y2=repmat(y2,1,model.B1);

y2=sort(y2);


Y=[y1;y1;y2;y2];

Y=Y*(1/model.B2);


GEO=cat(1,def,X,Y);


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5

%%%mesh and boundary creation

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%Mesh generation

dl=decsg(GEO);

if model.mesh_type==1

    [p,e,t]=initmesh(dl);

     not=length(t);

     while not <= model.number_of_triangles_overall

         [p,e,t]=refinemesh(dl, p, e, t);

         not=length(t);

     end

else

    wert=1/(sqrt(model.B1*model.B2*model.number_of_triangles_block/3));

    [p,e,t]=initmesh(dl, 'hmax',wert);

end


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%ATTENTION%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%following section will be moved to ellipt_compliant_detailed_simulation.m

%because of the parameter dependency of the boundary matrix


%

%

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% %%%%%create boundary matrix

%

% %set boundary condition

% %bottom+top: from left to right

% %left+right: from bottom to top

%

%

% noc_bm=length(dl);                          %number of columns

%

% check_left_region = dl(6,:);                %number of block left hand side

% check_right_region = dl(7,:);

% check_region = dl(6:7,:);

%

% zeros_left=find(check_left_region == 0);          %find 0's. this edge is an outer

% zeros_right=find(check_right_region == 0);        %boundary

%

% nonzeros_left=find(check_left_region);

% nonzeros_right=find(check_right_region);

%

% inner_boundary_index=intersect(nonzeros_left, nonzeros_right);  %index of non zero columns in dl

%

% zero_index=find(check_region==0);

% outer_boundary_index=ceil(zero_index/2);

% outer_boundary_index=outer_boundary_index';

%

% check_max_length = zeros(4,max(model.B1,model.B2));

% for i=1:model.B1

%     check_max_length(1,i)=length(model.gamma_bottom{i}{1})+length(model.gamma_bottom{i}{2})+4;          %bottom

%     check_max_length(2,i)=length(model.gamma_top{i}{3})+length(model.gamma_top{i}{4})+8;                %top

% end

% for l=1:model.B2

%     check_max_length(3,l)=length(model.gamma_left{l}{1})+length(model.gamma_left{l}{2})+4;              %left

%     check_max_length(4,l)=length(model.gamma_right{l}{1})+length(model.gamma_right{l}{2})+4;            %right

% end

%

% nor_bm=max(max(check_max_length));              %number of rows BM

%

%

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% %%%Boundary matrix

% BM=zeros(nor_bm, noc_bm);               %zeros BM Matrix

%

%

% %%%%% x1/x2, y1/y2 values

% check_y1=dl(4,:);

% check_y2=dl(5,:);

%

% check_x1=dl(2,:);

% check_x2=dl(3,:);

%

% %find indices of gamma top

% one_y1_index=find(check_y1 == 1);

% one_y2_index=find(check_y1 == 1);

%

% one_y_index = intersect(one_y1_index, one_y2_index);

%

% check_x1_one = check_x1(one_y_index);

%

% top_tmp=0;

% [tmp, top_tmp]=sort(check_x1_one);

%

% k=0;

% top_index=0;                            %

% for l=top_tmp

%     k=k+1;

%     top_index(k)=one_y_index(l);   %top_index: top indices from left to right

% end

%

% %find indices of gamma bottom

%

% zero_y1_index=find(check_y1 == 0);

% zero_y2_index=find(check_y2 == 0);

%

% zero_y_index= intersect(zero_y1_index, zero_y2_index);

%

% check_x1_zero=check_x1(zero_y_index);

%

% bottom_tmp=0;

% [tmp, bottom_tmp]=sort(check_x1_zero);

% k=0;

% bottom_index=0;                            %

% for l=bottom_tmp

%     k=k+1;

%     bottom_index(k)=zero_y_index(l);   %bottom index: bottom indices from left to right

% end

%

% %find left index

%

% zero_x1_index=find(check_x1 == 0);

% zero_x2_index=find(check_x2 == 0);

%

% zero_x_index=intersect(zero_x1_index, zero_x2_index);

%

% check_y1_zero=check_y1(zero_x_index);      %

%

% left_tmp=0;

% [tmp, left_tmp]=sort(check_y1_zero);

% k=0;

% left_index=0;                            %

% for l=left_tmp

%     k=k+1;

%     left_index(k)=zero_x_index(l);   %left index: left indices from left to right

% end

%

% %find right index

%

% one_x1_index=find(check_x1 == 1);

% one_x2_index=find(check_x2 == 1);

%

% one_x_index=intersect(one_x1_index, one_x2_index);

%

% check_y1_one=check_y1(one_x_index);      %

%

% right_tmp=0;

% [tmp, right_tmp]=sort(check_y1_one);

% k=0;

% right_index=0;                            %

% for l=right_tmp

%     k=k+1;

%     right_index(k)=one_x_index(l);   %right index: right indices from left to right

% end

%

%

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% %%%%%%final creation of boundaries

%

% %outer bound indices

% outer_bound=[top_index, bottom_index, right_index, left_index];

% outer_bound=sort(outer_bound);

%

% %inner bound indices

% tmp_noc_vec=1:noc_bm;

% for i=outer_bound

%     tmp_noc_vec(i)=-1;

% end

% inner_bound=find(tmp_noc_vec >=0);

%

%

%

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55

% %%%inner boundaries

% std_inner_boundary_char=['0','0','1','0'];        %double(...) -> 48 48 49 48

% std_inner_boundary_num=[0,repmat(1,1,5)];

%

% std_inner_boundary=[std_inner_boundary_num,double(std_inner_boundary_char)];

%

% if length(std_inner_boundary) < nor_bm              %fill up with zeros

%     std_inner_boundary(nor_bm);

% end

%

% %positioning of inner boundary condition

% for i=inner_bound

%     BM(:,i)=std_inner_boundary(:);

% end

%

%

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% %%%outer boundaries

%

% %top boundary

%

% %max length of vec

% %l_top_bound=max(check_max_length(2,:));

%

% %TOP=zeros(l_top_bound, model.B1);

% clear fronttail

% fronttail=repmat(1,1,4);        %definition of dim, diri, lq, lg

%

% for i=1:model.B1

%     lh=length(model.gamma_top{i}{3});        %number of chars for h

%     lr=length(model.gamma_top{i}{4});        %number of chars for r

%     outer_topboundary_num(i,:)=[fronttail,lh, lr, 48, 48];

% end

%

% outer_topboundary_num=outer_topboundary_num';

%

% tmp_char={};

% for i=1:model.B1

%     tmp_char3=model.gamma_top{i}{3};

%     tmp_char4=model.gamma_top{i}{4};

%

%     tmp_char{i}=[tmp_char3,tmp_char4];          %

% end

%

% for i=1:model.B1

%     outer_topboundary(:,i)=[outer_topboundary_num(:,i);double(tmp_char{i})'];

%     if length(outer_topboundary(:,i)) < nor_bm

%         outer_topboundary(nor_bm,i) = 0;

%     end

% end

%

%

% %positioning of top boundary conditions

% outer_topbound_loop=outer_topboundary;

% for i=top_index

%         BM(:,i)=outer_topbound_loop(:,1);

%         outer_topbound_loop(:,1)=[];

% end

%

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% %%%bottom boundaries

%

% clear fronttail

% fronttail=[1,0];        %definition of dim, diri, lq, lg

%

% for i=1:model.B1

%     lq=length(model.gamma_bottom{i}{1});        %number of chars for

%     lg=length(model.gamma_bottom{i}{2});        %number of chars for

%     outer_bottomboundary_num(i,:)=[fronttail,lq, lg];

% end

%

% outer_bottomboundary_num=outer_bottomboundary_num';

%

% tmp_char={};

% for i=1:model.B1

%     tmp_char1=model.gamma_bottom{i}{1};

%     tmp_char2=model.gamma_bottom{i}{2};

%

%     tmp_char{i}=[tmp_char1,tmp_char2];          %

% end

%

% a=[];

% for i=1:model.B1

%     a=[];

%     a=[outer_bottomboundary_num(:,i);double(tmp_char{i})'];

%     if length(a) < nor_bm

%         a(nor_bm);

%     end

%     outer_bottomboundary(:,i)=a;

%

% %     outer_bottomboundary(:,i)=[outer_bottomboundary_num(:,i);double(tmp_char{i})'];

% %     if length(outer_bottomboundary(:,i)) < nor_bm

% %         outer_bottomboundary(nor_bm,i) = 0;

% %     end

% end

%

% %positioning of bottom boundary conditions

% outer_bottombound_loop=outer_bottomboundary;

% for i=bottom_index

%         BM(:,i)=outer_bottombound_loop(:,1);

%         outer_bottombound_loop(:,1)=[];

% end

%

%

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% %%%%%%left boundaries

% clear fronttail

% fronttail=[1,0];        %definition of dim, diri, lq, lg

%

% for i=1:model.B2

%     lq=length(model.gamma_left{i}{1});        %number of chars for

%     lg=length(model.gamma_left{i}{2});        %number of chars for

%     outer_leftboundary_num(i,:)=[fronttail,lq, lg];

% end

%

% outer_leftboundary_num=outer_leftboundary_num';

%

% tmp_char={};

% for i=1:model.B2

%     tmp_char1=model.gamma_left{i}{1};

%     tmp_char2=model.gamma_left{i}{2};

%

%     tmp_char{i}=[tmp_char1,tmp_char2];          %

% end

%

% a=[];

% for i=1:model.B2

%     a=[];

%     a=[outer_leftboundary_num(:,i);double(tmp_char{i})'];

%     if length(a) < nor_bm

%         a(nor_bm);

%     end

%     outer_leftboundary(:,i)=a;

%

% %     outer_bottomboundary(:,i)=[outer_bottomboundary_num(:,i);double(tmp_char{i})'];

% %     if length(outer_bottomboundary(:,i)) < nor_bm

% %         outer_bottomboundary(nor_bm,i) = 0;

% %     end

% end

%

% %positioning of left boundary conditions

% outer_leftbound_loop=outer_leftboundary;

% for i=left_index

%         BM(:,i)=outer_leftbound_loop(:,1);

%         outer_leftbound_loop(:,1)=[];

% end

%

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% %%%%%%right boundaries

% clear fronttail

% fronttail=[1,0];        %definition of dim, diri, lq, lg

%

% for i=1:model.B2

%     lq=length(model.gamma_right{i}{1});        %number of chars for

%     lg=length(model.gamma_right{i}{2});        %number of chars for

%     outer_rightboundary_num(i,:)=[fronttail,lq, lg];

% end

%

% outer_rightboundary_num=outer_rightboundary_num';

%

% tmp_char={};

% for i=1:model.B2

%     tmp_char1=model.gamma_right{i}{1};

%     tmp_char2=model.gamma_right{i}{2};

%

%     tmp_char{i}=[tmp_char1,tmp_char2];          %

% end

%

% a=[];

% for i=1:model.B2

%     a=[];

%     a=[outer_rightboundary_num(:,i);double(tmp_char{i})'];

%     if length(a) < nor_bm

%         a(nor_bm);

%     end

%     outer_rightboundary(:,i)=a;

%

% %     outer_bottomboundary(:,i)=[outer_bottomboundary_num(:,i);double(tmp_char{i})'];

% %     if length(outer_bottomboundary(:,i)) < nor_bm

% %         outer_bottomboundary(nor_bm,i) = 0;

% %     end

% end

%

% %positioning of left boundary conditions

% outer_rightbound_loop=outer_rightboundary;

% for i=right_index

%         BM(:,i)=outer_rightbound_loop(:,1);

%         outer_rightbound_loop(:,1)=[];

% end


grid = [];

grid.p = p;

grid.e = e;

grid.t = t;


model_data.grid = grid;


model_data.dl=dl;

%model_data.dl = dl;

%model_data.boundary_matrix=BM;


%eof


% inner = [3, 7, 2, 9];

% all = 1:100;


% b = zeros(1,length(all));

% b(inner)= 1;

% outer = find(b==0);


% b = [1,5,3,8]

% setdiff(1:100,b)