private static Matrix ReadMatrixFromMatrixFile(BufferedFileReader reader, int networkId)
{
networkId = reader.JumpToNetworkId(networkId, true);
if (networkId == -1)
{
throw new Exception("Network ID does not exist!");
}
reader.ReadLine(); // Skip first line with the network id
string[] colLabels = reader.ReadLine().Split(',');
//Yushan
List <string> rowLabels = new List <string>();
//string[] colLabels = tempLabels.Split(',');
int rows = reader.CountLines(networkId) - 2; // Subtract off header columns
//int rows = reader.CountLinesAlt(tempIdStr) - 2;
int cols = colLabels.Length - 1;
Matrix matrix = new Matrix(rows, cols);
matrix.NetworkId = networkId;
matrix.NetworkIdStr = reader.GetNetworkRealId(networkId);
matrix.ColLabels.SetLabels(colLabels);
for (int row = 0; row < rows; ++row)
{
string temp = reader.ReadLine();
if (temp == null)
{
break;
}
//string[] parts = reader.ReadLine().Split(',');
string[] parts = temp.Split(',');
if (parts.Length > cols + 1) // one extra for header
{
throw new FileLoadException("Matrix file has too many entries for network id: " + matrix.NetworkIdStr + ", row " + parts[0]);
}
if (parts.Length == 0)
{
throw new FileLoadException("Matrix file has no entries for network id: " + matrix.NetworkIdStr);
}
//matrix.RowLabels[row] = parts[0];
rowLabels.Add(parts[0]);
for (int i = 1; i < parts.Length; ++i)
{
matrix[row, i - 1] = ExtractDouble(parts[i]);
}
}
matrix.RowLabels.SetLabels(rowLabels);
reader.closeStream();
reader.Dispose(); /* Change made 12/3/2010 - PM */
return(matrix);
}
public static List <Matrix> ReadMatrixFromMultipleDyadicFile(string filename, int networkId)
{
BufferedFileReader reader = BufferedFileTable.GetFile(filename);
lock (reader)
{
networkId = reader.JumpToNetworkId(networkId, true);
//Dictionary<string, int> labels = reader.GetDyadicLabels(networkId);
//int rows = labels.Count;
// initialize the matrices
List <Matrix> multipleMatrices = new List <Matrix>();
List <string> rowLabels = new List <string>();
List <string> colLabels = new List <string>();
//string[] topLabels = reader.TopLine.Split(',');
//string[] labels = new string[reader.CountVarsInDyadicFile()];
int varCnt = reader.CountVarsInDyadicFile();
List <double> flatMatrixList = new List <double>();
//for (int var = 0; var < reader.CountVarsInDyadicFile(); var++)
//{
// Matrix matrix = new Matrix(rows, rows);
// matrix.NetworkId = networkId;
// matrix.Name = topLabels[var + 3];
// matrix.RowLabels.SetLabels(labels.Keys);
// matrix.ColLabels.SetLabels(labels.Keys);
// multipleMatrices.Add(matrix);
//}
int totalLines = reader.CountLines(networkId);
for (int i = 0; i < totalLines; ++i)
{
string s = reader.ReadLine();
string[] parts = s.Split(',');
if (!rowLabels.Contains(parts[1]))
{
rowLabels.Add(parts[1]);
}
if (!colLabels.Contains(parts[2]))
{
colLabels.Add(parts[2]);
}
//if (parts.Length < 3 + reader.CountVarsInDyadicFile())
// throw new FileLoadException("Missing value for line: " + s);
for (int var = 0; var < varCnt; var++)
{
//multipleMatrices[var][labels[parts[1]], labels[parts[2]]] = ExtractDouble(parts[3 + var]);
flatMatrixList.Add(ExtractDouble(parts[3 + var]));
}
}
int rows = rowLabels.Count;
int cols = colLabels.Count;
for (int var = 0; var < varCnt; var++)
{
Matrix matrix = new Matrix(rows, cols);
matrix.RowLabels.SetLabels(rowLabels);
matrix.ColLabels.SetLabels(colLabels);
matrix.NetworkId = networkId;
matrix.NetworkIdStr = reader.GetNetworkRealId(networkId);
multipleMatrices.Add(matrix);
}
for (int i = 0; i < rows; i++)
{
for (int j = 0; j < cols; j++)
{
for (int var = 0; var < varCnt; var++)
{
multipleMatrices[var][i, j] = flatMatrixList[i * cols * varCnt + j * varCnt + var];
}
}
}
/*
* if (networkId < 1000)
* networkId = int.Parse("1" + networkId);
* else
* networkId = int.Parse("2" + networkId);
*/
return(multipleMatrices);
}
}
private static Matrix ReadMatrixFromDyadicFile(BufferedFileReader reader, int networkId, int dyadicVariable)
{
networkId = reader.JumpToNetworkId(networkId, true);
// Yushan
List <string> rowLabels = new List <string>();
List <string> colLabels = new List <string>();
List <double> flatMatrix = new List <double>();
//Dictionary<string, int> labels = reader.GetDyadicLabels(networkId);
//int rows = labels.Count;
//Matrix matrix = new Matrix(rows, rows);
//matrix.RowLabels.SetLabels(labels.Keys);
//matrix.ColLabels.SetLabels(labels.Keys);
int totalLines = reader.CountLines(networkId);
//Console.WriteLine("Total Lines: {0}", totalLines);
for (int i = 0; i < totalLines; ++i)
{
string s = reader.ReadLine();
string[] parts = s.Split(',');
if (parts.Length < 3 + dyadicVariable)
{
throw new FileLoadException("Missing value for line: " + s);
}
// Yushan
if (!rowLabels.Contains(parts[1]))
{
rowLabels.Add(parts[1]);
}
if (!colLabels.Contains(parts[2]))
{
colLabels.Add(parts[2]);
}
flatMatrix.Add(ExtractDouble(parts[3 + dyadicVariable]));
//matrix[labels[parts[1]], labels[parts[2]]] = ExtractDouble(parts[3 + dyadicVariable]);
}
/*
* if (networkId < 1000)
* networkId = int.Parse("1" + networkId);
* else
* networkId = int.Parse("2" + networkId);
*/
int rows = rowLabels.Count;
int cols = colLabels.Count;
Matrix matrix = new Matrix(rows, cols);
matrix.NetworkId = networkId;
matrix.NetworkIdStr = reader.GetNetworkRealId(networkId);
matrix.RowLabels.SetLabels(rowLabels);
matrix.ColLabels.SetLabels(colLabels);
for (int i = 0; i < rows; i++)
{
for (int j = 0; j < cols; j++)
{
matrix[i, j] = flatMatrix[i * cols + j];
}
}
return(matrix);
}
