armadillo

    0熱度

    1回答

    我嘗試一些R代碼裏面翻譯成RcppArmadillo,因此我也想執行以下操作: 假設有一個非負的矢量v和矩陣M,都與例如m行。我希望擺脫矩陣M中的所有行,只要在矢量v的對應行中存在零並且此後也除去矢量v中的零的所有條目。使用R這是單純只是執行以下操作: M = M[v>0,] v = v[v>0] 所以我的問題是,如果有一種方法在RcppArmadillo做到這一點。由於我對任何編程語言都

    1熱度

    1回答

    我試圖在C++中實現Cholesky分解,這是以前在lapack dpotrf_中實現的。 Cholesky分解:R' * R = A 代碼: #include <iostream> #include <armadillo> long my_chol( arma::mat &R, const arma::mat A, long lda ) {

    4熱度

    1回答

    我學會了如何使用Eigen找到矩陣的逆矩陣。但是,當我發現一個數組,它是的功能的輸出I得到一個錯誤 會員在「X」,「逆」,這是無級的請求類型的逆 「雙**」 請幫助我,在使用C++庫來查找矩陣的逆。 我寫的代碼是: #include <iostream> #include <armadillo> #include <cmath> #include <Eigen/Dense> using

    1熱度

    1回答

    我想在Armadillo包的矩陣中使用FFTW。我需要在2D矩陣上執行N個獨立的FFT。按照FFTW和其他資源在線手冊,我有以下代碼: #include <armadillo> #include <iostream> #include "fftw3.h" using namespace std; using namespace arma; fftw_plan fplan; int

    0熱度

    1回答

    我想在R包內編譯C++ 11下的sitmo prng。 The problematic code has been packaged and is available here.此R包的目標是提供sitmo header file以便其他包能夠使用說明中的LinkTo字段。作爲額外的獎勵,該軟件包計劃隨附Armadillo + OpenMP示例。還有一個包,mvnfast,它使用sitmo,但只在

    0熱度

    2回答

    我有一個arma::cube mycube(5,10,15);,我想置換其尺寸,一心只想在MATLAB做: mycube = ones(5,10,15); mycube = permute(mycube,[3 1 2]); size(mycube) % returns (15 5 10) 有沒有辦法做到這一點? 效率太低了嗎? 我實際上想做一個3D FFT,所以我想排列第一個和第三個維度,

    2熱度

    1回答

    我最近安裝了犰狳並嘗試了稀疏矩陣的特徵值問題。不幸的是,分解失敗是參數'N'(下面的代碼)太大e.q.我想知道這裏發生了什麼。矩陣不是很複雜 - 它具有對角結構。 UPDATE 數學也有這個矩陣的問題。它告訴我Arnoldi算法不會收斂。也許我需要在arnoldi arpack例程中手動指定一些參數以確保收斂? 這裏是我的代碼: #include <armadillo> int main()

    3熱度

    1回答

    我想採取arma::vec對象並將其重塑爲arma::cube對象。 例如: vec param(mm*n*g); param.randn(); cube LL = reshape(param,mm,n,g); // this line doesn't work 我能得到這個工作,最簡單的方法是: paramtemp = as<NumericVector>(wrap(param)); c

    0熱度

    1回答

    我有一個Visual Studio社區2015項目(x64解決方案平臺)使用Armadillo線性代數庫(版本6.500.5穩定),如果我嘗試接收構建錯誤使用C++ 11中的語法。該文檔指出「使用C++ 11標準時,Mat,Col,Row中的元素可通過初始化列表」設置。例如,如果我從犰狳文檔describing element initialization,取例如右 // C++11 mat A

    0熱度

    1回答

    我最近試圖在iOS上使用Armadillo來做一些矩陣計算。該應用程序在我的開發iPhone上工作,但Apple試圖在Appstore中發佈它時給了我錯誤消息。看來,犰狳稱爲內部的一些BLAS功能。我用信息搜索了網頁,但沒有發現任何有用的東西。我還發現用「cblas_」前綴調用BLAS函數,例如直接從我的代碼cblas_dgemv不會導致錯誤。但是,這使得使用犰狳毫無意義。 我想知道是否有人遇到同