27 template <typename TYPE, size_t C, size_t R>
32 template <typename TYPE, size_t C, size_t R>
33 mat<TYPE, C, R>& doAssign(  in doAssign()  argument
34         mat<TYPE, C, R>& lhs,  in doAssign()
36     for (size_t i=0 ; i<C ; i++)  in doAssign()
42 template <typename TYPE, size_t C, size_t R, size_t D>
43 mat<TYPE, C, R> PURE doMul(  in doMul()  argument
45         const mat<TYPE, C, D>& rhs)  in doMul()
47     mat<TYPE, C, R> res;  in doMul()
48     for (size_t c=0 ; c<C ; c++) {  in doMul()
76 template <typename TYPE, size_t C, size_t R>
77 mat<TYPE, C, R> PURE doMul(  in doMul()  argument
79         const mat<TYPE, C, 1>& rhs)  in doMul()
81     mat<TYPE, C, R> res;  in doMul()
82     for (size_t c=0 ; c<C ; c++) {  in doMul()
90 template <typename TYPE, size_t C, size_t R>
91 mat<TYPE, C, R> PURE doMul(  in doMul()  argument
92         const mat<TYPE, C, R>& rhs,  in doMul()
95     mat<TYPE, C, R> res;  in doMul()
96     for (size_t c=0 ; c<C ; c++) {  in doMul()
104 template <typename TYPE, size_t C, size_t R>
105 mat<TYPE, C, R> PURE doMul(  in doMul()  argument
107         const mat<TYPE, C, R>& rhs)  in doMul()
109     mat<TYPE, C, R> res;  in doMul()
110     for (size_t c=0 ; c<C ; c++) {  in doMul()
123 template <typename TYPE, size_t C, size_t R>
124 class mat : public vec< vec<TYPE, R>, C > {
126     typedef vec< vec<TYPE, R>, C > base;
133     size_type size() const { return R*C; }  in size()
134     enum { ROWS = R, COLS = C };
187             const mat<TYPE, C, D>& rhs) {
193             const mat& lhs, const vec<TYPE, C>& rhs) {
199             const vec<TYPE, R>& lhs, const mat<TYPE, C, 1>& rhs) {
241             const column_builder<C-2>& lhs,
243         lhs.matrix[C-1] = rhs;
289 template <typename TYPE, size_t C, size_t R>
290 mat<TYPE, R, C> PURE transpose(const mat<TYPE, C, R>& m) {  in transpose()  argument
291     mat<TYPE, R, C> r;  in transpose()
293         for (size_t j=0 ; j<C ; j++)  in transpose()
299 template <typename TYPE, size_t C> static TYPE trace(const mat<TYPE, C, C>& m) {  in trace()  argument
301     for (size_t i=0 ; i<C ; i++)  in trace()
307 template <typename TYPE, size_t C>
308 static bool isPositiveSemidefinite(const mat<TYPE, C, C>& m, TYPE tolerance) {  in isPositiveSemidefinite()  argument
309     for (size_t i=0 ; i<C ; i++)  in isPositiveSemidefinite()
313     for (size_t i=0 ; i<C ; i++)  in isPositiveSemidefinite()
314       for (size_t j=i+1 ; j<C ; j++)  in isPositiveSemidefinite()