Files
Rbm/source/Layer.cpp
T
2019-10-25 17:56:43 +00:00

124 lines
2.5 KiB
C++

/*
* To change this license header, choose License Headers in Project Properties.
* To change this template file, choose Tools | Templates
* and open the template in the editor.
*/
/*
* File: Layer.cpp
* Author: jens
*
* Created on 25. Oktober 2019, 08:13
*/
#include "Layer.hpp"
using namespace std;
Layer::Layer(const string &prjname, size_t id, size_t numVisibleX, size_t numVisibleY, size_t numHidden)
: upper(nullptr)
, lower(nullptr)
, m_prjname(prjname)
, m_id(id)
, m_numVisibleX(numVisibleX)
, m_numVisibleY(numVisibleY)
, m_rbm(m_rbm_params, numVisibleX*numVisibleY, numHidden)
{
m_weightsFile = m_prjname + string(".weights.") + to_string((int)m_id) + string(".dat");
}
Layer::Layer(const Layer& orig)
: upper(nullptr)
, lower(nullptr)
, m_prjname(orig.m_prjname)
, m_weightsFile(orig.m_weightsFile)
, m_id(orig.m_id)
, m_numVisibleX(orig.m_numVisibleX)
, m_numVisibleY(orig.m_numVisibleY)
, m_rbm_params(orig.m_rbm_params)
, m_rbm(orig.m_rbm)
{
}
Layer::~Layer()
{
}
void Layer::saveWeights()
{
FILE *pFile = fopen(m_weightsFile.c_str(), "w");
if (!pFile)
{
std::cout << "Could not open " << m_weightsFile << "!" << std::endl;
return;
}
const arma::mat &bv = rbm().bv();
const arma::mat &bh = rbm().bh();
const arma::mat &w = rbm().w();
size_t numHidden = bh.n_elem;
size_t numVisible = bv.n_elem;
fprintf(pFile, "%d %d %d\n", (int)m_numVisibleX, (int)m_numVisibleY, (int)numHidden);
uint32_t i, j;
for (i=0; i < numVisible; i++)
{
fprintf(pFile, "%3.6f\n", bv(i));
}
for (i=0; i < numHidden; i++)
{
fprintf(pFile, "%3.6f\n", bh(i));
}
for (i=0; i < numVisible; i++)
{
for (j=0; j < numHidden; j++)
{
fprintf(pFile, "%3.6f ", w(i,j));
}
fprintf(pFile, "\n");
}
fclose(pFile);
}
Rbm& Layer::rbm()
{
return m_rbm;
}
Json::Value Layer::toJson() const
{
Json::Value layer;
layer["name"] = string("Layer ") + to_string((int)m_id);
layer["weights_file"] = m_weightsFile;
layer["numVisibleX"] = to_string((int)m_numVisibleX);
layer["numVisibleY"] = to_string((int)m_numVisibleX);
layer["rbm"] = m_rbm.toJson();
return layer;
}
arma::mat Layer::up_pass(const arma::mat &hidden)
{
arma::mat reconstruction = rbm().toVisibleProbs(hidden);
if (upper)
{
return upper->up_pass(reconstruction);
}
return rbm().toHiddenProbs(reconstruction);
}
arma::mat Layer::down_pass(const arma::mat &visible)
{
arma::mat hidden = rbm().toHiddenProbs(visible);
if (lower)
{
return lower->down_pass(hidden);
}
return rbm().toVisibleProbs(hidden);
}