rbmatlab  1.16.09
 All Classes Namespaces Files Functions Variables Modules Pages
p_part_rb_simulation.m
1 function rb_sim_data = p_part_rb_simulation(model,reduced_data)
2 %function rb_sim_data = p_part_rb_simulation(model,reduced_data)
3 %
4 % function selecting correct reduced_data from the list and
5 % rb_simulation of base_model
6 
7 % Bernard Haasdonk 12.2.2010
8 % Markus Dihlmann 17.02.2010
9 
10 % search correct part number for current mu
11 
12 mu = get_mu(model);
13 %<<<<<<< HEAD:rbasis/p_part/p_part_rb_simulation.m
14 %disp('to be implemented: determine part-number for mu')
15 % dummy: simulate for part 1:
16 %p_part_number = 1;
17 % for this, in reduced_data, the ranges are available explicitly
18 %keyboard;
19 %=======
20 dim = length(mu);
21 %p_part_number = coord2leaf_element(reduced_data.pgrid, mu);
22 p_part_number=0;
23 for p=1:length(reduced_data.parts_reduced_data)
24  sum=0;
25  for i=1:dim
26  if((reduced_data.parts_reduced_data{p}.ranges{i}(1)<=mu(i))&&(reduced_data.parts_reduced_data{p}.ranges{i}(2)>=mu(i)))
27  sum=sum+1;
28  end;
29  end;
30  if sum==dim
31  p_part_number=p;
32  end;
33 end;
34 
35 
36 model.base_model.N=reduced_data.parts_reduced_data{p_part_number}.N;
37 rb_sim_data = ...
38  model.base_model.rb_simulation(...
39  model.base_model,...
40  reduced_data.parts_reduced_data{p_part_number});
41 
42 %| \docupdate