stokes_gen_model_data.m

function model_data = stokes_gen_model_data(model)
%function model_data = stokes_gen_model_data(model)
%
% generated fields of model_data:
% - grid
% - df_info
% - bc_info: structure holding information on the boundary conditions. One
% can specify boundary parts via model.bnd_rect_corner1 and
% model.bnd_rect_corner2. For boundary integrals one has to group the
% rectangles with numbers < -1 in model.bnd_rect_index. The boundaries not
% assigned in this way are used for Dirichlet conditions by default. It is
% also possible to define Dirichlet constraints for chosen dimensions of
% the field variable in model.bnd_dimrange_index. Then the rectangles have
% to be given explicitly.

% I. Maier, 22.01.2014


%% grid
grid = construct_grid(model);


%% discrete function info
switch model.df_type
    
    case 'lagrangian' % specify polynomial degree in model
        df_info = Fem.Lagrange.DefaultInfo(model.pdeg, grid);
        
        if model.dimrange > 1
            ids = cell(1, model.dimrange);
            for i = 1:model.dimrange
                ids{i} = ['dim', num2str(i)];
            end
            
            df_info = Fem.CompositeFunctionSpace(ids, repmat({df_info}, 1, model.dimrange));
        end
        
    case 'stokes_mini'
        velocity_comp = Fem.Lagrange.BubbleInfo([], grid);
        
        ids = {'velocity', 'pressure'};
        values = struct(...
            'velocity', Fem.CompositeFunctionSpace(velocity_comp, 'dim1', velocity_comp, 'dim2'), ...
            'pressure', Fem.Lagrange.DefaultInfo(1, grid));
        
        df_info = Fem.CompositeFunctionSpace(ids, values);
        
    case 'taylor_hood' % specify polynomial degree of velocity components in model
        velocity_comp = Fem.Lagrange.DefaultInfo(model.pdeg, grid);
        
        ids = {'velocity', 'pressure'};
        values = struct(...
            'velocity', Fem.CompositeFunctionSpace(velocity_comp, 'dim1', velocity_comp, 'dim2'), ...
            'pressure', Fem.Lagrange.DefaultInfo(model.pdeg - 1, grid));
        
        df_info = Fem.CompositeFunctionSpace(ids, values);
end


%% boundary condition info
bc_info = Fem.BcInfo(model, grid, df_info);


%% gather data
model_data.grid = grid;
model_data.df_info = df_info;
model_data.bc_info = bc_info;

model_data.operators = Fem.OperatorsDefault();
model_data.W = model.get_inner_product_matrix(model_data);

end