test_poisson_solver.m
46050b96d3admaster

test_poisson_solver.m
View Options

			r = linspace(1e-3,1e-2,1000);

			n0=1e17;
			rm=6e-3;
			rp=7e-3;
			n=zeros(size(r));
			n(r>=rm & r<=rp)=n0;

			%% Numerical solution
			[Er,phi]=poisson_solver(-n,r,0,0);

			%% Analytical solution of the potential
			[~,~,rana,phiana,~,~]=analytical_solver(1,0,0.17,1e-3,1e-2,6e-3,7e-3,1e17);

			%% Nonlinear solver solution
			resfile='/misc/nonlinear/results_m6_Np1e6_Nr512_Nthet1024.h5';
			M=nonlinhdf5(resfile);

			rnonlin=M.rgrid;
			phinonlin=M.phi(:,1,1);


			figure
			plot(r,phi,rana,phiana,rnonlin,phinonlin)
			figure
			plot(r,n)