| 682 | |
| 683 | |
| 684 | void Lzero (RowMatrix& M) |
| 685 | /* Zero strict lower triangle of Matrix |
| 686 | */ |
| 687 | { |
| 688 | std::size_t i,j; |
| 689 | const std::size_t n = M.size1(); |
| 690 | assert (n == M.size2()); |
| 691 | for (i = 1; i < n; ++i) |
| 692 | { |
| 693 | RowMatrix::Row Ui(M,i); |
| 694 | for (j = 0; j < i; ++j) |
| 695 | { |
| 696 | Ui[j] = 0; |
| 697 | } |
| 698 | } |
| 699 | } |
| 700 | |
| 701 | void Uzero (RowMatrix& M) |
| 702 | /* Zero strict upper triangle of Matrix |