|
2 | 2 | //
|
3 | 3 | // fastLm.cpp: Rcpp/Eigen example of a simple lm() alternative
|
4 | 4 | //
|
5 |
| -// Copyright (C) 2011 - 2015 Douglas Bates, Dirk Eddelbuettel and Romain Francois |
| 5 | +// Copyright (C) 2011 - 2022 Douglas Bates, Dirk Eddelbuettel and Romain Francois |
6 | 6 | //
|
7 | 7 | // This file is part of RcppEigen.
|
8 | 8 | //
|
@@ -144,11 +144,11 @@ namespace lmsol {
|
144 | 144 | if (m < n || S.size() != n || Vt.rows() != n || Vt.cols() != n)
|
145 | 145 | throw std::invalid_argument("dimension mismatch in gesvd"); // #nocov
|
146 | 146 | F77_CALL(dgesdd)("O", &m, &n, A.data(), &m, S.data(), A.data(),
|
147 |
| - &m, Vt.data(), &n, &wrk, &mone, &iwork[0], &info); |
| 147 | + &m, Vt.data(), &n, &wrk, &mone, &iwork[0], &info FCONE); |
148 | 148 | int lwork(wrk);
|
149 | 149 | std::vector<double> work(lwork);
|
150 | 150 | F77_CALL(dgesdd)("O", &m, &n, A.data(), &m, S.data(), A.data(),
|
151 |
| - &m, Vt.data(), &n, &work[0], &lwork, &iwork[0], &info); |
| 151 | + &m, Vt.data(), &n, &work[0], &lwork, &iwork[0], &info FCONE); |
152 | 152 | return info;
|
153 | 153 | }
|
154 | 154 |
|
|
0 commit comments