ChebyshevTools: A=expD(M,ord,mod); - 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.

A=expD(M,ord,mod);

function A=expD(M,ord,mod);

if nargin==1
    ord=1;
    mod=0;
elseif nargin==2
    mod=0;
end

% build the matrix operator A = Exp[D^ord] 

D=deriv(M)^ord;

A=zeros(M);

for j=1:ceil(M/ord)
   A=A+(D^(j-1))*(1/factorial(j-1));
end

%modify the first line
if (mod==1)
    A(1,:)=A(1,:)*2;
end

Highlights from ChebyshevTools

Highlights from
ChebyshevTools