Thanks for you fast reply. You are right: pardiso uses 1-base indexed arrays. After converting to 0-base indexed arrays ILU and gmres will work without problems. If have another question. Is it possible to use a parallel preconditioner and initialize the compressed_matrix with arrays from GPU memory, so I could save the copy time from GPU-CPU-GPU?
Hello, I will use the gmres solver for solving my linear equation system on GPU (CUDA). I have the 3 Arrays of the CSR Sparse format as pure C-Arrays (int, int, double) and try to initialize the vcl_compressed_matrix with the set function. typedef double ScalarType; viennacl::compressed_matrix<ScalarType> vcl_matrix; vcl_matrix.set(ia,ja,val,ld,ld,nnz); But after solving with gmres (using ILU preconditioner) I get a vector with only nan entries. Is there something wrong with my way of initializing?...
Hello, I will use the gmres solver for solving my linear equation system on GPU (CUDA). I have the 3 Arrays of the CSR Sparse format as pure C-Arrays (int, int, double) and try to initialize the vcl_compressed_matrix with the set function. typedef double ScalarType; viennacl::compressed_matrix<ScalarType> vcl_matrix; vcl_matrix.set(ia,ja,val,ld,ld,nnz); But after solving with gmres (using ILU preconditioner) I get a vector with only nan entries. Is there something wrong with my way of initializing?...
Hello, I will use the gmres solver for solving my linear equation system on GPU (CUDA). I have the 3 Arrays of the CSR Sparse format as pure C-Arrays (int, int, double) and try to initialize the vcl_compressed_matrix with the set function. viennacl::compressed_matrix<ScalarType> vcl_matrix; vcl_matrix.set(ia,ja,val,ld,ld,nnz); But after solving with gmres (using ILU preconditioner) I get a vector with only nan entries. Is there something wrong with my way of initializing? A few more information:...
Hello, I will use the gmres solver for solving my linear equation system on GPU (CUDA). I have the 3 Arrays of the CSR Sparse format as pure C-Arrays (int, int, double) and try to initialize the vcl_compressed_matrix with the set function. iennacl::compressed_matrix<ScalarType> vcl_matrix; vcl_matrix.set(ia,ja,val,ld,ld,nnz); But after solving with gmres (using ILU preconditioner) I get a vector with only nan entries. Is there something wrong with my way of initializing? A few more information: I...