Mercurial > jwe > octave

comparison liboctave/array/dRowVector.cc @ 21120:499b851fbfae

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Replace pattern if/err_XXX/else/code with if/err_XXX/ code. * schur.cc, ov-complex.h, ov-cx-mat.cc, ov-cx-sparse.cc, ov-float.h, ov-flt-complex.h, ov-flt-cx-mat.cc, ov-flt-re-mat.cc, ov-range.cc, ov-re-mat.cc, ov-re-sparse.cc, ov-scalar.h, ov-str-mat.cc, ops.h, pt-idx.cc, Array.cc, CColVector.cc, CMatrix.cc, CRowVector.cc, MSparse.cc, PermMatrix.cc, Sparse.cc, dColVector.cc, dMatrix.cc, dRowVector.cc, dSparse.cc, fCColVector.cc, fCMatrix.cc, fCRowVector.cc, fColVector.cc, fMatrix.cc, fRowVector.cc: Replace pattern if/err_XXX/else/code with if/err_XXX/ code.
author Rik <rik@octave.org>
date Wed, 20 Jan 2016 16:58:59 -0800
parents e39e05d90788
children 7cac4e7458f2
comparison
equal deleted inserted replaced
21119:90cd0f9442d5 21120:499b851fbfae
199 octave_idx_type a_nr = a.rows (); 199 octave_idx_type a_nr = a.rows ();
200 octave_idx_type a_nc = a.cols (); 200 octave_idx_type a_nc = a.cols ();
201 201
202 if (a_nr != len) 202 if (a_nr != len)
203 err_nonconformant ("operator *", 1, len, a_nr, a_nc); 203 err_nonconformant ("operator *", 1, len, a_nr, a_nc);
204
205 if (len == 0)
206 retval.resize (a_nc, 0.0);
204 else 207 else
205 { 208 {
206 if (len == 0) 209 // Transpose A to form A'*x == (x'*A)'
207 retval.resize (a_nc, 0.0); 210
208 else 211 octave_idx_type ld = a_nr;
209 { 212
210 // Transpose A to form A'*x == (x'*A)' 213 retval.resize (a_nc);
211 214 double *y = retval.fortran_vec ();
212 octave_idx_type ld = a_nr; 215
213 216 F77_XFCN (dgemv, DGEMV, (F77_CONST_CHAR_ARG2 ("T", 1),
214 retval.resize (a_nc); 217 a_nr, a_nc, 1.0, a.data (),
215 double *y = retval.fortran_vec (); 218 ld, v.data (), 1, 0.0, y, 1
216 219 F77_CHAR_ARG_LEN (1)));
217 F77_XFCN (dgemv, DGEMV, (F77_CONST_CHAR_ARG2 ("T", 1),
218 a_nr, a_nc, 1.0, a.data (),
219 ld, v.data (), 1, 0.0, y, 1
220 F77_CHAR_ARG_LEN (1)));
221 }
222 } 220 }
223 221
224 return retval; 222 return retval;
225 } 223 }
226 224
322 320
323 octave_idx_type a_len = a.numel (); 321 octave_idx_type a_len = a.numel ();
324 322
325 if (len != a_len) 323 if (len != a_len)
326 err_nonconformant ("operator *", len, a_len); 324 err_nonconformant ("operator *", len, a_len);
327 else if (len != 0) 325
326 if (len != 0)
328 F77_FUNC (xddot, XDDOT) (len, v.data (), 1, a.data (), 1, retval); 327 F77_FUNC (xddot, XDDOT) (len, v.data (), 1, a.data (), 1, retval);
329 328
330 return retval; 329 return retval;
331 } 330 }
332 331