tridia-solver moved to seperate module; the tridiagonal matrix coefficients of array tri are calculated only once at the beginning
