f.tridiagonal.solver: Tridiagonal solver

In TESTgroup-BNL/LeafGasExchange: A package containing functions for fitting photosynthetic response curves, simulat leaf and canopy photosynthesis

Tridiagonal solver

Description

of length N, D is a vector of length N, and R is an N x N tridiagonal matrix defined by the vectors A, B, C each of length N. A(1) and C(N) are undefined and are not referenced.

|B(1) C(1) ... ... ... | |A(2) B(2) C(2) ... ... | R = | A(3) B(3) C(3) ... | | ... A(N-1) B(N-1) C(N-1)| | ... ... A(N) B(N) |

The system of equations is written as:

A_i * U_i-1 + B_i * U_i + C_i * U_i+1 = D_i

for i = 1 to N. The solution is found by rewriting the equations so that:

U_i = F_i - E_i * U_i+1

Usage

f.tridiagonal.solver(a, b, c, d, n)

Arguments

a	See description.
b	See description.
c	See description.
d	See description.
n	See description.

Details

Tridiagonal solver Converted into a R code from the original code of Gordon Bonan: Bonan, G. (2019). Climate Change and Terrestrial Ecosystem Modeling. Cambridge: Cambridge University Press. doi:10.1017/9781107339217

Value

Solution U

TESTgroup-BNL/LeafGasExchange documentation built on May 17, 2024, 8:11 p.m.