Diana
0.8.3
|
00001 /* ----------------------------------------------------------------------------- 00002 * Diana process modelling, simulation and analysis software 00003 * Copyright (c) 2005, Michael Krasnyk 00004 * All rights reserved. 00005 * 00006 * This file is a part of Diana process modelling, simulation and analysis software 00007 * 00008 * Diana is free software; you can redistribute it and/or modify it 00009 * under the terms of the GNU General Public License as published 00010 * by the Free Software Foundation (see accompanying file LICENSE) 00011 * ----------------------------------------------------------------------------- 00012 * $Id: DianaSparseSolver.hpp 8160 2009-03-16 18:36:30Z miha $ 00013 * ----------------------------------------------------------------------------- 00014 * Description: 00015 */ 00016 00017 #ifndef DIANA_SPARSE_SOLVER_HPP 00018 #define DIANA_SPARSE_SOLVER_HPP 00019 00020 #include <vector> 00021 #include <Basic/DianaSparseArray.hpp> 00022 #include "IDianaLinearAlgebraSolver.hpp" 00023 00024 namespace Diana { 00025 00027 class DianaSparseSolver : virtual public IDianaLinearAlgebraSolver { 00028 public: 00029 DianaSparseSolver(); 00030 virtual ~DianaSparseSolver(); 00031 00032 DianaSparseSolver(const DianaSparseArray& spar) 00033 throw (Common::Error::ECapeUnknown, Common::Error::ECapeBadArgument, 00034 Common::Error::ECapeOutOfBounds, Common::Error::ECapeSolvingError, 00035 Common::Error::ECapeNoMemory); 00036 00037 DianaSparseSolver& operator=(const DianaSparseArray& spar) 00038 throw (Common::Error::ECapeUnknown, Common::Error::ECapeBadArgument, 00039 Common::Error::ECapeOutOfBounds, Common::Error::ECapeSolvingError, 00040 Common::Error::ECapeNoMemory); 00041 00042 void decompose() 00043 throw (Common::Error::ECapeUnknown, Common::Error::ECapeBadArgument, 00044 Common::Error::ECapeOutOfBounds, Common::Error::ECapeSolvingError, 00045 Common::Error::ECapeNoMemory); 00046 00047 void solve(const Common::Types::CapeArrayDouble& rhs, 00048 Common::Types::CapeArrayDouble& sol, 00049 Common::Types::CapeBoolean transpose=false) 00050 throw (Common::Error::ECapeUnknown, Common::Error::ECapeBadArgument, 00051 Common::Error::ECapeOutOfBounds, Common::Error::ECapeSolvingError, 00052 Common::Error::ECapeNoMemory); 00053 00054 void solve(int n, const double* rhs, double* sol, 00055 Common::Types::CapeBoolean transpose=false) 00056 throw (Common::Error::ECapeUnknown, Common::Error::ECapeBadArgument, 00057 Common::Error::ECapeOutOfBounds, Common::Error::ECapeSolvingError, 00058 Common::Error::ECapeNoMemory); 00059 00060 void mulvec(const Common::Types::CapeArrayDouble& vec, 00061 Common::Types::CapeArrayDouble& res, 00062 Common::Types::CapeBoolean transpose=false) 00063 throw (Common::Error::ECapeUnknown, Common::Error::ECapeBadArgument, 00064 Common::Error::ECapeOutOfBounds, Common::Error::ECapeSolvingError, 00065 Common::Error::ECapeNoMemory); 00066 00067 void mulvec(int nvec, const double* vec, int nres, double* res, 00068 Common::Types::CapeBoolean transpose=false) 00069 throw (Common::Error::ECapeUnknown, Common::Error::ECapeBadArgument, 00070 Common::Error::ECapeOutOfBounds, Common::Error::ECapeSolvingError, 00071 Common::Error::ECapeNoMemory); 00072 private: 00073 void init_matrix(const DianaSparseArray& spar) 00074 throw (Common::Error::ECapeUnknown, Common::Error::ECapeBadArgument, 00075 Common::Error::ECapeOutOfBounds, Common::Error::ECapeSolvingError, 00076 Common::Error::ECapeNoMemory); 00077 00078 void lusymb() 00079 throw (Common::Error::ECapeUnknown, Common::Error::ECapeBadArgument, 00080 Common::Error::ECapeOutOfBounds,Common::Error::ECapeSolvingError, 00081 Common::Error::ECapeNoMemory); 00082 00083 private: 00084 Common::Types::CapeBoolean bNewStructure; 00085 Common::Types::CapeBoolean bNewValues; 00086 DianaSparseArray::crc_type crc32sa; 00087 00088 void* Symbolic; 00089 void* Numeric; 00090 long nRows; 00091 long nCols; 00092 long NZ; 00093 std::vector<double> Info; 00094 std::vector<double> Control; 00095 std::vector<double> WorkD; 00096 std::vector<double> Ax; 00097 std::vector<long> WorkI; 00098 std::vector<long> Ti; 00099 std::vector<long> Tj; 00100 std::vector<long> Ap; 00101 std::vector<long> Ai; 00102 std::vector<long> Mapping; 00103 }; 00104 }; 00105 00106 #endif // DIANA_SPARSE_SOLVER_HPP