%diff_testproblem    sets up Reference Examples 1.1 to 1.4
%   IFISS scriptfile: DJS, HCE; 18 March 2005, PC version. 
% Copyright (c) 2005 D.J. Silvester, H.C. Elman, A. Ramage 
gohome
clear 
fprintf('\nspecification of reference Poisson problem.')
fprintf('\nchoose specific example');
fprintf('\n     1  Square domain, constant source')
fprintf('\n     2  L-shaped domain, constant source')
fprintf('\n     3  Square domain, analytic solution')
fprintf('\n     4  L-shaped domain, analytic solution\n');
sn = default('',1);
%
if sn==1
   !copy .\diffusion\test_problems\unit_rhs.m .\diffusion\specific_rhs.m
   !copy .\diffusion\test_problems\zero_bc.m .\diffusion\specific_bc.m
   square_diff
elseif sn==2
   !copy .\diffusion\test_problems\unit_rhs.m .\diffusion\specific_rhs.m
   !copy .\diffusion\test_problems\zero_bc.m .\diffusion\specific_bc.m
   ell_diff
elseif sn==3
   !copy .\diffusion\test_problems\zero_rhs.m .\diffusion\specific_rhs.m
   !copy .\diffusion\test_problems\solution3_bc.m .\diffusion\specific_bc.m
   square_diff
elseif sn==4
   !copy .\diffusion\test_problems\zero_rhs.m .\diffusion\specific_rhs.m
   !copy .\diffusion\test_problems\solution4_bc.m .\diffusion\specific_bc.m
   ell_diff
else
   error('reference problem datafile not found!')
end