1 function U_der2 = lin_evol_opt_fd_fd_second_derivative(model, model_data, mu)
2 %
function U_der = lin_evol_opt_fd_fd_second_derivative(model, model_data, mu)
4 %
function calculating the second derivative solution via PDE and finite difference
5 % procedure
for the parameter mu.
7 % Markus Dihlmann 03.02.2011
13 old_mu = get_mu(model);
17 compute_derivative_old = model.compute_derivative_info;
19 model.compute_derivative_info =1;
21 %sim_data = detailed_simulation(model, model_data);
25 U_der_fix=lin_evol_opt_fd_derivative(model, model_data,mu);
26 U_der2 = cell(1,length(mu));
31 %is mu+h inside mu_range?
32 if (mu(i)+h>=model.mu_ranges{i}(1))&&(mu(i)+h<=model.mu_ranges{i}(2))
34 model = set_mu(model,mu);
35 U_der = lin_evol_opt_fd_derivative(model, model_data,mu);
37 model = set_mu(model,mu);
39 U_der2{i} = (U_der{i}-U_der_fix{i})./h;
42 model = set_mu(model,mu);
43 U_der = lin_evol_opt_fd_derivative(model, model_data,mu);
45 model = set_mu(model,mu);
47 U_der2{i} = (-U_der{i}+U_der_fix{i})./h;
52 model = set_mu(model, old_mu);
53 model.compute_derivative_info = compute_derivative_old;