1.13.10/doc/test__p__part_8m_source.html

function OK = test_p_part

% Test of p-partition functionality

%

% Consitency check of p-partition algorithm is included in

% p_part_gen_detailed_data, hence no consisstency check neccessary in this

% test-file.

%

% Markus Dihlmann 30.04.2010


OK=1;

disp('Test of p-partition');

params.coarse_factor = 8;

model = advection_fv_output_opt_model(params);

mu = [1,1,1];

model.RB_numintervals = [2,2,2];


model = set_mu(model,mu);


model.verbose=8;


model.RB_generation_mode = 'greedy_uniform_fixed';

model.RB_error_indicator = 'estimator';

model.RB_extension_algorithm = @RB_extension_PCA_fixspace;


model.RB_stop_epsilon=2;

model.RB_stop_timeout=100*60;

model.RB_stop_Nmax=4;


model.p_part_max_refinement_level=1;


%params.p_part_generation_mode = 'adaptive';

params.p_part_generation_mode = 'adaptive_epsilon_fixed';

params.p_part_early_refinement = 1;

params.p_part_numintervals = [1,1,1];


model = p_part_model(model,params);


try

    model_data = gen_model_data(model);

    detailed_data = gen_detailed_data(model,model_data);

    detailed_data = clean_up_part_detailed_data(detailed_data);

catch ME

    disp('Error in p-partition');

    OK=0;

end