%%% UserParameter.m 
%%% Oct 4, 2019.

%%% General definitions
yes = 1;
no = 0;

%%% parameters about the geometry
params_geom.shape = 2;
% 1 = disk
% 2 = rectangle

params_geom.height = 5;
params_geom.width = 3;
params_geom.radius = 20;

params_geom.plot = yes; 

%%% parameters about the finite elements mesh
params_mesh.hmax = 1;

params_mesh.plot = yes;

%%% the type of PDE
params_pde.type = 1;
%%% heat = 1;  please set d = 1, m = 0;
%%% wave = 2;  please set d = 0, m = 1;

parms_bc.type = 1;
%%% neumann = 1;
%%% dirichlet = 2;

%%% the PDE coefficients
%%% the coefficients are of the form
%%% gg+cc*exp(-((x-x0)^2+(y-y0)^2)/aa)
params_ccoeff.gg = 1;
params_ccoeff.cc = 0;
params_ccoeff.aa = 1; 
params_ccoeff.x0 = 0; 
params_ccoeff.y0 = 0;

params_acoeff.gg = 0;
params_acoeff.cc = 0;
params_acoeff.aa = 1; 
params_acoeff.x0 = 0; 
params_acoeff.y0 = 0;


params_fcoeff.gg = 0;
params_fcoeff.cc = 1;
params_fcoeff.aa = 1;
params_fcoeff.x0 = 0;
params_fcoeff.y0 = 0;

%%% the BC coefficients
%%% the coefficients are of the form
%%% coeff = gg
%%% for neumann
params_qcoeff.gg = 0;
params_gcoeff.gg = 1;

%%% for dirichlet
params_hcoeff.gg = 1;
params_rcoeff.gg = 1;

%%% the Initial Condition coefficients
%%% the coefficients are of the form
%%% gg+cc*exp(-((x-x0)^2+(y-y0)^2)/aa)
params_ic.gg = 0;
params_ic.cc = 1;
params_ic.aa = 1;
params_ic.x0 = 0;
params_ic.y0 = 0;

%%% ODE parameters
params_ode.method = 1;
%%% 1 = ode23t
%%% 2 = ode15s
params_ode.tlist = linspace(0,1,21);

params_ode.rtol = 1e-4;
params_ode.atol = 1e-6;

%%% PDE solution
params_solution.plot = yes;


