| 988 | |
| 989 | template<typename typ> |
| 990 | static inline void |
| 991 | triu_matrix(typ *matrix, size_t n) |
| 992 | { |
| 993 | size_t i, j; |
| 994 | matrix += n; |
| 995 | for (i = 1; i < n; ++i) { |
| 996 | for (j = 0; j < i; ++j) { |
| 997 | matrix[j] = numeric_limits<typ>::zero; |
| 998 | } |
| 999 | matrix += n; |
| 1000 | } |
| 1001 | } |
| 1002 | |
| 1003 | |
| 1004 | /* -------------------------------------------------------------------------- */ |