Cusp中的實數比例稀疏矩陣向量乘法？

[英]Real scaled Sparse matrix vector multiplication in Cusp?

在cusp中，有一個乘法來計算spmv（稀疏矩陣向量乘法），它需要進行歸約和合並： 從接口看來，對於任何矩陣/矢量乘法，自定義合並和歸約應該是可能的。 我認為cusp除了使用乘法和加號來計算spmv外，還支持使用其他在推力/functional.h中定義的組合和歸約函數。 例如，我可 ...

找不到cusp :: extract_diagonal

[英]cusp::extract_diagonal not found

我將CUSP用於GPU加速代碼的一部分。 我需要從尖峰矩陣C中提取對角線。應在下面的代碼中將其放入對角線。 但是在編譯時，我得到“尖峰”沒有成員“ extract_diagonal” 如果我有： 如在示例中，我找不到format_utils.h。 如果我插入詳細信息以完 ...

CUDA CUSTP線性求解器中的內存使用

[英]Memory use in cuda cusp linear solver

我正在使用cusp :: bicgstab求解線性系統Ax = b ，其中A是MxNxP網格上的3D泊松， x是未知數， b是RHS。 我有一個擁有12GB內存的K40m Tesla 。 我用M = 2000，N = 2000，P = 20 （8000萬個未知數）進行了測試，變量類型為d ...

如何從傳遞的數組中正確構造CUSP coo矩陣

[英]How to properly construct CUSP coo matrix from passed arrays

我正在嘗試將CUSP集成到現有的Fortran代碼中。 現在，我只是想為推力/ CUSP提供基本設置，以便從Fortran中饋入數組，並使用它們來構建CUSP矩陣（當前為coo格式）。 由於有了這個線程，我能夠獲得一個像C例程這樣的包裝器，以編譯成一個庫並將其與Fortran代碼鏈接： un ...

從傳遞的FORTRAN數組中生成CUSP coo_matrix

[英]Generating CUSP coo_matrix from passed FORTRAN arrays

我正在將CUSP求解器集成到現有的FORTRAN代碼中。 第一步，我只是嘗試傳遞一對整數數組和一個來自FORTRAN的浮點數（在FORTRAN中為real * 4），該浮點數將用於構造並打印COO格式的CUSP矩陣。 到目前為止，我已經能夠跟蹤該線程並進行所有編譯和鏈接： 使用IFORT ...

cuSPARSE乘法函數輸出錯誤的值

[英]cuSPARSE multiplication function outputting incorrect value

我正在嘗試使用cuSPARSE庫實現稀疏矩陣乘法。 我已經使用了文檔中的大多數代碼，也就是這里的代碼。 即使我獲得了正確的輸出行指針，列指針，也得到了錯誤的輸出值。 我已按照以下步驟操作：1.生成coo格式矩陣。 2.將其轉換為csr格式3.在輸出矩陣中找到非零元素矢量以及非零元素的 ...

結合使用IFORT和nvcc和CUSP的未解析引用

[英]Unresolved references using IFORT with nvcc and CUSP

我有一個像這樣編譯的程序： 一切工作正常，直到我在其中具有此包裝程序（bicgstab.cu）的地方添加了CUSP支持： Nvcc編譯並生成對象，但是在最后一條命令中，由於鏈接導致我將所有錯誤鏈接在一起時，出現了以下錯誤： 我相信這是因為ifort正在添加或刪除下划線，添 ...

CUSP稀疏庫給出了奇怪的結果

[英]CUSP sparse library giving strange results

我正在嘗試在CUSP中執行以下操作： 是的 另一方面，我從CUSP得到一個奇怪的答案 結果看起來像 最后一個看起來不太正確。 誰知道我做錯了什么？ 我使用的代碼是錯誤的還是我們應該有一個非常好的猜測解決方案+使用預處理器？ ...

帶有OpenMP的CUSP主機功能並行化

[英]CUSP host function parallelization with openMP

CUSP主機功能是否也與openMP並行實現？ 我有點困惑，因為在github討論中它說“ openMP中增加了功能”，但是風口浪尖上沒有明確聲明它已完全用openMP實現。 ...

如何從尖峰庫矩陣格式獲取原始指針

[英]How to get raw pointer from cusp library matrix format

我需要從尖峰庫矩陣格式獲取原始指針。 例如： 如何獲取指向row_indices，column_indices和values數組的原始指針？ 我需要將它們傳遞給內核，並且如果可能的話，我希望避免不必要的數據復制。 ...

