ChebyshevTools: B=pseudo_eye(M,P) - File Exchange

No BSD License

Highlights from
ChebyshevTools

A=expD(M,ord,mod);
A_res=restrictA(A,m) this function extracts the submatrix A(1+m:end,1:end-m)
B=inv_deriv1(M); This function returns one matrix, for solving
B=inv_derivP(M,P) this function returns the Pth quasi inverse for Chebycheb differentiation
B=pseudo_eye(M,P) this function returns the pseudo-identity matrix I^(P), where there are
D=deriv(M); returns the first order chebyshev spectral differentiation matrix
D=deriv2(M);
P=cheb2poly(U) this function converts the spectral representation of a function in Chebyshev polynomials
S=stencil_mat(M,typ) this function returns the stencil matrix associated with the related
Y=cheb2phys2(X);
Y=cheb2phys3(X); 3d transform chebyshev to physical space
Y=phys2cheb2(X);
Y=phys2cheb3(X); 3d transform physical space to chebyshev
[BIG_OP,BIG_rhs]=bc_2d_new(P,... Generate the tau line boundary conditions and the appropriate rhs vector
[BIG_OP,BIG_rhs]=bc_2d_new_to... Generate the tau line boundary conditions and the appropriate rhs vector
[OP,rhs]=bc_1d(M,bc_type,bc_v...
f=galerkin2cheb(g,typ) this function transforms from the galerkin basis to the chebyshev basis
f=galerkin2cheb2(g,typ) this function transforms from the galerkin basis to the chebyshev basis in
f=galerkin2cheb3(g,typ) this function transforms from the galerkin basis to the chebyshev basis in
g=cheb2galerkin(f,typ) this function transforms from the chebyshev basis to the galerkin basis
g=cheb2galerkin2(f,typ) this function transforms from the galerkin basis to the chebyshev basis in
g=cheb2galerkin3(f,typ) this function transforms from the galerkin basis to the chebyshev basis in
test_problem(M,N,tc) This function is run from the command line taking argumetns
v=getdiag(A,n); extract the nth diag from matrix A
x=cheb_grid(N); this function returns a chebyshev grid on the interval [-1,1]
y=cheb2phys(F); Written by Mike Watson, November 29
y=definite_integral_cart(F,X1... Calculate the definite integral of a function defined by spectral
y=myC(p)
y=phys2cheb(f); Written by Mike Watson, November 29
deriv3.m
View all files

from ChebyshevTools by Michael Watson
Tool box for solving ODE/PDEs using spectral Chebyshev differentiation matrices.

B=pseudo_eye(M,P)

function B=pseudo_eye(M,P)

%this function returns the pseudo-identity matrix I^(P), where there are
%ones on the main diagonal, except the top P entries which are zero

B=eye(M);

if P>0
   B(1:P,:)=0;
else
   B(end-P+1:end,:)=0;
end
B=sparse(B);

Highlights from ChebyshevTools

Highlights from
ChebyshevTools