45 typedef complex<double>
Cd_t;
46 typedef valarray<double>
Vd_t;
56 D4.
coef << 1.0, 0.0, 0.0, 0.0, 0.0;
57 D3.coef << 1.0, 0.0, 0.0, 0.0;
58 D2.coef << 1.0, 0.0, 0.0;
63 Lmat << D4/16.0 +
Vd_t(2.0*
y), 1.0,
75 bcs.
eval << -1.0, 0.0,
84 eigs.
compute(Lmat, Mmat, 6, bcs);
85 std::cout <<
"Eigenvalues: \n" <<
'\n';
BcMat will hold general Boundary conditions as LinopMats at evealuation points, as given by operator ...
Eigen::Matrix< T, Eigen::Dynamic, Eigen::Dynamic > eval
Linop This class creates a Linear operator to solve TPBVPs.
Eigen::Matrix< std::complex< T >, Eigen::Dynamic, 1 > eigenvalues
This class represents a block matrix operator. It is a matrix of operators.
Eigen::Matrix< T, Eigen::Dynamic, 1 > coef
Stores the coefficients in the differential equation.
void compute(Linop< T > L, Linop< T > M, int num_vals)
Call this with an input Linear operator to solve for eigenvalues and vectors. The number of Eigen val...
int N
Specifies number of Chebyshev polynomials, default N = 31.
This class will solve the generalized eigenvalue problem for two linear operators. One of them can be singular.
std::valarray< SIS_TYPE > y(N+1)
This class holds a matrix of Chebfuns.