使用尖點庫將2個稀疏矩陣相乘

[英]Multiply 2 sparse matrices using cusp library

我剛開始使用cusp庫作為cuda。 我正在嘗試為CUDA實現修訂的單純形算法。 為此，我需要將2個稀疏矩陣相乘以更新基本矩陣。 所以問題是-如何使用尖峰庫將2個稀疏矩陣相乘（實際上什么格式都不重要）？ 還有一種方法可以知道結果矩陣要包含多少個非零元素（出於內存分配的目的）？ ...

在Visual Studio C ++項目中使用cusp

[英]Using cusp in a Visual Studio C++ project

我正在VS2012 Ultimate上使用cusp v.0.4.0和cuda V5.5。 我使用新的項目向導創建了一個CUDA項目，並將尖點路徑添加到其項目properties \\ VC ++ Directories \\ Include Directories 。 我在VS2012生成的 ...

CUSP CG融合

[英]CUSP CG convergence

我使用CUSP共軛梯度法來求解對稱稀疏矩陣。 而且我不知道為什么它不能收斂。 我使用的矩陣尺寸並不大（1K至100K）。 MKL可以很容易地解決相同的線性系統，因此矩陣不會病態。 但是我嘗試添加預處理器，但沒有結果： 對角預處理器和AINV（不完整的Cholesky）使殘留物無限增 ...

從設備內核訪問尖點變量元素

[英]Accessing cusp variable element from device kernel

我在從設備/全局內核訪問和分配具有cusp array1d類型的變量時遇到問題。 隨附的代碼給出了錯誤 下面的代碼 ...

CUSP函數生成具有隨機值的矩陣

[英]CUSP function to generate a matrix with random values

我想知道CUSP庫是否提供了一種功能，該功能可以創建具有特定數量的列，行和任何隨機值的矩陣？ 我找到了poisson5pt函數，但是沒有返回我指定尺寸的矩陣！ 提前致謝 ...

cusp繼承了Class和* .cu設備代碼編譯

[英]cusp inherited Class and *.cu device code compilation

請原諒我，如果這是非平凡的。 我有一個小類，包圍cusp::csr_matrix<int,float,cusp::device_memory>如下所示： 到目前為止，我沒有做任何我想象的需要為設備編譯的東西。 但是我計划在這個特定的類中編寫代碼，作為成員/函數，應該在設備 ...

無法理解CUSP中CSR表示的輸出

[英]Not able to understand output of CSR Representation in CUSP

我正在嘗試使用CUSP庫。 我正在讀取.txt文件，這些文件基本上是稀疏的COO表示。 我正在使用CUSP轉換為CSR格式。 當我用cusp::print()打印矩陣時，它會打印出COO表示的正確結果。 然而，當我將矩陣轉換為CSR時，我已經編寫了自己的打印功能，但結果並不是我想要的 ...

如何從MEX中的CUSP :: csr_matrix返回稀疏矩陣到Matlab？

[英]How to return sparse matrix from CUSP::csr_matrix in MEX to Matlab?

我正在開發用於CUDA的稀疏矩陣計算的MEX文件。 我正在使用CUSP庫。 我不知道如何將cusp :: csr_matrix返回到Matlab。 例如，我有 因此，它是CSR格式的At矩陣，可以說，我已經計算過了。 現在，我必須將其返回給Matlab。 像這樣： 但是 ...

如何將row_indices從COO矩陣復制到CUSP中的推力向量

[英]How do I copy the row_indices from a COO matrix to a thrust vector in CUSP

我一直在嘗試將行索引，列索引和COO矩陣的值復制到單獨的推力向量中，但我發現我無法這樣做。 下面是代碼 sim_row是一個推力向量指針，它已經在以前的代碼中分配了內存，在這里不相關。 代碼編譯，但在運行時失敗並出現錯誤： 在拋出'thrust :: system ...

cusp :: io :: write_matrix_market_file（）提供訪問沖突異常

[英]cusp::io::write_matrix_market_file() gives Access Violation Exception

我正在嘗試使用CUSP將Matlab的272 x 544雙矩陣讀取為ELLPACK格式的稀疏矩陣。 到目前為止，我唯一能找到的方法是將數據讀入mxArray，將mxArray讀入double *數組（使用mxGetPr（）），然后將值復制到cusp :: array2d。 到達那里后，我需要 ...