partquantity.m
e4669700ed0dgeom_testmumps

partquantity.m
View Options

			classdef partquantity < h5quantity

			properties

			end

			methods
			function obj=partquantity(filename, dataset,scale)
			if nargin<3
			scale=1;
			end
			obj=obj@h5quantity(filename, dataset, 0, 0, scale);

			end

			function quantity=val(obj,indices)
			if indices{1}==':'
			r=1:obj.nr;
			else
			r=indices{1};
			end
			if indices{2}==':'
			z=1:obj.nz;
			else
			z=indices{2};
			end
			if indices{3}==':'
			t=1:obj.nt;
			else
			t=indices{3};
			end
			if obj.centered
			quantity=zeros((obj.nz-1)*(obj.nr-1),length(t));
			else
			quantity=zeros((obj.nz)*(obj.nr),length(t));
			end
			for i=1:length(t)
			quantity(:,i) = h5read(obj.filename, obj.dataset,[1 t(i)],[Inf 1]) ;
			end
			if obj.centered
			quantity=reshape(quantity,obj.nz-1,obj.nr-1,[]);
			else
			quantity=reshape(quantity,obj.nz,obj.nr,[]);
			end
			quantity=permute(quantity,[2,1,3]);

			if obj.centered
			quantity=padarray(quantity,[1,1],0,'post');
			end
			quantity=quantity(r,z,:)*obj.scale;
			end
			end
			end