octave: liboctave/operators/Sparse-perm-op-defs.h annotate

annotate liboctave/operators/Sparse-perm-op-defs.h @ 33248:7f73e4805a1f

replace calls to assert with liboctave_panic functions liboctave * lo-error.h (liboctave_panic_impossible, liboctave_panic_if, liboctave_panic_unless): New macros. Replace all calls to assert with panic_if, panic_impossible, or panic_unless as appropriate. Include "lo-error.h" as needed. Don't include <cassert>. Affected files: Array-base.cc, Array-util.cc, Array.h, CSparse.cc, DiagArray2.cc, 1 DiagArray2.h, Sparse.cc, Sparse.h, dim-vector.h, idx-vector.cc, idx-vector.h, CollocWt.cc, Quad.cc, oct-rand.cc, qrp.cc, svd.cc, Sparse-perm-op-defs.h, and oct-sort.cc.

author	John W. Eaton <jwe@octave.org>
date	Sun, 24 Mar 2024 18:12:34 -0400
parents	2e484f9f1f18
children

rev	line source
27923 bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	1 ////////////////////////////////////////////////////////////////////////
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	2 //
32632 2e484f9f1f18 maint: update Octave Project Developers copyright for the new year John W. Eaton <jwe@octave.org> parents: 31706 diff changeset	3 // Copyright (C) 2009-2024 The Octave Project Developers
27923 bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	4 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	5 // See the file COPYRIGHT.md in the top-level directory of this
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	6 // distribution or <https://octave.org/copyright/>.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	7 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	8 // This file is part of Octave.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	9 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	10 // Octave is free software: you can redistribute it and/or modify it
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	11 // under the terms of the GNU General Public License as published by
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	12 // the Free Software Foundation, either version 3 of the License, or
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	13 // (at your option) any later version.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	14 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	15 // Octave is distributed in the hope that it will be useful, but
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	16 // WITHOUT ANY WARRANTY; without even the implied warranty of
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	17 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	18 // GNU General Public License for more details.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	19 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	20 // You should have received a copy of the GNU General Public License
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	21 // along with Octave; see the file COPYING. If not, see
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	22 // <https://www.gnu.org/licenses/>.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	23 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	24 ////////////////////////////////////////////////////////////////////////
8968 91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	25
20791 f7084eae3318 maint: Use Octave coding conventions for #if statements. Rik <rik@octave.org> parents: 19697 diff changeset	26 #if ! defined (octave_Sparse_perm_op_defs_h)
17822 ebb3ef964372 maint: Use common #define syntax "octave_filename_h" in h_files. Rik <rik@octave.org> parents: 17744 diff changeset	27 #define octave_Sparse_perm_op_defs_h 1
8968 91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	28
21244 1473547f50f5 include octave-config.h in public header files John W. Eaton <jwe@octave.org> parents: 21129 diff changeset	29 #include "octave-config.h"
1473547f50f5 include octave-config.h in public header files John W. Eaton <jwe@octave.org> parents: 21129 diff changeset	30
23442 53f5f8231c37 allow most header files to be compiled separately John W. Eaton <jwe@octave.org> parents: 23220 diff changeset	31 #include "PermMatrix.h"
53f5f8231c37 allow most header files to be compiled separately John W. Eaton <jwe@octave.org> parents: 23220 diff changeset	32 #include "lo-array-errwarn.h"
53f5f8231c37 allow most header files to be compiled separately John W. Eaton <jwe@octave.org> parents: 23220 diff changeset	33 #include "oct-locbuf.h"
53f5f8231c37 allow most header files to be compiled separately John W. Eaton <jwe@octave.org> parents: 23220 diff changeset	34 #include "oct-sort.h"
53f5f8231c37 allow most header files to be compiled separately John W. Eaton <jwe@octave.org> parents: 23220 diff changeset	35 #include "quit.h"
53f5f8231c37 allow most header files to be compiled separately John W. Eaton <jwe@octave.org> parents: 23220 diff changeset	36
8968 91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	37 // Matrix multiplication
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	38
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	39 template <typename SM>
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	40 SM octinternal_do_mul_colpm_sm (const octave_idx_type *pcol, const SM& a)
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	41 // Relabel the rows according to pcol.
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	42 {
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	43 const octave_idx_type nr = a.rows ();
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	44 const octave_idx_type nc = a.cols ();
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	45 const octave_idx_type nent = a.nnz ();
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	46 SM r (nr, nc, nent);
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	47
12424 afb65a7bc065 update NEWS John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	48 octave_sort<octave_idx_type> sort;
afb65a7bc065 update NEWS John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	49
8968 91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	50 for (octave_idx_type j = 0; j <= nc; ++j)
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	51 r.xcidx (j) = a.cidx (j);
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	52
12424 afb65a7bc065 update NEWS John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	53 for (octave_idx_type j = 0; j < nc; j++)
afb65a7bc065 update NEWS John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	54 {
afb65a7bc065 update NEWS John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	55 octave_quit ();
13141 e81ddf9cacd5 maint: untabify and remove trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 12424 diff changeset	56
15018 3d8ace26c5b4 maint: Use Octave coding conventions for cuddled parentheses in liboctave/. Rik <rik@octave.org> parents: 14846 diff changeset	57 OCTAVE_LOCAL_BUFFER (octave_idx_type, sidx, r.xcidx (j+1) - r.xcidx (j));
3d8ace26c5b4 maint: Use Octave coding conventions for cuddled parentheses in liboctave/. Rik <rik@octave.org> parents: 14846 diff changeset	58 for (octave_idx_type i = r.xcidx (j), ii = 0; i < r.xcidx (j+1); i++)
12424 afb65a7bc065 update NEWS John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	59 {
afb65a7bc065 update NEWS John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	60 sidx[ii++]=i;
afb65a7bc065 update NEWS John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	61 r.xridx (i) = pcol[a.ridx (i)];
afb65a7bc065 update NEWS John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	62 }
15018 3d8ace26c5b4 maint: Use Octave coding conventions for cuddled parentheses in liboctave/. Rik <rik@octave.org> parents: 14846 diff changeset	63 sort.sort (r.xridx () + r.xcidx (j), sidx, r.xcidx (j+1) - r.xcidx (j));
3d8ace26c5b4 maint: Use Octave coding conventions for cuddled parentheses in liboctave/. Rik <rik@octave.org> parents: 14846 diff changeset	64 for (octave_idx_type i = r.xcidx (j), ii = 0; i < r.xcidx (j+1); i++)
3d8ace26c5b4 maint: Use Octave coding conventions for cuddled parentheses in liboctave/. Rik <rik@octave.org> parents: 14846 diff changeset	65 r.xdata (i) = a.data (sidx[ii++]);
12424 afb65a7bc065 update NEWS John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	66 }
afb65a7bc065 update NEWS John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	67
8968 91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	68 return r;
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	69 }
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	70
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	71 template <typename SM>
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	72 SM octinternal_do_mul_pm_sm (const PermMatrix& p, const SM& a)
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	73 {
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	74 const octave_idx_type nr = a.rows ();
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	75 if (p.cols () != nr)
22402 4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22327 diff changeset	76 octave::err_nonconformant ("operator *",
4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22327 diff changeset	77 p.rows (), p.cols (), a.rows (), a.cols ());
8968 91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	78
18849 aa9ca67f09fb make all permutation matrices column permutations (bug #42418) David Spies <dnspies@gmail.com> parents: 17822 diff changeset	79 return octinternal_do_mul_colpm_sm (p.col_perm_vec ().data (), a);
8968 91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	80 }
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	81
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	82 template <typename SM>
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	83 SM octinternal_do_mul_sm_rowpm (const SM& a, const octave_idx_type *prow)
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	84 // For a row permutation, iterate across the source a and stuff the
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	85 // results into the correct destination column in r.
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	86 {
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	87 const octave_idx_type nr = a.rows ();
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	88 const octave_idx_type nc = a.cols ();
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	89 const octave_idx_type nent = a.nnz ();
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	90 SM r (nr, nc, nent);
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	91
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	92 for (octave_idx_type j_src = 0; j_src < nc; ++j_src)
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	93 r.xcidx (prow[j_src]) = a.cidx (j_src+1) - a.cidx (j_src);
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	94 octave_idx_type k = 0;
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	95 for (octave_idx_type j = 0; j < nc; ++j)
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	96 {
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	97 const octave_idx_type tmp = r.xcidx (j);
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	98 r.xcidx (j) = k;
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	99 k += tmp;
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	100 }
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	101 r.xcidx (nc) = nent;
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	102
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	103 octave_idx_type k_src = 0;
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	104 for (octave_idx_type j_src = 0; j_src < nc; ++j_src)
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	105 {
10142 829e69ec3110 make OCTAVE_QUIT a function Jaroslav Hajek <highegg@gmail.com> parents: 8968 diff changeset	106 octave_quit ();
8968 91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	107 const octave_idx_type j = prow[j_src];
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	108 const octave_idx_type kend_src = a.cidx (j_src + 1);
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	109 for (k = r.xcidx (j); k_src < kend_src; ++k, ++k_src)
10312 cbc402e64d83 untabify liboctave header files John W. Eaton <jwe@octave.org> parents: 10142 diff changeset	110 {
cbc402e64d83 untabify liboctave header files John W. Eaton <jwe@octave.org> parents: 10142 diff changeset	111 r.xridx (k) = a.ridx (k_src);
cbc402e64d83 untabify liboctave header files John W. Eaton <jwe@octave.org> parents: 10142 diff changeset	112 r.xdata (k) = a.data (k_src);
cbc402e64d83 untabify liboctave header files John W. Eaton <jwe@octave.org> parents: 10142 diff changeset	113 }
8968 91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	114 }
33248 7f73e4805a1f replace calls to assert with liboctave_panic functions liboctave John W. Eaton <jwe@octave.org> parents: 32632 diff changeset	115 liboctave_panic_unless (k_src == nent);
8968 91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	116
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	117 return r;
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	118 }
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	119
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	120 template <typename SM>
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	121 SM octinternal_do_mul_sm_colpm (const SM& a, const octave_idx_type *pcol)
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	122 // For a column permutation, iterate across the destination r and pull
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	123 // data from the correct column of a.
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	124 {
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	125 const octave_idx_type nr = a.rows ();
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	126 const octave_idx_type nc = a.cols ();
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	127 const octave_idx_type nent = a.nnz ();
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	128 SM r (nr, nc, nent);
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	129
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	130 for (octave_idx_type j = 0; j < nc; ++j)
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	131 {
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	132 const octave_idx_type j_src = pcol[j];
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	133 r.xcidx (j+1) = r.xcidx (j) + (a.cidx (j_src+1) - a.cidx (j_src));
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	134 }
33248 7f73e4805a1f replace calls to assert with liboctave_panic functions liboctave John W. Eaton <jwe@octave.org> parents: 32632 diff changeset	135 liboctave_panic_unless (r.xcidx (nc) == nent);
8968 91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	136
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	137 octave_idx_type k = 0;
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	138 for (octave_idx_type j = 0; j < nc; ++j)
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	139 {
10142 829e69ec3110 make OCTAVE_QUIT a function Jaroslav Hajek <highegg@gmail.com> parents: 8968 diff changeset	140 octave_quit ();
8968 91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	141 const octave_idx_type j_src = pcol[j];
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	142 octave_idx_type k_src;
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	143 const octave_idx_type kend_src = a.cidx (j_src + 1);
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	144 for (k_src = a.cidx (j_src); k_src < kend_src; ++k_src, ++k)
10312 cbc402e64d83 untabify liboctave header files John W. Eaton <jwe@octave.org> parents: 10142 diff changeset	145 {
cbc402e64d83 untabify liboctave header files John W. Eaton <jwe@octave.org> parents: 10142 diff changeset	146 r.xridx (k) = a.ridx (k_src);
cbc402e64d83 untabify liboctave header files John W. Eaton <jwe@octave.org> parents: 10142 diff changeset	147 r.xdata (k) = a.data (k_src);
cbc402e64d83 untabify liboctave header files John W. Eaton <jwe@octave.org> parents: 10142 diff changeset	148 }
8968 91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	149 }
33248 7f73e4805a1f replace calls to assert with liboctave_panic functions liboctave John W. Eaton <jwe@octave.org> parents: 32632 diff changeset	150 liboctave_panic_unless (k == nent);
8968 91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	151
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	152 return r;
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	153 }
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	154
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	155 template <typename SM>
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	156 SM octinternal_do_mul_sm_pm (const SM& a, const PermMatrix& p)
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	157 {
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	158 const octave_idx_type nc = a.cols ();
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	159 if (p.rows () != nc)
22402 4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22327 diff changeset	160 octave::err_nonconformant ("operator *",
4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22327 diff changeset	161 a.rows (), a.cols (), p.rows (), p.cols ());
8968 91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	162
18849 aa9ca67f09fb make all permutation matrices column permutations (bug #42418) David Spies <dnspies@gmail.com> parents: 17822 diff changeset	163 return octinternal_do_mul_sm_colpm (a, p.col_perm_vec ().data ());
8968 91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	164 }
91d53dc37f79 Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations. Jason Riedy <jason@acm.org> parents: diff changeset	165
20945 9d9270e2f98f eliminate comments after preprocessor conditionals John W. Eaton <jwe@octave.org> parents: 20791 diff changeset	166 #endif

Mercurial > octave

annotate liboctave/operators/Sparse-perm-op-defs.h @ 33248:7f73e4805a1f