15 typedef complex<double>
Cd_t;
16 typedef valarray<double>
Vd_t;
26 D4.
coef << 1.0, 0.0, 0.0, 0.0, 0.0;
27 D3.coef << 1.0, 0.0, 0.0, 0.0;
28 D2.coef << 1.0, 0.0, 0.0;
33 Lmat << D4 - (8.0 * D2) + (16.0 * D0);
34 Mmat << D2 - (4.0 * D0);
47 eigs.
compute(Lmat, Mmat, 6, bcs);
48 std::cout <<
"Eigenvalues: \n" <<
'\n';
BcMat will hold general Boundary conditions as LinopMats at evealuation points, as given by operator ...
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.
This class holds a matrix of Chebfuns.