Mercurial > hg > octave-lyh

comparison liboctave/floatQR.cc @ 8811:20dfb885f877

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
int -> octave_idx fixes
author John W. Eaton <jwe@octave.org>
date Wed, 18 Feb 2009 23:34:03 -0500
parents c7864bb74914
children b03953732530
comparison
equal deleted inserted replaced
8810:c9e1db15035b 8811:20dfb885f877
106 // workspace query. 106 // workspace query.
107 float rlwork; 107 float rlwork;
108 F77_XFCN (sgeqrf, SGEQRF, (m, n, afact.fortran_vec (), m, tau, &rlwork, -1, info)); 108 F77_XFCN (sgeqrf, SGEQRF, (m, n, afact.fortran_vec (), m, tau, &rlwork, -1, info));
109 109
110 // allocate buffer and do the job. 110 // allocate buffer and do the job.
111 octave_idx_type lwork = rlwork; lwork = std::max (lwork, 1); 111 octave_idx_type lwork = rlwork;
112 lwork = std::max (lwork, static_cast<octave_idx_type> (1));
112 OCTAVE_LOCAL_BUFFER (float, work, lwork); 113 OCTAVE_LOCAL_BUFFER (float, work, lwork);
113 F77_XFCN (sgeqrf, SGEQRF, (m, n, afact.fortran_vec (), m, tau, work, lwork, info)); 114 F77_XFCN (sgeqrf, SGEQRF, (m, n, afact.fortran_vec (), m, tau, work, lwork, info));
114 } 115 }
115 116
116 form (n, afact, tau, qr_type); 117 form (n, afact, tau, qr_type);
173 float rlwork; 174 float rlwork;
174 F77_XFCN (sorgqr, SORGQR, (m, k, min_mn, q.fortran_vec (), m, tau, 175 F77_XFCN (sorgqr, SORGQR, (m, k, min_mn, q.fortran_vec (), m, tau,
175 &rlwork, -1, info)); 176 &rlwork, -1, info));
176 177
177 // allocate buffer and do the job. 178 // allocate buffer and do the job.
178 octave_idx_type lwork = rlwork; lwork = std::max (lwork, 1); 179 octave_idx_type lwork = rlwork;
180 lwork = std::max (lwork, static_cast<octave_idx_type> (1));
179 OCTAVE_LOCAL_BUFFER (float, work, lwork); 181 OCTAVE_LOCAL_BUFFER (float, work, lwork);
180 F77_XFCN (sorgqr, SORGQR, (m, k, min_mn, q.fortran_vec (), m, tau, 182 F77_XFCN (sorgqr, SORGQR, (m, k, min_mn, q.fortran_vec (), m, tau,
181 work, lwork, info)); 183 work, lwork, info));
182 } 184 }
183 } 185 }