porsche_gen_model_data.m

function model_data=porsche_gen_model_data(model)
%function model_data=porsche_gen_model_data(model)
%
% function generating a macro-triangulation
%
% needed fields of model:
% pmacro: n-by-2-vector of all the macro-points
% tmacro: m-by-3-vector of the points belonging to one macro-triangle
% cp:     k-by-2-vector of all the carpoints (the model is described as
%         polygonial of k points)
%
% Oliver Zeeb, 31.01.11

tmacro = model.tmacro;
pmacro = model.pmacro;
cp=model.cp;


% Create macro-grid
% =================
macrogrid = triagrid(pmacro', tmacro',[]);
disp('created macrogrid...')


% Create micro-grid: (either with pdetoolbox or with presaved data)
% ==================
if model.use_saved_pet_for_model_data % In case of pre-saved p, e, t - data: 
                                      % microgrid can be obtained with
                                      % these data.
    load porsche_turbo_p_e_t_microgrid;
    model_data.pdetoolbox_mesh.p = p;
    model_data.pdetoolbox_mesh.e = e;
    model_data.pdetoolbox_mesh.t = t;
    disp('p, e, t, data loaded from porsche_turbo_p_e_t_microgrid.mat')
else % p, e, t data generated by pdetoolbox:
    % initialize pdetoolbox
    [pde_fig,ax]=pdeinit;
    pdetool('appl_cb',1);
    set(ax,'DataAspectRatio',[1 1 1]);
    set(ax,'PlotBoxAspectRatio',[1.5 1 1]);
    set(ax,'XLim',[-50 750]);
    set(ax,'YLim',[-50 350]);
    set(ax,'XTickMode','auto');
    set(ax,'YTickMode','auto');

    %create the macro-triangles as defined in the model
    for k=1:size(tmacro,1)
        tria_pts = pmacro(tmacro(k,:),:); % points of the triangle
        pdepoly(tria_pts(:,1), tria_pts(:,2), ['T',num2str(k)])
    end

    % whole car: 
    pdepoly(cp(:,1),cp(:,2),'PorscheTurbo')

    %set area to: sum(triangles) - PorscheTurbo
    set(findobj(get(pde_fig,'Children'),'Tag','PDEEval'),'String',...
        '(T1+T2+T3+T4+T5+T6+T7+T8+T9+T10+T11+T12+T13+T14+T15+T16+T17+T18+T19+T20+T21+T22)-PorscheTurbo')

    % Mesh generation:
    setappdata(pde_fig,'Hgrad',1.3);
    setappdata(pde_fig,'refinemethod','regular');
    setappdata(pde_fig,'jiggle',char('on','mean',''));
    pdetool('initmesh')

    % Export grid from PDE-Toolbox by hand
    disp('====================================================================')
    disp('Please export the grid information from the pdetoolbox.')
    disp('In the PDE Toolbox choose "Mesh" --> "Export Mesh..." and press "OK"')
    disp('then continue by typing "dbcont"')
    disp('====================================================================')
    keyboard

    %PDE Toolbox exports mesh to the workspace, so usage of 'evalin' is needed
    p = evalin('base', 'p');
    e = evalin('base', 'e');
    t = evalin('base', 't');

    model_data.pdetoolbox_mesh.p = p;
    model_data.pdetoolbox_mesh.e = e;
    model_data.pdetoolbox_mesh.t = t;

    evalin('base', 'clear p e t');
    disp('p, e, t, data created and exported from pdetoolbox')
end

microgrid = triagrid(p,t,[]);
microgrid = set_boundary_types(microgrid, model);
disp('created microgrid...')
df_info = feminfo(model,microgrid);


% Setting grid and dfinfo in the model_data
% =========================================
model_data.grid=microgrid;
model_data.macrogrid=macrogrid;
model_data.df_info=df_info;


% Information about the front-part of the car
% ===========================================
%find the elements belonging to the boundary with corresponding
%element-index and edge-index
bnd_ind = find(model_data.grid.NBI(:)<0);
[el_ind_bnd,local_edge_ind_bnd] =  ind2sub(size(model_data.grid.NBI),bnd_ind); %lokaler Edge-Index (Kante 1, 2 oder 3 vom Element)

[x_car_min, x_car_min_ind] = min(model.cp(:,1));
[y_car_min, y_car_min_ind] = min(model.cp(:,2));
[y_car_max, y_car_max_ind] = max(model.cp(:,2));
x_car_max = model.cp(y_car_max_ind,1); % x-value of the car point with maximum y_value

%find points on the front of the car
ind_front=find(model_data.grid.ECX(bnd_ind) > x_car_min-eps & ...
       model_data.grid.ECY(bnd_ind) > y_car_min+eps & ... %damit werden punkte des unterbodens ausgeschlossen
       model_data.grid.ECY(bnd_ind) < y_car_max+eps & ...
       model_data.grid.ECX(bnd_ind) < x_car_max+eps);

el_ind_front = el_ind_bnd(ind_front);
local_edge_ind_front = local_edge_ind_bnd(ind_front);
linear_ind_front=sub2ind(size(model_data.grid.VI), el_ind_front, local_edge_ind_front); %single indices are needed for getting the points in model_data.grid.VI

%%%%%%%%%%%%%%%%%
edge_length=model_data.grid.EL(linear_ind_front); %length of the edges belonging to the front

%get the tangent vectors: rotate normal-vector by 90 degree!
% rot_angle = pi/2;
% rot_matrix=[cos(rot_angle), -sin(rot_angle); sin(rot_angle), cos(rot_angle)];
nx_front=model_data.grid.NX(linear_ind_front); %normal vectors of the edges belonging to the front  
ny_front=model_data.grid.NY(linear_ind_front);
rot_matrix=[0, -1 ; 1 0];
tangent_vectors=rot_matrix * [nx_front'; ny_front'];
tx_front = tangent_vectors(1,:)';
ty_front = tangent_vectors(2,:)';
edge_midpts_x=model_data.grid.ECX(linear_ind_front);
edge_midpts_y=model_data.grid.ECY(linear_ind_front);
%%%%%%%%%%%%%%%%%%%%%%%%

%Adding the information to model_data
model_data.front_ref_domain.el_ind_front = el_ind_front;
model_data.front_ref_domain.local_edge_ind_front = local_edge_ind_front;
model_data.front_ref_domain.linear_ind_front = linear_ind_front;
model_data.front_ref_domain.edge_length_front=edge_length;
model_data.front_ref_domain.normal_front=[nx_front, ny_front];
model_data.front_ref_domain.tangent_front=[tx_front, ty_front];
model_data.front_ref_domain.edge_midpts_front=[edge_midpts_x, edge_midpts_y];