displaydirichletweight.m
4ebefa948a9cmaster

displaydirichletweight.m
View Options

			function displaydirichletweight(obj)
			%DISPLAYGTILDE Displays the dirichlet weight for the geometry
			%% dirichlet weight function
			f=figure();
			geomw=obj.geomweight(:,:,1);
			%gtilde=(Mtest.gtilde(:,:,1)./(1-Mtest.geomweight(:,:,1))-2).*(1-Mtest.geomweight(:,:,1));%/Mtest.phinorm-1;
			gtilde=(obj.dirichletweight(:,:,1));
			gtilde(obj.geomweight(:,:,1)<0)=NaN;
			contourf(obj.zgrid1e3, obj.rgrid1e3,gtilde,20)
			hold on
			contour(obj.zgrid1e3,obj.rgrid1e3,geomw,[0 0],'r-','linewidth',1.5)
			%title('\tilde{g}')
			xlabel('z [mm]')
			ylabel('r [mm]')
			xlim([obj.zgrid(1) obj.zgrid(end)]*1e3)
			ylim([obj.rgrid(1) obj.rgrid(end)]*1e3)
			c = colorbar;
			c.Label.String='w(r,z) [a.u.]';
			axis equal
			obj.savegraph(f,sprintf('%s/%s_gtilde',obj.folder,obj.name));
			end