Mercurial > octave-nkf

comparison src/DLD-FUNCTIONS/kron.cc @ 14854:5ae9f0f77635

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
maint: Use Octave coding conventions for coddling parenthis is DLD-FUNCTIONS directory * __fltk_uigetfile__.cc, __glpk__.cc, __init_fltk__.cc, __magick_read__.cc, besselj.cc, bsxfun.cc, ccolamd.cc, cellfun.cc, chol.cc, colamd.cc, daspk.cc, dasrt.cc, dassl.cc, dmperm.cc, fft.cc, filter.cc, find.cc, gcd.cc, kron.cc, lsode.cc, lu.cc, luinc.cc, quad.cc, quadcc.cc, rand.cc, regexp.cc, schur.cc, str2double.cc, symbfact.cc, symrcm.cc, tril.cc, urlwrite.cc: Use Octave coding conventions for coddling parenthis is DLD-FUNCTIONS directory.
author Rik <octave@nomad.inbox5.com>
date Mon, 09 Jul 2012 13:01:49 -0700
parents 460a3c6d8bf1
children
comparison
equal deleted inserted replaced
14853:72b8b39e12be 14854:5ae9f0f77635
86 86
87 for (octave_idx_type ja = 0; ja < dla; ja++) 87 for (octave_idx_type ja = 0; ja < dla; ja++)
88 for (octave_idx_type jb = 0; jb < ncb; jb++) 88 for (octave_idx_type jb = 0; jb < ncb; jb++)
89 { 89 {
90 octave_quit (); 90 octave_quit ();
91 mx_inline_mul (nrb, &c.xelem(ja*nrb, ja*ncb + jb), a.dgelem (ja), b.data () + nrb*jb); 91 mx_inline_mul (nrb, &c.xelem (ja*nrb, ja*ncb + jb), a.dgelem (ja), b.data () + nrb*jb);
92 } 92 }
93 93
94 return c; 94 return c;
95 } 95 }
96 96
108 for (octave_idx_type Bj = 0; Bj < B.columns (); Bj++) 108 for (octave_idx_type Bj = 0; Bj < B.columns (); Bj++)
109 { 109 {
110 octave_quit (); 110 octave_quit ();
111 for (octave_idx_type Ai = A.cidx (Aj); Ai < A.cidx (Aj+1); Ai++) 111 for (octave_idx_type Ai = A.cidx (Aj); Ai < A.cidx (Aj+1); Ai++)
112 { 112 {
113 octave_idx_type Ci = A.ridx(Ai) * B.rows (); 113 octave_idx_type Ci = A.ridx (Ai) * B.rows ();
114 const T v = A.data (Ai); 114 const T v = A.data (Ai);
115 115
116 for (octave_idx_type Bi = B.cidx (Bj); Bi < B.cidx (Bj+1); Bi++) 116 for (octave_idx_type Bi = B.cidx (Bj); Bi < B.cidx (Bj+1); Bi++)
117 { 117 {
118 C.data (idx) = v * B.data (Bi); 118 C.data (idx) = v * B.data (Bi);