% moje testovaci uloha cislo
x_l = -5; x_u = 15; y_l = -5; y_u = 15; z_l = -5; z_u = 15; T = 10;
a_advec = 1; b_advec = 1; c_advec = 1;
RHS = @(x,y,z,t) 0 * x + 0 * y + 0 * z + 0 * t;
b = 1; IC = @(x,y,z) exp(-log(2) / (b ^ 2) * (x .^ 2 + y .^ 2 + z .^ 2));
BC = @(x,y,z,t) IC(x - a_advec * t,y - b_advec * t,z - c_advec * t);

%[P,E] = MshToMatlab3D(['Advection3D_16elem.msh']);
mesh_size_parametr = 2.00;
[P,E] = MshToMatlab3D(['Advection3D_h=',sprintf('%0.2f',mesh_size_parametr),'.msh']);
%[P,E,E_partitions,E_partitions_neighbours] = MshToMatlab3D_with_partitions(['Advection3D_h=',sprintf('%0.2f',mesh_size_parametr),'_2partitions.msh']);

x_axis_limits = [x_l, x_u]; y_axis_limits = [y_l, y_u]; z_axis_limits = [z_l, z_u];


% moje testovaci uloha cislo - mensi na debug
% x_l = -5; x_u = 10; y_l = -5; y_u = 10; z_l = -5; z_u = 10; T = 5;
% a_advec = 1; b_advec = 1; c_advec = 1;
% RHS = @(x,y,z,t) 0 * x + 0 * y + 0 * z + 0 * t;
% b = 1; IC = @(x,y,z) exp(-log(2) / (b ^ 2) * (x .^ 2 + y .^ 2 + z .^ 2));
% BC = @(x,y,z,t) IC(x - a_advec * t,y - b_advec * t,z - c_advec * t);
% 
% mesh_size_parametr = 2.00;
% [P,E] = MshToMatlab3D(['Advection3D_smaller_h=',sprintf('%0.2f',mesh_size_parametr),'.msh']);
% 
% x_axis_limits = [x_l, x_u]; y_axis_limits = [y_l, y_u]; z_axis_limits = [z_l, z_u];