0001 function error_on_tiling = s2let_check_tiling(psi, phi, L, spin, J)
0002
0003
0004
0005
0006
0007
0008
0009
0010 identity = zeros(1,L);
0011 for l=abs(spin):L-1
0012 identity(1,l+1) = identity(1,l+1) + 4*pi/(2*l+1) * phi(l+1) * conj(phi(l+1));
0013 end
0014
0015 for j=0:J
0016 ind = spin*spin + 1;
0017 for l=abs(spin):L-1
0018 for m=-l:l
0019 identity(1,l+1) = identity(1,l+1) + 8*pi^2/(2*l+1) * psi(ind, j+1) * conj(psi(ind, j+1));
0020 ind = ind + 1;
0021 end
0022 end
0023 end
0024
0025 error_on_tiling = 0;
0026 for l=abs(spin):L-1
0027 error_on_tiling = error_on_tiling + identity(1,l+1) - 1.0;
0028 end
0029
0030 end