12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273 |
- goog.provide('goog.math.tdma');
- goog.math.tdma.solve = function(
- subDiag, mainDiag, supDiag, vecRight, opt_result) {
-
- mainDiag = mainDiag.slice();
- vecRight = vecRight.slice();
-
- var nDim = mainDiag.length;
-
-
- for (var i = 1; i < nDim; ++i) {
- var m = subDiag[i - 1] / mainDiag[i - 1];
- mainDiag[i] = mainDiag[i] - m * supDiag[i - 1];
- vecRight[i] = vecRight[i] - m * vecRight[i - 1];
- }
-
- var result = opt_result || new Array(vecRight.length);
- result[nDim - 1] = vecRight[nDim - 1] / mainDiag[nDim - 1];
- for (i = nDim - 2; i >= 0; --i) {
- result[i] = (vecRight[i] - supDiag[i] * result[i + 1]) / mainDiag[i];
- }
- return result;
- };
|