Hi, Karl, thank you for the response! The stripped down version fails for me as well: compressed_matrix of size (3, 3) with 3 nonzeros: (0, 0) 1 (1, 1) 1 (2, 2) 1 [3](1,2,3) [3](nan,nan,nan) However with a little modification to my original code, where I omit the conversions to viennacl types and simply do this call: arma_result = viennacl::linalg::solve(arma_sparsemat, arma_rhs, viennacl::linalg::cg_tag()); I get the neccessary result. Going ahead and adding another element does not work so well...
Hi, Karl, thank you for the response! The stripped down version fails for me as well: compressed_matrix of size (3, 3) with 3 nonzeros: (0, 0) 1 (1, 1) 1 (2, 2) 1 [3](1,2,3) [3](nan,nan,nan) However with a little modification to my original code, where I omit the conversions to viennacl types and simply do this call: arma_result = viennacl::linalg::solve(arma_sparsemat, arma_rhs, viennacl::linalg::cg_tag()); I get the neccessary result. Going ahead and adding another element does not work so well...
Hi, Karl, thank you for the response! The stripped down version fails for me as well: compressed_matrix of size (3, 3) with 3 nonzeros: (0, 0) 1 (1, 1) 1 (2, 2) 1 [3](1,2,3) [3](nan,nan,nan) However with a little modification to my original code, where I omit the conversions to viennacl types and simply do this call: arma_result = viennacl::linalg::solve(arma_sparsemat, arma_rhs, viennacl::linalg::cg_tag()); I get the result: sparsemat [matrix size: 3x3; n_nonzero: 4; density: 44.44%] (0, 0) 1.0000...
Hi, Karl, thank you for the response! The stripped down version runs just fine, however with a little modification it does not: with one element added so that the solver has to iterate, the solution is wrong (and the solution actually diverges with more allowed iterations). Example output: sparsemat [matrix size: 3x3; n_nonzero: 4; density: 44.44%] (0, 0) 1.0000 (0, 1) 0.5000 (1, 1) 1.0000 (2, 2) 1.0000 rhs 1.0000 2.0000 3.0000 arma result 0.0216 2.0175 3.0262 My OS: Windows 7 Ultimate x64 (version...
Hi, Karl, thank you for the response! The stripped down version runs just fine, however with a little modification it does not: with one element added so that the solver has to iterate, the solution is wrong (and the solution actually diverges with more allowed iterations). Example output: sparsemat [matrix size: 3x3; n_nonzero: 4; density: 44.44%] (0, 0) 1.0000 (0, 1) 0.5000 (1, 1) 1.0000 (2, 2) 1.0000 rhs 1.0000 2.0000 3.0000 arma result 0.0216 2.0175 3.0262 My OS: Windows 7 Ultimate x64 (version...
Hi, Karl, thank you for the response! The stripped down version runs just fine, however with a little modification it does not: with one element added so that the solver has to iterate, the solution is wrong (and the solution actually diverges with more allowed iterations). Example output: sparsemat [matrix size: 3x3; n_nonzero: 4; density: 44.44%] (0, 0) 1.0000 (0, 1) 0.5000 (1, 1) 1.0000 (2, 2) 1.0000 rhs 1.0000 2.0000 3.0000 arma result 0.0216 2.0175 3.0262 My OS: Windows 7 Ultimate x64 (version...
Hi, Karl, thank you for the response! The stripped down version runs just fine, however with a little modification it does not: with one element added so that the solver has to iterate, the solution is wrong (and the solution actually diverges with more allowed iterations). Example output: sparsemat [matrix size: 3x3; n_nonzero: 4; density: 44.44%] (0, 0) 1.0000 (0, 1) 0.5000 (1, 1) 1.0000 (2, 2) 1.0000 rhs 1.0000 2.0000 3.0000 arma result 0.0216 2.0175 3.0262 My OS: Windows 7 Ultimate x64 (version...
Hi! I am trying to solve a simple linear system of linear equations with a sparse matrix: #include <iostream> #define VIENNACL_WITH_OPENCL #define ENABLE_OPENCL // Armadillo headers (disable BLAS and LAPACK to avoid linking issues) #define ARMA_DONT_USE_BLAS #define ARMA_DONT_USE_LAPACK #include <armadillo> // IMPORTANT: Must be set prior to any ViennaCL includes if you want to use ViennaCL algorithms on Armadillo objects #define VIENNACL_WITH_ARMADILLO 1 // ViennaCL includes #include "viennacl/vector.hpp"...