|
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