Method
NumCosmoMathSBesselOdeSolversolve_dense
Declaration [src]
NcmVector*
ncm_sbessel_ode_solver_solve_dense (
NcmSBesselOdeSolver* solver,
const gdouble a,
const gdouble b,
guint ell,
NcmVector* rhs,
gint nrows
)
Description [src]
Solves the ODE using a dense matrix representation with standard linear algebra. This method is useful for validation and testing, providing a reference solution to compare against the adaptive QR method.
The system is truncated to nrows equations (including 2 boundary conditions),
and solved using LAPACK’s general linear solver (LU decomposition).
Parameters
a-
Type:
const gdoubleLeft endpoint.
b-
Type:
const gdoubleRight endpoint.
ell-
Type:
guint$\ell$ multipole order.
rhs-
Type:
NcmVectorRight-hand side vector (Chebyshev coefficients of f(x)).
The data is owned by the caller of the method. nrows-
Type:
gintSize of the truncated system to solve.
Return value
Type: NcmVector
Solution vector (Chebyshev coefficients).
| The caller of the method takes ownership of the returned data, and is responsible for freeing it. |