Mercurial > hg > octave-jordi

comparison src/DLD-FUNCTIONS/symbfact.cc @ 11586:12df7854fa7c

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
strip trailing whitespace from source files
author John W. Eaton <jwe@octave.org>
date Thu, 20 Jan 2011 17:24:59 -0500
parents 01f703952eff
children 990762e784fe
comparison
equal deleted inserted replaced
11585:1473d0cf86d2 11586:12df7854fa7c
269 269
270 // count the total number of entries in L 270 // count the total number of entries in L
271 octave_idx_type lnz = 0 ; 271 octave_idx_type lnz = 0 ;
272 for (octave_idx_type j = 0 ; j < n ; j++) 272 for (octave_idx_type j = 0 ; j < n ; j++)
273 lnz += ColCount [j] ; 273 lnz += ColCount [j] ;
274 274
275 275
276 // allocate the output matrix L (pattern-only) 276 // allocate the output matrix L (pattern-only)
277 SparseBoolMatrix L (n, n, lnz); 277 SparseBoolMatrix L (n, n, lnz);
278 278
279 // initialize column pointers 279 // initialize column pointers
290 octave_idx_type *W = First; 290 octave_idx_type *W = First;
291 for (octave_idx_type j = 0 ; j < n ; j++) 291 for (octave_idx_type j = 0 ; j < n ; j++)
292 W [j] = L.xcidx(j); 292 W [j] = L.xcidx(j);
293 293
294 // get workspace for computing one row of L 294 // get workspace for computing one row of L
295 cholmod_sparse *R = cholmod_allocate_sparse (n, 1, n, false, true, 295 cholmod_sparse *R = cholmod_allocate_sparse (n, 1, n, false, true,
296 0, CHOLMOD_PATTERN, cm); 296 0, CHOLMOD_PATTERN, cm);
297 octave_idx_type *Rp = static_cast<octave_idx_type *>(R->p); 297 octave_idx_type *Rp = static_cast<octave_idx_type *>(R->p);
298 octave_idx_type *Ri = static_cast<octave_idx_type *>(R->i); 298 octave_idx_type *Ri = static_cast<octave_idx_type *>(R->i);
299 299
300 // compute L one row at a time 300 // compute L one row at a time