1 # ifndef CPPAD_SPEED_MAT_SUM_SQ_HPP
2 # define CPPAD_SPEED_MAT_SUM_SQ_HPP
133 template <
class Vector>
137 for(i = 0; i < n; i++)
138 {
for(j = 0; j < n; j++)
140 for(k = 0; k < n; k++)
141 y[i * n + j] += x[i * n + k] * x[k * n + j];
145 for(i = 0; i < n; i++)
146 {
for(j = 0; j < n; j++)
147 z[0] += y[i * n + j];
void mat_sum_sq(size_t n, Vector &x, Vector &y, Vector &z)