mlosch/interp_llc/mdsiocompact.m

function b=mdsiocompact(a,flag);
  % hack to convert to and from mdsio compact format
  
  if nargin < 2
    flag = 1;
  end
  if flag
    % convert from "nice" to compact format
    [nx,ny,nz,nt]=size(a);
    if nx/4~=round(nx/4);
      error('not a llc field');
    end

    n4=nx/4;
    m=(ny-n4)/n4;
    m4=m*n4;
    for kt = 1:nt
      for kz = 1:nz
        for k=1:4
          sides{k} = a([1:n4]+(k-1)*n4,1:m4,kz,kt);
        end
        cap = a([1:n4]+1*n4,m4+1:end,kz,kt);
        % reformat
        btmp = [sides{1} ...
                sides{2} ...
                cap ...
                reshape(rot90(sides{3},1),[n4 m4]) ...
                reshape(rot90(sides{4},1),[n4 m4]) ...
                cap*0];
        b(:,:,kz,kt)=btmp;
      end
    end
  else
    % convert from compact to "nice" format
    [nx,m,nz,nt]=size(a);
    if m/nx~=round(m/nx);
      error('not a compact format llc field');
    end

    fac = 1;
    ntiles = m/nx;
    if mod(ntiles,2); % odd number of tiles
      tilesPerFace = (ntiles-1)/4;
    else
      tilesPerFace = (ntiles-2)/4;
    end
    n4=tilesPerFace*nx;
    for kt = 1:nt
      for kz = 1:nz
        for k=1:4
          jj = [1:n4]+(k-1)*n4;
          if k>2; jj = jj + nx; end
          sides{k} = a(:,jj,kz,kt);
        end
        cap = rot90(a(:,[1:nx]+2*n4,kz,kt),1);
        % reformat
        btmp = [sides{1}; ...
                sides{2}; ...
                rot90(reshape(sides{3}(:),[n4 nx]),-1); ...
                rot90(reshape(sides{4}(:),[n4 nx]),-1) ];
        b(:,:,kz,kt)=[btmp ...
                      [cap; rot90(cap,-1)*abs(fac); ...
                       rot90(cap,-2)*abs(fac); ...
                       rot90(cap,-3)*abs(fac)]];
      end
    end
  end
  
  return
1	mlosch	1.2	function b=mdsiocompact(a,flag);
2			% hack to convert to and from mdsio compact format
3	mlosch	1.1
4	mlosch	1.2	if nargin < 2
5			flag = 1;
6	mlosch	1.1	end
7	mlosch	1.2	if flag
8			% convert from "nice" to compact format
9			[nx,ny,nz,nt]=size(a);
10			if nx/4~=round(nx/4);
11			error('not a llc field');
12			end
13
14			n4=nx/4;
15			m=(ny-n4)/n4;
16			m4=m*n4;
17			for kt = 1:nt
18			for kz = 1:nz
19			for k=1:4
20			sides{k} = a([1:n4]+(k-1)*n4,1:m4,kz,kt);
21			end
22			cap = a([1:n4]+1*n4,m4+1:end,kz,kt);
23			% reformat
24			btmp = [sides{1} ...
25			sides{2} ...
26			cap ...
27			reshape(rot90(sides{3},1),[n4 m4]) ...
28			reshape(rot90(sides{4},1),[n4 m4]) ...
29			cap*0];
30			b(:,:,kz,kt)=btmp;
31			end
32			end
33			else
34			% convert from compact to "nice" format
35			[nx,m,nz,nt]=size(a);
36			if m/nx~=round(m/nx);
37			error('not a compact format llc field');
38			end
39	mlosch	1.1
40	mlosch	1.2	fac = 1;
41			ntiles = m/nx;
42			if mod(ntiles,2); % odd number of tiles
43			tilesPerFace = (ntiles-1)/4;
44			else
45			tilesPerFace = (ntiles-2)/4;
46			end
47			n4=tilesPerFace*nx;
48			for kt = 1:nt
49			for kz = 1:nz
50			for k=1:4
51			jj = [1:n4]+(k-1)*n4;
52			if k>2; jj = jj + nx; end
53			sides{k} = a(:,jj,kz,kt);
54			end
55			cap = rot90(a(:,[1:nx]+2*n4,kz,kt),1);
56			% reformat
57			btmp = [sides{1}; ...
58			sides{2}; ...
59			rot90(reshape(sides{3}(:),[n4 nx]),-1); ...
60			rot90(reshape(sides{4}(:),[n4 nx]),-1) ];
61			b(:,:,kz,kt)=[btmp ...
62			[cap; rot90(cap,-1)*abs(fac); ...
63			rot90(cap,-2)*abs(fac); ...
64			rot90(cap,-3)*abs(fac)]];
65	mlosch	1.1	end
66			end
67			end
68
69			return