Mercurial > jwe > octave

comparison liboctave/array/fCRowVector.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 228b65504557
comparison
equal deleted inserted replaced
21119:90cd0f9442d5 21120:499b851fbfae
325 octave_idx_type a_nr = a.rows (); 325 octave_idx_type a_nr = a.rows ();
326 octave_idx_type a_nc = a.cols (); 326 octave_idx_type a_nc = a.cols ();
327 327
328 if (a_nr != len) 328 if (a_nr != len)
329 err_nonconformant ("operator *", 1, len, a_nr, a_nc); 329 err_nonconformant ("operator *", 1, len, a_nr, a_nc);
330
331 if (len == 0)
332 retval.resize (a_nc, 0.0);
330 else 333 else
331 { 334 {
332 if (len == 0) 335 // Transpose A to form A'*x == (x'*A)'
333 retval.resize (a_nc, 0.0); 336
334 else 337 octave_idx_type ld = a_nr;
335 { 338
336 // Transpose A to form A'*x == (x'*A)' 339 retval.resize (a_nc);
337 340 FloatComplex *y = retval.fortran_vec ();
338 octave_idx_type ld = a_nr; 341
339 342 F77_XFCN (cgemv, CGEMV, (F77_CONST_CHAR_ARG2 ("T", 1),
340 retval.resize (a_nc); 343 a_nr, a_nc, 1.0, a.data (),
341 FloatComplex *y = retval.fortran_vec (); 344 ld, v.data (), 1, 0.0, y, 1
342 345 F77_CHAR_ARG_LEN (1)));
343 F77_XFCN (cgemv, CGEMV, (F77_CONST_CHAR_ARG2 ("T", 1),
344 a_nr, a_nc, 1.0, a.data (),
345 ld, v.data (), 1, 0.0, y, 1
346 F77_CHAR_ARG_LEN (1)));
347 }
348 } 346 }
349 347
350 return retval; 348 return retval;
351 } 349 }
352 350
450 448
451 octave_idx_type a_len = a.numel (); 449 octave_idx_type a_len = a.numel ();
452 450
453 if (len != a_len) 451 if (len != a_len)
454 err_nonconformant ("operator *", len, a_len); 452 err_nonconformant ("operator *", len, a_len);
455 else if (len != 0) 453
454 if (len != 0)
456 F77_FUNC (xcdotu, XCDOTU) (len, v.data (), 1, a.data (), 1, retval); 455 F77_FUNC (xcdotu, XCDOTU) (len, v.data (), 1, a.data (), 1, retval);
457 456
458 return retval; 457 return retval;
459 } 458 }
460 459