Mercurial > octave
diff libinterp/corefcn/__ilu__.cc @ 31189:820d2c802247 stable
fortran_vec: Fix instances of method chaining that may cause dangling pointer.
* libinterp/corefcn/quandcc.cc (Fquadcc): Replace method chaining with a
separate line for NDArray object to avoid going out of scope.
* libinterp/corefcn/__ilu__.cc (F__iluc__): Obtain RowVector and ColumnVector
objects in separate line to remain in scope.
* libinterp/corefcn/__ichol__.cc (F__icholt__): Obtain RowVector object in
separate line to remain in scope.
For more details: https://octave.discourse.group/t/3132
author | Maged Rifaat <magedrifaat@gmail.com> |
---|---|
date | Mon, 15 Aug 2022 00:22:55 +0200 |
parents | 796f54d4ddbf |
children | e88a07dec498 |
line wrap: on
line diff
--- a/libinterp/corefcn/__ilu__.cc Fri Aug 12 15:56:50 2022 +0200 +++ b/libinterp/corefcn/__ilu__.cc Mon Aug 15 00:22:55 2022 +0200 @@ -483,9 +483,11 @@ SparseMatrix sm_l = Ftril (ovl (sm, -1))(0).sparse_matrix_value (); SparseMatrix U, L; + RowVector sm_col_norms = xcolnorms (sm); + ColumnVector sm_row_norms = xrownorms (sm); ilu_crout <SparseMatrix, double> (sm_l, sm_u, L, U, - xcolnorms (sm).fortran_vec (), - xrownorms (sm).fortran_vec (), + sm_col_norms.fortran_vec (), + sm_row_norms.fortran_vec (), droptol, milu); SparseMatrix speye (DiagMatrix (L.cols (), L.cols (), 1.0));