# HG changeset patch # User Rik # Date 1630292306 25200 # Node ID 7215f2324c8aa91c4beef4d844ef16268453e15e # Parent 2fbb98fb1b25c18cb73d752e22a211b9e5788f96 maint: use "m_" prefix for member variables in class sparse_lu. * sparse-lu.cc, sparse-lu.h: use "m_" prefix for member variables in class sparse_lu. diff -r 2fbb98fb1b25 -r 7215f2324c8a liboctave/numeric/sparse-lu.cc --- a/liboctave/numeric/sparse-lu.cc Sun Aug 29 19:51:50 2021 -0700 +++ b/liboctave/numeric/sparse-lu.cc Sun Aug 29 19:58:26 2021 -0700 @@ -407,7 +407,7 @@ template sparse_lu::sparse_lu (const lu_type& a, const Matrix& piv_thres, bool scale) - : Lfact (), Ufact (), Rfact (), cond (0), P (), Q () + : m_Lfact (), m_Ufact (), m_Rfact (), m_cond (0), m_P (), m_Q () { #if defined (HAVE_UMFPACK) octave_idx_type nr = a.rows (); @@ -489,7 +489,7 @@ &Numeric, control, info); umfpack_free_symbolic (&Symbolic); - cond = Info (UMFPACK_RCOND); + m_cond = Info (UMFPACK_RCOND); if (status < 0) { @@ -523,39 +523,39 @@ octave_idx_type n_inner = (nr < nc ? nr : nc); if (lnz < 1) - Lfact = lu_type (n_inner, nr, - static_cast (1)); + m_Lfact = lu_type (n_inner, nr, + static_cast (1)); else - Lfact = lu_type (n_inner, nr, lnz); + m_Lfact = lu_type (n_inner, nr, lnz); - octave_idx_type *Ltp = Lfact.cidx (); - octave_idx_type *Ltj = Lfact.ridx (); - lu_elt_type *Ltx = Lfact.data (); + octave_idx_type *Ltp = m_Lfact.cidx (); + octave_idx_type *Ltj = m_Lfact.ridx (); + lu_elt_type *Ltx = m_Lfact.data (); if (unz < 1) - Ufact = lu_type (n_inner, nc, - static_cast (1)); + m_Ufact = lu_type (n_inner, nc, + static_cast (1)); else - Ufact = lu_type (n_inner, nc, unz); + m_Ufact = lu_type (n_inner, nc, unz); - octave_idx_type *Up = Ufact.cidx (); - octave_idx_type *Uj = Ufact.ridx (); - lu_elt_type *Ux = Ufact.data (); + octave_idx_type *Up = m_Ufact.cidx (); + octave_idx_type *Uj = m_Ufact.ridx (); + lu_elt_type *Ux = m_Ufact.data (); - Rfact = SparseMatrix (nr, nr, nr); + m_Rfact = SparseMatrix (nr, nr, nr); for (octave_idx_type i = 0; i < nr; i++) { - Rfact.xridx (i) = i; - Rfact.xcidx (i) = i; + m_Rfact.xridx (i) = i; + m_Rfact.xcidx (i) = i; } - Rfact.xcidx (nr) = nr; - double *Rx = Rfact.data (); + m_Rfact.xcidx (nr) = nr; + double *Rx = m_Rfact.data (); - P.resize (dim_vector (nr, 1)); - octave_idx_type *p = P.fortran_vec (); + m_P.resize (dim_vector (nr, 1)); + octave_idx_type *p = m_P.fortran_vec (); - Q.resize (dim_vector (nc, 1)); - octave_idx_type *q = Q.fortran_vec (); + m_Q.resize (dim_vector (nc, 1)); + octave_idx_type *q = m_Q.fortran_vec (); octave_idx_type do_recip; status = umfpack_get_numeric (Ltp, Ltj, Ltx, @@ -575,22 +575,22 @@ } else { - Lfact = Lfact.transpose (); + m_Lfact = m_Lfact.transpose (); if (do_recip) for (octave_idx_type i = 0; i < nr; i++) Rx[i] = 1.0 / Rx[i]; umfpack_report_matrix (nr, n_inner, - Lfact.cidx (), - Lfact.ridx (), - Lfact.data (), + m_Lfact.cidx (), + m_Lfact.ridx (), + m_Lfact.data (), static_cast (1), control); umfpack_report_matrix (n_inner, nc, - Ufact.cidx (), - Ufact.ridx (), - Ufact.data (), + m_Ufact.cidx (), + m_Ufact.ridx (), + m_Ufact.data (), static_cast (1), control); umfpack_report_perm (nr, p, control); @@ -620,7 +620,7 @@ const Matrix& piv_thres, bool scale, bool FixedQ, double droptol, bool milu, bool udiag) - : Lfact (), Ufact (), Rfact (), cond (0), P (), Q () + : m_Lfact (), m_Ufact (), m_Rfact (), m_cond (0), m_P (), m_Q () { #if defined (HAVE_UMFPACK) @@ -727,7 +727,7 @@ &Numeric, control, info); umfpack_free_symbolic (&Symbolic); - cond = Info (UMFPACK_RCOND); + m_cond = Info (UMFPACK_RCOND); if (status < 0) { @@ -761,39 +761,39 @@ octave_idx_type n_inner = (nr < nc ? nr : nc); if (lnz < 1) - Lfact = lu_type (n_inner, nr, - static_cast (1)); + m_Lfact = lu_type (n_inner, nr, + static_cast (1)); else - Lfact = lu_type (n_inner, nr, lnz); + m_Lfact = lu_type (n_inner, nr, lnz); - octave_idx_type *Ltp = Lfact.cidx (); - octave_idx_type *Ltj = Lfact.ridx (); - lu_elt_type *Ltx = Lfact.data (); + octave_idx_type *Ltp = m_Lfact.cidx (); + octave_idx_type *Ltj = m_Lfact.ridx (); + lu_elt_type *Ltx = m_Lfact.data (); if (unz < 1) - Ufact = lu_type (n_inner, nc, - static_cast (1)); + m_Ufact = lu_type (n_inner, nc, + static_cast (1)); else - Ufact = lu_type (n_inner, nc, unz); + m_Ufact = lu_type (n_inner, nc, unz); - octave_idx_type *Up = Ufact.cidx (); - octave_idx_type *Uj = Ufact.ridx (); - lu_elt_type *Ux = Ufact.data (); + octave_idx_type *Up = m_Ufact.cidx (); + octave_idx_type *Uj = m_Ufact.ridx (); + lu_elt_type *Ux = m_Ufact.data (); - Rfact = SparseMatrix (nr, nr, nr); + m_Rfact = SparseMatrix (nr, nr, nr); for (octave_idx_type i = 0; i < nr; i++) { - Rfact.xridx (i) = i; - Rfact.xcidx (i) = i; + m_Rfact.xridx (i) = i; + m_Rfact.xcidx (i) = i; } - Rfact.xcidx (nr) = nr; - double *Rx = Rfact.data (); + m_Rfact.xcidx (nr) = nr; + double *Rx = m_Rfact.data (); - P.resize (dim_vector (nr, 1)); - octave_idx_type *p = P.fortran_vec (); + m_P.resize (dim_vector (nr, 1)); + octave_idx_type *p = m_P.fortran_vec (); - Q.resize (dim_vector (nc, 1)); - octave_idx_type *q = Q.fortran_vec (); + m_Q.resize (dim_vector (nc, 1)); + octave_idx_type *q = m_Q.fortran_vec (); octave_idx_type do_recip; status = umfpack_get_numeric (Ltp, Ltj, Ltx, @@ -813,22 +813,22 @@ } else { - Lfact = Lfact.transpose (); + m_Lfact = m_Lfact.transpose (); if (do_recip) for (octave_idx_type i = 0; i < nr; i++) Rx[i] = 1.0 / Rx[i]; umfpack_report_matrix (nr, n_inner, - Lfact.cidx (), - Lfact.ridx (), - Lfact.data (), + m_Lfact.cidx (), + m_Lfact.ridx (), + m_Lfact.data (), static_cast (1), control); umfpack_report_matrix (n_inner, nc, - Ufact.cidx (), - Ufact.ridx (), - Ufact.data (), + m_Ufact.cidx (), + m_Ufact.ridx (), + m_Ufact.data (), static_cast (1), control); umfpack_report_perm (nr, p, control); @@ -865,30 +865,30 @@ lu_type sparse_lu::Y (void) const { - octave_idx_type nr = Lfact.rows (); - octave_idx_type nz = Lfact.cols (); - octave_idx_type nc = Ufact.cols (); + octave_idx_type nr = m_Lfact.rows (); + octave_idx_type nz = m_Lfact.cols (); + octave_idx_type nc = m_Ufact.cols (); - lu_type Yout (nr, nc, Lfact.nnz () + Ufact.nnz () - (nr::Pr (void) const { - octave_idx_type nr = Lfact.rows (); + octave_idx_type nr = m_Lfact.rows (); SparseMatrix Pout (nr, nr, nr); for (octave_idx_type i = 0; i < nr; i++) { Pout.cidx (i) = i; - Pout.ridx (P (i)) = i; + Pout.ridx (m_P(i)) = i; Pout.data (i) = 1; } @@ -922,12 +922,12 @@ ColumnVector sparse_lu::Pr_vec (void) const { - octave_idx_type nr = Lfact.rows (); + octave_idx_type nr = m_Lfact.rows (); ColumnVector Pout (nr); for (octave_idx_type i = 0; i < nr; i++) - Pout.xelem (i) = static_cast (P(i) + 1); + Pout.xelem (i) = static_cast (m_P(i) + 1); return Pout; } @@ -936,21 +936,21 @@ PermMatrix sparse_lu::Pr_mat (void) const { - return PermMatrix (P, false); + return PermMatrix (m_P, false); } template SparseMatrix sparse_lu::Pc (void) const { - octave_idx_type nc = Ufact.cols (); + octave_idx_type nc = m_Ufact.cols (); SparseMatrix Pout (nc, nc, nc); for (octave_idx_type i = 0; i < nc; i++) { Pout.cidx (i) = i; - Pout.ridx (i) = Q (i); + Pout.ridx (i) = m_Q(i); Pout.data (i) = 1; } @@ -963,12 +963,12 @@ ColumnVector sparse_lu::Pc_vec (void) const { - octave_idx_type nc = Ufact.cols (); + octave_idx_type nc = m_Ufact.cols (); ColumnVector Pout (nc); for (octave_idx_type i = 0; i < nc; i++) - Pout.xelem (i) = static_cast (Q(i) + 1); + Pout.xelem (i) = static_cast (m_Q(i) + 1); return Pout; } @@ -977,7 +977,7 @@ PermMatrix sparse_lu::Pc_mat (void) const { - return PermMatrix (Q, true); + return PermMatrix (m_Q, true); } // Instantiations we need. diff -r 2fbb98fb1b25 -r 7215f2324c8a liboctave/numeric/sparse-lu.h --- a/liboctave/numeric/sparse-lu.h Sun Aug 29 19:51:50 2021 -0700 +++ b/liboctave/numeric/sparse-lu.h Sun Aug 29 19:58:26 2021 -0700 @@ -53,7 +53,7 @@ typedef typename lu_type::element_type lu_elt_type; sparse_lu (void) - : Lfact (), Ufact (), Rfact (), cond (0), P (), Q () { } + : m_Lfact (), m_Ufact (), m_Rfact (), m_cond (0), m_P (), m_Q () { } OCTAVE_API sparse_lu (const lu_type& a, const Matrix& piv_thres = Matrix (), @@ -66,19 +66,19 @@ bool milu = false, bool udiag = false); sparse_lu (const sparse_lu& a) - : Lfact (a.Lfact), Ufact (a.Ufact), Rfact (), cond (a.cond), - P (a.P), Q (a.Q) + : m_Lfact (a.m_Lfact), m_Ufact (a.m_Ufact), m_Rfact (), + m_cond (a.m_cond), m_P (a.m_P), m_Q (a.m_Q) { } sparse_lu& operator = (const sparse_lu& a) { if (this != &a) { - Lfact = a.Lfact; - Ufact = a.Ufact; - cond = a.cond; - P = a.P; - Q = a.Q; + m_Lfact = a.m_Lfact; + m_Ufact = a.m_Ufact; + m_cond = a.m_cond; + m_P = a.m_P; + m_Q = a.m_Q; } return *this; @@ -86,11 +86,11 @@ virtual ~sparse_lu (void) = default; - lu_type L (void) const { return Lfact; } + lu_type L (void) const { return m_Lfact; } - lu_type U (void) const { return Ufact; } + lu_type U (void) const { return m_Ufact; } - SparseMatrix R (void) const { return Rfact; } + SparseMatrix R (void) const { return m_Rfact; } OCTAVE_API lu_type Y (void) const; @@ -106,22 +106,22 @@ OCTAVE_API PermMatrix Pr_mat (void) const; - const octave_idx_type * row_perm (void) const { return P.data (); } + const octave_idx_type * row_perm (void) const { return m_P.data (); } - const octave_idx_type * col_perm (void) const { return Q.data (); } + const octave_idx_type * col_perm (void) const { return m_Q.data (); } - double rcond (void) const { return cond; } + double rcond (void) const { return m_cond; } protected: - lu_type Lfact; - lu_type Ufact; - SparseMatrix Rfact; + lu_type m_Lfact; + lu_type m_Ufact; + SparseMatrix m_Rfact; - double cond; + double m_cond; - MArray P; - MArray Q; + MArray m_P; + MArray m_Q; }; } }