1 // SPDX-License-Identifier: Apache-2.0 2 // 3 // Copyright 2008-2016 Conrad Sanderson (http://conradsanderson.id.au) 4 // Copyright 2008-2016 National ICT Australia (NICTA) 5 // 6 // Licensed under the Apache License, Version 2.0 (the "License"); 7 // you may not use this file except in compliance with the License. 8 // You may obtain a copy of the License at 9 // http://www.apache.org/licenses/LICENSE-2.0 10 // 11 // Unless required by applicable law or agreed to in writing, software 12 // distributed under the License is distributed on an "AS IS" BASIS, 13 // WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. 14 // See the License for the specific language governing permissions and 15 // limitations under the License. 16 // ------------------------------------------------------------------------ 17 18 19 namespace newarp 20 { 21 22 23 //! Define matrix operations on existing matrix objects 24 template<typename eT> 25 class SparseGenRealShiftSolve 26 { 27 private: 28 29 #if defined(ARMA_USE_SUPERLU) 30 // The following objects are read-only in perform_op() 31 mutable superlu_supermatrix_wrangler l; 32 mutable superlu_supermatrix_wrangler u; 33 mutable superlu_array_wrangler<int> perm_c; 34 mutable superlu_array_wrangler<int> perm_r; 35 #endif 36 37 38 public: 39 40 bool valid = false; 41 42 const uword n_rows; // number of rows of the underlying matrix 43 const uword n_cols; // number of columns of the underlying matrix 44 45 inline SparseGenRealShiftSolve(const SpMat<eT>& mat_obj, const eT shift); 46 47 inline void perform_op(eT* x_in, eT* y_out) const; 48 }; 49 50 51 } // namespace newarp 52