Skip to content

Instantly share code, notes, and snippets.

/* Fills solution into x. Warning: will modify c and d! */
// this is the old signature:
//void TridiagonalSolve(const double *a, const double *b, double *c, double *d, double *x, unsigned int n){
// this is the new signature signature:
// - we don't use arrays, instead we use the matrix and vector objects introduced in earlier pechstein exercises
// - the "unsigned int n", which denotes the size of the matrix/vectors, can be dropped because it is
// stored inside the matrix/vector objects and is redundant information.
void TridiagonalSolve(const Matrix &Kh, const Vector &fh, const Vector &uh){
int i;