Method

NumCosmoMathSBesselOdeSolverget_operator_matrix_colmajor

Declaration [src]

NcmMatrix*
ncm_sbessel_ode_solver_get_operator_matrix_colmajor (
  NcmSBesselOdeSolver* solver,
  const gdouble a,
  const gdouble b,
  guint ell,
  gint nrows
)

Description [src]

Builds and returns a dense matrix representation of the differential operator truncated to nrows rows in column-major format (suitable for LAPACK/BLAS routines).

The matrix includes: - Row 0: boundary condition u(-1) = 0 - Row 1: boundary condition u(+1) = 0 - Rows 2 to nrows-1: differential operator rows

The matrix is square (nrows x nrows) and stored in column-major order as expected by LAPACK routines like dgesv.

Parameters

a

Type: const gdouble

Left endpoint.

b

Type: const gdouble

Right endpoint.

ell

Type: guint

$\ell$ multipole order.

nrows

Type: gint

Number of rows to extract (including 2 boundary condition rows).

Return value

Type: NcmMatrix

Dense matrix representation of the operator (column-major).

The caller of the method takes ownership of the returned data, and is responsible for freeing it.