44 ivec
find(
const bvec &invector);
63 mat
ones(
int rows,
int cols);
67 imat
ones_i(
int rows,
int cols);
69 cmat
ones_c(
int rows,
int cols);
81 mat
zeros(
int rows,
int cols);
85 imat
zeros_i(
int rows,
int cols);
87 cmat
zeros_c(
int rows,
int cols);
199 template <
typename Num_T>
202 int n_rows = c.
size();
203 int n_cols = r.
size();
205 for (
int i = 0; i < n_rows; ++i) {
206 int j_limit =
std::min(n_cols, n_rows - i);
207 for (
int j = 0; j < j_limit; ++j) {
208 output(i + j, j) = c(i);
211 for (
int j = 1; j < n_cols; ++j) {
212 int i_limit =
std::min(n_rows, n_cols - j);
213 for (
int i = 0; i < i_limit; ++i) {
214 output(i, i + j) = r(j);
221 template <
typename Num_T>
226 for (
int i = 0; i < s; ++i) {
227 for (
int j = 0; j < s - i; ++j) {
228 output(i + j, j) = c(i);
231 for (
int j = 1; j < s; ++j) {
232 for (
int i = 0; i < s - j; ++i) {
233 output(i, i + j) = c(j);
256 void house(
const vec &x, vec &v,
double &beta);
262 void givens(
double a,
double b,
double &c,
double &s);
268 void givens(
double a,
double b, mat &m);
274 mat
givens(
double a,
double b);
280 void givens_t(
double a,
double b, mat &m);
473 #endif // #ifndef SPECMAT_H