autoTableOfReal FFNet_extractWeights (FFNet me, long layer) {
try {
FFNet_checkLayerNumber (me, layer);
long numberOfUnitsFrom = my nUnitsInLayer[layer - 1] + 1;
long numberOfUnitsTo = my nUnitsInLayer[layer];
autoTableOfReal thee = TableOfReal_create (numberOfUnitsFrom, numberOfUnitsTo);
char32 label[40];
for (long i = 1; i <= numberOfUnitsFrom - 1; i++) {
Melder_sprint (label,40, U"L", layer - 1, U"-", i);
TableOfReal_setRowLabel (thee.peek(), i, label);
}
TableOfReal_setRowLabel (thee.peek(), numberOfUnitsFrom, U"Bias");
for (long i = 1; i <= numberOfUnitsTo; i++) {
Melder_sprint (label,40, U"L", layer, U"-", i);
TableOfReal_setColumnLabel (thee.peek(), i, label);
}
long node = 1;
for (long i = 0; i < layer; i++) {
node += my nUnitsInLayer[i] + 1;
}
for (long i = 1; i <= numberOfUnitsTo; i++, node++) {
long k = 1;
for (long j = my wFirst[node]; j <= my wLast[node]; j++) {
thy data[k++][i] = my w[j];
}
}
return thee;
} catch (MelderError) {
Melder_throw (me, U": no TableOfReal created.");
}
}
FFNet FFNet_and_TabelOfReal_to_FFNet (FFNet me, TableOfReal him, long layer)
{
FFNet thee = NULL; long i, node = 1;
if (! FFNet_checkLayerNumber (me, layer)) return NULL;
if ((my nUnitsInLayer[layer] != his numberOfColumns) ||
(my nUnitsInLayer[layer] == his numberOfColumns && my nUnitsInLayer[layer-1]+1 == his numberOfRows))
{
long i, try[3], rows[3], cols[3], ntry = my nLayers > 3 ? 3 : my nLayers, ok = 0;
if (my nLayers > 3) return Melder_errorp1 (L"Dimensions don't fit.");
for (i = 1; i <= ntry; i++)
{
cols[i] = my nUnitsInLayer[i] == his numberOfColumns;
rows[i] = my nUnitsInLayer[i - 1] + 1 == his numberOfRows;
try[i] = rows[i] && cols[i];
if (try[i]) ok ++;
}
if (! rows[layer])
(void) Melder_error3 (L"The number of rows in the TableOfReal does not equal \n"
"the number of units in the layer that connect to layer ", Melder_integer (layer), L".\n");
else
(void) Melder_error3 (L"The number of columns in the TableOfReal does not equal \n"
"the number of units in layer ", Melder_integer (layer), L".\n");
if (ok == 0)
return Melder_errorp1 (L"Please quit, there is no appropriate layer in the FFNet for this TableOfReal.");
else
{
if (ok == 1)
return Melder_errorp3 (L"Please try again with layer number ",
Melder_integer (try[1] ? try[1] : (try[2] ? try[2] : try[3])), L".");
else
return Melder_errorp1 (L"Please try again with one of the other two layer numbers.");
}
}
thee = Data_copy (me);
if (thee == NULL) return NULL;
for (i = 0; i < layer; i++)
{
node += thy nUnitsInLayer[i] + 1;
}
for (i = 1; i <= thy nUnitsInLayer[layer]; i++, node++)
{
long j, k = 1;
for (j = thy wFirst[node]; j <= thy wLast[node]; j++, k++)
{
thy w[j] = his data[k][i];
}
}
return thee;
}
autoFFNet FFNet_and_TabelOfReal_to_FFNet (FFNet me, TableOfReal him, long layer) {
try {
FFNet_checkLayerNumber (me, layer);
if ((my nUnitsInLayer[layer] != his numberOfColumns) ||
(my nUnitsInLayer[layer] == his numberOfColumns && my nUnitsInLayer[layer - 1] + 1 == his numberOfRows)) {
long trys[1+3], rows[1+3], cols[1+3], ntry = my nLayers > 3 ? 3 : my nLayers, ok = 0;
if (my nLayers > 3) {
Melder_throw (U"Dimensions don't fit.");
}
for (long i = 1; i <= ntry; i++) {
cols[i] = my nUnitsInLayer[i] == his numberOfColumns;
rows[i] = my nUnitsInLayer[i - 1] + 1 == his numberOfRows;
trys[i] = rows[i] && cols[i];
if (trys[i]) {
ok ++;
}
}
if (! rows[layer]){
Melder_throw (U"The number of rows in the TableOfReal does not equal \n"
U"the number of units in the layer that connect to layer ", layer, U".");
} else {
Melder_throw (U"The number of columns in the TableOfReal does not equal \n"
U"the number of units in layer ", layer, U".");
}
if (ok == 0) {
Melder_throw (U"Please quit, there is no appropriate layer in the FFNet for this TableOfReal.");
} else {
if (ok == 1) {
Melder_throw (U"Please try again with layer number ", trys[1] ? trys[1] : (trys[2] ? trys[2] : trys[3]), U".");
} else {
Melder_throw (U"Please try again with one of the other two layer numbers.");
}
}
}
autoFFNet thee = Data_copy (me);
long node = 1;
for (long i = 0; i < layer; i++) {
node += thy nUnitsInLayer[i] + 1;
}
for (long i = 1; i <= thy nUnitsInLayer[layer]; i++, node++) {
long k = 1;
for (long j = thy wFirst[node]; j <= thy wLast[node]; j++, k++) {
thy w[j] = his data[k][i];
}
}
return thee;
} catch (MelderError) {
Melder_throw (me, U": no FFNet created.");
}
}
TableOfReal FFNet_extractWeights (FFNet me, long layer)
{
TableOfReal thee = NULL;
long i, node = 1, numberOfUnitsFrom, numberOfUnitsTo;
wchar_t label[20];
if (! FFNet_checkLayerNumber (me, layer)) return NULL;
numberOfUnitsFrom = my nUnitsInLayer[layer-1] + 1;
numberOfUnitsTo = my nUnitsInLayer[layer];
thee = TableOfReal_create (numberOfUnitsFrom, numberOfUnitsTo);
if (thee == NULL) return NULL;
for (i = 1; i <= numberOfUnitsFrom - 1; i++)
{
swprintf (label,20,L"L%ld-%ld", layer-1, i);
TableOfReal_setRowLabel (thee, i, label);
}
TableOfReal_setRowLabel (thee, numberOfUnitsFrom, L"Bias");
for (i = 1; i <= numberOfUnitsTo; i++)
{
swprintf (label,20,L"L%ld-%ld", layer, i);
TableOfReal_setColumnLabel (thee, i, label);
}
for (i = 0; i < layer; i++) node += my nUnitsInLayer[i] + 1;
for (i = 1; i <= numberOfUnitsTo; i++, node++)
{
long j, k = 1;
for (j = my wFirst[node]; j <= my wLast[node]; j++)
{
thy data[k++][i] = my w[j];
}
}
return thee;
}
