follicle_rect_experiments.m File Reference

experiments with the model of the human follicle growth More...

Go to the source code of this file.

Functions
function res =	follicle_rect_experiments (step, model, model_data, detailed_data)
	experiments with the model of the human follicle growth More...
function J =	folliclerectexperiments>my_opt_target (mu, Ytarget, model, reduced_data)
function J =	folliclerectexperiments>my_opt_target_single_param (mu1, Ytarget, model, reduced_data)

Detailed Description

experiments with the model of the human follicle growth

Definition in file follicle_rect_experiments.m.

Function Documentation

function res = follicle_rect_experiments	(		step,
			model,
			model_data,
			detailed_data
	)

experiments with the model of the human follicle growth

model range is n=150, noutput = 50, nt = 500, time = 19.02 years implicit Euler time-discretization, accept relatively large error at initial time, for large overall time interval and few timesteps and "interactivity"

For the results of the paper, simply run step 3.4, 5, 6, 7, 8 in this order.

(step = 1 simulation and visualization of detailed dynamical system, using filecaching => do not call directly) step = 1.1, n=150, initial data distributed over larger domain init_range = 20 and different output functional expectation_M n smaller, e.g. 120: decay of output!!!! (step = 2 reduced simulation and visualization of approximation => do not call directly) (step = 3 detailed and reduced simulation and visualization of error and estimator => do not call directly) step = 3.1, n=150, n_output = 50, lower states as reduced basis expected large error, but test of functionality. step = 3.2, n=150, n_output = 50, generate POD basis of single trajectory basis is saved for use in further steps and plotted step = 3.3, n=150, n_output , generate BT basis of single trajectory basis is saved for use in further steps and plotted step = 3.4, n=150, generate PODGreedy basis of 81 trajectories basis is saved for use in further steps and plotted step = 3.5 experiments with a part of POD-basis from 3.4 step = 3.8 generate system matrices and export to file for remote BT generation. System A = 20301x20301 => ss = sys(A,B,C,D) => 3 GB system ss!!!!! balreal takes about ... step 4.0: demo_rb_gui of POD basis from 3.4 step 5.0: Optimization with RB model with POD basis from 3.4 generation of optimization plots for initial submission of paper step 5.1: Optimization with RB model with POD basis from 3.4 generation of optimization plots for revised paper step 6.0: Generation of trajectory plots for paper step 7.0: Generation of basis vector plots for paper step 8.0: Generation of POD-Greedy plots for paper

some steps explicitly require a model and model_data

Todo

generate BT basis compare quality of BT to POD on single trajectory compare quality of BT to POD on parameter corners

Parameters

step	step
model	model
model_data	model data
detailed_data	detailed data

Return values

res

res

Required fields of model:

• get_mu —  get mu
• debug —  flag indicating wether debug output shall be turned on
• set_mu —  set mu
• dim_x —  dim x
• plot_sim_data_state —  plot sim data state
• enable_error_estimator —  enable error estimator
• dim_y —  dim y
• nt —  number of time steps for evolution discretizations
• plot_detailed_data —  plot detailed data
• A_function_ptr —  A function ptr
• B_function_ptr —  B function ptr
• C_function_ptr —  C function ptr
• D_function_ptr —  D function ptr
• transformed_colormap —  transformed colormap
• plot_slice —  plot slice

Required fields of detailed_data:

• G —  G
• V —  V
• RB_info.max_err_sequence —  RB info.max err sequence
• RB_info.mu_sequence —  RB info.mu sequence

Required fields of model_data:

• Ns —  Ns
• Ms —  Ms
• G —  G
• RB —  RB

Definition at line 17 of file follicle_rect_experiments.m.

Here is the call graph for this function: