NaphtaliSandholmMethodSolveBlockTridiagonal Method

Solves a block-tridiagonal linear system J*dx = rhs using the Thomas algorithm for block matrices. O(nBlocks * blockSize^3) instead of O(n^3) for dense solve.

private static double[] SolveBlockTridiagonal(
	double[,] J,
	double[] rhs,
	int blockSize,
	int nBlocks
)

Private Shared Function SolveBlockTridiagonal ( 
	J As Double(,),
	rhs As Double(),
	blockSize As Integer,
	nBlocks As Integer
) As Double()