dune_laplace.m

if ~(exist('steps', 'var'))
  steps = [0, 2];
end

%setenv('DUNERBHOME', '/home/mdroh_01/projects/dune/modules/dune-rb/build_2.1/dune/rb/matlabcomm/client')

for i = 1:length(steps)
  step = steps(i);
  switch (step)
    case 0
      descr = laplace_dune_descr;
      [dmodel, rmodel] = gen_models(descr);
      model_data = gen_model_data(dmodel);
      set_mu(dmodel, [1 0]);

    case 1
      sim_data = detailed_simulation(dmodel, model_data);
      plot_sim_data(dmodel, model_data, sim_data, []);
    case 2
      tic
      detailed_data = gen_detailed_data(rmodel, model_data);
      t_offline = toc;
      reduced_data = gen_reduced_data(rmodel, detailed_data);
    case 3
      nmus = 100;
      mus = rand_uniform(nmus, rmodel.mu_ranges);
      dtimes = zeros(1, nmus);
      rtimes = zeros(1, nmus);
      rrtimes = zeros(1, nmus);
      errs = zeros(1, nmus);
      for mui = 1:size(mus,2)
        disp(['computing parameter ', num2str(mui)]);
        mu = mus(:,mui);
        set_mu(dmodel, mu);
        set_mu(rmodel, mu);
        tic;
        detailed_simulation(dmodel, model_data);
        dtimes(mui) = toc;
        tic;
        rb_sim_data = rb_simulation(rmodel, reduced_data);
        rtimes(mui) = toc;
        tic;
        rb_sim_data = rb_reconstruction(rmodel, detailed_data, rb_sim_data);
        rrtimes(mui) = toc;
        errs(mui) = rb_sim_data.Delta;
      end
      maxerr = max(errs);
      avgrt  = mean(rtimes);
      avgrrt = mean(rrtimes);
      avgdt  = mean(dtimes);
    case 4
      mus = {[1 0], [1 0.1], [1 0.2], [10 0], [10 0.1], [10 02]};
      nmus = length(mus);
      for mui = 1:nmus
        mu = mus{mui};
        set_mu(dmodel, mu);
        sim_data = detailed_simulation(dmodel, model_data);
        p = plot_sim_data(dmodel, model_data, sim_data, []);
        pause(10);
        dest = ['/home/mdrohma/projects/diss/the_diss/images/laplace_', num2str(mui)];
        movefile('/home/mdrohma/projects/dune/modules/dune-rb/build_2.1/dune/rb/examples/pdelab/laplaceccfv/vtk/out.png', dest);
        close(p);
      end
    otherwise
  end
end