1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
//!
//! # BVP_sci - Advanced Boundary Value Problem Solver Module
//!
//! This module provides a comprehensive suite of high-performance solvers for Boundary Value Problems (BVP)
//! in ordinary differential equations. It implements state-of-the-art numerical methods with both symbolic
//! and direct numerical approaches.
//!
//! ## Key Features
//! - **Multiple Linear Algebra Backends**: Support for both `faer` and `nalgebra` crates
//! - **Symbolic Integration**: Automatic conversion from symbolic expressions to numerical functions
//! - **Analytical Jacobians**: Symbolic differentiation for faster convergence
//! - **Adaptive Mesh Refinement**: Automatic grid adaptation for optimal accuracy
//! - **Collocation Methods**: 4th-order collocation algorithms with residual control
//! - **Sparse Matrix Support**: Efficient handling of large, sparse systems
//! - **Variable Bounds**: Numerical stability through variable constraints
//! - **Comprehensive Testing**: Extensive test suites for reliability
//!
//! ## Module Structure
//! - `BVP_sci_faer`: Core BVP solver using faer linear algebra (recommended for performance)
//! - `BVP_sci_numerical`: Pure numerical no-symbolic user-facing API on top of the faer core
//! - `BVP_sci_nalgebra`: Alternative implementation using nalgebra
//! - `BVP_sci_symb`: High-level symbolic wrapper for easy problem setup
//! - `BVP_sci_symbolic_functions`: Symbolic-to-numerical conversion utilities
//! - `BVP_sci_utils`: Common utilities and helper functions
//!
//! ## Supported Problem Types
//! - Linear and nonlinear BVPs
//! - Systems of coupled ODEs
//! - Problems with parameters
//! - Stiff and non-stiff equations
//! - Multi-point boundary conditions
//!
//! ## Performance Characteristics
//! - Parallel symbolic differentiation
//! - Sparse matrix optimizations
//! - Bandwidth-aware Jacobian storage
//! - Memory-efficient mesh handling
//! - Adaptive error control
//!
///BVP solver using faer crate for matrix&vector operations
/// BVP solver using nalgebra crate for matrix&vector operations