Code: Select all

V = zeros(S,1);
pi = zeros(S,tau,'int8');

Q = zeros(S,A);

for t=1:tau
    disp(t)
    for a=1:A
        Q(:,a) = MDP{a}.R + MDP{a}.T*V;
    end
    [V pi(:,t)] = max(Q,[],2);
end


Code: Select all

void MDP::plan() {
   V = MatrixXf(S,1);
   V.setConstant(0.0);
   policy = MatrixXf(S,tau);
   Q = MatrixXf(S,A);
   Q.setConstant(0.0);

   for (int t = 0; t < tau; t++) {
      for (int a = 0; a < A; a++) {
         Q.col(a) = R[a] + T[a]*(V);
      }
   }
}


Code: Select all


S = largeNumber (~640k);

T = new SparseMatrix<double>[A];
R = new SparseVector<double>[A];

...

T[a] = SparseMatrix<double>(S,S);
R[a] = SparseVector<double>(S);


g++ -I/usr/include/eigen3/ -O0 -g3 -Wall -c -fmessage-length=0 -MMD -MP -MF"MDP.d" -MT"MDP.d" -o "MDP.o" "../MDP.cpp" -std=c++0x -lm
In file included from /usr/include/eigen3/Eigen/Core:313:0,
from /usr/include/eigen3/Eigen/Dense:1,
from ../MDP.h:11,
from ../MDP.cpp:8:
/usr/include/eigen3/Eigen/src/Core/ProductBase.h: In instantiation of ‘struct Eigen::internal::traits<Eigen::ProductBase<Eigen::SparseTimeDenseProduct<Eigen::SparseMatrix<double>, Eigen::Matrix<float, -1, -1> >, Eigen::SparseMatrix<double>, Eigen::Matrix<float, -1, -1> > >’:
/usr/include/eigen3/Eigen/src/SparseCore/SparseDenseProduct.h:133:8: required from ‘struct Eigen::internal::traits<Eigen::SparseTimeDenseProduct<Eigen::SparseMatrix<double>, Eigen::Matrix<float, -1, -1> > >’
/usr/include/eigen3/Eigen/src/Core/DenseBase.h:31:34: required from ‘class Eigen::DenseBase<Eigen::SparseTimeDenseProduct<Eigen::SparseMatrix<double>, Eigen::Matrix<float, -1, -1> > >’
/usr/include/eigen3/Eigen/src/Core/MatrixBase.h:48:34: required from ‘class Eigen::MatrixBase<Eigen::SparseTimeDenseProduct<Eigen::SparseMatrix<double>, Eigen::Matrix<float, -1, -1> > >’
/usr/include/eigen3/Eigen/src/Core/ProductBase.h:63:7: required from ‘class Eigen::ProductBase<Eigen::SparseTimeDenseProduct<Eigen::SparseMatrix<double>, Eigen::Matrix<float, -1, -1> >, Eigen::SparseMatrix<double>, Eigen::Matrix<float, -1, -1> >’
/usr/include/eigen3/Eigen/src/SparseCore/SparseDenseProduct.h:238:7: required from ‘class Eigen::SparseTimeDenseProduct<Eigen::SparseMatrix<double>, Eigen::Matrix<float, -1, -1> >’
../MDP.cpp:162:35: required from here
/usr/include/eigen3/Eigen/src/Core/ProductBase.h:27:98: error: invalid use of incomplete type ‘struct Eigen::internal::scalar_product_traits<double, float>’
In file included from /usr/include/eigen3/Eigen/Core:245:0,
from /usr/include/eigen3/Eigen/Dense:1,
from ../MDP.h:11,
from ../MDP.cpp:8:
/usr/include/eigen3/Eigen/src/Core/util/Meta.h:189:41: error: declaration of ‘struct Eigen::internal::scalar_product_traits<double, float>’
In file included from /usr/include/eigen3/Eigen/Core:270:0,
from /usr/include/eigen3/Eigen/Dense:1,
from ../MDP.h:11,
from ../MDP.cpp:8: