Я использую компилятор mingw в кодовых блоках. Когда я компилирую свой код, он компилируется отлично. Однако, пока я запускаю его, я получаю следующую ошибку. Хотя я изменил файл config.hpp, чтобы разрешить использование Lapck, я все еще получаю эту ошибку. Также использование #define ARMA_DONT_USE_WRAPPER дает ту же ошибку.
ошибка: svd(): необходимо включить использование LAPACK
завершение вызывается после создания экземпляра 'std::logic_error'
what(): svd(): необходимо включить использование LAPACKЭто приложение запросило у среды выполнения необычное завершение его работы. Пожалуйста, свяжитесь со службой поддержки приложения для получения дополнительной информации.
Процесс вернул 3 (0x3) время выполнения: 4,227 с Нажмите любую клавишу, чтобы продолжить.
Ниже приведен мой код:
#include<fstream>
#include<sstream>
#include<iostream>
#include<conio.h>
#include<stdio.h>
#include<process.h>
#include<string>
#include<iomanip>
#include <random>
#include <math.h>
#define ARMA_DONT_USE_WRAPPER
#include <armadillo>
#include <cstdio>
#include <cstdlib>
using namespace std;
using namespace arma;
//function to solve the nmf
double** pnmf(double** Xget, int row, int column)
{ int rows=row;
int columns=column;
int rinit=2;
double X[rows][columns];
for(int i = 0; i < rows; i++) {
for(int j = 0; j < columns; j++){ X[i][j] = Xget[i][j]; }// sample set value;
}
// finding the transpose
double Xt[columns][rows];
for(int j=0; j<columns; j++) {
for(int i=0; i<rows; i++){
Xt[j][i]=X[i][j];
}
}
int k=0;
double XX[rows][columns];
double r1[rows][rinit];
double r2[rows][rinit];
double r3[rinit][rinit];
double r4[rows][rinit];
double r5[rinit][rinit];
double r6[rows][rinit];
double r7[rows][rinit];
double r8[rows][rinit];
double r9[rows][rinit];
double r10[rows][rinit];
for(int i=0;i< rows;i++)
{
for(int j=0; j < columns ;j++)
{
XX[i][j] = 0;
for(k=0;k< columns;k++)
{
XX[i][j] = XX[i][j] + X[i][k] * Xt[k][j];
}
} // end of j sub loop
} // end of i main loop
// initializing W;
double W[rows][rinit];
double Wt[rinit][rows];
double W_Wold[rinit][rows];
double Wold[rows][rinit];
for(int i = 0; i < rows; ++i){
for(int j = 0; j < rinit; j++){
W[i][j] = (double)rand()/(double)RAND_MAX;
}
}
int max_iter =100;
double tol=0.00001;
double sigma[rinit];
double sigma_diag[rinit][rinit];
double bsx[rows][rinit];
for(int iter=0; iter<max_iter; iter++){
// W_old
for(int i = 0; i < rows; ++i){
for(int j = 0; j < rinit; ++j){
Wold[i][j] = W[i][j];
}
}
mat B(rows,rinit);
// calculate transpose of W
for(int j=0; j<rinit; j++) {
for(int i=0; i<rows; i++){
Wt[j][i]=W[i][j];
}
}
// calculate bsxfun
for(int i = 0; i < rows; ++i){
for(int j = 0; j < rinit; ++j){
bsx[i][j] = W[i][j]*W[i][j];
}
}
// calculate sigma
for(int j = 0; j < rinit; ++j){
sigma[j]=0;
for(int i = 0; i < rows; ++i){
sigma[j] = sigma[j]+bsx[i][j];
}
}
//creating diagonal matrix of sigma
for(int i = 0; i < rinit; ++i){
for(int j = 0; j < rinit; ++j){
if (i==j) {sigma_diag[i][j] = sigma[j];}
else {sigma_diag[i][j] = 0; }
}
}
// creating terms to update W
for(int i=0;i< rows;i++)
{
for(int j=0; j < rinit ;j++)
{
r1[i][j] = 0;
for(k=0;k< rows;k++)
{
r1[i][j] = r1[i][j] + XX[i][k] * W[k][j];
}
} // end of j sub loop
}
for(int i=0;i< rows;i++)
{
for(int j=0; j < rinit ;j++)
{
r2[i][j] = 0;
for(k=0;k< rinit;k++)
{
r2[i][j] = r2[i][j] + r1[i][k] * sigma_diag[k][j];
}
} // end of j sub loop
}
for(int i=0;i< rinit; i++)
{
for(int j=0; j < rinit ;j++)
{
r3[i][j] = 0;
for(k=0;k< rows;k++)
{
r3[i][j] = r3[i][j] + Wt[i][k] * r1[k][j];
}
} // end of j sub loop
}
for(int i=0;i< rows;i++)
{
for(int j=0; j < rinit ;j++)
{
r4[i][j] = 0;
for(k=0;k< rinit; k++)
{
r4[i][j] = r4[i][j] + W[i][k] * r3[k][j];
}
} // end of j sub loop
}
for(int i=0;i< rinit;i++)
{
for(int j=0; j < rinit ;j++)
{
r5[i][j] = 0;
for(k=0;k< rows;k++)
{
r5[i][j] = r5[i][j] + Wt[i][k] * W[k][j];
}
} // end of j sub loop
}
for(int i=0;i< rows;i++)
{
for(int j=0; j < rinit ;j++)
{
r6[i][j] = 0;
for(k=0;k< rinit;k++)
{
r6[i][j] = r6[i][j] + r1[i][k] * r5[k][j];
}
} // end of j sub loop
}
for(int i=0;i< rows;i++)
{
for(int j=0; j < rinit ;j++)
{
r7[i][j]=r4[i][j]+r6[i][j];
} // end of j sub loop
}
for(int i=0;i< rows;i++)
{
for(int j=0; j < rinit ;j++)
{
r8[i][j] = 0;
for(k=0;k< rinit;k++)
{
r8[i][j] = r8[i][j] + r7[i][k] * sigma_diag[k][j];
}
} // end of j sub loop
}
for(int i=0;i< rows;i++)
{
for(int j=0; j < rinit ;j++)
{
r9[i][j]=W[i][j]+r8[i][j];
} // end of j sub loop
}
for(int i=0;i< rows;i++)
{
for(int j=0; j < rinit ;j++)
{
r10[i][j]=r2[i][j]/r9[i][j];
} // end of j sub loop
}
for(int i=0;i< rows;i++)
{
for(int j=0; j < rinit ;j++)
{
W[i][j]=W[i][j]*r10[i][j];
} // end of j sub loop
}
// finding norm of W
//assigning W to A
mat A(rows,rinit);
for(int i=0;i< rows;i++)
{
for(int j=0; j < rinit ;j++)
{
A(i,j) = W[i][j];
} // end of j sub loop
}
double norm_W = norm(A);
// final W
for(int i=0;i< rows;i++)
{
for(int j=0; j < rinit ;j++)
{
W[i][j]=W[i][j]/norm_W;
} // end of j sub loop
}
//check convergence
for(int i=0;i< rows;i++)
{
for(int j=0; j < rinit ;j++)
{
W_Wold[i][j]=Wold[i][j] - W[i][j];
} // end of j sub loop
}
// assigning W_Wold to B and W_old to A
for(int i=0;i< rows;i++)
{
for(int j=0; j < rinit ;j++)
{
A(i,j) = Wold[i][j];
} // end of j sub loop
}
for(int i=0;i< rows;i++)
{
for(int j=0; j < rinit ;j++)
{
B(i,j) = W_Wold[i][j];
} // end of j sub loop
}
double diffW= norm(B, "fro")/ norm(A, "fro");
if (diffW<tol)
{
cout<<"\nconverged after"<<iter<<"steps";
break;
}
}
double** table = new double*[rows];
for(int i = 0; i < rows; i++) {
table[i] = new double[rinit];
for(int j = 0; j < rinit; j++){ table[i][j] = W[i][j]; }// sample set value;
}
return table;
}
